Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,183,345:1 +T 19.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.T19.3% 565.1 / 163.9 197intergenic (+22/+468)zupT/ribBzinc transporter/3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base . (85/74);  new base T (26/12);  total (111/86)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

Only 100 of 205 total aligned reads displayed.

CCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGC  >  NC_000913/3183246‑3183444
                                                                                                    |                                                                                                   
cccAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCa                                                                                                      <  2:2705076/100‑1 (MQ=255)
 ccAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑g                                                                                                    >  2:71713/1‑100 (MQ=255)
 ccAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCACTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑g                                                                                                    <  1:3088486/100‑1 (MQ=255)
    gCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑gcgc                                                                                                 <  2:3784103/100‑1 (MQ=255)
       aTGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATg                                                                                              <  2:838462/100‑1 (MQ=255)
       aTGGCGTACTATATGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATg                                                                                              >  2:2894791/1‑100 (MQ=255)
        tGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATgg                                                                                             <  2:4144153/100‑1 (MQ=255)
        tGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATgg                                                                                             >  1:2506852/1‑100 (MQ=255)
         ggCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGa                                                                                            <  2:1578237/100‑1 (MQ=255)
          gCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGAc                                                                                           >  1:789365/1‑100 (MQ=255)
            gTACTATGTGGAATGTCAGTGATGGAATTCACTTTAGTACTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAAAAGTTATACA‑GCGCATGGACaa                                                                                         >  1:3479598/1‑100 (MQ=255)
                tATCTGGAATGTCAGTGATGAGATGCAGTTTAGCGCTGCTACAAACGGCGGGAACTGGGTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAACt                                                                                       <  1:1209185/98‑1 (MQ=25)
                 aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCt                                                                                    >  1:4286585/1‑100 (MQ=255)
                 aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCt                                                                                    <  2:3401031/100‑1 (MQ=255)
                   gtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAAAAACATCAGAGGTTGTCCA‑GCGCATGGACAACTTATGc                                                                                  >  1:2589002/1‑100 (MQ=255)
                    tgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGc                                                                                  >  1:3286103/1‑99 (MQ=255)
                     ggAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAActgctgat                                                                                 <  1:2090242/99‑3 (MQ=38)
                     ggAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAACTTCTgatg                                                                                <  2:3690350/100‑4 (MQ=38)
                      gAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAACTTCTgatgt                                                                               <  2:3192283/100‑5 (MQ=37)
                        aTGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCt                                                                             <  2:3051969/100‑1 (MQ=255)
                         tGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGCAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCtt                                                                            >  1:2751942/1‑100 (MQ=255)
                         tGTCAGTGATGGGATTCAGTTTAGTGCTGCCACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTACTGCCAGCtt                                                                            >  1:2892837/1‑100 (MQ=255)
                          gTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTc                                                                           >  2:168496/1‑100 (MQ=255)
                           tCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCa                                                                          >  2:2127474/1‑100 (MQ=255)
                             aGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACg                                                                        >  2:4020561/1‑100 (MQ=255)
                             aGTGATGGAATTCAGTTTAGTGCTGCTACAAACGGCGCGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCTTTTTGACAACTTCTGCCAGCTTCACg                                                                        <  1:2050127/100‑1 (MQ=255)
                              gTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGa                                                                       >  2:3846952/1‑100 (MQ=255)
                              gTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGa                                                                       >  2:1071569/1‑100 (MQ=255)
                                 aTGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAg                                                                    <  2:1028236/100‑1 (MQ=255)
                                   gggATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAAGAACATCAGGAGTTGTCCA‑GCGCATGGACCACTTCTGCCAGCTTCACGAAAGcc                                                                  >  1:1781675/1‑100 (MQ=255)
                                   gggATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGcc                                                                  <  2:3288049/100‑1 (MQ=255)
                                     gATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCtt                                                                >  1:4354721/1‑100 (MQ=255)
                                     gATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCtt                                                                >  2:143853/1‑100 (MQ=255)
                                       ttCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTAt                                                              <  1:1517819/100‑1 (MQ=255)
                                        tCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTAtt                                                             <  1:1650586/100‑1 (MQ=255)
                                         cAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTc                                                            <  1:1470133/100‑1 (MQ=255)
                                           gTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTa                                                          <  2:4362126/100‑1 (MQ=255)
                                           gTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTa                                                          >  2:2854471/1‑100 (MQ=255)
                                           gTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTa                                                          >  1:3025046/1‑100 (MQ=255)
                                            tttAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTAc                                                         <  2:1155544/100‑1 (MQ=255)
                                             ttAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTAcc                                                        <  2:3042525/100‑1 (MQ=255)
                                             ttAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTAcc                                                        <  2:789366/100‑1 (MQ=255)
                                                   ctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCg                                                  <  2:4136481/100‑1 (MQ=255)
                                                    tgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCgg                                                 >  2:2243196/1‑100 (MQ=255)
                                                    tgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCgg                                                 >  2:2768374/1‑100 (MQ=255)
                                                     gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                <  1:3750096/100‑1 (MQ=255)
                                                     gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                <  2:501377/100‑1 (MQ=255)
                                                       tACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTg                                              <  1:2127467/100‑1 (MQ=255)
                                                            aCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctct                                         <  2:2547734/100‑1 (MQ=255)
                                                                cggGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttctc                                     >  1:447798/1‑100 (MQ=255)
                                                                cggGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttctc                                     >  1:1728302/1‑100 (MQ=255)
                                                                cggGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttct                                      >  2:3018336/1‑99 (MQ=255)
                                                                    aaTTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGTATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGt                                 <  1:3551353/100‑1 (MQ=255)
                                                                    aaTTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGt                                 >  1:2238590/1‑100 (MQ=255)
                                                                       tGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACt                              >  2:2479282/1‑100 (MQ=255)
                                                                        ggTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAACCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTg                             >  1:1629371/1‑100 (MQ=255)
                                                                         gTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGt                            >  1:925603/1‑100 (MQ=255)
                                                                         gTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGt                            >  1:1653183/1‑100 (MQ=255)
                                                                         gTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGt                            >  1:4075988/1‑100 (MQ=255)
                                                                         gTTAAAAACAACATCAGAAGCTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGt                            >  2:4151183/1‑100 (MQ=255)
                                                                           tAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGttt                          >  2:1278761/1‑100 (MQ=255)
                                                                           tAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGttt                          >  1:1280003/1‑100 (MQ=255)
                                                                           tAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGGCTCTTCTCTCGTACTGttt                          <  1:1805368/100‑1 (MQ=255)
                                                                            atctacaacaTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTa                         <  1:527186/96‑1 (MQ=255)
                                                                              aaaCAGCTTCAGAAGTTGTCCC‑GCGACTGGACAACTTCTGCCAGCTTCATGAAAGCCTTATTCTACCGTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAg                       <  1:3689934/100‑1 (MQ=255)
                                                                              aaaCAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAg                       >  1:1105074/1‑100 (MQ=255)
                                                                              aaaCAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAg                       <  1:3886966/100‑1 (MQ=255)
                                                                               aacaacATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGt                      <  2:655794/100‑1 (MQ=255)
                                                                               aacaacATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGt                      <  1:2660543/100‑1 (MQ=255)
                                                                                acaacaTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAgat                     >  2:527187/1‑98 (MQ=255)
                                                                                acaacaTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTa                     >  1:3598084/1‑100 (MQ=255)
                                                                                 caacaTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAt                    >  2:3219483/1‑100 (MQ=255)
                                                                                 caacaTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAt                    <  1:4184354/100‑1 (MQ=255)
                                                                                  aacaTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAtt                   <  1:2808905/100‑1 (MQ=255)
                                                                                    caTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTg                 <  2:30948/100‑1 (MQ=255)
                                                                                    caTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTg                 <  1:2939841/100‑1 (MQ=255)
                                                                                    caTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTg                 >  1:984144/1‑100 (MQ=255)
                                                                                     ggcAGAAGTTGTGCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTCTCCTACATTTTCGGCCGTCTCTTCTCTTGTATCGTTCACGTATTTg                 >  1:3046477/3‑99 (MQ=16)
                                                                                     ggcAGAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt                >  2:1890081/3‑100 (MQ=38)
                                                                                     aTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt                >  2:462110/1‑100 (MQ=255)
                                                                                     aTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt                >  2:3912236/1‑100 (MQ=255)
                                                                                      tcaggggTTGTCGA‑GCGCATGGACAACTTCGGCCGGCTTCACGAAAGCCTTATTCTGTCTTTTCGGCTGTCTCTTCTCTCGTGCTGTGTAAGTATTGGtt               <  1:3658047/94‑1 (MQ=255)
                                                                                      gcAGAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGtt               >  1:1913560/2‑100 (MQ=38)
                                                                                       cAGAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTc              >  2:3636207/1‑100 (MQ=38)
                                                                                       cAGAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTc              >  1:3298160/1‑100 (MQ=38)
                                                                                       cAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTc              >  2:2458922/1‑100 (MQ=255)
                                                                                        aGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTcc             >  2:3571619/1‑100 (MQ=255)
                                                                                         gAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  1:4094438/1‑100 (MQ=39)
                                                                                         gAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  2:758203/1‑100 (MQ=39)
                                                                                         gAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  1:630940/1‑100 (MQ=39)
                                                                                         gAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCGCGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  1:2365671/1‑100 (MQ=255)
                                                                                         gAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  1:3923181/1‑100 (MQ=255)
                                                                                          aaGTTGTCCA‑GCGCATGGACATCTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgc           <  1:2012652/100‑1 (MQ=255)
                                                                                          aaGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgc           <  1:3004217/100‑1 (MQ=255)
                                                                                           aGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgcg          >  1:3204609/1‑100 (MQ=39)
                                                                                           aGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgcg          >  2:4223480/1‑100 (MQ=255)
                                                                                           aGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTCTAAGTATTTGTTCcgcg          >  1:2402651/1‑100 (MQ=255)
                                                                                             ttGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCt        <  2:2774708/100‑1 (MQ=255)
                                                                                               gTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTc      >  1:2266873/1‑100 (MQ=255)
                                                                                                  cA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAg   >  2:3356828/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
CCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGC  >  NC_000913/3183246‑3183444

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