Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,183,345:1 +T 15.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.T15.2% 599.8 / 117.6 184intergenic (+22/+468)zupT/ribBzinc transporter/3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base . (92/64);  new base T (14/14);  total (106/78)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

Only 100 of 195 total aligned reads displayed.

CCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGC  >  NC_000913/3183246‑3183444
                                                                                                    |                                                                                                   
   aGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAgcg                                                                                                  >  1:2018802/1‑100 (MQ=255)
    gCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAgcgc                                                                                                 <  2:2256596/100‑1 (MQ=255)
         ggCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGa                                                                                            >  1:2462439/1‑100 (MQ=255)
         ggCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGa                                                                                            >  1:1436101/1‑100 (MQ=255)
              aCTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACt                                                                                       <  1:2701537/100‑1 (MQ=39)
                tATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTc                                                                                     <  1:911526/100‑1 (MQ=39)
                tATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTc                                                                                     <  1:2733050/100‑1 (MQ=39)
                 aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCt                                                                                    <  1:2989962/100‑1 (MQ=38)
                 aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCt                                                                                    <  1:1138543/100‑1 (MQ=38)
                 aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACctct                                                                                    >  1:3499292/1‑100 (MQ=255)
                 aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCt                                                                                    >  1:989022/1‑100 (MQ=255)
                 aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTACTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCt                                                                                    >  2:2507537/1‑100 (MQ=255)
                  tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTg                                                                                   <  2:533322/100‑1 (MQ=255)
                  tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTg                                                                                   <  2:1084120/100‑1 (MQ=255)
                  tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTg                                                                                   <  1:3253714/100‑1 (MQ=255)
                  tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTg                                                                                   <  2:816841/100‑1 (MQ=255)
                   gtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCTga                                                                                  <  1:637995/100‑2 (MQ=38)
                   gtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGc                                                                                  >  1:2933675/1‑100 (MQ=255)
                   gtgGAATGTCAGTGATGGGATTCAGTATAGTTCTGCTACAGACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCTg                                                                                   <  2:1119390/99‑1 (MQ=37)
                    tgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCTgat                                                                                 <  1:1718137/100‑3 (MQ=38)
                    tgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGcc                                                                                 >  1:1189753/1‑100 (MQ=255)
                    tagAATGTCAGTGATGGGATTCAGTATAGTGCTGCTACCAACGGCGGGAATTGGTTACAATCAACATCAGATGTTGTCCAGCGTATGGACTTCTTCTgtc                                                                                 <  2:1267271/98‑3 (MQ=255)
                      gtaggTCGGTGATGGGGTTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCTgatgt                                                                               <  2:1988066/96‑5 (MQ=25)
                       aaTGTCAGGGATAGATTTCAGTTTAGTGCTGCTACTAACAGCGGGAATTGGTTAAAAACCACATCTGAAGTTGTCCAGCGCATAGACACCTTCTg                                                                                   <  2:2830174/95‑1 (MQ=255)
                        aTGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCt                                                                             <  1:1307027/100‑1 (MQ=255)
                           tCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCa                                                                          <  1:3190405/100‑1 (MQ=255)
                           tCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCa                                                                          >  1:1616125/1‑100 (MQ=255)
                           tCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGCTTAAGAACATCATCAGAAGTTGTCCAGCGCATGGACAATTTCTGCCCGc                                                                              <  2:1026046/96‑1 (MQ=255)
                            cAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTGCAGCGCATGGACAACTTCTGCCAGCTTCAc                                                                         >  1:816677/1‑100 (MQ=255)
                            cAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCAc                                                                         >  2:1222184/1‑100 (MQ=255)
                                  tGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGc                                                                   >  1:2308831/1‑100 (MQ=255)
                                  tGGGATTCAGTTTAGTGCTGCTAAAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGGCAACTTCTGCCAGCTTCACGa                                                                       <  2:262020/96‑1 (MQ=255)
                                      aTTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTa                                                               <  1:1878702/100‑1 (MQ=255)
                                       ttCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTGCACGAAAGCCTTAt                                                              <  2:265705/100‑1 (MQ=255)
                                        tCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTAtt                                                             <  2:3411908/100‑1 (MQ=255)
                                           gTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTa                                                          >  1:1579895/1‑100 (MQ=255)
                                           gTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTa                                                          >  1:425667/1‑100 (MQ=255)
                                            tttAGTGCTGCTACAAACGGCGCGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTAc                                                         <  2:891873/100‑1 (MQ=255)
                                             ttAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTAcc                                                        <  1:2164492/100‑1 (MQ=255)
                                             ttAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTAcc                                                        <  1:2808381/100‑1 (MQ=255)
                                             ttAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTAcc                                                        <  2:875943/100‑1 (MQ=255)
                                               aGTGCTGGTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGCCAACTTCTGCCCGCTTCACGAAAGCCTTATCCTACCtt                                                      >  2:3761768/1‑100 (MQ=255)
                                               aGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCtt                                                      >  2:2867892/1‑100 (MQ=255)
                                                 tgctgcTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCtttt                                                    >  2:1661728/1‑100 (MQ=255)
                                                  gctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTc                                                   <  1:2150626/100‑1 (MQ=255)
                                                  gctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTGATTCTACCGTTTc                                                   >  1:2473515/1‑100 (MQ=255)
                                                   ctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGCAAGCCTTATTCTACCTTTTCg                                                  >  2:1114211/1‑100 (MQ=255)
                                                   ctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTACAGCGCATGGACAACTTCTGCCAGATTCACGAAAGCCTTATTCTACCtttt                                                    >  2:1582487/1‑98 (MQ=255)
                                                     gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCgg                                                 >  2:897990/1‑99 (MQ=255)
                                                     gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                >  1:351508/1‑100 (MQ=255)
                                                     gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                <  1:1409227/100‑1 (MQ=255)
                                                        aCAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGt                                             >  1:2538942/1‑100 (MQ=255)
                                                           aaCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctc                                          >  2:2096582/1‑100 (MQ=255)
                                                              ggcggGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttc                                       >  1:2437125/1‑100 (MQ=255)
                                                              ggcggGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttc                                       >  2:2465797/1‑100 (MQ=255)
                                                                cggGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttctc                                     <  2:328205/100‑1 (MQ=255)
                                                                  ggAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTtctctc                                   <  1:1243001/100‑1 (MQ=255)
                                                                    tccgatctAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGt                                 <  2:1614704/93‑1 (MQ=255)
                                                                    aaTTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTCTCGGCTGTCTCTTCTCTCGt                                 >  2:3219414/1‑100 (MQ=255)
                                                                      ttGGTTAAAAACACCATCAGAAGTGGTCCAGCGCATGGACAACTTCGGCCAGCTGCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTAc                               <  2:3236291/100‑1 (MQ=255)
                                                                        ggTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTg                             >  2:644777/1‑100 (MQ=255)
                                                                          ttAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGtt                           >  1:1625291/1‑100 (MQ=255)
                                                                            aaaaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTa                         >  2:1280756/1‑100 (MQ=255)
                                                                            aaaaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTa                         <  2:907756/100‑1 (MQ=255)
                                                                              aaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAg                       <  1:3291307/100‑1 (MQ=255)
                                                                               aacaacATCAGAAGTTGTCCAGCGCATGGGCAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTGAAg                       <  2:2605946/99‑1 (MQ=255)
                                                                                acaacaTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTa                     >  1:290127/1‑100 (MQ=255)
                                                                                 caacaTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGATTAAGTAt                    >  2:1897261/1‑100 (MQ=255)
                                                                                   ctggcAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAttt                  >  1:752790/5‑100 (MQ=37)
                                                                                   acaTCAGAAGTTGTCCAGCGCATTGACAACTTCTGCCAGGTTCACGAAAGCCTTATTCTATCTTTTCGGCTGTCTCTGCTCTCGTACTGTTTAAGTAttt                  <  2:2825257/100‑1 (MQ=255)
                                                                                   acaTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAttt                  >  1:1488045/1‑100 (MQ=255)
                                                                                    tggcAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTg                 >  1:1788401/4‑100 (MQ=38)
                                                                                    caTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTg                 >  1:1086572/1‑100 (MQ=255)
                                                                                    caTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTg                 >  1:3387910/1‑100 (MQ=255)
                                                                                     ggcAGAGGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt                >  2:2755920/3‑100 (MQ=38)
                                                                                     aTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt                >  1:890713/1‑100 (MQ=255)
                                                                                      tCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGtt               <  1:1945542/100‑1 (MQ=255)
                                                                                      tCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGtt               >  1:2190229/1‑100 (MQ=255)
                                                                                       ccgAAGTTGTCCAGCGCATGGACAACTTTTGCCTGCTTCGCGAATGCCTTATTCTAGCTTTTCGGCTCTCTCTTCTCTCGTACTGTTTAAGTATTTGtt               <  2:3138335/97‑1 (MQ=255)
                                                                                       cAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTc              >  1:2857210/1‑100 (MQ=255)
                                                                                        aGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTcc             >  2:2998770/1‑100 (MQ=38)
                                                                                        aGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTcc             <  1:484380/100‑1 (MQ=255)
                                                                                        aGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTcc             >  2:2944289/1‑100 (MQ=255)
                                                                                         gAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  1:3645369/1‑100 (MQ=39)
                                                                                         gAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  1:2524557/1‑100 (MQ=39)
                                                                                         gAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  1:3350155/1‑100 (MQ=39)
                                                                                         gAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  1:3808655/1‑100 (MQ=39)
                                                                                         gAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  2:1330636/1‑100 (MQ=39)
                                                                                         gAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            <  1:3166295/100‑1 (MQ=255)
                                                                                          aaGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgc           >  2:177987/1‑100 (MQ=39)
                                                                                           aGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgcg          >  2:2745432/1‑100 (MQ=39)
                                                                                            gTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCgcgc         >  1:2226265/1‑100 (MQ=255)
                                                                                            gTTGTCCAGCGCAGGGACAACTGCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCGCTCGTACTGTTTAAGTATTTGTTCCgcgc         <  2:1522107/100‑1 (MQ=255)
                                                                                              tGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGttagatcgga       >  1:3272895/1‑92 (MQ=255)
                                                                                              tGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCtt       >  1:1032295/1‑100 (MQ=255)
                                                                                                 ccAGCGCATGGACGACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCaa    >  1:2594894/1‑100 (MQ=255)
                                                                                                  cAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAg   >  1:3743467/1‑100 (MQ=255)
                                                                                                   aGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGc  <  1:3087396/100‑1 (MQ=255)
                                                                                                   aGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGc  <  1:962532/100‑1 (MQ=255)
                                                                                                   aGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGc  <  1:1543269/100‑1 (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: