Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,183,345:1 +T 18.3% 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.T18.3% 566.7 / 147.2 186intergenic (+22/+468)zupT/ribBzinc transporter/3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base . (84/68);  new base T (19/15);  total (103/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.91e-01

Only 100 of 198 total aligned reads displayed.

CCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGC  >  NC_000913/3183246‑3183444
                                                                                                    |                                                                                                   
 ccAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCat                                                                                                     >  2:1014867/1‑99 (MQ=255)
    gCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAgcgc                                                                                                 <  1:2272912/100‑1 (MQ=255)
       aTGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATg                                                                                              >  1:2239348/1‑100 (MQ=255)
       aTGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATg                                                                                              >  2:404960/1‑100 (MQ=255)
       aTGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATg                                                                                              <  2:2656600/100‑1 (MQ=255)
        tGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATgg                                                                                             <  2:562147/100‑1 (MQ=255)
        tGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATgg                                                                                             >  1:1274477/1‑100 (MQ=255)
        tGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATgg                                                                                             <  2:3613064/100‑1 (MQ=255)
        tGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATgg                                                                                             >  2:3600138/1‑100 (MQ=255)
         ggCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGa                                                                                            >  1:3072613/1‑100 (MQ=255)
         ggCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGa                                                                                            >  1:3532635/1‑100 (MQ=255)
          gCGTACTATGTGGAATGTCAGTGATGGGATTCATTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAATATCAGAAGTTGTCCAGCGCATGGAc                                                                                           <  2:828441/100‑1 (MQ=255)
          gCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGAc                                                                                           >  1:644349/1‑100 (MQ=255)
           cGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACa                                                                                          >  1:910186/1‑100 (MQ=255)
           cGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACa                                                                                          >  2:2965099/1‑100 (MQ=255)
              aCTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACt                                                                                       <  1:3101625/100‑1 (MQ=39)
              aCTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACt                                                                                       <  1:2763571/100‑1 (MQ=255)
              aCTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACt                                                                                       <  2:3283056/100‑1 (MQ=255)
                tATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACt                                                                                       <  2:3176020/98‑1 (MQ=39)
                  tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTg                                                                                   <  2:881683/100‑1 (MQ=255)
                   gtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCTga                                                                                  <  1:1613838/100‑2 (MQ=38)
                    tgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGTCGGGAATTGGTTAAAAACATCATCAGAAGTTGTCCATGCGCTGGACAACTTCTgat                                                                                 <  2:2079067/100‑3 (MQ=37)
                    tgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGcc                                                                                 <  1:1585026/100‑1 (MQ=255)
                    tgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGcc                                                                                 <  1:404962/100‑1 (MQ=255)
                     ggAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCa                                                                                >  2:2267968/1‑100 (MQ=255)
                     ggAATAGCAGTTATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACACCATCAGAAGTTGTCCCGCGCATGGACAACTTCt                                                                                    <  2:3342072/96‑1 (MQ=255)
                       aaTGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGc                                                                              <  2:1271763/100‑1 (MQ=255)
                       aaTGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGc                                                                              <  2:1969136/100‑1 (MQ=255)
                          gTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTc                                                                           >  1:2519713/1‑100 (MQ=255)
                          gTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTc                                                                           >  1:2606063/1‑100 (MQ=255)
                          gTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTc                                                                           >  1:958164/1‑100 (MQ=255)
                             aGTGATGGGATTCAGATGAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGCAGTTGTCCAGCGTGTGGACAACTTTTGCCAGCTTc                                                                           <  2:1695930/97‑1 (MQ=255)
                              gTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGa                                                                       >  1:3288846/1‑100 (MQ=255)
                               tGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGaa                                                                      >  1:3410304/1‑100 (MQ=255)
                                gATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGaaa                                                                     >  2:1968138/1‑100 (MQ=255)
                                gATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGaaa                                                                     >  1:2097215/1‑100 (MQ=255)
                                  tGGGATGCAGTCTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAGCTTCTGCCAGTTTCACGAAAg                                                                    <  2:1777415/99‑1 (MQ=255)
                                   gggATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGcc                                                                  <  1:1039115/100‑1 (MQ=255)
                                   gggATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGcc                                                                  <  1:1039084/100‑1 (MQ=255)
                                    ggATTCAGTTTAGTGCTGCTAGAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCt                                                                 <  2:2898977/100‑1 (MQ=255)
                                    ggATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCt                                                                 >  1:3335668/1‑100 (MQ=255)
                                     gATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCtt                                                                >  1:1689242/1‑100 (MQ=255)
                                         cAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAACCAACACCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTc                                                            <  2:1923943/100‑1 (MQ=255)
                                          aGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCt                                                           >  2:1198133/1‑100 (MQ=255)
                                          aGTTTAGTGCTGCCACAAACGGCGGGAATTGGTTGAAGACGACATCAGAAGTTGTCTAGCGCATGGACAACTTCTGCCAGGTTCACGGATGCCTTATTCt                                                           <  2:155409/100‑1 (MQ=255)
                                           gTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTa                                                          >  1:2682789/1‑100 (MQ=255)
                                              tAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCt                                                       >  1:3060098/1‑100 (MQ=255)
                                                  gctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTc                                                   >  1:2169676/1‑100 (MQ=255)
                                                   ctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCTCGAAAGCCTTATTCTACCTTTTCg                                                  <  2:2600961/100‑1 (MQ=255)
                                                   ctgctACAAACGGCGAGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTAACTAAAGCCTTATTCTACATTTTc                                                   >  2:2133575/1‑99 (MQ=255)
                                                    tgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCgg                                                 >  1:3618574/1‑100 (MQ=255)
                                                       tACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTg                                              <  1:1721565/100‑1 (MQ=255)
                                                       tACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTg                                              <  1:656206/100‑1 (MQ=255)
                                                       tACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTg                                              <  1:861424/100‑1 (MQ=255)
                                                        aCAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGt                                             <  1:3371187/100‑1 (MQ=255)
                                                         cAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtc                                            >  1:3309994/1‑100 (MQ=255)
                                                           aaCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctc                                          <  1:1302089/100‑1 (MQ=255)
                                                               gcggGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttct                                      >  2:2934154/1‑100 (MQ=255)
                                                               gcggGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttct                                      >  2:3435964/1‑100 (MQ=255)
                                                                 gggAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttctct                                    <  1:1828589/100‑1 (MQ=255)
                                                                  ggAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTtctctc                                   >  2:889166/1‑100 (MQ=255)
                                                                   gAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCg                                  <  2:1599731/100‑1 (MQ=255)
                                                                    aaTTGGTTAAAAACATCATCAGAAGTTAGCCAGCGCATTGACTTCTTCTGCCAGCTTCACTAAAGCCTTATCCTACCTTTTCGGCTGTCTCTTCTCTCGt                                 <  2:1111452/100‑1 (MQ=255)
                                                                    aaTTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGt                                 <  2:541872/100‑1 (MQ=255)
                                                                     aTTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTa                                >  1:2969435/1‑100 (MQ=255)
                                                                           tAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGttt                          >  1:1198674/1‑100 (MQ=255)
                                                                           tAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAAATTCTTCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTGTCTTCTCTCGTACTGt                            >  2:3091143/1‑98 (MQ=255)
                                                                              aaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAg                       <  1:2783113/100‑1 (MQ=255)
                                                                              aaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAg                       <  1:1545220/100‑1 (MQ=255)
                                                                                acaacaTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTa                     >  1:1449993/1‑100 (MQ=255)
                                                                                acaacaTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTa                     <  1:2699465/100‑1 (MQ=255)
                                                                                 caacaTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAt                    >  1:2379300/1‑100 (MQ=255)
                                                                                  aacaTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAtt                   <  1:2063998/100‑1 (MQ=255)
                                                                                   ctggcAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAttt                  >  1:467120/5‑100 (MQ=37)
                                                                                    caTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTg                 <  2:3546988/100‑1 (MQ=255)
                                                                                    caTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTg                 <  2:2760536/100‑1 (MQ=255)
                                                                                     aTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt                >  1:2176923/1‑100 (MQ=255)
                                                                                     aTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTATACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt                >  1:2722643/1‑100 (MQ=255)
                                                                                      gcAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGtt               >  1:187131/2‑100 (MQ=38)
                                                                                       cAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTc              >  1:1486358/1‑100 (MQ=38)
                                                                                       cAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTc              >  1:1126134/1‑100 (MQ=38)
                                                                                        aGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTcc             <  1:245824/100‑1 (MQ=255)
                                                                                         gAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  1:611480/1‑100 (MQ=39)
                                                                                         gAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  1:1779150/1‑100 (MQ=39)
                                                                                         gAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  2:786423/1‑100 (MQ=39)
                                                                                         gAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  2:3217666/1‑100 (MQ=39)
                                                                                         gAAGTTGTACATGCGCTGGACAAATTCTGCCAGCTTCACGAAAGCCTTATTCGACCTTTTCGGCTGTCTCTTCTCACGTACTGTTTAAGTATTTGTTcc             >  2:692742/1‑99 (MQ=255)
                                                                                          aaGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgc           >  1:160224/1‑100 (MQ=39)
                                                                                           aGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgcg          >  2:3281496/1‑100 (MQ=39)
                                                                                           aGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgcg          >  1:742009/1‑100 (MQ=39)
                                                                                           aGTTGTCCAGCGCATGGACAACTTCTGCCAGCTGCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgc           <  2:1484750/99‑1 (MQ=255)
                                                                                            gTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTCTTCTACCTTTTATGCTGTGTCTTCTCTCGTCCTGTTTAAGTATTTGTTc              >  2:1804659/1‑95 (MQ=37)
                                                                                              tGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCtt       >  1:55846/1‑100 (MQ=255)
                                                                                              tGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCtt       <  2:613705/100‑1 (MQ=255)
                                                                                               gTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTc      <  1:2988532/100‑1 (MQ=255)
                                                                                               gTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTc      <  2:2820248/100‑1 (MQ=255)
                                                                                                   aGCGCATGGACATCTTCTGCCAGCTTCACGAAAGCGTTCTTCTTCCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGc  <  2:1883172/100‑1 (MQ=255)
                                                                                                   aGCGCATGGACAACTTTTGCCAGCTTCACGAAAGCCTTATTCTATCTTTTCGGCTGTCTCTTCTCGCGTACTGTTTAAGTATTTGTTCCGCGATTCa     <  2:3337553/97‑1 (MQ=255)
                                                                                                   aGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGc  <  1:3600165/100‑1 (MQ=255)
                                                                                                   aGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGc  >  1:2088297/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
CCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGC  >  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: