Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,183,345:1 +T 20.2% 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.T20.2% 333.0 / 101.8 114intergenic (+22/+468)zupT/ribBzinc transporter/3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base . (56/35);  new base T (12/11);  total (68/46)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-02

Only 100 of 119 total aligned reads displayed.

CCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCA  >  NC_000913/3183247‑3183441
                                                                                                   |                                                                                                
ccAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAg                                                                                                 <  1:1658388/100‑1 (MQ=255)
ccAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAg                                                                                                 <  1:1587861/100‑1 (MQ=255)
  aacTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAgcg                                                                                               >  1:1789904/3‑100 (MQ=255)
   gCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAgcgc                                                                                              >  2:2386399/1‑100 (MQ=255)
   gCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAgcgc                                                                                              >  2:921469/1‑100 (MQ=255)
       tGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATgg                                                                                          >  2:2603417/1‑100 (MQ=255)
       tGGCGTACTATGTAGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATgg                                                                                          <  2:1799949/100‑1 (MQ=255)
        ggCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGa                                                                                         >  1:1142816/1‑100 (MQ=255)
        ggCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGa                                                                                         >  2:424724/1‑100 (MQ=255)
            tACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAAc                                                                                     <  2:1384627/100‑1 (MQ=255)
            tACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAAc                                                                                     <  2:1333914/100‑1 (MQ=255)
             aCTATGTGGAATGTCAGTGATGTTATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACt                                                                                    <  1:323559/100‑1 (MQ=38)
              cTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACtt                                                                                   >  1:2175091/1‑100 (MQ=255)
               tATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCATAAGTTGTCCATGCGCTGGACAACTTc                                                                                  <  2:1281038/100‑1 (MQ=39)
               tATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTc                                                                                  <  1:2231070/100‑1 (MQ=39)
               tATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTc                                                                                  <  1:1127405/100‑1 (MQ=39)
               tATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACtt                                                                                   <  2:851308/99‑1 (MQ=255)
               tATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTc                                                                                  <  1:811823/100‑1 (MQ=255)
                ttgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCt                                                                                 <  2:1087141/99‑1 (MQ=38)
                aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGGAGTTGTCCAGCGCATGGACAACTTCt                                                                                 >  1:1977985/1‑100 (MQ=255)
                aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCt                                                                                 >  1:1880777/1‑100 (MQ=255)
                 tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCta                                                                                >  1:1087151/1‑99 (MQ=38)
                 tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTg                                                                                >  1:1232746/1‑100 (MQ=255)
                  gtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGc                                                                               >  1:2520786/1‑100 (MQ=255)
                  gtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGc                                                                               >  1:2159382/1‑100 (MQ=255)
                   tgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCTgat                                                                              <  2:1147597/100‑3 (MQ=38)
                   tgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCTgat                                                                              <  1:1747293/100‑3 (MQ=38)
                   tgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGcc                                                                              >  1:720590/1‑100 (MQ=255)
                    ggAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCa                                                                             <  2:1306802/100‑1 (MQ=255)
                    ggAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCa                                                                             >  1:2121249/1‑100 (MQ=255)
                     gAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCa                                                                             <  2:1819511/99‑1 (MQ=255)
                       aTGTTAGTGATTGGATTCAGTTTAGTGCTGCTACAAACGGCTGGAATTGGTTGAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCt                                                                          <  1:2183702/100‑1 (MQ=255)
                          tCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCTgatgttgttt                                                                       <  1:356108/100‑10 (MQ=25)
                           cAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCAc                                                                      >  2:597106/1‑100 (MQ=255)
                            aGTGATGGGATTCAGTTTAGTGCTGTTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACg                                                                     <  1:464012/100‑1 (MQ=255)
                             gTGATGGGGTTCAGTTTAGTGCTGATACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAATTGTCAAGCGCATGGACAACtcatacca                                                                             >  2:2544106/1‑84 (MQ=255)
                               gATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGaaa                                                                  >  1:247232/1‑100 (MQ=255)
                                   ggATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCt                                                              >  1:1547235/1‑100 (MQ=255)
                                           tttAGTGCTGCTACAATCGTCGTTAATTGGTTAAAAACAATATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTAc                                                      <  1:2344188/100‑1 (MQ=255)
                                           tttAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTAc                                                      >  1:2370895/1‑100 (MQ=255)
                                              aGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGTAGTTGTCCAACGCATGGACAACTTCTGCCAGCTACACAAAAGCCGGATTCGACCtt                                                   >  1:2386180/1‑100 (MQ=255)
                                                 gctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGATTCACGAAAGCCTTATTCTACCTTTTc                                                >  1:2265365/1‑100 (MQ=255)
                                                  ctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCg                                               <  1:1030331/100‑1 (MQ=255)
                                                  ctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCg                                               <  1:57978/100‑1 (MQ=255)
                                                  ctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCg                                               <  1:1164233/100‑1 (MQ=255)
                                                    gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                             <  2:1166988/100‑1 (MQ=255)
                                                    gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                             >  2:1683193/1‑100 (MQ=255)
                                                      tACAAACGGCGGGAATTTGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTg                                           <  2:2265337/100‑1 (MQ=255)
                                                      tACAAACGGCGGGAATTGGTTAAAAACAACATCAGGAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTg                                           <  2:1977967/100‑1 (MQ=255)
                                                      tACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTg                                           <  2:1587252/100‑1 (MQ=255)
                                                      tACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTg                                           >  1:1206690/1‑100 (MQ=255)
                                                        cAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAACCTTATTCTACCTTTTCGGCTGtc                                         >  1:1193241/1‑100 (MQ=255)
                                                          aaCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGAAAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctc                                       >  2:2342196/1‑100 (MQ=255)
                                                           actgcgtGAATTTGTTAAAAACAACATCAGAAGTTTTCCAGCGCATTGACAACTTCTGCCAGCTTCACAAAAGCCTTATTCTACCTTTTCGGCTGtctc                                       <  2:549224/96‑1 (MQ=255)
                                                             ggcgtGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttc                                    <  1:1874513/100‑1 (MQ=255)
                                                             ggcggGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCTAGCTTCACGAAAGCATTATTCTACATTTTCGCCTGtctct                                      >  1:1543133/1‑98 (MQ=255)
                                                             ggcggGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttc                                    >  2:1930420/1‑100 (MQ=255)
                                                             ggcggGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttc                                    <  1:921478/100‑1 (MQ=255)
                                                             ggcggGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttc                                    <  1:2553583/100‑1 (MQ=255)
                                                               cggGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGATAACTTCTGCCAGCTTCTCTAAAGCCTTATTCTTCTTTTTCGGCTGtctgttctc                                  <  2:1264957/100‑1 (MQ=255)
                                                               cggGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttctc                                  >  2:1118022/1‑100 (MQ=255)
                                                                gggAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttctct                                 <  1:597110/100‑1 (MQ=255)
                                                                 ggAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTtctctc                                >  2:535165/1‑100 (MQ=255)
                                                                 ggAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTtctctc                                >  1:1380435/1‑100 (MQ=255)
                                                                  gAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCg                               <  1:1809943/100‑1 (MQ=255)
                                                                     ttGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTAc                            >  2:2092495/1‑100 (MQ=255)
                                                                       ggTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTg                          >  2:957329/1‑100 (MQ=255)
                                                                       ggTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTg                          >  2:398075/1‑100 (MQ=255)
                                                                       ggTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTACTGCCCGATTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACt                           >  2:2366101/1‑99 (MQ=255)
                                                                        gTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGt                         >  2:975487/1‑100 (MQ=255)
                                                                         ttAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGtt                        >  2:600603/1‑100 (MQ=255)
                                                                            gatctgctggcAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTaa                     <  2:409550/90‑1 (MQ=25)
                                                                            gatctgctggcAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTaa                     <  2:789995/90‑1 (MQ=25)
                                                                            aaaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTaa                     >  1:2305006/1‑100 (MQ=255)
                                                                            aaaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTaa                     <  1:1850073/100‑1 (MQ=255)
                                                                             aaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAg                    <  2:4286/100‑1 (MQ=255)
                                                                              gaagctggcAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGt                   >  1:240369/9‑100 (MQ=25)
                                                                              aacaacATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTATACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGt                   >  1:152621/1‑100 (MQ=255)
                                                                                caacaTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTATACCTTTTCGGCTGTATATTCTCTCGTACTGTTTAAGTAt                 >  1:539108/1‑100 (MQ=255)
                                                                                 gctggcAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTTACGAAAGCCTTATTCTACCTTTTCGTCTGTCTCTTCTCTCGTACTGTTTaaagatc                >  1:790003/6‑95 (MQ=25)
                                                                                  ctggcAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAttt               >  2:1135713/5‑100 (MQ=37)
                                                                                   tggcAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTg              >  1:842975/4‑100 (MQ=38)
                                                                                    ggcAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt             >  1:1664019/3‑100 (MQ=38)
                                                                                    ggcAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt             >  1:2100523/3‑100 (MQ=38)
                                                                                    aTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt             >  1:1708609/1‑100 (MQ=255)
                                                                                    aTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt             >  2:438003/1‑100 (MQ=255)
                                                                                    aTCAGAAGTTGTCAAGCGCATGGACAACTTCTGCCAGCGTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt             >  2:1637349/1‑100 (MQ=255)
                                                                                     tCAGAAGTTGTCCAGTGCATGGACAACTTCTGCTAGCTTCCCGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt             <  1:1550985/99‑1 (MQ=255)
                                                                                     tCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGtt            >  2:1021775/1‑100 (MQ=255)
                                                                                      cAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTc           >  1:2092691/1‑100 (MQ=255)
                                                                                      cAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTc           >  1:2494488/1‑100 (MQ=255)
                                                                                       aGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTcc          <  2:247230/100‑1 (MQ=255)
                                                                                        gAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg         <  2:53818/100‑1 (MQ=255)
                                                                                          aGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgcg       >  1:1927412/1‑100 (MQ=39)
                                                                                          aGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgcg       >  1:1741722/1‑100 (MQ=39)
                                                                                            ttGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCt     <  1:1227433/100‑1 (MQ=255)
                                                                                             tGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCtt    >  1:166336/1‑100 (MQ=255)
                                                                                              gTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTc   >  2:240969/1‑100 (MQ=255)
                                                                                               tCCAGCGCATTGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTGCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCa  <  1:106705/100‑1 (MQ=255)
                                                                                               tCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCa  >  2:116140/1‑100 (MQ=255)
                                                                                                   |                                                                                                
CCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCA  >  NC_000913/3183247‑3183441

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