Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,978 A→C 22.6% intergenic (‑188/+120) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9780AC22.6% 676.1 / 187.5 314intergenic (‑188/+120)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base A (107/120);  new base C (59/2);  total (192/122)
Fisher's exact test for biased strand distribution p-value = 2.90e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-02

Only 100 of 320 total aligned reads displayed.

GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCA  >  NC_000913/3705879‑3706076
                                                                                                   |                                                                                                       
gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                         >  1:2796947/1‑100 (MQ=255)
gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                         >  2:3772715/1‑100 (MQ=255)
 aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                        >  1:3458102/1‑100 (MQ=255)
        gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGAGGCGTTAACGCCTTa                                                                                                 >  2:4205702/1‑100 (MQ=255)
        gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                 >  1:2252936/1‑100 (MQ=255)
          tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                               >  2:1922185/1‑100 (MQ=255)
           gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                              >  1:815914/1‑100 (MQ=255)
             cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                            >  2:2086148/1‑100 (MQ=255)
              ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                           >  2:2285025/1‑100 (MQ=255)
               gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                          <  2:3968923/100‑2 (MQ=37)
               gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGGTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                          <  1:3149483/100‑1 (MQ=255)
                aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                         <  2:557144/100‑1 (MQ=255)
                aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                         <  2:646497/100‑3 (MQ=37)
                    cgttGCACGCAGAGGTACAGCGCCCCCGGGTATATAGATTGTTGTCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                     <  1:3858042/96‑1 (MQ=21)
                    cgcTGTTCGCAGAGGCACAGCGCTACCGGGTATACAGATTGGTGACAAAGTGCGCTTTGTTTATTCCGGATAAGGCGTTATCGCCGCATCCGGCATaaac                                                                                     <  1:1002361/100‑4 (MQ=21)
                     gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                                   >  1:1268005/1‑100 (MQ=255)
                       tGCAGGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                                 <  1:2319343/100‑1 (MQ=255)
                       tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                                 <  2:1124123/100‑1 (MQ=255)
                          aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                              >  2:1196977/1‑100 (MQ=255)
                           cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                             <  2:1050710/100‑1 (MQ=255)
                           cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                             <  1:401646/100‑1 (MQ=255)
                           cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                             <  1:1867245/100‑1 (MQ=255)
                            ggagagGCACAGCACCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                            >  2:2876369/3‑100 (MQ=255)
                            gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                            <  2:1994244/100‑1 (MQ=255)
                                 ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                       >  1:2217753/1‑100 (MQ=255)
                                 gacacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                       >  2:3654339/3‑100 (MQ=255)
                                   cacgGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                     >  2:1513226/1‑100 (MQ=255)
                                    acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCaa                                                                    <  1:2889123/100‑1 (MQ=255)
                                     caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                   <  1:505650/100‑1 (MQ=255)
                                     caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                   <  1:2854196/100‑1 (MQ=255)
                                           cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                             >  2:1687148/1‑100 (MQ=255)
                                           cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                             <  2:2369967/100‑1 (MQ=255)
                                            aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                            <  2:2502954/100‑1 (MQ=255)
                                              cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑g                                                         <  1:3031665/100‑1 (MQ=255)
                                              cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑g                                                         <  1:3330410/100‑1 (MQ=255)
                                                      cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCg                                              <  2:2109352/100‑1 (MQ=255)
                                                        gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTa                                            <  2:1379094/100‑1 (MQ=255)
                                                          atcttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAgg                                          <  1:1686241/97‑1 (MQ=255)
                                                            gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGcc                                        >  1:727115/1‑100 (MQ=255)
                                                             ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCt                                       >  1:1478996/1‑100 (MQ=255)
                                                               gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGa                                     <  1:3536186/100‑1 (MQ=255)
                                                               gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGa                                     >  2:2911931/1‑100 (MQ=255)
                                                                 cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATa                                   >  2:3252168/1‑100 (MQ=255)
                                                                     gTGCGCTTTGTTTATGCCGGATGTGGCGTTAACTCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCg                               <  2:1842329/100‑1 (MQ=255)
                                                                     gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCg                               >  2:2175300/1‑100 (MQ=255)
                                                                      tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCgg                              >  1:2823462/1‑100 (MQ=255)
                                                                      tGCGCGTTGTTTATGCCGGATTCGGAGTTATCGGCTTATCTGGCCTAC‑AAAATCGTGCAAACTCAAGACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCgg                              <  1:2624705/100‑1 (MQ=18)
                                                                         gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAgc                           <  1:1567754/100‑1 (MQ=255)
                                                                         gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAgc                           <  1:1911322/100‑1 (MQ=255)
                                                                          cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAgcg                          <  1:1531116/100‑1 (MQ=255)
                                                                          cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAgcg                          >  2:289564/1‑100 (MQ=255)
                                                                          cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAgcg                          <  1:2047092/100‑1 (MQ=255)
                                                                          cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAgcg                          <  2:1782461/100‑1 (MQ=255)
                                                                            ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCa                        >  2:3497768/1‑100 (MQ=255)
                                                                            atttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑AT‑CCCTTGTAGGCCTGATAAGCGTAGCGCAt                       >  1:3903465/9‑100 (MQ=18)
                                                                             tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGCA‑AT‑TC‑T‑CTTGTAGGCCTGATAAGCGTAGCGCAt                       >  2:27686/8‑100 (MQ=18)
                                                                             tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCAt                       >  1:1600225/1‑100 (MQ=255)
                                                                             tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCAt                       >  1:2553322/1‑100 (MQ=255)
                                                                              gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATc                      >  2:2683847/1‑100 (MQ=255)
                                                                               tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GA‑GA‑T‑TATGTAGGCCTGATAAGCGTAGCGCATCa                     >  1:639892/6‑100 (MQ=17)
                                                                                gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCC‑CTTGTAGGCCTGATAAGCGTAGCGCATCAgg                   >  2:679164/5‑100 (MQ=18)
                                                                                 taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAgg                   >  2:251099/4‑100 (MQ=18)
                                                                                 taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGAT‑C‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAgg                   >  2:271092/4‑100 (MQ=17)
                                                                                 tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAGATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAgg                   >  1:2997110/1‑100 (MQ=255)
                                                                                  aggCCGGATAAGGCGTTTACGCCGCACCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑T‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGc                  >  2:2596725/3‑100 (MQ=17)
                                                                                  aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                  >  2:384718/3‑100 (MQ=18)
                                                                                  aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑AA‑A‑CACGTAGGCCTGATAAGCGAAGCGCATCAGGc                  >  1:955708/3‑100 (MQ=18)
                                                                                  aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGc                  >  2:3035489/3‑100 (MQ=18)
                                                                                  aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGc                  <  1:1432903/100‑1 (MQ=255)
                                                                                  aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGc                  <  1:1191859/100‑1 (MQ=255)
                                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑C‑CTTCTAGGCCTGATAAGCGTAGCGCATCAGGCa                 >  1:2564641/2‑100 (MQ=17)
                                                                                   cgctggataaggcgTTCACGCAGCATCCAGCCTAC‑AAAATTGTGCAAAATCAATGGATTGCA‑CA‑TT‑T‑AGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                 >  1:364030/11‑100 (MQ=18)
                                                                                    gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                 >  1:3434330/1‑99 (MQ=255)
                                                                                    gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt                >  2:3743114/1‑100 (MQ=255)
                                                                                    gCCGGATAAGGCGTTTACGCCTTATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑C‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCag                >  2:4055131/1‑99 (MQ=18)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑C‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCag                >  1:2953754/1‑99 (MQ=17)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑A‑TTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaa                >  2:3164096/1‑99 (MQ=25)
                                                                                     ccggataaggcgTTTACGCCGCATCCGGTCTAC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CT‑T‑CTGGTAGGCCTGATAAGCATCGCGCATCAGGCaag               >  2:15831/9‑98 (MQ=18)
                                                                                     ccggataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt               >  2:1889770/9‑100 (MQ=18)
                                                                                      cggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑T‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt              >  1:1786424/8‑100 (MQ=17)
                                                                                      cggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtt              >  2:3410148/8‑97 (MQ=18)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt             >  1:1127275/7‑100 (MQ=18)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt             >  1:3477226/7‑100 (MQ=18)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑C‑CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAattc             >  2:1327207/7‑99 (MQ=21)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑C‑CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAattc             >  1:3066090/7‑99 (MQ=21)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGAT‑T‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt             >  2:494128/7‑100 (MQ=17)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GA‑GA‑T‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCAattt             >  1:1320148/7‑100 (MQ=17)
                                                                                       ggATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt             <  1:1960470/100‑1 (MQ=255)
                                                                                       ggATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt             >  2:3306813/1‑100 (MQ=255)
                                                                                        gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCAGGG‑AC‑T‑AT‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctta            >  2:3466834/6‑95 (MQ=21)
                                                                                        gataacgcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCC‑GA‑AA‑C‑CCTGTAGGCCTCAACAGCGGCGCGGATCACGc                  >  1:2085452/7‑94 (MQ=9)
                                                                                         aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg           <  2:4112468/100‑1 (MQ=255)
                                                                                         aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg           <  2:2629454/100‑1 (MQ=255)
                                                                                          taaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑T‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgc          >  1:4078602/4‑100 (MQ=17)
                                                                                          tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAGGCGGAGCGCATCAGGCATTtttgc          >  2:1977453/1‑100 (MQ=255)
                                                                                          tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc          >  2:2425261/1‑100 (MQ=255)
                                                                                           gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACACTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc         <  1:3745252/100‑1 (MQ=255)
                                                                                                gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc    >  2:771379/1‑100 (MQ=255)
                                                                                                 ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc    <  2:871169/99‑1 (MQ=255)
                                                                                                  tAACGCCTTATCCGGCCTAC‑AAAATCGTACAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa  >  2:4273223/1‑100 (MQ=255)
                                                                                                   |                                                                                                       
GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCA  >  NC_000913/3705879‑3706076

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