Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,048 G→T 22.8% 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,0480GT22.8% 523.8 / 128.4 215intergenic (‑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 (79/82);  new base T (47/2);  total (131/84)
Fisher's exact test for biased strand distribution p-value = 2.88e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 220 total aligned reads displayed.

GCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705950‑3706146
                                                                                                        |                                                                                                  
gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGa                                                                                                   >  1:1695818/1‑100 (MQ=255)
  gcTTTGTTTATGCCGTATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCATGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGa                                                                                                   <  1:845356/98‑1 (MQ=255)
  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                 <  2:1180077/100‑1 (MQ=255)
   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                <  2:1907433/100‑1 (MQ=255)
   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                 <  1:1219329/99‑1 (MQ=255)
      tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCAt                                                                                             <  1:3753949/100‑1 (MQ=255)
      tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCAt                                                                                             <  2:2250984/100‑1 (MQ=255)
       ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AACC‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATc                                                                                            >  1:1135662/7‑100 (MQ=18)
         ttATGCCGTATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                                                          <  2:615806/100‑1 (MQ=255)
         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑ACCA‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                          >  1:1949991/5‑100 (MQ=18)
         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AGTT‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                         >  2:3120682/5‑100 (MQ=21)
          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                                                         >  1:1550646/1‑100 (MQ=255)
           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCA‑TG‑GGCC‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                        >  1:3201060/3‑100 (MQ=17)
           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGATCA‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                        >  2:3365943/3‑100 (MQ=17)
           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGCA‑AT‑CTCC‑GT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                        >  1:722526/3‑100 (MQ=18)
           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑AAAC‑AC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGc                                                                                        >  2:1029086/3‑100 (MQ=18)
           aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                                                        >  2:2151747/3‑100 (MQ=18)
            ggcaggataaggcgTTCACGCCGCATCCTGCCTAC‑AAAAGCACGCAAATTCAATATATTGCA‑GA‑GATG‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:1258959/11‑100 (MQ=14)
            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCATGG‑GCCA‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  1:471172/2‑100 (MQ=17)
            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑ATCC‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:1784483/2‑100 (MQ=21)
             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑ATCC‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                      >  1:3347711/1‑99 (MQ=21)
             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGATCA‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                      >  2:3302709/1‑99 (MQ=17)
              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑ATCC‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                      >  2:939175/9‑98 (MQ=18)
              ccGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑ACCA‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                                     >  1:929404/1‑98 (MQ=18)
               cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑ATCT‑CCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                                                    >  2:335103/8‑97 (MQ=18)
               cGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                                                                                    >  2:2844637/1‑100 (MQ=255)
                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AGCC‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  1:2524882/7‑100 (MQ=17)
                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑ATCC‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCGGGCaattc                                                                                   >  2:635156/7‑96 (MQ=21)
                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGATCA‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                                                   >  2:3181850/7‑100 (MQ=17)
                ggataaggcgTTCACGCCGCATCCAGCCTAC‑AAAATTGTGCAAAATCAATGGATTGCA‑CA‑GCTA‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagatt                                                                                   >  1:2434254/7‑96 (MQ=25)
                ggATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                                                                                   >  1:1090170/1‑100 (MQ=255)
                 taggcaaacatcacttcgaaaacgtctttatcgaataattaaaagccatacaagactgatggcaaaggcaaaaatGCCTGATAAGCGGAGCGCATCAGGc                                                                                        >  1:1390709‑M1/76‑100 (MQ=255)
                 gataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgcg                                                                                 >  1:3655984/6‑97 (MQ=21)
                 gataaggcgTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CTTC‑TG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaagctt                                                                                  >  2:3568092/6‑95 (MQ=21)
                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  1:1032340/100‑1 (MQ=255)
                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  2:70289/100‑1 (MQ=255)
                   taaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GTTC‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgc                                                                                >  2:2685643/4‑100 (MQ=17)
                     gctcgTTAACGGCGTATGCGGCCTAC‑AAAATCTTGCACAGTCAATACATTGCA‑GA‑GTTT‑TC‑‑GTATGCCTGATAAGCGGAGCGCATCAGGCATGtttgc                                                                                <  1:3087845/95‑1 (MQ=17)
                     cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                              >  2:151445/1‑100 (MQ=255)
                       gcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt                                                                            <  2:1597117/100‑1 (MQ=255)
                       aacatcacttcgaaaacgtctctatcgaataattaaaagccatacaagactgatggcaaaggcaaaaatGCCTGATAAGCGGAGCGCATCAGGCAttttt                                                                                  >  1:179055‑M1/70‑100 (MQ=255)
                        cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                           >  2:3285306/1‑100 (MQ=255)
                        acatcacttcgaaaacgtctctatcgaataattaaaagccatacaagactgatggcaaaggcaaaaatGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 >  1:3021037‑M1/69‑100 (MQ=255)
                        acatcacttcgaaaacgtctctatcgaataattaaaagccatacaagactgatggcaaaggcaaaaatGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 >  1:2595677‑M1/69‑100 (MQ=255)
                          ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                                         <  2:2016149/100‑1 (MQ=255)
                             aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                      <  2:2523203/100‑1 (MQ=255)
                             aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                      <  2:3495213/100‑1 (MQ=255)
                              cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                      >  1:3679804/1‑99 (MQ=255)
                                ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATCT‑TC‑‑GTGGCCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                                                   >  1:164981/1‑100 (MQ=255)
                                   tATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt                                                                >  2:194285/1‑100 (MQ=255)
                                   tATCCGGACTAC‑AAAATCGTGCAAATTCAATCCATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                                                 <  1:535698/99‑1 (MQ=255)
                                        ggCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                                           >  1:396405/1‑100 (MQ=255)
                                        ggCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                                           <  1:2333133/100‑1 (MQ=255)
                                        ggCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                                           >  2:1337342/1‑100 (MQ=255)
                                          ccTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                         <  2:2366167/100‑1 (MQ=255)
                                            tAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                                       <  2:1478849/100‑1 (MQ=255)
                                             ac‑acAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCTTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                                       <  1:2758435/99‑1 (MQ=255)
                                             aC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                                      >  2:984993/1‑100 (MQ=255)
                                                 aaaTCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                                                    >  1:3335552/1‑99 (MQ=255)
                                                    tCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatta                                                <  2:2307533/100‑1 (MQ=255)
                                                       tGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  2:1067822/100‑1 (MQ=255)
                                                       tGCAAATGCAATACATTGCA‑GA‑ATTT‑TA‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTACCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  1:2252360/100‑1 (MQ=255)
                                                        gCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                            <  1:3550230/100‑1 (MQ=255)
                                                            aTTCAATACATTGTA‑GA‑ATTT‑TC‑‑GTAGCTCTGATAAGCGGATCGTATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGTCTTTTAATTATTc                                             <  1:2350371/95‑1 (MQ=37)
                                                             ttCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga                                       >  2:1874643/1‑100 (MQ=255)
                                                             ttCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga                                       <  2:882728/100‑1 (MQ=255)
                                                             ttCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga                                       <  1:1297889/100‑1 (MQ=255)
                                                                    cATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGttt                                >  1:2555635/1‑100 (MQ=255)
                                                                      ttGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc                              <  1:265547/100‑1 (MQ=255)
                                                                          a‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                          <  1:1335982/100‑1 (MQ=255)
                                                                            gA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                         >  2:881959/1‑100 (MQ=255)
                                                                            gA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                         <  1:198666/100‑1 (MQ=255)
                                                                               aTTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGa                       <  1:66391/100‑1 (MQ=255)
                                                                                    tC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGtt                   <  1:2844620/100‑1 (MQ=255)
                                                                                    tC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGtt                   >  2:3051554/1‑100 (MQ=255)
                                                                                     c‑‑GTAGGCTTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGttt                  <  1:555104/100‑1 (MQ=255)
                                                                                        gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                 <  1:3483059/100‑1 (MQ=255)
                                                                                          aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  1:1400031/100‑1 (MQ=255)
                                                                                           ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt              >  2:969394/1‑100 (MQ=255)
                                                                                           ggCCTGATAAGCGGAGCGCATCAGGCATTTATGACTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt              >  1:2501592/1‑100 (MQ=255)
                                                                                             ccTGATAAGCGGAGCGCATCAGGCATTTTTTTTTTTTTTTTCAGTCTTGTATGGCTTTTAATTATTCTTTAGAGAAGTTTTCGAAGTGATGTTTGCCt              >  1:770004/1‑98 (MQ=25)
                                                                                             CCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTaa            <  2:2873811‑M1/100‑1 (MQ=255)
                                                                                              CTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAAt           <  2:1638078‑M1/100‑1 (MQ=255)
                                                                                              CTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAAt           >  2:3732197‑M1/1‑100 (MQ=255)
                                                                                              CTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAAt           >  2:608271‑M1/1‑100 (MQ=255)
                                                                                               TGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTGGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg          <  1:1866607‑M1/100‑1 (MQ=255)
                                                                                               TGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg          >  1:3452776‑M1/1‑100 (MQ=255)
                                                                                                gATAAGCGGAGCGCATCAGGCATTTTTTTTTTTTTTATCAGTCTTGTATGGCTTTTAATTATTTTATAGAGACGTTTTCGAAGTGATGTTTGCCTAAt           >  1:1230176/1‑98 (MQ=255)
                                                                                                GATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATgg         >  2:1401882‑M1/1‑100 (MQ=255)
                                                                                                GATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATgg         >  1:2198281‑M1/1‑100 (MQ=255)
                                                                                                GATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATgg         <  2:772033‑M1/100‑1 (MQ=255)
                                                                                                GATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATgg         >  2:2816142‑M1/1‑100 (MQ=255)
                                                                                                  TAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGAt       <  2:1611590‑M1/100‑1 (MQ=255)
                                                                                                   aaGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATc      <  2:1695825/100‑1 (MQ=255)
                                                                                                    aGCGGAGCGCATCAGGCAGTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGGTAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa     >  1:749820/1‑100 (MQ=255)
                                                                                                      cGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAc   <  1:2167308/100‑1 (MQ=255)
                                                                                                      cGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAc   <  2:396407/100‑1 (MQ=255)
                                                                                                      cGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAc   >  2:3526130/1‑100 (MQ=255)
                                                                                                       ggAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc  <  2:1163944/100‑1 (MQ=255)
                                                                                                       ggAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc  <  1:2823039/100‑1 (MQ=255)
                                                                                                        |                                                                                                  
GCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705950‑3706146

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