Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,048 G→T 27.0% intergenic (‑258/+50) 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,0480GT27.0% 573.4 / 186.3 263intergenic (‑258/+50)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (90/97);  new base T (64/3);  total (163/100)
Fisher's exact test for biased strand distribution p-value = 2.11e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

Only 100 of 266 total aligned reads displayed.

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705949‑3706146
                                                                                                         |                                                                                                  
 gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGa                                                                                                   >  1:2010450/1‑100 (MQ=255)
  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAg                                                                                                  <  1:2709829/100‑1 (MQ=255)
   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                 >  2:2063265/1‑100 (MQ=255)
    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                <  1:1419714/100‑1 (MQ=255)
      ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                                                              >  2:3356334/1‑100 (MQ=255)
       tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                                                                                             >  1:1198358/1‑100 (MQ=255)
        ttgtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGCA‑GA‑AA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATc                                                                                            >  2:2307935/7‑100 (MQ=18)
        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                                                            <  2:1256428/100‑1 (MQ=255)
         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                                                                                           >  2:2341635/1‑100 (MQ=255)
           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                        >  1:2069677/4‑100 (MQ=18)
            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGCA‑AT‑TC‑TCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                         >  2:1063249/3‑99 (MQ=17)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑AA‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGc                                                                                        >  2:172828/3‑100 (MQ=18)
            aTGCCGGATGCGGCGTTAACGCCGTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                                                         <  1:1470/99‑1 (MQ=255)
             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                                       >  1:1945101/1‑100 (MQ=255)
             ggtcggataaggcgTTCACGCCGCATCCGGCCCAC‑AAAGTCTTGCAAAATCAATAAATTGCAGGA‑AC‑GAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:2767944/11‑100 (MQ=18)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  1:386876/2‑100 (MQ=14)
              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                      >  1:2885851/1‑100 (MQ=255)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GG‑CT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCTCATCAGGCag                                                                                      >  1:2497197/1‑99 (MQ=18)
               ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                                     >  2:1005942/9‑98 (MQ=18)
               ccggataaggcgTTCACGCCGCACCCGGCCTAC‑AAATCCTTGCGAATTCAATATATTGTATGG‑AA‑ATTC‑‑‑‑AGGCCTGATAAGCGTAGCGCATCAGGCttt                                                                                     >  1:131964/9‑99 (MQ=25)
                cggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCCTT‑‑CTAGGCCTGATAAGCGTAGCGCATCAGGCAtttg                                                                                   >  1:3469626/8‑99 (MQ=18)
                cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGAT‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                                                    >  1:1803314/8‑97 (MQ=17)
                 ggataaggcgTTTACGCCTTATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCATGG‑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:1032296/7‑98 (MQ=18)
                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                                                  >  1:1381035/7‑98 (MQ=18)
                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                                                   >  1:368198/7‑96 (MQ=25)
                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGAT‑CG‑AT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  2:557841/7‑100 (MQ=21)
                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattc                                                                                   >  2:2193868/7‑96 (MQ=21)
                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattc                                                                                   >  2:1211845/7‑96 (MQ=21)
                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GG‑AG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  1:2177461/7‑100 (MQ=17)
                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCC‑GG‑GA‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                                                   >  1:3378774/7‑96 (MQ=18)
                 ggataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCTGAGCGCATCAGGCAtt                                                                                     >  1:3122850/7‑98 (MQ=17)
                  gataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagttct                                                                                  >  1:2310443/6‑95 (MQ=18)
                   ataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttg                                                                                 <  2:2145583/96‑1 (MQ=17)
                   aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  2:2232917/100‑1 (MQ=255)
                   aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  1:3025232/100‑1 (MQ=255)
                   aTGCGGCGTTAACGCCTTATCCGGACTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  1:3610125/100‑1 (MQ=255)
                    tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                                                <  2:1952107/100‑1 (MQ=255)
                       ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                                                             >  2:3596102/1‑100 (MQ=255)
                         cgTTTACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAACATATTGC‑‑CA‑CT‑CACCCAGTAGGCCTGATAAGCGCAGCGCATCAGGCAattttacattt                                                                            >  2:2714125/1‑94 (MQ=25)
                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                           >  2:3673206/1‑100 (MQ=255)
                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                           <  2:2958997/100‑1 (MQ=255)
                           ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcn                                                                         <  1:281638/100‑2 (MQ=255)
                            tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                                        <  1:874478/100‑1 (MQ=255)
                            tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                                        >  2:1786395/1‑100 (MQ=255)
                             aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                                                       >  1:1838723/1‑100 (MQ=255)
                              aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                      >  2:1623522/1‑100 (MQ=255)
                              aCGCCTTATCCGGACTAC‑AAAATCGTGAAAATTAAATACATTGCA‑GC‑AT‑TTTA‑‑GTAGGAGTGATAAGAGGAGCGCATCAGGCATTtttgcctttggtatc                                                                      >  1:2364465/1‑95 (MQ=18)
                               cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGACTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                                                     >  1:3549058/1‑100 (MQ=255)
                                 ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                                                   >  1:837584/1‑100 (MQ=255)
                                   ttATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                                                 <  2:597113/100‑1 (MQ=255)
                                     aTCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAtacaagact                                                               <  2:469193/100‑9 (MQ=255)
                                       ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGTGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                             <  2:2692026/100‑1 (MQ=255)
                                       ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                             >  2:896758/1‑100 (MQ=255)
                                           ccTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                         <  1:2799469/100‑1 (MQ=255)
                                            cTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                                                        >  1:3082533/1‑100 (MQ=255)
                                             tAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                                       <  1:2642341/100‑1 (MQ=255)
                                             tAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                                       >  1:3267130/1‑100 (MQ=255)
                                                 aaaaTCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                                                    <  1:853209/100‑1 (MQ=255)
                                                  aaaTCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTaa                                                   <  1:1230150/100‑1 (MQ=255)
                                                    aTCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatt                                                 <  2:721508/100‑1 (MQ=255)
                                                      cGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                                               >  1:221091/1‑100 (MQ=255)
                                                        tGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  2:1109801/100‑1 (MQ=255)
                                                         gCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                            >  1:27021/1‑100 (MQ=255)
                                                         gCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                            <  2:2477566/100‑1 (MQ=255)
                                                            aaTTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa                                         >  2:2385970/1‑100 (MQ=255)
                                                             aTTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag                                        <  2:3628064/100‑1 (MQ=255)
                                                               tCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag                                      <  1:1971069/100‑1 (MQ=255)
                                                                cAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga                                     >  1:2422744/1‑100 (MQ=255)
                                                                 aaTACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGAc                                    >  2:2521314/1‑100 (MQ=255)
                                                                 aaTACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGAc                                    >  2:402391/1‑100 (MQ=255)
                                                                     cATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGttt                                >  2:3212731/1‑100 (MQ=255)
                                                                         gCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            <  2:30052/100‑1 (MQ=255)
                                                                         gCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            >  2:2886699/1‑100 (MQ=255)
                                                                          cA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGaa                           >  2:19145/1‑100 (MQ=255)
                                                                           a‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGTGACGTTTTCGAAg                          <  2:2588363/100‑1 (MQ=255)
                                                                             gA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTCTGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                          >  2:2673547/1‑99 (MQ=255)
                                                                             gA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                         >  2:801042/1‑100 (MQ=255)
                                                                              a‑aT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                         >  1:529011/1‑99 (MQ=255)
                                                                              a‑aT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg                        <  2:2797521/100‑1 (MQ=255)
                                                                                aT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGa                       >  2:1528450/1‑100 (MQ=255)
                                                                                 t‑tttC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGATTCGTTTTCGAAGTGAt                      >  1:1577528/1‑100 (MQ=255)
                                                                                     tC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGtt                   <  1:580638/100‑1 (MQ=255)
                                                                                         gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGTCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                 <  2:949173/100‑1 (MQ=255)
                                                                                         gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                 <  2:13740/100‑1 (MQ=255)
                                                                                           aGGCCTGATAAGCGGAGCGCATCAGGTATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  1:1905825/100‑1 (MQ=255)
                                                                                            ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt              >  2:3262512/1‑100 (MQ=255)
                                                                                            ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt              >  1:2767700/1‑100 (MQ=255)
                                                                                             gCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTa             <  1:1528441/100‑1 (MQ=255)
                                                                                              ccTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTaa            >  2:3016328/1‑100 (MQ=255)
                                                                                                   tAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGAt       >  1:3754554/1‑100 (MQ=255)
                                                                                                   tAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGAt       <  1:2759358/100‑1 (MQ=255)
                                                                                                   tAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGAt       <  1:107219/100‑1 (MQ=255)
                                                                                                   tAAGCGGAGCGCACCAGGTATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGAt       <  1:445300/100‑1 (MQ=255)
                                                                                                    aaGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATTAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATc      >  2:3277672/1‑100 (MQ=255)
                                                                                                     aGCGGAGCGCATCAGGCATTTTTTTTTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa     >  1:2588637/1‑100 (MQ=255)
                                                                                                     aGCGGAGCGCATCAGGCATTTTTTCTTTTTTTTTCCGTCTAGTATGGCTTATAATTATTCAATAGAAAAGTTTTCGTAGTGATGTTTGCCTCTTGGATCa     >  1:799969/1‑100 (MQ=255)
                                                                                                      gCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCaa    >  2:1717800/1‑100 (MQ=255)
                                                                                                      gCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCaa    >  2:1448122/1‑100 (MQ=255)
                                                                                                         |                                                                                                  
TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705949‑3706146

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