Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 35,858 T→A 42.6% intergenic (+89/+19) USA300HOU_0025 → / ← USA300HOU_0026 hypothetical protein/IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073035,8580TA42.6% 80.8 / 92.9 76intergenic (+89/+19)USA300HOU_0025/USA300HOU_0026hypothetical protein/IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (21/21);  new base A (8/24);  total (29/47)
Fisher's exact test for biased strand distribution p-value = 3.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

AAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTT  >  CP000730/35719‑35994
                                                                                                                                           |                                                                                                                                         
aaaGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAAttt                                                                                                                                          <  2:118183/141‑1 (MQ=255)
aaaGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTGGGCGAAAAATGATTACTTGAACAAGTGCTAAAATAAAATTGTGATCTATTAACAATTGTAGGTTTTGATAAAAAAAAAAAAAATATAGCTAACCACTAAttt                                                                                                                                          <  2:14841/141‑1 (MQ=255)
aaaGTATTTAATGTAGAAGCCATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGAACTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAAttt                                                                                                                                          <  2:178454/141‑1 (MQ=255)
 aaGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTa                                                                                                                                         >  1:177984/1‑141 (MQ=255)
      tttAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATg                                                                                                                                    <  1:323201/141‑1 (MQ=255)
      tttAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATg                                                                                                                                    >  2:134912/1‑141 (MQ=255)
      tttAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATg                                                                                                                                    >  2:289885/1‑141 (MQ=255)
       ttAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGt                                                                                                                                   <  2:44254/141‑1 (MQ=255)
          aTGTACAAGCCACTCAGTATAAAGTTTTTGGCCAAAAATGAGTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGTTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAg                                                                                                                                <  1:48866/141‑1 (MQ=37)
                   cTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCtt                                                                                                                       <  1:207222/141‑1 (MQ=255)
                   cTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCtt                                                                                                                       <  1:176168/141‑1 (MQ=255)
                        cAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACtt                                                                                                                  >  1:61403/1‑141 (MQ=255)
                          ggataAAGTTTTTGGCGGAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGc                                                                                                                <  2:23588/139‑1 (MQ=255)
                          ggataAAGTTTATGGAGAAAAATGATTACATGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTATGAAGAAAAAAAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGc                                                                                                                <  2:182102/139‑1 (MQ=25)
                          ataggaaGTTTTTGGCGGAAAATGATGAATTGCAGTAGGGCCAAATTAAAATTGAAAAATAAAAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGc                                                                                                                <  2:208359/136‑1 (MQ=17)
                           tataAAGTTTTTGGCGAAAAATGATTACTTGAACACGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAACAATATCGCTAACCACTAATTTTTCATGTCAGTGTTCGCTTAACTTGCt                                                                                                               <  2:280519/141‑1 (MQ=255)
                           tataAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCt                                                                                                               <  1:99099/141‑1 (MQ=255)
                               aaGTTTTTGGCGAAACATGCTTACTTGAACAAGAGCTAAAATAAAATTGTGATCTACTAAAACCAGAAGGCTCTGTTGCACAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTCTTCGCTTAACTTGCTaccc                                                                                                           >  2:349042/1‑138 (MQ=21)
                               aaGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAACAATAGAAGGTTCTGTTGCAAAGTAAAAAAATGTAGCTAACCACTAATTTATCATGTCAGTGTTCGCTGAACTTGCTAGCa                                                                                                           >  2:168274/1‑141 (MQ=255)
                                    tttGGCGAAAAATGATAACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatgatg                                                                                                      >  1:195544/1‑141 (MQ=255)
                                    caggggggtagaTGATTACTATAACAAGAGCTAAAATAAAATTGTGATCTAATAACAATAGAAGGTTCTGTTGAAAAGTAAAAAAATATAGCTAACCACTAATGTAGCATGTCAGTGTTCGCTTAACTTGCTAGCatgatg                                                                                                      <  2:94336/130‑1 (MQ=18)
                                     ttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGc                                                                                                     <  1:343172/141‑1 (MQ=255)
                                     taggaaaaaaaTGACTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTAAAAAAAAAAAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGc                                                                                                     <  2:99470/135‑1 (MQ=21)
                                          gAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGCCAGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTGGCTAGCATGATGCTAAtt                                                                                                >  1:72633/1‑141 (MQ=255)
                                          gAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGACTTCGCGTACCTTGCTAGCATGATGCTAAtt                                                                                                >  2:239565/1‑141 (MQ=255)
                                            aaaaTGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTc                                                                                              >  2:50591/1‑141 (MQ=255)
                                                         gaacaaagagtaaaataaaaTTGTAATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGaa                                                                                 <  2:303541/131‑1 (MQ=255)
                                                              agagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATcc                                                                            >  1:320738/1‑141 (MQ=255)
                                                              agagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGCCACTCGTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATcc                                                                            >  1:358312/1‑141 (MQ=255)
                                                               gagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAGCTTGCTAGCATGATGCTAATTTCGTGGCCTGGCGAAAATCCg                                                                           >  1:2429/1‑141 (MQ=255)
                                                               gagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCg                                                                           >  2:259048/1‑141 (MQ=255)
                                                                agCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGt                                                                          <  2:26462/141‑1 (MQ=255)
                                                                       acaaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                                                                 <  1:3457/127‑1 (MQ=11)
                                                                       acaaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                                                                 >  1:28203/15‑141 (MQ=11)
                                                                        taaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCGTGGCGAAAATCCGTACATCTGa                                                                  >  1:256736/1‑141 (MQ=255)
                                                                        caaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                                                                <  1:54206/128‑1 (MQ=11)
                                                                        caaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATAACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                                                                >  1:230165/14‑141 (MQ=11)
                                                                          aaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAagag                                                              <  1:31791/130‑1 (MQ=11)
                                                                           aagagtattttaataaAAATGATGGTGCTTGTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAgaga                                                             <  1:362686/131‑1 (MQ=11)
                                                                            agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATTTATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                           <  1:50098/132‑1 (MQ=11)
                                                                            agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                           <  1:60447/132‑1 (MQ=11)
                                                                            agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                           <  1:136768/132‑1 (MQ=11)
                                                                            agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                           <  1:236434/132‑1 (MQ=11)
                                                                            agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                           <  2:53954/132‑1 (MQ=11)
                                                                            agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGACCTGAAGAGAc                                                           >  1:288450/10‑141 (MQ=11)
                                                                             gagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                                                          <  1:361002/133‑1 (MQ=11)
                                                                               ttgATTTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                                                          <  2:37991/140‑1 (MQ=255)
                                                                               gtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                                        <  1:243031/135‑1 (MQ=11)
                                                                                tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTCGATCTGAAGAGACCTGc                                                       <  2:192795/136‑1 (MQ=11)
                                                                                tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                                       <  2:196476/136‑1 (MQ=11)
                                                                                tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                                       <  1:366229/136‑1 (MQ=11)
                                                                                 taTCTAATAAAAATAGAAGGGTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATGTATCATGTCAGTGTTCGCTTAACTTGCTAGCATTATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                                        <  1:364514/140‑1 (MQ=255)
                                                                                 attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                                      >  2:312142/5‑141 (MQ=11)
                                                                                 attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                                      >  2:258630/5‑141 (MQ=11)
                                                                                  ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                                     <  1:245200/138‑1 (MQ=11)
                                                                                  ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                                     <  1:81635/138‑1 (MQ=11)
                                                                                  ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCGAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                                     <  2:291873/138‑1 (MQ=11)
                                                                                  aTCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                                       >  1:216593/1‑141 (MQ=255)
                                                                                   tttaataaAAATTATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTCATATATCATGTCAGTGTTCGTTTAACCTGCTAGTTTGATGCTAATTTCGTGGTATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                                    <  1:174426/139‑1 (MQ=11)
                                                                                   tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                                    <  2:313926/139‑1 (MQ=11)
                                                                                   tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                                    >  2:312786/3‑141 (MQ=11)
                                                                                   tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTCCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                                    <  1:162065/139‑1 (MQ=11)
                                                                                   tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCAATAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGAGCTGAAGAGACCTGCGGt                                                    >  1:75182/3‑141 (MQ=11)
                                                                                   tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATACAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                                    <  2:208841/139‑1 (MQ=11)
                                                                                   tCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTAAagtggcgcgcg                                                       >  2:57802/1‑132 (MQ=255)
                                                                                    ttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGtt                                                   <  2:68014/140‑1 (MQ=11)
                                                                                     taataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTc                                                  <  2:244793/141‑1 (MQ=11)
                                                                                      aataaAAATGGAAGTTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTGATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGtt                                                   <  2:182580/141‑1 (MQ=255)
                                                                                       ataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGTAGAGACCTGCGGTTCtt                                                >  2:257491/1‑141 (MQ=11)
                                                                                        taaAAATAGAAGGTTCTGTTGCAAAGTAAACAAATATAGCTAACCACTAATTTATCACGTCCGTGTTCGCTTAACTTGCTAGCAGGAAGCTAATTTCGTGGCAGGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGGTCt                                                 >  2:39149/1‑141 (MQ=11)
                                                                                        taaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATTTAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCt                                                 >  2:47803/1‑141 (MQ=38)
                                                                                        taaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCt                                                 >  1:277386/1‑141 (MQ=38)
                                                                                           aaaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGCGTTCGCTCAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtttt                                              <  1:182349/141‑1 (MQ=25)
                                                                                             aTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAGCTGCGGTTCTTTTta                                            <  1:284926/141‑1 (MQ=18)
                                                                                              ccgAAGTTTCTTTTGCAAAGTAAAAAATATAGCTACTCACTAATATATCATGTCAGTGTTCGTTTAGCTTGCTAGCAGGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtata                                          <  2:284489/139‑1 (MQ=11)
                                                                                               aGAAGGTTCTGTGGCAAAGTAAAAAAATATAGCTAACAACTGATTTATCATATCAGTGTTCGATTAACTTACTAGCATGATGCTAATTCCGTGGCATGGCGAAAGTCCGTAGATCTGAAGAGACCTGCGGTGCttttgaga                                          >  2:95131/1‑137 (MQ=11)
                                                                                                 aaGTTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCCAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtatata                                        <  2:256738/141‑1 (MQ=11)
                                                                                                 aaGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGTCCTGCGGTTCTTTTtatata                                        <  1:379500/141‑1 (MQ=11)
                                                                                                                               ccaatcACTAATATTGCATGGCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAg          <  1:10667/136‑1 (MQ=1)
                                                                                                                                       tactttCTCATGTAAGTATTCGCTTAACTTGCGAGCATGATGCTAATTTCGTGGGATTGCGAAAAACCGTAGAACTAAAAGATCCTTCAATTCTTTGTATATCCTGCGTAAATACACTCAATACCTTTTTAATTATTCttt  <  2:174431/138‑1 (MQ=0)
                                                                                                                                           |                                                                                                                                         
AAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTT  >  CP000730/35719‑35994

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: