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.

CCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGC  >  NC_000913/3183246‑3183444
                                                                                                    |                                                                                                   
cccAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCa                                                                                                      <  2:2705076/100‑1 (MQ=255)
 ccAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAg                                                                                                    >  2:71713/1‑100 (MQ=255)
 ccAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCACTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAg                                                                                                    <  1:3088486/100‑1 (MQ=255)
    gCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAgcgc                                                                                                 <  2:3784103/100‑1 (MQ=255)
       aTGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATg                                                                                              <  2:838462/100‑1 (MQ=255)
       aTGGCGTACTATATGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATg                                                                                              >  2:2894791/1‑100 (MQ=255)
        tGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATgg                                                                                             <  2:4144153/100‑1 (MQ=255)
        tGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATgg                                                                                             >  1:2506852/1‑100 (MQ=255)
         ggCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGa                                                                                            <  2:1578237/100‑1 (MQ=255)
          gCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGAc                                                                                           >  1:789365/1‑100 (MQ=255)
            gTACTATGTGGAATGTCAGTGATGGAATTCACTTTAGTACTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAAAAGTTATACAGCGCATGGACaa                                                                                         >  1:3479598/1‑100 (MQ=255)
                tATCTGGAATGTCAGTGATGAGATGCAGTTTAGCGCTGCTACAAACGGCGGGAACTGGGTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACt                                                                                       <  1:1209185/98‑1 (MQ=25)
                 aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCt                                                                                    >  1:4286585/1‑100 (MQ=255)
                 aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCt                                                                                    <  2:3401031/100‑1 (MQ=255)
                   gtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAAAAACATCAGAGGTTGTCCAGCGCATGGACAACTTATGc                                                                                  >  1:2589002/1‑100 (MQ=255)
                    tgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGc                                                                                  >  1:3286103/1‑99 (MQ=255)
                     ggAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAActgctgat                                                                                 <  1:2090242/99‑3 (MQ=38)
                     ggAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCTgatg                                                                                <  2:3690350/100‑4 (MQ=38)
                      gAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCTgatgt                                                                               <  2:3192283/100‑5 (MQ=37)
                        aTGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCt                                                                             <  2:3051969/100‑1 (MQ=255)
                         tGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGCAGTTGTCCAGCGCATGGACAACTTCTGCCAGCtt                                                                            >  1:2751942/1‑100 (MQ=255)
                         tGTCAGTGATGGGATTCAGTTTAGTGCTGCCACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTACTGCCAGCtt                                                                            >  1:2892837/1‑100 (MQ=255)
                          gTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTc                                                                           >  2:168496/1‑100 (MQ=255)
                           tCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCa                                                                          >  2:2127474/1‑100 (MQ=255)
                             aGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACg                                                                        >  2:4020561/1‑100 (MQ=255)
                             aGTGATGGAATTCAGTTTAGTGCTGCTACAAACGGCGCGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCTTTTTGACAACTTCTGCCAGCTTCACg                                                                        <  1:2050127/100‑1 (MQ=255)
                              gTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGa                                                                       >  2:3846952/1‑100 (MQ=255)
                              gTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGa                                                                       >  2:1071569/1‑100 (MQ=255)
                                 aTGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAg                                                                    <  2:1028236/100‑1 (MQ=255)
                                   gggATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAAGAACATCAGGAGTTGTCCAGCGCATGGACCACTTCTGCCAGCTTCACGAAAGcc                                                                  >  1:1781675/1‑100 (MQ=255)
                                   gggATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGcc                                                                  <  2:3288049/100‑1 (MQ=255)
                                     gATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCtt                                                                >  1:4354721/1‑100 (MQ=255)
                                     gATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCtt                                                                >  2:143853/1‑100 (MQ=255)
                                       ttCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTAt                                                              <  1:1517819/100‑1 (MQ=255)
                                        tCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTAtt                                                             <  1:1650586/100‑1 (MQ=255)
                                         cAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTc                                                            <  1:1470133/100‑1 (MQ=255)
                                           gTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTa                                                          <  2:4362126/100‑1 (MQ=255)
                                           gTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTa                                                          >  2:2854471/1‑100 (MQ=255)
                                           gTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTa                                                          >  1:3025046/1‑100 (MQ=255)
                                            tttAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTAc                                                         <  2:1155544/100‑1 (MQ=255)
                                             ttAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTAcc                                                        <  2:3042525/100‑1 (MQ=255)
                                             ttAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTAcc                                                        <  2:789366/100‑1 (MQ=255)
                                                   ctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCg                                                  <  2:4136481/100‑1 (MQ=255)
                                                    tgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCgg                                                 >  2:2243196/1‑100 (MQ=255)
                                                    tgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCgg                                                 >  2:2768374/1‑100 (MQ=255)
                                                     gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                <  1:3750096/100‑1 (MQ=255)
                                                     gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                <  2:501377/100‑1 (MQ=255)
                                                       tACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTg                                              <  1:2127467/100‑1 (MQ=255)
                                                            aCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctct                                         <  2:2547734/100‑1 (MQ=255)
                                                                cggGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttctc                                     >  1:447798/1‑100 (MQ=255)
                                                                cggGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttctc                                     >  1:1728302/1‑100 (MQ=255)
                                                                cggGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttct                                      >  2:3018336/1‑99 (MQ=255)
                                                                    aaTTGGTTAAAAACAACATCAGAAGTTGTCCAGCGTATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGt                                 <  1:3551353/100‑1 (MQ=255)
                                                                    aaTTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGt                                 >  1:2238590/1‑100 (MQ=255)
                                                                       tGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACt                              >  2:2479282/1‑100 (MQ=255)
                                                                        ggTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAACCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTg                             >  1:1629371/1‑100 (MQ=255)
                                                                         gTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGt                            >  1:925603/1‑100 (MQ=255)
                                                                         gTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGt                            >  1:1653183/1‑100 (MQ=255)
                                                                         gTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGt                            >  1:4075988/1‑100 (MQ=255)
                                                                         gTTAAAAACAACATCAGAAGCTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGt                            >  2:4151183/1‑100 (MQ=255)
                                                                           tAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGttt                          >  2:1278761/1‑100 (MQ=255)
                                                                           tAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGttt                          >  1:1280003/1‑100 (MQ=255)
                                                                           tAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGGCTCTTCTCTCGTACTGttt                          <  1:1805368/100‑1 (MQ=255)
                                                                            atctacaacaTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTa                         <  1:527186/96‑1 (MQ=255)
                                                                              aaaCAGCTTCAGAAGTTGTCCCGCGACTGGACAACTTCTGCCAGCTTCATGAAAGCCTTATTCTACCGTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAg                       <  1:3689934/100‑1 (MQ=255)
                                                                              aaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAg                       >  1:1105074/1‑100 (MQ=255)
                                                                              aaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAg                       <  1:3886966/100‑1 (MQ=255)
                                                                               aacaacATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGt                      <  2:655794/100‑1 (MQ=255)
                                                                               aacaacATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGt                      <  1:2660543/100‑1 (MQ=255)
                                                                                acaacaTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAgat                     >  2:527187/1‑98 (MQ=255)
                                                                                acaacaTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTa                     >  1:3598084/1‑100 (MQ=255)
                                                                                 caacaTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAt                    >  2:3219483/1‑100 (MQ=255)
                                                                                 caacaTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAt                    <  1:4184354/100‑1 (MQ=255)
                                                                                  aacaTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAtt                   <  1:2808905/100‑1 (MQ=255)
                                                                                    caTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTg                 <  2:30948/100‑1 (MQ=255)
                                                                                    caTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTg                 <  1:2939841/100‑1 (MQ=255)
                                                                                    caTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTg                 >  1:984144/1‑100 (MQ=255)
                                                                                     ggcAGAAGTTGTGCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTCTCCTACATTTTCGGCCGTCTCTTCTCTTGTATCGTTCACGTATTTg                 >  1:3046477/3‑99 (MQ=16)
                                                                                     ggcAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt                >  2:1890081/3‑100 (MQ=38)
                                                                                     aTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt                >  2:462110/1‑100 (MQ=255)
                                                                                     aTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt                >  2:3912236/1‑100 (MQ=255)
                                                                                      tcaggggTTGTCGAGCGCATGGACAACTTCGGCCGGCTTCACGAAAGCCTTATTCTGTCTTTTCGGCTGTCTCTTCTCTCGTGCTGTGTAAGTATTGGtt               <  1:3658047/94‑1 (MQ=255)
                                                                                      gcAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGtt               >  1:1913560/2‑100 (MQ=38)
                                                                                       cAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTc              >  2:3636207/1‑100 (MQ=38)
                                                                                       cAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTc              >  1:3298160/1‑100 (MQ=38)
                                                                                       cAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTc              >  2:2458922/1‑100 (MQ=255)
                                                                                        aGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTcc             >  2:3571619/1‑100 (MQ=255)
                                                                                         gAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  1:4094438/1‑100 (MQ=39)
                                                                                         gAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  2:758203/1‑100 (MQ=39)
                                                                                         gAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  1:630940/1‑100 (MQ=39)
                                                                                         gAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCGCGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  1:2365671/1‑100 (MQ=255)
                                                                                         gAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  1:3923181/1‑100 (MQ=255)
                                                                                          aaGTTGTCCAGCGCATGGACATCTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgc           <  1:2012652/100‑1 (MQ=255)
                                                                                          aaGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgc           <  1:3004217/100‑1 (MQ=255)
                                                                                           aGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgcg          >  1:3204609/1‑100 (MQ=39)
                                                                                           aGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgcg          >  2:4223480/1‑100 (MQ=255)
                                                                                           aGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTCTAAGTATTTGTTCcgcg          >  1:2402651/1‑100 (MQ=255)
                                                                                             ttGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCt        <  2:2774708/100‑1 (MQ=255)
                                                                                               gTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTc      >  1:2266873/1‑100 (MQ=255)
                                                                                                  cAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAg   >  2:3356828/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
CCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGC  >  NC_000913/3183246‑3183444

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