Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 972,118 T→A 6.7% L136Q (CTG→CAG)  secD → protein translocase subunit SecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947972,1180TA 6.7% 409.3 / 17.0 134L136Q (CTG→CAG) secDprotein translocase subunit SecD
Reads supporting (aligned to +/- strand):  ref base T (65/59);  new base A (4/5);  total (69/65)
Fisher's exact test for biased strand distribution p-value = 7.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 146 total aligned reads displayed.

GCTGCCAGCCAAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGC  >  NC_002947/971969‑972265
                                                                                                                                                     |                                                                                                                                                   
gCTGCCAGCCAAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCct                                                                                                                                                     >  2:3058150/1‑150 (MQ=255)
  tGCCAGCCAAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCctgc                                                                                                                                                   >  1:2053426/1‑150 (MQ=255)
  tGCCAGCCAAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCctgc                                                                                                                                                   <  2:1393339/150‑1 (MQ=255)
  tGCCAGCCAAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCctgc                                                                                                                                                   <  1:2548964/150‑1 (MQ=255)
  tGCCAGCCAAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCctgc                                                                                                                                                   <  2:597184/150‑1 (MQ=255)
   gccagccaAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCctgct                                                                                                                                                  <  1:1017428/150‑1 (MQ=255)
    ccagccaAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCctgctg                                                                                                                                                 <  1:1872363/150‑1 (MQ=255)
       gccaAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGaa                                                                                                                                              <  1:2187875/150‑1 (MQ=255)
       gccaAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGaa                                                                                                                                              >  2:2159082/1‑150 (MQ=255)
         caAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGt                                                                                                                                            >  2:1693211/1‑150 (MQ=255)
          aaGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCctgc                                                                                                                                                   <  2:1001047/142‑1 (MQ=255)
          aaGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCctgc                                                                                                                                                   >  1:1001047/1‑142 (MQ=255)
          aaGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtg                                                                                                                                           >  2:2453119/1‑150 (MQ=255)
           aGGATGTAGTGCGAAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtgg                                                                                                                                          <  2:21216/150‑1 (MQ=255)
            ggATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGGCCAGACTACCCCGCACTGGCTGCGCACCCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGAGGTGTGCCCTTCCTGCTGCaatagaa                                                                                                                                         >  2:2322438/1‑145 (MQ=255)
            ggATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTCCACTTCCAGCAGGAAGtgca                                                                                                                                         <  1:364150/150‑3 (MQ=255)
             gATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCAGCAGGAAGtgcac                                                                                                                                        <  1:392781/150‑1 (MQ=255)
                  aGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggacatgga                                                                                                                                   <  2:764435/150‑1 (MQ=255)
                     gcgcAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACaa                                                                                                                                <  1:1104251/150‑1 (MQ=255)
                        cAAGGCACTGGGCGATGGTTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCAGCAGGAAGtgcaca                                                                                                                                       >  2:1576575/1‑140 (MQ=255)
                        cAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCAGCAGGAAGtgcaca                                                                                                                                       <  1:1576575/140‑1 (MQ=255)
                          aGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCa                                                                                                                           >  1:1778059/1‑150 (MQ=255)
                          aGGCACTGGGCGATGATTACGTCGTAGCGCCGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCa                                                                                                                           <  2:2688893/150‑1 (MQ=255)
                           ggCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCAt                                                                                                                          <  1:378799/150‑1 (MQ=255)
                           ggCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCAt                                                                                                                          >  2:1127355/1‑150 (MQ=255)
                            gCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATg                                                                                                                         <  2:1825851/150‑1 (MQ=255)
                             cACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGTCATGa                                                                                                                        <  2:2279248/150‑1 (MQ=255)
                                    cGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCctgc                                                                                                                                                   >  1:1286304/1‑116 (MQ=255)
                                    cGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCg                                                                                                                 >  1:2257204/1‑150 (MQ=255)
                                       tgatTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCAt                                                                                                              <  2:2053426/150‑1 (MQ=255)
                                           tACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGCCATGGACAAGGCCATGACTGCCCGCATGaaa                                                                                                          <  2:1892749/150‑1 (MQ=255)
                                                 gtGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTAc                                                                                                    <  2:1759058/150‑1 (MQ=255)
                                                      ccTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGTCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggacat                                                                                                                                      >  2:1799352/1‑111 (MQ=255)
                                                      ccTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGGAGTGGACATGGACAAGGCCATGACTGCCCGCATGaaa                                                                                                          >  1:536284/1‑139 (MQ=255)
                                                      ccTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTc                                                                                                       <  1:2317931/142‑1 (MQ=255)
                                                       cTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGc                                                                                              >  2:2016/1‑150 (MQ=255)
                                                       cTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGc                                                                                              >  2:2878022/1‑150 (MQ=255)
                                                        tGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggacatggac                                                                                                                                  >  2:789552/1‑113 (MQ=255)
                                                         gAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTACGGCGGGAATTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGa                                                                                            <  1:1911253/150‑1 (MQ=255)
                                                            ccTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggacat                                                                                                                                      <  1:3012040/105‑1 (MQ=255)
                                                            ccTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggacat                                                                                                                                      >  2:3012040/1‑105 (MQ=255)
                                                              tGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCa                                                                                       >  1:1719181/1‑150 (MQ=255)
                                                                gCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGt                                                                                         <  1:1750311/146‑1 (MQ=255)
                                                                  ccAGACTACCCCGCAATGGCTGCGCCACCTGGGCGCAAGCCCGCTGAAACTGGGCCTGGACCTGTCCGGTGGTGGGCACTTCATGCTGGAAGTGGACATGGACAAAGCCATGACTGCCCGCATGAAAGTCTACGAAATCGAGGTCAAGAc                                                                                   >  2:941057/1‑150 (MQ=255)
                                                                  ccAGACTACCCCGCAATGGCTGCGCAAGCTGGGCGCAAGCCCGATGAAACTGGGCCGGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGAc                                                                                   >  2:468469/1‑150 (MQ=255)
                                                                  ccAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGa                                                                                    <  1:732637/149‑1 (MQ=255)
                                                                  ccAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGAc                                                                                   >  1:1016256/1‑150 (MQ=255)
                                                                  ccAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGAc                                                                                   <  1:1240687/150‑1 (MQ=255)
                                                                        tACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTtgct                                                                             >  2:2755670/1‑150 (MQ=255)
                                                                           cccGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTgcg                                                                          >  2:3016749/1‑150 (MQ=255)
                                                                            ccGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTgcgc                                                                         >  2:2435541/1‑150 (MQ=255)
                                                                             cGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACg                                                                                                   <  1:1198670/123‑1 (MQ=255)
                                                                             cGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACg                                                                                                   >  2:1198670/1‑123 (MQ=255)
                                                                                aaTGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAg                                                                                           <  2:960111/128‑1 (MQ=255)
                                                                                  tGGCTGCGCAACCTGGGCGAAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAagag                                                                   <  1:1693211/150‑1 (MQ=255)
                                                                                    gCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTtg                                                                               <  2:2292280/136‑1 (MQ=255)
                                                                                     cTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAgcgc                                                                <  1:2414292/150‑1 (MQ=255)
                                                                                        cgcAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTc                                                             >  1:1475624/1‑150 (MQ=255)
                                                                                         gcAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTcc                                                            <  1:2755670/150‑1 (MQ=255)
                                                                                         gcAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTCCACTTCCAGCAGGAAGtgca                                                                                                                                         >  2:389882/1‑71 (MQ=37)
                                                                                           aaCCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGc                                                          <  1:1970687/150‑1 (MQ=255)
                                                                                             ccTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTa                                                        >  1:2591005/1‑150 (MQ=255)
                                                                                             ccTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTa                                                        >  1:2845170/1‑150 (MQ=255)
                                                                                                      aaGCCCGATGAAACTGGGCCTGGACCCGTGCGGTGGCGGGCACTTCCTGCTGGAAGCCCAAATGGCCAAGGCCATAACTGACCGCATGAAAGTCTACGAAGGCGACGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACAGCAgcct                                               >  2:500651/1‑150 (MQ=255)
                                                                                                        gCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAgcctgc                                             >  1:2711842/1‑150 (MQ=255)
                                                                                                        gCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAgcctgc                                             <  2:450980/150‑1 (MQ=255)
                                                                                                         cccGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAgcctgcc                                            >  1:866865/1‑150 (MQ=255)
                                                                                                          ccGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAgcctgcct                                           >  1:296181/1‑150 (MQ=255)
                                                                                                            gATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTca                                         <  1:2865955/150‑1 (MQ=255)
                                                                                                            gATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTca                                         >  2:25034/1‑150 (MQ=255)
                                                                                                            gATGAAACTGGGCCTGGACCTGTCCGGTGGTGTCCACTTCCAGCAGGAAGtg                                                                                                                                           >  1:15358/1‑52 (MQ=25)
                                                                                                            gATGAAACTGGGCCTGGACCTGTCCGGTGGTGTCCACTTCCAGCAGGAAGtg                                                                                                                                           <  2:15358/52‑1 (MQ=25)
                                                                                                               gAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTcagca                                      <  1:292766/150‑1 (MQ=255)
                                                                                                                    tGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATg                                 >  2:1678814/1‑150 (MQ=255)
                                                                                                                    tGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATg                                 >  2:1255567/1‑150 (MQ=255)
                                                                                                                      ggCCTGGACCTGTCCGGTGGTGTGAACTTCCTGCTGGAAGTGGACATGGACAAGAACCCGGCGGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATggc                               <  1:1571834/150‑1 (MQ=255)
                                                                                                                        ccTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggac                                                                                                                                        <  1:451409/43‑1 (MQ=255)
                                                                                                                        ccTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATggcgg                             >  1:926357/1‑150 (MQ=255)
                                                                                                                         cTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGt                                                                                         <  2:2326604/89‑1 (MQ=255)
                                                                                                                         cTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGt                                                                                         >  1:2326604/1‑89 (MQ=255)
                                                                                                                          tGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGca                           <  1:3016749/150‑1 (MQ=255)
                                                                                                                               cTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGaaa                                                                                                          >  2:1255248/1‑66 (MQ=255)
                                                                                                                               cTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGaaa                                                                                                          <  1:1255248/66‑1 (MQ=255)
                                                                                                                               cTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAgc                                                                  <  2:698698/106‑1 (MQ=255)
                                                                                                                               cTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAgc                                                                  >  1:698698/1‑106 (MQ=255)
                                                                                                                               cTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATg                      >  2:2901699/1‑150 (MQ=255)
                                                                                                                                  tCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTg                   >  2:171316/1‑150 (MQ=255)
                                                                                                                                      gtggtgTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGa                                                                                                            >  1:1085849/1‑57 (MQ=255)
                                                                                                                                        ggtgTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTc             <  2:1719181/150‑1 (MQ=255)
                                                                                                                                         gtgtGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCg            >  1:293196/1‑150 (MQ=255)
                                                                                                                                           gtgCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCt          >  1:626456/1‑150 (MQ=255)
                                                                                                                                           gtccACTTCCAGCAGGAAGTGGACATGGACAAGTCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCt          >  2:2940034/4‑150 (MQ=255)
                                                                                                                                            tgCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATggcg                              >  1:98531/1‑129 (MQ=255)
                                                                                                                                             gCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTg                   >  1:1557971/1‑139 (MQ=255)
                                                                                                                                             gCACTTCCTGCTGGAAGTGGACATGGACAAGGCAAGGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTACTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGa        <  1:2322438/150‑1 (MQ=255)
                                                                                                                                              cACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCCTGCTGGGCTTCGCTGAc       <  2:293196/150‑1 (MQ=255)
                                                                                                                                                cTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTggg                 <  2:2444054/138‑1 (MQ=255)
                                                                                                                                                cTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTggg                 >  1:2444054/1‑138 (MQ=255)
                                                                                                                                                 ttCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGa        <  2:2394989/146‑1 (MQ=255)
                                                                                                                                                   ccTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGc  <  2:1016256/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCTGCCAGCCAAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGC  >  NC_002947/971969‑972265

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