Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,183,345:1 +T 15.4% 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.T15.4% 336.9 / 67.5 104intergenic (+22/+468)zupT/ribBzinc transporter/3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base . (47/41);  new base T (10/6);  total (57/47)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

Only 100 of 116 total aligned reads displayed.

CAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGC  >  NC_000913/3183248‑3183444
                                                                                                  |                                                                                                   
 aGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑gcg                                                                                                  <  1:1246733/100‑1 (MQ=255)
  gCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑gcgc                                                                                                 <  2:591827/100‑1 (MQ=255)
  gCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑gcgc                                                                                                 >  1:2132884/1‑100 (MQ=255)
  gCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑gcgc                                                                                                 >  1:2120914/1‑100 (MQ=255)
  gCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑gcgc                                                                                                 >  2:1369292/1‑100 (MQ=255)
      tGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCatgcgctgg                                                                                              <  1:408013/100‑9 (MQ=255)
      tGGCGTACTATGTGGAATGGCAGTGATGGTATTCAGTCTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCatgcgctgg                                                                                              <  2:1504859/100‑9 (MQ=255)
        gCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGAc                                                                                           >  1:692383/1‑100 (MQ=255)
          tttcTATCTGCAATGTCAGTGATGGGATTCAGTCTAGTTCGGCTACAAACGGCGGGTATTGGTTAAAATCATCATCAGAAGTTGTCCA‑GCGCATGGACa                                                                                          <  2:2257156/96‑1 (MQ=255)
            aCTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACt                                                                                       <  1:488503/100‑1 (MQ=255)
             cTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAACtt                                                                                      <  1:998178/100‑1 (MQ=39)
              tATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACGACATCAGAAGTTGTCCA‑GCGCATGGACAACTTc                                                                                     <  1:1084208/100‑1 (MQ=255)
              tATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAACTTc                                                                                     <  2:7185/100‑1 (MQ=39)
              tATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAACTTc                                                                                     <  1:1522384/100‑1 (MQ=39)
               aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCt                                                                                    <  1:2275079/100‑1 (MQ=255)
                tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAACTTCt                                                                                    <  1:8800/99‑1 (MQ=38)
                tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTg                                                                                   <  1:660781/100‑1 (MQ=255)
                tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTg                                                                                   >  2:1647304/1‑100 (MQ=255)
                 gtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGc                                                                                  >  2:1888195/1‑100 (MQ=255)
                   ggAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAACTTCTgatg                                                                                <  1:1844073/100‑4 (MQ=38)
                   ggAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCa                                                                                <  1:1120912/100‑1 (MQ=255)
                     aaTGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGc                                                                              <  2:2199192/100‑1 (MQ=255)
                      aTGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGCAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCt                                                                             <  2:2264901/100‑1 (MQ=255)
                        gTCAGTTATGTTATTCAGTTTAGTGCTGCTACAAACGGCGCGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTc                                                                           <  2:735443/100‑1 (MQ=255)
                         tCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCa                                                                          >  1:291260/1‑100 (MQ=255)
                          cAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAACTTCTgatgttgtttt                                                                         <  1:1289188/100‑11 (MQ=25)
                          cAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCAc                                                                         >  2:1542018/1‑100 (MQ=255)
                          cAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCAc                                                                         >  2:977054/1‑100 (MQ=255)
                          cAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAAAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCAc                                                                         <  1:1887924/100‑1 (MQ=255)
                            gTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGa                                                                       >  1:2258673/1‑100 (MQ=255)
                             tGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGaa                                                                      >  2:1843856/1‑100 (MQ=255)
                                tGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGc                                                                   <  2:626101/100‑1 (MQ=255)
                                tGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGc                                                                   >  2:1345108/1‑100 (MQ=255)
                                   gATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCtt                                                                >  2:1969136/1‑100 (MQ=255)
                                   gATTCAGTGTAGTGCTACTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCCTGGACAACTTCTGCCAGCTTCACGAAAGCCtt                                                                >  2:764142/1‑100 (MQ=255)
                                         gTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTa                                                          >  1:2305392/1‑100 (MQ=255)
                                         gTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTa                                                          >  1:2357614/1‑100 (MQ=255)
                                         gTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCATTATTCTa                                                          >  2:1064391/1‑100 (MQ=255)
                                          tttAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTAc                                                         <  2:1361381/100‑1 (MQ=255)
                                            tAGTGCCGCTACAACCGAAGGAAATAGATTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCt                                                       <  2:1215671/100‑1 (MQ=255)
                                              gTGCTGCTACAACCGGCGGGAATTGGGTAAAAACAACATCAGAACTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCATGAACGGCTTAGTCTACCttt                                                     >  2:1514514/1‑100 (MQ=255)
                                              gTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCttt                                                     >  1:517919/1‑100 (MQ=255)
                                               tgctgcTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCtttt                                                    <  1:1261631/100‑1 (MQ=255)
                                                gctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAAAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTc                                                   >  2:1357915/1‑100 (MQ=255)
                                                 ctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCg                                                  >  1:1868609/1‑100 (MQ=255)
                                                 ctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCg                                                  <  2:796237/100‑1 (MQ=255)
                                                   gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                <  1:395112/100‑1 (MQ=255)
                                                   gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                <  2:1713710/100‑1 (MQ=255)
                                                   gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                >  1:956494/1‑100 (MQ=255)
                                                    ctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACCACTTCTGCCAGCTTCCCGAAAGCCTTATTCTACCTTTTCGGCt                                               <  1:2011357/100‑1 (MQ=255)
                                                     tACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTg                                              >  1:506267/1‑100 (MQ=255)
                                                     tACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTg                                              >  2:1948491/1‑100 (MQ=255)
                                                       cAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtc                                            >  2:308221/1‑100 (MQ=255)
                                                       cAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtc                                            >  2:9098/1‑100 (MQ=255)
                                                        aaaCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtct                                           <  2:378198/100‑1 (MQ=255)
                                                         aaCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctc                                          >  1:1510521/1‑100 (MQ=255)
                                                            ggcggGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttc                                       >  2:1537974/1‑100 (MQ=255)
                                                            ggcggGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttc                                       >  2:1935654/1‑100 (MQ=255)
                                                            ggcggGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttc                                       <  1:1283217/100‑1 (MQ=255)
                                                            ggcggGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttc                                       <  2:1766122/100‑1 (MQ=255)
                                                             gcggGAATTGGTTAAAAACAACATCAGATGTTTTCCA‑GCGCATGGACAACTTCTGCCATCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttc                                       <  2:96689/99‑1 (MQ=255)
                                                             gcggGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttcn                                      <  2:1331976/100‑2 (MQ=255)
                                                              cggGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttctc                                     <  1:328648/100‑1 (MQ=255)
                                                                ggAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTtctctc                                   >  1:850950/1‑100 (MQ=255)
                                                                 gAATTGGTTAAAAACAACATCAGAAGTTGTTCA‑GAGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACGTTTTCGGCTGTCTCTTCTCTCg                                  >  2:1081219/1‑100 (MQ=255)
                                                                  aaTTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGt                                 <  1:759184/100‑1 (MQ=255)
                                                                   aTTGGTTAAAACCAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCCCGAAAGCCTTATTCTACCTTGTCGGCTGTCTCTTCTCTCg                                  <  2:1008474/98‑1 (MQ=255)
                                                                    ttGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTAc                               <  1:992912/100‑1 (MQ=255)
                                                                     tGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACt                              >  1:1886513/1‑100 (MQ=255)
                                                                      ggTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTg                             >  1:1122234/1‑100 (MQ=255)
                                                                      ggTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTg                             >  1:614328/1‑100 (MQ=255)
                                                                       gTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGt                            >  1:547216/1‑100 (MQ=255)
                                                                          aaaaaCAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTa                         <  1:37932/100‑1 (MQ=255)
                                                                           aaaaCAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTaa                        >  1:1644642/1‑100 (MQ=255)
                                                                            aaaCAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAg                       <  1:1827507/100‑1 (MQ=255)
                                                                             aacaacATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCTGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAg                       <  2:1905811/99‑1 (MQ=255)
                                                                             aacaacATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGt                      >  2:1302072/1‑100 (MQ=255)
                                                                              acaacaTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTa                     <  1:9098/100‑1 (MQ=255)
                                                                                aacaTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAtt                   <  2:2002801/100‑1 (MQ=255)
                                                                                aacaTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAtt                   >  1:224870/1‑100 (MQ=255)
                                                                                   aTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt                >  1:587773/1‑100 (MQ=255)
                                                                                   aTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt                >  1:2333386/1‑100 (MQ=255)
                                                                                    tCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGtt               >  1:2110572/1‑100 (MQ=255)
                                                                                    tCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGtt               >  2:681232/1‑100 (MQ=255)
                                                                                     cAGAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGtt               >  2:1000951/1‑99 (MQ=38)
                                                                                     cAGAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTc              >  1:1106422/1‑100 (MQ=38)
                                                                                     cAGAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTc              >  2:1159447/1‑100 (MQ=38)
                                                                                     cAGAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTc              >  1:829477/1‑100 (MQ=38)
                                                                                     cAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTc              >  2:1078203/1‑100 (MQ=255)
                                                                                       gAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  1:2297708/1‑100 (MQ=39)
                                                                                       gAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  1:2244696/1‑100 (MQ=39)
                                                                                       gAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  1:1336640/1‑100 (MQ=39)
                                                                                        aaGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  1:1208109/1‑99 (MQ=39)
                                                                                          gTTGTCCA‑GCTCATGGACAGCTTCTGCCAGCTTCACGAATGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCgcgc         <  2:1080459/100‑1 (MQ=255)
                                                                                          gTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCgcgc         >  2:906767/1‑100 (MQ=255)
                                                                                             gTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTc      >  1:2282804/1‑100 (MQ=255)
                                                                                             gTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTc      <  1:1302079/100‑1 (MQ=255)
                                                                                                cA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAg   >  1:2098918/1‑100 (MQ=255)
                                                                                                cA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATGCTACCTTCTCGGCTGTCTCTTCTCTCGTACTGTTTAAGGATTTGTTCCGCTTTTCa     <  2:888772/98‑1 (MQ=255)
                                                                                                 a‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGc  <  1:977057/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
CAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGC  >  NC_000913/3183248‑3183444

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