Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,028 T→C 23.4% intergenic (‑238/+70) 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,706,0280TC23.4% 480.5 / 143.2 248intergenic (‑238/+70)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (94/90);  new base C (54/2);  total (156/92)
Fisher's exact test for biased strand distribution p-value = 1.96e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 252 total aligned reads displayed.

TACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGAT  >  NC_000913/3705931‑3706126
                                                                                                       |                                                                                                    
taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑At‑ttt                                                                                                    >  1:2480225/1‑100 (MQ=255)
 aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTc                                                                                                   <  1:3467121/100‑1 (MQ=255)
  cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑g                                                                                                >  1:1271249/1‑100 (MQ=255)
  cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑g                                                                                                >  1:3224683/1‑100 (MQ=255)
  cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑g                                                                                                <  1:2800359/100‑1 (MQ=255)
  cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑g                                                                                                >  2:3218815/1‑100 (MQ=255)
      ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAgg                                                                                            >  1:3510846/1‑100 (MQ=255)
       tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGc                                                                                           <  2:2242086/100‑1 (MQ=255)
         ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCt                                                                                         >  1:1226680/1‑100 (MQ=255)
          tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTg                                                                                        >  1:3438451/1‑100 (MQ=255)
          tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTg                                                                                        >  1:2712378/1‑100 (MQ=255)
                  tGCGCGTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCgg                                                                                <  1:1307701/100‑1 (MQ=255)
                     gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGa                                                                               <  1:2325865/98‑1 (MQ=255)
                       cttttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAgcgc                                                                           >  2:2007855/10‑100 (MQ=17)
                        ttttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCa                                                                          >  2:1505621/9‑100 (MQ=18)
                        ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                                          >  2:3356334/1‑100 (MQ=255)
                        ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                                          <  1:373835/100‑1 (MQ=255)
                        ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                                          <  2:1240099/100‑1 (MQ=255)
                         cttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TG‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                                         >  1:3669143/8‑100 (MQ=17)
                          ttgtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGC‑‑A‑GA‑AA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATc                                                                        >  2:2307935/7‑100 (MQ=18)
                          ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATc                                                                        >  2:1720692/7‑100 (MQ=18)
                          gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                                        <  2:1256428/100‑1 (MQ=255)
                           tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                                                                       >  2:2341635/1‑100 (MQ=255)
                            gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                     >  1:1548395/5‑100 (MQ=18)
                            gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                      >  2:1563847/5‑100 (MQ=21)
                            gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑‑A‑GG‑AG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                      >  2:2838790/5‑100 (MQ=18)
                              nggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑ATACCGCGCAAATTCAATATATTGCGGA‑GA‑AA‑AT‑‑‑‑GTAGGCCTGATAAGCGTAGCGCGTCAGGct                                                                   >  1:2043341/3‑99 (MQ=18)
                              aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGC‑‑A‑AT‑TC‑TCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                     >  2:1063249/3‑99 (MQ=17)
                              aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                     >  1:1429532/3‑99 (MQ=18)
                               tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                   <  2:406035/100‑1 (MQ=255)
                               ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                   >  2:1737619/2‑100 (MQ=17)
                               ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                   >  1:386876/2‑100 (MQ=14)
                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑A‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                   >  2:1389603/2‑100 (MQ=18)
                                gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                  >  1:1362251/1‑99 (MQ=18)
                                gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑A‑‑A‑AT‑CAACATGTAGGCCGGATAAGCGCAGCGCATCAGGCa                                                                   >  1:2571710/1‑100 (MQ=18)
                                 ctggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGC‑‑A‑GAGAT‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                                 <  2:2915479/92‑3 (MQ=17)
                                 ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TG‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                 >  1:3558594/9‑98 (MQ=17)
                                 ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                                 >  1:1152417/9‑98 (MQ=21)
                                 ccggataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                                 >  2:1799584/9‑100 (MQ=18)
                                  cGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGATCGCATCAGGCAtt                                                                 >  2:1094665/1‑99 (MQ=255)
                                   ggataaggcgTTTACGCCTTATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGC‑‑ATGG‑GC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                               >  2:148673/7‑100 (MQ=18)
                                   ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                              >  2:1381571/7‑98 (MQ=18)
                                   ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                              >  1:3778619/7‑98 (MQ=21)
                                   ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                              >  1:1032296/7‑98 (MQ=18)
                                   ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑A‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                               >  1:2235449/7‑100 (MQ=18)
                                   ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑A‑GGGAC‑TAT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                               >  1:793883/7‑96 (MQ=18)
                                   ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑A‑GG‑GA‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                               >  1:220340/7‑96 (MQ=18)
                                   ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGC‑‑A‑GAGAT‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatag                                                               >  1:2915465/7‑96 (MQ=17)
                                   ggataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCTGAGCGCATCAGGCAtt                                                                 >  1:3122850/7‑98 (MQ=17)
                                     ataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TG‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttg                                                             <  2:2145583/96‑1 (MQ=17)
                                     aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                             <  1:3025232/100‑1 (MQ=255)
                                      tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                            <  2:1952107/100‑1 (MQ=255)
                                         ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                                         >  2:3596102/1‑100 (MQ=255)
                                           cgTTTACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAACATATTGC‑‑‑‑CA‑CT‑CACCCAGTAGGCCTGATAAGCGCAGCGCATCAGGCAattttacattt                                                        >  2:2714125/1‑94 (MQ=25)
                                           cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                       <  2:2958997/100‑1 (MQ=255)
                                           cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                       <  2:2084588/100‑1 (MQ=255)
                                           cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                       >  2:3673206/1‑100 (MQ=255)
                                            gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                      <  1:1401484/100‑1 (MQ=255)
                                             ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcn                                                     <  1:281638/100‑2 (MQ=255)
                                             ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                     <  1:833241/100‑1 (MQ=255)
                                             ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                     <  2:2447027/100‑1 (MQ=255)
                                              tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                    <  1:874478/100‑1 (MQ=255)
                                               aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                                   >  1:1838723/1‑100 (MQ=255)
                                               aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                                   >  1:3092567/1‑100 (MQ=255)
                                                aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                  >  2:1623522/1‑100 (MQ=255)
                                                aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                  <  1:1818766/100‑1 (MQ=255)
                                                aCGCCTTATCCGGACTAC‑AAAATCGTGAAAATTAAATACATTGC‑‑A‑GC‑AT‑TTTA‑‑GTAGGAGTGATAAGAGGAGCGCATCAGGCATTtttgcctttggtatc                                                  >  1:2364465/1‑95 (MQ=18)
                                                 ccccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCCTCa                                                 >  2:2576738/3‑100 (MQ=255)
                                                 cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                                 >  1:2086127/1‑100 (MQ=255)
                                                   ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGG‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGACTTTGCCATCAGt                                               >  1:949001/1‑100 (MQ=255)
                                                   ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                               >  2:2224573/1‑100 (MQ=255)
                                                     ttATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                             <  2:597113/100‑1 (MQ=255)
                                                       aTCCGGCCTAC‑AGAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                           <  1:2446146/100‑1 (MQ=255)
                                                       aTCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAtacaagact                                           <  2:469193/100‑9 (MQ=255)
                                                       aTCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                           <  2:537967/100‑1 (MQ=255)
                                                         ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGTGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                         <  2:2692026/100‑1 (MQ=255)
                                                         ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGTCTTCAGTCTTGTa                                         >  1:2703393/1‑100 (MQ=255)
                                                         ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                         >  2:926617/1‑100 (MQ=255)
                                                           ggCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                       <  1:1672870/100‑1 (MQ=255)
                                                               tAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                   >  1:3328773/1‑100 (MQ=255)
                                                                   aaaaTCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                                <  1:853209/100‑1 (MQ=255)
                                                                      aTCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatt                             <  2:721508/100‑1 (MQ=255)
                                                                          tGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTATTAATTATTc                         >  1:1967836/1‑100 (MQ=255)
                                                                           gCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAATTTTGCCTTTTCCATCAGGCTTGTATGGCTTTTAATTATTCg                        >  1:2670376/1‑100 (MQ=255)
                                                                            cAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa                       >  1:1935694/1‑100 (MQ=255)
                                                                               aTTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag                    <  1:458449/100‑1 (MQ=255)
                                                                               aTTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag                    <  2:3265389/100‑1 (MQ=255)
                                                                               aTTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag                    <  2:3628064/100‑1 (MQ=255)
                                                                                   aaTACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGAc                >  2:2521314/1‑100 (MQ=255)
                                                                                       cATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGttt            <  2:3011312/100‑1 (MQ=255)
                                                                                           gC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa        >  2:1839928/1‑100 (MQ=255)
                                                                                           gC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa        >  2:2886699/1‑100 (MQ=255)
                                                                                            c‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTTTTTTTTTTTTCAGTCTTGTATGGCTTTTAATTATTCGCTTGAGActgttt           >  1:276568/1‑91 (MQ=21)
                                                                                                 gA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTCTGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg      >  2:2673547/1‑99 (MQ=255)
                                                                                                 gA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt     <  1:19144/100‑1 (MQ=255)
                                                                                                 gA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt     <  1:197850/100‑1 (MQ=255)
                                                                                                 gA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt     >  2:801042/1‑100 (MQ=255)
                                                                                                  a‑aT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt     >  1:529011/1‑99 (MQ=255)
                                                                                                    aT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGa   <  2:3119193/100‑1 (MQ=255)
                                                                                                     t‑tttC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGAt  <  2:2737722/100‑1 (MQ=255)
                                                                                                       |                                                                                                    
TACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGAT  >  NC_000913/3705931‑3706126

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: ‑