Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,183,345:1 +T 17.1% intergenic (+22/+468) zupT → / ← ribB zinc transporter/3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,183,3451.T17.1% 396.6 / 91.1 123intergenic (+22/+468)zupT/ribBzinc transporter/3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base . (54/48);  new base T (10/11);  total (64/59)
Fisher's exact test for biased strand distribution p-value = 8.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.77e-01

Only 100 of 135 total aligned reads displayed.

CCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGC  >  NC_000913/3183246‑3183444
                                                                                                    |                                                                                                   
cccAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCa                                                                                                      >  2:1991756/1‑100 (MQ=255)
atctgCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCa                                                                                                      <  1:1887037/96‑1 (MQ=255)
atctgCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCa                                                                                                      <  1:1886971/96‑1 (MQ=255)
 ccAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑g                                                                                                    >  2:1294496/1‑100 (MQ=255)
 ccAGCTATGGCGTACTATGTGGAATGTCAGTGATGCGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCa                                                                                                      <  2:503613/99‑1 (MQ=255)
    gCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCaagat                                                                                                  >  2:1886963/1‑96 (MQ=255)
    gCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCaagat                                                                                                  >  2:1887029/1‑96 (MQ=255)
    gCTATGGCGTACTATGTGGAATGCCAGTGATGGGCTTCAGTTTAGTGCTGCTACAAACGGCGGGAATTCGATAAACACAACATCAGAAGTTGTCCA‑gcgc                                                                                                 >  2:1831142/1‑100 (MQ=255)
     cTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCa                                                                                                <  2:1023488/100‑1 (MQ=255)
       aTGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATg                                                                                              <  1:2233766/100‑1 (MQ=255)
         ggCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGa                                                                                            <  1:1769685/100‑1 (MQ=255)
          gCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGAc                                                                                           <  2:905734/100‑1 (MQ=255)
           cGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACa                                                                                          >  2:602000/1‑100 (MQ=255)
             tACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAAc                                                                                        <  2:956536/100‑1 (MQ=39)
             tACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAAc                                                                                        <  1:292138/100‑1 (MQ=39)
             tACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAAc                                                                                        <  2:2517011/100‑1 (MQ=255)
              aCTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACt                                                                                       >  2:607617/1‑100 (MQ=255)
                tgtgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAACTTc                                                                                     <  2:1612966/98‑1 (MQ=38)
                tATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTc                                                                                     <  1:683555/100‑1 (MQ=255)
                 aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAACTTCt                                                                                    <  1:2174823/100‑1 (MQ=38)
                 aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAACTTCt                                                                                    <  2:18771/100‑1 (MQ=38)
                 aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCt                                                                                    >  1:736294/1‑100 (MQ=255)
                  tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTg                                                                                   <  1:458165/100‑1 (MQ=255)
                  tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTg                                                                                   <  1:779520/100‑1 (MQ=255)
                  tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTg                                                                                   <  2:1151818/100‑1 (MQ=255)
                  tgtgGAATGTCAGTGATGGGATTCAGTATAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGGCCA‑GCGCATGGACAACTTCTg                                                                                   <  2:2467528/100‑1 (MQ=255)
                  tgtgGAATATCAGTGATGGGATTCAGTTTAGGGCTGCTACAATCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTg                                                                                   <  2:1905610/100‑1 (MQ=255)
                   gtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGc                                                                                  <  1:1950560/100‑1 (MQ=255)
                   gtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGc                                                                                  >  1:756962/1‑100 (MQ=255)
                   gtgGAATGTCAGTGATGGGATTCAGTTTAGGGCTGCCACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACATCTTCTg                                                                                   <  2:415037/99‑1 (MQ=255)
                    tgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGcc                                                                                 <  1:359113/100‑1 (MQ=255)
                    tgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGcc                                                                                 <  1:1402941/100‑1 (MQ=255)
                    tgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGcc                                                                                 <  2:2717158/100‑1 (MQ=255)
                     ggAATGTCAGTGGTGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGcc                                                                                 <  2:2241239/99‑1 (MQ=255)
                     ggAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCa                                                                                <  1:1018347/100‑1 (MQ=255)
                     ggAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCa                                                                                >  1:1122971/1‑100 (MQ=255)
                        aTGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCt                                                                             <  1:2530187/100‑1 (MQ=255)
                          gTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTc                                                                           >  1:958377/1‑100 (MQ=255)
                            cAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCAc                                                                         <  1:2516205/100‑1 (MQ=255)
                            cAGTGAGGGGATTCAGTTTAGTGCTGCTACAAACGGTGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAACTTCTgatgttgtttt                                                                         <  2:917247/100‑11 (MQ=21)
                             aGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCCGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACg                                                                        <  2:36356/100‑1 (MQ=255)
                               tGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCCTCACGaa                                                                      >  2:85060/1‑100 (MQ=255)
                                  tGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGc                                                                   >  1:1966805/1‑100 (MQ=255)
                                         cAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTAtt                                                             <  2:2616168/99‑1 (MQ=255)
                                         cAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTc                                                            <  1:179839/100‑1 (MQ=255)
                                           gTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTa                                                          >  1:246918/1‑100 (MQ=255)
                                             ttAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTAcc                                                        <  1:1054724/100‑1 (MQ=255)
                                                 tgctgcTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCtttt                                                    >  2:811072/1‑100 (MQ=255)
                                                  gctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCGTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTc                                                   <  2:2251251/100‑1 (MQ=255)
                                                  gctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTc                                                   >  1:1792207/1‑100 (MQ=255)
                                                   ctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCg                                                  >  1:915708/1‑100 (MQ=255)
                                                   ctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCg                                                  >  1:212472/1‑100 (MQ=255)
                                                   ctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCg                                                  <  2:2072333/100‑1 (MQ=255)
                                                    tgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCgg                                                 <  2:1327407/100‑1 (MQ=255)
                                                     gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCgga                                                >  1:1327412/1‑99 (MQ=255)
                                                     gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                <  1:2319234/100‑1 (MQ=255)
                                                     gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                >  2:349785/1‑100 (MQ=255)
                                                     gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCATCACGAAAGCCGCAGTCTACCTTTTCGGc                                                >  2:1648820/1‑100 (MQ=255)
                                                      ctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCt                                               <  1:1209323/100‑1 (MQ=255)
                                                      ctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCt                                               >  1:944770/1‑100 (MQ=255)
                                                       tACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTg                                              <  1:1101402/100‑1 (MQ=255)
                                                       tACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTg                                              <  1:747323/100‑1 (MQ=255)
                                                         cAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtc                                            >  2:1921241/1‑100 (MQ=255)
                                                         cAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtc                                            <  2:614933/100‑1 (MQ=255)
                                                         cAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCGGtc                                            >  2:80951/1‑100 (MQ=255)
                                                           aaCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctc                                          >  2:1257210/1‑100 (MQ=255)
                                                              ggcggGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttc                                       >  1:1977894/1‑100 (MQ=255)
                                                                   gAATTGGTTAACAACAACATCCGAAGTTGCCCA‑GCGCATGGACAACTTCTGACAGCTTCACAAAACCCTTCTTCTACCTTTTCGGCTGTCTCTtctctc                                   <  2:950753/99‑1 (MQ=255)
                                                                    aaTTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGCTAACTTCTGCCCGCTTCACGCAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCg                                  <  2:1919404/99‑1 (MQ=255)
                                                                      ttGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTAc                               <  1:2146674/100‑1 (MQ=255)
                                                                        ggTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTg                             >  2:122916/1‑100 (MQ=255)
                                                                         gTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGt                            >  1:1445040/1‑100 (MQ=255)
                                                                           tAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAATTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGttt                          >  2:1466575/1‑100 (MQ=255)
                                                                             aaaaCAACATCAGAAGTTGTCCA‑GCGAATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTaa                        >  1:1677587/1‑100 (MQ=255)
                                                                             aaaaCAACATCAGAAGTCGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCGTTTCGGCTGTCTCTTCTCTCGTACTGTTTa                         <  2:1596155/99‑1 (MQ=255)
                                                                                   acaTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTTCCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAttt                  <  2:117368/100‑1 (MQ=255)
                                                                                   acaTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAttt                  <  1:692070/100‑1 (MQ=255)
                                                                                   acaTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAttt                  >  2:2095322/1‑100 (MQ=255)
                                                                                    tggcAGAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTg                 >  1:1726816/4‑100 (MQ=38)
                                                                                    tggcAGAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTg                 >  2:1831761/4‑100 (MQ=38)
                                                                                    caTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTg                 >  1:1014826/1‑100 (MQ=255)
                                                                                    caTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTg                 >  1:948132/1‑100 (MQ=255)
                                                                                     ggcAGAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt                >  1:1641405/3‑100 (MQ=38)
                                                                                       cAGAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTc              >  2:794931/1‑100 (MQ=38)
                                                                                        aGAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTcc             >  1:446019/1‑100 (MQ=38)
                                                                                         gAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  2:387119/1‑100 (MQ=39)
                                                                                         gAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  1:697083/1‑100 (MQ=39)
                                                                                          aaGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgc           >  1:1000763/1‑100 (MQ=255)
                                                                                           aGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCTTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTGCcgcg          <  2:702137/100‑1 (MQ=255)
                                                                                           aGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgcg          <  1:2127015/100‑1 (MQ=255)
                                                                                           aGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgcg          >  1:624248/1‑100 (MQ=255)
                                                                                            gTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCCCGTACTGTTTAAGTATTTGTTCCgcgc         >  2:690480/1‑100 (MQ=255)
                                                                                               gTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTc      >  1:780886/1‑100 (MQ=255)
                                                                                                tCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCa     <  2:720056/100‑1 (MQ=255)
                                                                                                 ccA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCaa    >  1:2749670/1‑100 (MQ=255)
                                                                                                  cA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAg   >  2:1181916/1‑100 (MQ=255)
                                                                                                  cA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAg   <  1:1266499/100‑1 (MQ=255)
                                                                                                   a‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGc  >  1:956849/1‑100 (MQ=255)
                                                                                                   a‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGc  <  1:1181921/100‑1 (MQ=255)
                                                                                                   a‑GCGCATGGACAACTTCTCCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGc  <  2:1198308/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
CCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGC  >  NC_000913/3183246‑3183444

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