Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,048 G→T 23.4% 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,0480GT23.4% 670.4 / 177.8 286intergenic (‑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 (82/133);  new base T (55/10);  total (143/143)
Fisher's exact test for biased strand distribution p-value = 1.76e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

Only 100 of 289 total aligned reads displayed.

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705949‑3706146
                                                                                                       |                                                                                                  
tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCgg                                                                                                    >  2:75980/1‑100 (MQ=255)
     attttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAATCGTGCTAATTCAATATATTGC‑A‑GAAA‑CCTTGTAGGCCTGATAAGCGTAgcgc                                                                                               >  1:1406379/10‑100 (MQ=17)
     aatttgtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGC‑A‑GAAA‑TCATGTAGGCCTGATAAGCGTAgcgc                                                                                               >  1:3320330/10‑100 (MQ=18)
        ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑A‑TGAG‑CCACGTAGGCCTGATAAGCGTAGCGCa                                                                                              >  1:1294540/7‑98 (MQ=17)
         tgtaggCCGGATAAGGCGTTCACTCCGCATCCGGCCTAC‑AAAAGCACGCAAATTCAATATATTGC‑A‑GAGA‑TTGCGTAGGCCTGATAAGCGTAGCGCATc                                                                                            >  1:3992390/6‑99 (MQ=11)
          gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAATCGTGCTAATTCAATATATTGC‑A‑GAAA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCa                                                                                           >  2:4062215/5‑99 (MQ=17)
          gtaggCCGGATAAGGCGCTTACGCCGCATCCGGCCTAC‑AAGAACGTGCAAATTCAACAAATTGTGATGCAT‑CCT‑GTAGGCCTGATAAGCGTAGCGCATCa                                                                                           >  1:1023285/5‑100 (MQ=18)
           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAgg                                                                                         >  1:2918897/1‑100 (MQ=255)
           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTCC‑A‑GAAT‑TTTCGTAGACCTGATAAGAGGAGCGCATCAgg                                                                                         <  2:2567147/100‑1 (MQ=255)
            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑A‑GAAA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                        >  1:3633099/3‑100 (MQ=17)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGC‑A‑‑ATA‑TGACGTAGGCCTGATAAGCGTAGCGCATCAGGcg                                                                                       >  1:638026/3‑99 (MQ=21)
            aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAgg                                                                                         >  2:3368986/1‑99 (MQ=255)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑G‑GAAA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  1:2546909/2‑100 (MQ=17)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGT‑G‑GGCT‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:1516650/2‑100 (MQ=21)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑A‑GGAG‑CTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  1:3830939/2‑100 (MQ=18)
              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                      >  1:1754180/1‑100 (MQ=255)
              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                      >  2:1639382/1‑100 (MQ=255)
              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑GAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                                                     >  1:743393/1‑100 (MQ=18)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGA‑T‑CGAT‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                      >  1:3517283/1‑99 (MQ=21)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGC‑A‑GGAA‑ACACGTAGGCCTGATAAGCGAAGCGCATCAGGCag                                                                                      >  1:2784636/1‑99 (MQ=18)
               ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑A‑GGAC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                                     >  2:2910422/9‑98 (MQ=18)
                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑GAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                                                  >  2:228456/7‑98 (MQ=18)
                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑A‑GAAA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  1:303077/7‑100 (MQ=18)
                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑A‑GGAG‑CTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  1:535835/7‑100 (MQ=17)
                 ggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGC‑AGGAAC‑GAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                                                                                   >  1:3390970/7‑96 (MQ=25)
                 ggataaggcgTTCACGCCGCATCCAGCCTAC‑AAAATTGTGCAAAATCAATGGATTGC‑A‑CAGC‑TAGCGTAGGCCTGATAAGCGTAGCGCATCAGGCagatt                                                                                   >  2:2796848/7‑96 (MQ=25)
                  gataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑A‑GGAC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                                                                                  >  2:2828621/6‑100 (MQ=18)
                   aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  2:2382480/100‑1 (MQ=255)
                   aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  1:3438729/100‑1 (MQ=255)
                     gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                                               <  2:3608705/100‑1 (MQ=255)
                      cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                              >  2:3730951/1‑100 (MQ=255)
                       ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                                                             <  2:11452/100‑1 (MQ=255)
                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                           <  2:320447/100‑1 (MQ=255)
                         cgTTAACGCCTTATCCGACCAAC‑‑AAATTATGCAACTCCACAATATTGC‑A‑GAAA‑TTAAGTAGGACTTATAAGCATAACGCATAAAACAtttattcgttt                                                                            >  2:329578/1‑90 (MQ=255)
                          gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                                          <  1:2669803/100‑1 (MQ=255)
                          gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                           <  2:474306/99‑1 (MQ=255)
                           ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                                         <  1:899379/100‑1 (MQ=255)
                            tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                                        >  1:4091024/1‑100 (MQ=255)
                              aCGCCTTATCCGGCTTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                      <  1:3690947/100‑1 (MQ=255)
                              aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                      <  1:384417/100‑1 (MQ=255)
                                gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                                                    >  1:259524/1‑100 (MQ=255)
                                 ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                                                   >  1:3105986/1‑100 (MQ=255)
                                    tATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt                                                                <  1:1890908/100‑1 (MQ=255)
                                    tATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCACCAGTCtt                                                                >  1:3390128/1‑100 (MQ=255)
                                      tCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGt                                                              >  1:3224180/1‑100 (MQ=255)
                                       ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                             <  1:1582880/100‑1 (MQ=255)
                                       ccGGCATAC‑AAAATCGTGCTAATTCGATACATTGC‑A‑GACT‑TTTCGTAGGCTTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCGATCAGTCGTGTa                                                             <  2:2425323/100‑1 (MQ=18)
                                          gCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATgg                                                          >  2:2502602/1‑100 (MQ=255)
                                          gCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATgg                                                          <  1:4043368/100‑1 (MQ=255)
                                           ccTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                         <  2:3057963/100‑1 (MQ=255)
                                             tAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                                       <  2:1963119/100‑1 (MQ=255)
                                             tAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                                       <  1:2730040/100‑1 (MQ=255)
                                             tAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                                       <  1:606860/100‑1 (MQ=255)
                                              aC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                                      >  1:222163/1‑100 (MQ=255)
                                                 aaaaTCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                                                    <  1:3150127/100‑1 (MQ=255)
                                                   aaTCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                                                  <  1:3430994/100‑1 (MQ=255)
                                                   aaTCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                                                  <  2:3522212/100‑1 (MQ=255)
                                                   aaTCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCAACAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                                                  >  2:2663186/1‑100 (MQ=255)
                                                    aTCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                                                  <  2:3699686/99‑1 (MQ=255)
                                                        tGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             >  1:2980534/1‑100 (MQ=255)
                                                        tGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  1:1676166/100‑1 (MQ=255)
                                                         gCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                            <  2:1735904/100‑1 (MQ=255)
                                                          cAAATTCAATACATTGC‑A‑GAAA‑TATGGTAGGCGTGATAAGATAAGCGCATTAGGCATTTTTACATTAATAATAAGACTAGTATGGCTCGGAATTATTc                                             >  2:942887/1‑98 (MQ=9)
                                                           aaaTTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt                                          <  2:3390235/100‑1 (MQ=255)
                                                            aaTTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa                                         >  2:1261108/1‑100 (MQ=255)
                                                                 aaTACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGAc                                    >  1:3564456/1‑100 (MQ=255)
                                                                  aTACATTGC‑A‑GAAT‑TTTCGTAGGACTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCGTGTCTGGCTTTTAATTATTCGATAGAGACg                                   >  2:3988596/1‑100 (MQ=255)
                                                                   tACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGt                                  >  1:3498028/1‑100 (MQ=255)
                                                                     cATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGttt                                >  2:911901/1‑100 (MQ=255)
                                                                     cATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGttt                                >  1:11339/1‑100 (MQ=255)
                                                                            a‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                          <  2:1079983/100‑1 (MQ=255)
                                                                            a‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                          <  2:1370910/100‑1 (MQ=255)
                                                                              gAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                         >  1:1844958/1‑100 (MQ=255)
                                                                                 t‑tttCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGAt                      <  1:1160293/100‑1 (MQ=255)
                                                                                 t‑tttCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGCTTTCGCAGTGAt                      >  1:3916876/1‑100 (MQ=255)
                                                                                      cGTAGGCCTGATAAGCGTAGCGCATCAGGCctacgaaaatcttgcaatgtattgaattttcatgattttgtaggccggataaggcgttcacgccgcatcc                  <  1:1534681‑M2/100‑71 (MQ=255)
                                                                                       gTAGGCCTGATAAGCGTAGCGCATCAGGCctacgaaaatcttgcaatgtattgaattttcatgattttgtaggccggataaggcgttcacgccgcatccg                 <  1:1483842‑M2/100‑72 (MQ=255)
                                                                                        tAGGCCTGATAAGCGCAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTATATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGc                >  2:3938122/1‑100 (MQ=255)
                                                                                         aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTTTTTTTTTTTTCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGc                >  2:315705/1‑99 (MQ=37)
                                                                                         aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  1:2166590/100‑1 (MQ=255)
                                                                                         aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  2:1934575/100‑1 (MQ=255)
                                                                                         aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  1:2577527/100‑1 (MQ=255)
                                                                                          ggCCTGATAAGCGTAGCGCATCAGGCctacgaaaatcttgcaatgtattgaattttcatgattttgtaggccggataaggcgttcacgccgcatccggca              <  1:211633‑M2/100‑75 (MQ=255)
                                                                                          ggCCTGATAAGCGGAGCGCATCAGGCctacgaaaattctgcaatgtattgaatttgcacgattttgtaggccggataaggcgttaacgccgcatccggca              <  1:1447241‑M2/100‑75 (MQ=255)
                                                                                           gCCTGATAAGCGTAGCGCATCAGGCctacgaaaatcttgcaatgtattgaattttcatgattttgtaggccggataaggcgttcacgccgcatccggcaa             <  1:3874311‑M2/100‑76 (MQ=255)
                                                                                             cTGATAAGCGGAGCGCATCAGGAATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAAt           <  2:4142988/100‑1 (MQ=255)
                                                                                              tGATAAGCGTAGCGCATCAGGCctacgaaaatcttgcaatgtattgaattttcatgattttgtaggccggataaggcgttcacgccgcatccggcaagaa          <  2:2101865‑M2/100‑79 (MQ=255)
                                                                                              tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg          <  1:1089222/100‑1 (MQ=255)
                                                                                              tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg          <  2:1046475/100‑1 (MQ=255)
                                                                                              tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg          <  1:1331146/100‑1 (MQ=255)
                                                                                              tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg          <  1:2378466/100‑1 (MQ=255)
                                                                                               gATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATgg         >  2:205706/1‑100 (MQ=255)
                                                                                               gATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATgg         <  1:2293630/100‑1 (MQ=255)
                                                                                                  aaGCGGAGCGCATCAGGCATTTTTGCCTTTACCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATc      <  1:1206315/100‑1 (MQ=255)
                                                                                                   aGCGGAGCGCATCAGGCATTTTTTTTTTTTTTATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa     >  2:1920843/1‑100 (MQ=255)
                                                                                                   aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa     <  1:2478949/100‑1 (MQ=255)
                                                                                                    gCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCaa    <  2:871727/100‑1 (MQ=255)
                                                                                                     cGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAc   >  1:3777612/1‑100 (MQ=255)
                                                                                                      ggAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc  <  2:3848214/100‑1 (MQ=255)
                                                                                                      ggAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc  <  1:2004929/100‑1 (MQ=255)
                                                                                                       |                                                                                                  
TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705949‑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: ‑