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.

CCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGC  >  NC_000913/3183247‑3183444
                                                                                                   |                                                                                                   
     tATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCAt                                                                                               <  1:3860657/100‑1 (MQ=255)
     tATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCAt                                                                                               <  2:1001030/100‑1 (MQ=255)
      aTGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATg                                                                                              <  1:77745/100‑1 (MQ=255)
        ggCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGa                                                                                            >  2:3792221/1‑100 (MQ=255)
        ggCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGa                                                                                            >  2:2250422/1‑100 (MQ=255)
         gCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGAc                                                                                           >  2:3922125/1‑100 (MQ=255)
          cgttctagatggaatggaaGTGATGGGAGTCAGTTGAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCGA‑GCGCAGGGACa                                                                                          <  1:4121766/91‑1 (MQ=255)
          cGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACa                                                                                          >  2:1079153/1‑100 (MQ=255)
            tACTATGTGGAATGTCGGTGGTGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAGGTTGTCCA‑GCGCATGGACAAc                                                                                        <  1:976116/100‑1 (MQ=255)
             cttATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACAGCAGAAGTTGTCCATGCGC‑TGGACAACt                                                                                       <  1:450727/98‑1 (MQ=38)
             aCTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAACt                                                                                       <  2:2197032/100‑1 (MQ=39)
              cTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACtt                                                                                      >  2:2458346/1‑100 (MQ=255)
              cTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACtt                                                                                      <  1:2385691/100‑1 (MQ=255)
               tATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTc                                                                                     >  1:4194670/1‑100 (MQ=255)
               tATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACAGGAGAAGTTGTCCATGCGC‑TGGACAACtag                                                                                     >  2:450729/1‑98 (MQ=38)
                aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAACTTc                                                                                     <  1:2842544/99‑1 (MQ=38)
                aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAACTTCt                                                                                    <  2:4062044/100‑1 (MQ=38)
                aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCt                                                                                    <  1:3538190/100‑1 (MQ=255)
                aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCt                                                                                    <  2:4035469/100‑1 (MQ=255)
                 tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAGCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTg                                                                                   <  1:595520/100‑1 (MQ=255)
                 tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAACTTCt                                                                                    <  1:650593/99‑1 (MQ=38)
                 tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTg                                                                                   <  2:2952747/100‑1 (MQ=255)
                   tgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGcc                                                                                 <  2:1753728/100‑1 (MQ=255)
                   tgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGcc                                                                                 <  1:2303373/100‑1 (MQ=255)
                    ggAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCa                                                                                >  2:1612570/1‑100 (MQ=255)
                     gAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTCCTGCCAg                                                                               <  1:2042056/100‑1 (MQ=255)
                       aTGTCAGTGATGGGATTCAGTTTAGTGCTGCTAGAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGc                                                                              <  1:4030824/99‑1 (MQ=255)
                          tCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCa                                                                          <  2:1168159/100‑1 (MQ=255)
                           cAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAACTTCTgatgttgtttt                                                                         <  2:3036062/100‑11 (MQ=25)
                            aGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACg                                                                        <  2:3029955/100‑1 (MQ=255)
                            aGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACg                                                                        <  1:4221744/100‑1 (MQ=255)
                               gATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGaaa                                                                     >  2:3735310/1‑100 (MQ=255)
                                aTGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATTGACAACTTCTGCCAGCTTCACGAAAg                                                                    <  2:2680780/100‑1 (MQ=255)
                                 tGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGc                                                                   <  2:3618893/100‑1 (MQ=255)
                                 tGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGc                                                                   >  2:919078/1‑100 (MQ=255)
                                  gggATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGcc                                                                  >  2:1251912/1‑100 (MQ=255)
                                       tCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTAt                                                              <  1:523179/99‑1 (MQ=255)
                                       tCAGTTTAGTGCTGCTACAAACGGCGGGAATAACATAAAAACAACATCAGACGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGCAAGCCTTAtt                                                             <  1:3508024/100‑1 (MQ=255)
                                        cAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTc                                                            >  1:3181254/1‑100 (MQ=255)
                                         aGTTTAGTGGTGCTACAAACGGCGGGAATTGGGTAAAAACAACATCAGAAGTTGGCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAGGCCTTATTc                                                            <  1:1038361/99‑1 (MQ=255)
                                          gTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTa                                                          >  1:1241505/1‑100 (MQ=255)
                                            ttAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACCACGTCAGACGTTGTCCA‑GCGCATGGACAACTTGTGCCAGCTTCACGAAAGCCTTATTCTccc                                                        <  1:201290/100‑1 (MQ=255)
                                              agggctgctACAAACGGCGGGACGTGGTTAAAAACAACATCAGAAGTGGTGCA‑GCGCATGGACAACTTCTGCCAGCTTCGCGAAAGCCTTATTCTAcc                                                        <  1:1605250/95‑1 (MQ=255)
                                              aGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCtt                                                      <  2:4355060/100‑1 (MQ=255)
                                              aGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCtt                                                      <  1:2172791/100‑1 (MQ=255)
                                                tgctgcTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCtttt                                                    >  1:1337923/1‑100 (MQ=255)
                                                 gctgctACAAATGGCGGGAATTGGTTAAAAACAACATTAGAAGTTGTCAA‑GCGCATGGACAAATTATGCCAGCTTCACGAAAGCCTGATTCTACATTATc                                                   >  1:4160892/1‑100 (MQ=255)
                                                   tgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCgg                                                 >  2:2167990/1‑100 (MQ=255)
                                                    gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACACCTTCTGCCAGCCTCACGAAAGCCTTATTCTACCTTTTCGGc                                                <  1:427648/100‑1 (MQ=255)
                                                    gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                >  1:244280/1‑100 (MQ=255)
                                                    gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                >  2:2379167/1‑100 (MQ=255)
                                                    gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                >  2:1920402/1‑100 (MQ=255)
                                                    gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                >  2:1198643/1‑100 (MQ=255)
                                                      tACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTg                                              <  2:1759411/100‑1 (MQ=255)
                                                      tACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTg                                              <  2:2671416/100‑1 (MQ=255)
                                                       aCAAACGGCGGGAATTGGTTAAAAACAACACCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAACGCCTTATTCTACCTTGTGGGCTg                                              <  1:2789275/99‑1 (MQ=255)
                                                         aaaCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtct                                           <  1:4199411/100‑1 (MQ=255)
                                                          aaCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctc                                          <  1:2279329/100‑1 (MQ=255)
                                                          aaCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctc                                          >  2:2402501/1‑100 (MQ=255)
                                                           aCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctct                                         >  2:648794/1‑100 (MQ=255)
                                                            cggcggGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctctt                                        <  2:60021/100‑1 (MQ=255)
                                                               cggGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttct                                      >  2:3842476/1‑99 (MQ=255)
                                                                gggAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttctct                                    <  1:4000465/100‑1 (MQ=255)
                                                                  gAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCg                                  >  2:1302014/1‑100 (MQ=255)
                                                                   aaTTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGt                                 <  2:2800281/100‑1 (MQ=255)
                                                                    aTTGGTTAAAAACAACATCATATGTTGGCCA‑GCGCATGGACATCTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGt                                 <  1:1999077/99‑1 (MQ=255)
                                                                     ttGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACACCTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTAc                               <  2:4285719/100‑1 (MQ=255)
                                                                       ggTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTg                             <  2:716854/100‑1 (MQ=255)
                                                                       ggTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTg                             >  2:1424471/1‑100 (MQ=255)
                                                                        gTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGt                            >  2:16946/1‑100 (MQ=255)
                                                                        gTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGt                            >  2:3599492/1‑100 (MQ=255)
                                                                          tAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACGTCTGCCAGCTTCACGAAAGCCTTATTCGACCTTTTCTGCGGTCTCGTCTCTCGTACGGtgt                          >  1:3668342/1‑98 (MQ=255)
                                                                           aaaaaCAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTa                         >  2:2104767/1‑100 (MQ=255)
                                                                           aaaaaCAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTa                         <  1:2645088/100‑1 (MQ=255)
                                                                            aaaaCAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTaa                        >  1:2504089/1‑100 (MQ=255)
                                                                             aaaCAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAg                       <  2:807051/100‑1 (MQ=255)
                                                                             aaaCAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAg                       <  2:3520488/100‑1 (MQ=255)
                                                                               acaacaTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTa                     >  1:1546802/1‑100 (MQ=255)
                                                                                  acaTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAttt                  <  1:1211122/100‑1 (MQ=255)
                                                                                   tggcAGAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTg                 >  2:2527182/4‑100 (MQ=38)
                                                                                    ggcAGAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt                >  2:1579304/3‑100 (MQ=38)
                                                                                    ggcAGAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt                >  1:1481909/3‑100 (MQ=38)
                                                                                     tCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGtt               >  2:66597/1‑100 (MQ=255)
                                                                                      cAGAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGtt               >  2:1761114/1‑99 (MQ=38)
                                                                                      cAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGtt               >  2:1447716/1‑99 (MQ=255)
                                                                                      cAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTc              >  1:2343115/1‑100 (MQ=255)
                                                                                        gAAGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg            >  2:523135/1‑100 (MQ=39)
                                                                                         aaGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTGCcgc           >  2:3404535/1‑100 (MQ=38)
                                                                                          aGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgc           >  2:1523503/1‑99 (MQ=39)
                                                                                          aGTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTATAGCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTCAGTATTTGTTCcgcg          >  2:395523/1‑100 (MQ=38)
                                                                                          aGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgcg          <  2:155816/100‑1 (MQ=255)
                                                                                           gTTGTCCA‑GCGTATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTTCCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgcg          <  1:1952411/99‑1 (MQ=255)
                                                                                           gTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCgcgc         >  1:3253942/1‑100 (MQ=255)
                                                                                           gTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCgcgc         >  2:1512726/1‑100 (MQ=255)
                                                                                           gTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCgcgc         >  2:447578/1‑100 (MQ=255)
                                                                                              gTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTc      >  2:2631201/1‑100 (MQ=255)
                                                                                                 cA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAg   >  1:1058337/1‑100 (MQ=255)
                                                                                                 cA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAg   <  2:3170892/100‑1 (MQ=255)
                                                                                                  a‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGc  >  2:4191400/1‑100 (MQ=255)
                                                                                                  a‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGc  >  2:1478739/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGC  >  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: ‑