Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,183,345:1 +T 16.5% 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.T16.5% 619.8 / 136.8 194intergenic (+22/+468)zupT/ribBzinc transporter/3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base . (74/88);  new base T (22/10);  total (96/98)
Fisher's exact test for biased strand distribution p-value = 2.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

Only 100 of 199 total aligned reads displayed.

CCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGC  >  NC_000913/3183247‑3183444
                                                                                                   |                                                                                                   
     tATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCAt                                                                                               <  1:3860657/100‑1 (MQ=255)
     tATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCAt                                                                                               <  2:1001030/100‑1 (MQ=255)
      aTGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATg                                                                                              <  1:77745/100‑1 (MQ=255)
        ggCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGa                                                                                            >  2:3792221/1‑100 (MQ=255)
        ggCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGa                                                                                            >  2:2250422/1‑100 (MQ=255)
         gCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGAc                                                                                           >  2:3922125/1‑100 (MQ=255)
          cgttctagatggaatggaaGTGATGGGAGTCAGTTGAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCGAGCGCAGGGACa                                                                                          <  1:4121766/91‑1 (MQ=255)
          cGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACa                                                                                          >  2:1079153/1‑100 (MQ=255)
            tACTATGTGGAATGTCGGTGGTGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAGGTTGTCCAGCGCATGGACAAc                                                                                        <  1:976116/100‑1 (MQ=255)
             cttATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACAGCAGAAGTTGTCCATGCGCTGGACAACt                                                                                       <  1:450727/98‑1 (MQ=38)
             aCTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACt                                                                                       <  2:2197032/100‑1 (MQ=39)
              cTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACtt                                                                                      >  2:2458346/1‑100 (MQ=255)
              cTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACtt                                                                                      <  1:2385691/100‑1 (MQ=255)
               tATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTc                                                                                     >  1:4194670/1‑100 (MQ=255)
               tATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACAGGAGAAGTTGTCCATGCGCTGGACAACtag                                                                                     >  2:450729/1‑98 (MQ=38)
                aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTc                                                                                     <  1:2842544/99‑1 (MQ=38)
                aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCt                                                                                    <  2:4062044/100‑1 (MQ=38)
                aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCt                                                                                    <  1:3538190/100‑1 (MQ=255)
                aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCt                                                                                    <  2:4035469/100‑1 (MQ=255)
                 tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAGCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTg                                                                                   <  1:595520/100‑1 (MQ=255)
                 tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCt                                                                                    <  1:650593/99‑1 (MQ=38)
                 tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTg                                                                                   <  2:2952747/100‑1 (MQ=255)
                   tgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGcc                                                                                 <  2:1753728/100‑1 (MQ=255)
                   tgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGcc                                                                                 <  1:2303373/100‑1 (MQ=255)
                    ggAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCa                                                                                >  2:1612570/1‑100 (MQ=255)
                     gAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTCCTGCCAg                                                                               <  1:2042056/100‑1 (MQ=255)
                       aTGTCAGTGATGGGATTCAGTTTAGTGCTGCTAGAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGc                                                                              <  1:4030824/99‑1 (MQ=255)
                          tCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCa                                                                          <  2:1168159/100‑1 (MQ=255)
                           cAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCTgatgttgtttt                                                                         <  2:3036062/100‑11 (MQ=25)
                            aGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACg                                                                        <  2:3029955/100‑1 (MQ=255)
                            aGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACg                                                                        <  1:4221744/100‑1 (MQ=255)
                               gATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGaaa                                                                     >  2:3735310/1‑100 (MQ=255)
                                aTGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATTGACAACTTCTGCCAGCTTCACGAAAg                                                                    <  2:2680780/100‑1 (MQ=255)
                                 tGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGc                                                                   <  2:3618893/100‑1 (MQ=255)
                                 tGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGc                                                                   >  2:919078/1‑100 (MQ=255)
                                  gggATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGcc                                                                  >  2:1251912/1‑100 (MQ=255)
                                       tCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTAt                                                              <  1:523179/99‑1 (MQ=255)
                                       tCAGTTTAGTGCTGCTACAAACGGCGGGAATAACATAAAAACAACATCAGACGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGCAAGCCTTAtt                                                             <  1:3508024/100‑1 (MQ=255)
                                        cAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTc                                                            >  1:3181254/1‑100 (MQ=255)
                                         aGTTTAGTGGTGCTACAAACGGCGGGAATTGGGTAAAAACAACATCAGAAGTTGGCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAGGCCTTATTc                                                            <  1:1038361/99‑1 (MQ=255)
                                          gTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTa                                                          >  1:1241505/1‑100 (MQ=255)
                                            ttAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACCACGTCAGACGTTGTCCAGCGCATGGACAACTTGTGCCAGCTTCACGAAAGCCTTATTCTccc                                                        <  1:201290/100‑1 (MQ=255)
                                              agggctgctACAAACGGCGGGACGTGGTTAAAAACAACATCAGAAGTGGTGCAGCGCATGGACAACTTCTGCCAGCTTCGCGAAAGCCTTATTCTAcc                                                        <  1:1605250/95‑1 (MQ=255)
                                              aGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCtt                                                      <  2:4355060/100‑1 (MQ=255)
                                              aGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCtt                                                      <  1:2172791/100‑1 (MQ=255)
                                                tgctgcTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCtttt                                                    >  1:1337923/1‑100 (MQ=255)
                                                 gctgctACAAATGGCGGGAATTGGTTAAAAACAACATTAGAAGTTGTCAAGCGCATGGACAAATTATGCCAGCTTCACGAAAGCCTGATTCTACATTATc                                                   >  1:4160892/1‑100 (MQ=255)
                                                   tgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCgg                                                 >  2:2167990/1‑100 (MQ=255)
                                                    gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACACCTTCTGCCAGCCTCACGAAAGCCTTATTCTACCTTTTCGGc                                                <  1:427648/100‑1 (MQ=255)
                                                    gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                >  1:244280/1‑100 (MQ=255)
                                                    gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                >  2:2379167/1‑100 (MQ=255)
                                                    gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                >  2:1920402/1‑100 (MQ=255)
                                                    gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                >  2:1198643/1‑100 (MQ=255)
                                                      tACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTg                                              <  2:1759411/100‑1 (MQ=255)
                                                      tACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTg                                              <  2:2671416/100‑1 (MQ=255)
                                                       aCAAACGGCGGGAATTGGTTAAAAACAACACCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAACGCCTTATTCTACCTTGTGGGCTg                                              <  1:2789275/99‑1 (MQ=255)
                                                         aaaCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtct                                           <  1:4199411/100‑1 (MQ=255)
                                                          aaCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctc                                          <  1:2279329/100‑1 (MQ=255)
                                                          aaCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctc                                          >  2:2402501/1‑100 (MQ=255)
                                                           aCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctct                                         >  2:648794/1‑100 (MQ=255)
                                                            cggcggGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctctt                                        <  2:60021/100‑1 (MQ=255)
                                                               cggGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttct                                      >  2:3842476/1‑99 (MQ=255)
                                                                gggAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttctct                                    <  1:4000465/100‑1 (MQ=255)
                                                                  gAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCg                                  >  2:1302014/1‑100 (MQ=255)
                                                                   aaTTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGt                                 <  2:2800281/100‑1 (MQ=255)
                                                                    aTTGGTTAAAAACAACATCATATGTTGGCCAGCGCATGGACATCTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGt                                 <  1:1999077/99‑1 (MQ=255)
                                                                     ttGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACACCTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTAc                               <  2:4285719/100‑1 (MQ=255)
                                                                       ggTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTg                             <  2:716854/100‑1 (MQ=255)
                                                                       ggTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTg                             >  2:1424471/1‑100 (MQ=255)
                                                                        gTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGt                            >  2:16946/1‑100 (MQ=255)
                                                                        gTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGt                            >  2:3599492/1‑100 (MQ=255)
                                                                          tAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACGTCTGCCAGCTTCACGAAAGCCTTATTCGACCTTTTCTGCGGTCTCGTCTCTCGTACGGtgt                          >  1:3668342/1‑98 (MQ=255)
                                                                           aaaaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTa                         >  2:2104767/1‑100 (MQ=255)
                                                                           aaaaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTa                         <  1:2645088/100‑1 (MQ=255)
                                                                            aaaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTaa                        >  1:2504089/1‑100 (MQ=255)
                                                                             aaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAg                       <  2:807051/100‑1 (MQ=255)
                                                                             aaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAg                       <  2:3520488/100‑1 (MQ=255)
                                                                               acaacaTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTa                     >  1:1546802/1‑100 (MQ=255)
                                                                                  acaTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAttt                  <  1:1211122/100‑1 (MQ=255)
                                                                                   tggcAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTg                 >  2:2527182/4‑100 (MQ=38)
                                                                                    ggcAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt                >  2:1579304/3‑100 (MQ=38)
                                                                                    ggcAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt                >  1:1481909/3‑100 (MQ=38)
                                                                                     tCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGtt               >  2:66597/1‑100 (MQ=255)
                                                                                      cAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGtt               >  2:1761114/1‑99 (MQ=38)
                                                                                      cAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGtt               >  2:1447716/1‑99 (MQ=255)
                                                                                      cAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTc              >  1:2343115/1‑100 (MQ=255)
                                                                                        gAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  2:523135/1‑100 (MQ=39)
                                                                                         aaGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTGCcgc           >  2:3404535/1‑100 (MQ=38)
                                                                                          aGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgc           >  2:1523503/1‑99 (MQ=39)
                                                                                          aGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTATAGCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTCAGTATTTGTTCcgcg          >  2:395523/1‑100 (MQ=38)
                                                                                          aGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgcg          <  2:155816/100‑1 (MQ=255)
                                                                                           gTTGTCCAGCGTATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTTCCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgcg          <  1:1952411/99‑1 (MQ=255)
                                                                                           gTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCgcgc         >  1:3253942/1‑100 (MQ=255)
                                                                                           gTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCgcgc         >  2:1512726/1‑100 (MQ=255)
                                                                                           gTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCgcgc         >  2:447578/1‑100 (MQ=255)
                                                                                              gTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTc      >  2:2631201/1‑100 (MQ=255)
                                                                                                 cAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAg   >  1:1058337/1‑100 (MQ=255)
                                                                                                 cAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAg   <  2:3170892/100‑1 (MQ=255)
                                                                                                  aGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGc  >  2:4191400/1‑100 (MQ=255)
                                                                                                  aGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGc  >  2:1478739/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGC  >  NC_000913/3183247‑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: