Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,877,026 T→A 5.5% D211V (GAC→GTC)  ptsP ← phosphoenolpyruvate‑dependent regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,877,0260TA 5.5% 931.9 / 22.6 273D211V (GAC→GTC) ptsPphosphoenolpyruvate‑dependent regulator
Reads supporting (aligned to +/- strand):  ref base T (128/127);  new base A (9/7);  total (139/134)
Fisher's exact test for biased strand distribution p-value = 7.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 273 total aligned reads displayed.

GGCTGCGTCTTCAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAGGCCACGGATCGA  >  NC_002947/5876878‑5877174
                                                                                                                                                    |                                                                                                                                                    
  cTGCGTCTTCAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCAt                                                                                                                                                   <  1:712470/150‑1 (MQ=255)
  cTGCGTCTTCAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCAt                                                                                                                                                   <  2:4267184/150‑1 (MQ=255)
  cTGCGTCTTCAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGCTCGGCGTCGATGTCAt                                                                                                                                                   >  2:3547456/1‑150 (MQ=255)
    gCGTCTTCAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCg                                                                                                                                                 >  2:3911831/1‑150 (MQ=255)
     cGTCTTCAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGa                                                                                                                                                <  2:1346868/150‑1 (MQ=255)
        cTTCAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACgg                                                                                                                                             <  1:1177352/150‑1 (MQ=255)
        cTTCAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGACATCGACgc                                                                                                                                             <  1:1407773/150‑2 (MQ=255)
        cTTCAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGACATCGACgc                                                                                                                                             <  1:1675010/150‑2 (MQ=255)
         ttCAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTGTGGAACAGCTTGAGTTCGGCGTCGATGCCATCGACGGt                                                                                                                                            >  2:1758850/1‑150 (MQ=255)
         ttCAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGt                                                                                                                                            >  1:5247876/1‑150 (MQ=255)
         ttCAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGt                                                                                                                                            >  2:4171442/1‑150 (MQ=255)
         ttCAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGACATCGACgcc                                                                                                                                            <  1:1677025/150‑3 (MQ=255)
           cAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCt                                                                                                                                          >  2:2207818/1‑150 (MQ=255)
           cAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCt                                                                                                                                          >  1:871316/1‑150 (MQ=255)
            aaGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCtt                                                                                                                                         >  1:4525271/1‑150 (MQ=255)
            aaGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCtt                                                                                                                                         >  2:1201857/1‑150 (MQ=255)
            aaGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGACATCGACgccgaa                                                                                                                                         <  2:5462411/150‑6 (MQ=255)
              gCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGt                                                                                                                                       >  2:5147419/1‑150 (MQ=255)
                atcatcAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTcc                                                                                                                                     <  2:802213/150‑1 (MQ=255)
                   atcaGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCggg                                                                                                                                  <  2:5238543/150‑1 (MQ=255)
                    tcaGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGa                                                                                                                                 <  2:1480755/150‑1 (MQ=255)
                    tcaGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGa                                                                                                                                 <  2:2471681/150‑1 (MQ=255)
                    tcaGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGa                                                                                                                                 >  2:1117851/1‑150 (MQ=255)
                      aGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGAcc                                                                                                                               >  2:5011190/1‑150 (MQ=255)
                        aTAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCAc                                                                                                                             >  1:3757926/1‑150 (MQ=255)
                           gACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTc                                                                                                                          >  2:3450782/1‑150 (MQ=255)
                             cGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCa                                                                                                                        <  2:1965748/150‑1 (MQ=255)
                             cGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCa                                                                                                                        <  1:756396/150‑1 (MQ=255)
                               tCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAgg                                                                                                                      >  1:3409301/1‑150 (MQ=255)
                               tCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAgg                                                                                                                      >  1:3727460/1‑150 (MQ=255)
                                  aaCAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTcg                                                                                                                   >  2:2888583/1‑150 (MQ=255)
                                  aaCAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTcg                                                                                                                   <  1:2475333/150‑1 (MQ=255)
                                  aaCAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTcg                                                                                                                   <  1:4097554/150‑1 (MQ=255)
                                        gcACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTcggccggc                                                                                                             >  1:544557/1‑150 (MQ=255)
                                            gCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCa                                                                                                         >  2:1917497/1‑150 (MQ=255)
                                              tctTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAAc                                                                                                       >  1:2753789/1‑150 (MQ=255)
                                                    gggCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATcacc                                                                                                 <  1:4019233/150‑1 (MQ=255)
                                                                 tGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTcc                                                                                                                         >  1:4863490/1‑113 (MQ=255)
                                                                 tGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCa                                                                                    >  1:3327407/1‑150 (MQ=255)
                                                                  ggCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCAc                                                                                   <  1:933646/150‑1 (MQ=255)
                                                                    ccAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACcg                                                                                 >  1:587956/1‑150 (MQ=255)
                                                                     cAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACcgc                                                                                <  2:3971310/150‑1 (MQ=255)
                                                                     cAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACcgc                                                                                <  1:3079248/150‑1 (MQ=255)
                                                                     cAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACcgc                                                                                >  1:1316305/1‑150 (MQ=255)
                                                                       gTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCgg                                                                              <  1:4948225/150‑1 (MQ=255)
                                                                          tGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCAc                                                                           <  2:71530/150‑1 (MQ=255)
                                                                            gCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTcggcc                                                                                                                >  2:4210593/1‑111 (MQ=255)
                                                                            gCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTCGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCa                                                                         <  1:4879849/150‑1 (MQ=255)
                                                                             ccGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAg                                                                        <  1:1808801/150‑1 (MQ=255)
                                                                              cGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAgg                                                                       <  2:5312444/150‑1 (MQ=255)
                                                                               gACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTCCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGc                                                                      <  1:2630162/150‑1 (MQ=255)
                                                                                    cTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAAcc                                                                 >  1:5393094/1‑150 (MQ=255)
                                                                                      tGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAAcccn                                                               <  1:1088333/150‑2 (MQ=255)
                                                                                       gcgcATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCgg                                                              <  1:5602180/150‑1 (MQ=255)
                                                                                       gcgcATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCgg                                                              >  1:5247537/1‑150 (MQ=255)
                                                                                       gcgcATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCgg                                                              <  2:1488575/150‑1 (MQ=255)
                                                                                       gcgcATGTCGTCGCGCACACCGTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCgg                                                              <  2:4314772/150‑1 (MQ=255)
                                                                                        cgcATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGc                                                             <  2:2079546/150‑1 (MQ=255)
                                                                                          cATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGcac                                                           >  2:3579228/1‑150 (MQ=255)
                                                                                           aTGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGcaca                                                          >  1:2714572/1‑150 (MQ=255)
                                                                                              tcgtcgCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACAccc                                                       <  1:5064081/150‑1 (MQ=255)
                                                                                                 tcgCGCACACCTTCCAGGGCGTTCTGGAACAGCGTGAGTTCGGCGTCGATGGCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACg                                                    >  2:903871/1‑150 (MQ=255)
                                                                                                   gcgcACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGaa                                                  <  1:1901622/150‑1 (MQ=255)
                                                                                                   gcgcACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGaa                                                  >  2:1307394/1‑150 (MQ=255)
                                                                                                    cgcACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAg                                                 <  2:5422530/150‑1 (MQ=255)
                                                                                                    cgcACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAg                                                 >  2:3406739/1‑150 (MQ=255)
                                                                                                     gcACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAgc                                                >  1:1457380/1‑150 (MQ=255)
                                                                                                       acacCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAgcgc                                              >  1:2238205/1‑150 (MQ=255)
                                                                                                       acacCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAgcgc                                              >  1:3231829/1‑150 (MQ=255)
                                                                                                        cacCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAgcgcg                                             >  2:4403908/1‑150 (MQ=255)
                                                                                                           cTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCtt                                          <  1:3896433/150‑1 (MQ=255)
                                                                                                              ccAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCt                                       >  1:4217750/1‑150 (MQ=255)
                                                                                                              ccAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCt                                       >  1:1565922/1‑150 (MQ=255)
                                                                                                              ccAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCt                                       <  1:4752770/150‑1 (MQ=255)
                                                                                                               cAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTg                                      <  1:5323575/150‑1 (MQ=255)
                                                                                                               cAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTg                                      <  2:3308016/150‑1 (MQ=255)
                                                                                                                 gggCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCAGGAAGCGCGCTTCGTggg                                    <  1:5121141/150‑3 (MQ=255)
                                                                                                                 gggCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGaa                                    >  2:812188/1‑150 (MQ=255)
                                                                                                                   gCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATa                                  >  2:4451381/1‑150 (MQ=255)
                                                                                                                      ttCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCt                               >  1:2484689/1‑150 (MQ=255)
                                                                                                                        cTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCttt                             <  1:2298951/150‑1 (MQ=255)
                                                                                                                        cTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCttt                             <  1:4584639/150‑1 (MQ=255)
                                                                                                                        cTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCttt                             <  2:536290/150‑1 (MQ=255)
                                                                                                                         tGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTg                            >  1:3520027/1‑150 (MQ=255)
                                                                                                                         tGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTg                            <  1:3592271/150‑1 (MQ=255)
                                                                                                                           gAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGcc                          <  2:3791789/150‑1 (MQ=255)
                                                                                                                           gAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGcc                          <  2:860739/150‑1 (MQ=255)
                                                                                                                             aCAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCt                        <  1:4304477/150‑1 (MQ=255)
                                                                                                                                gCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGAc                     <  2:3693050/150‑1 (MQ=255)
                                                                                                                                 cTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCGTCCTGAATACCTTTGCCCTGACg                    <  1:2207818/150‑1 (MQ=255)
                                                                                                                                      gTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCa               <  1:2886670/150‑1 (MQ=255)
                                                                                                                                       ttCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAg              <  1:4180159/150‑1 (MQ=255)
                                                                                                                                       gtCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAg              <  2:4630841/149‑1 (MQ=255)
                                                                                                                                          accggcGATGACATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAGGcc           >  2:2427160/6‑150 (MQ=255)
                                                                                                                                          accGTCGATGACATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAGGcc           >  1:3995158/3‑150 (MQ=255)
                                                                                                                                           gCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAGGCCa          >  2:2764961/1‑150 (MQ=255)
                                                                                                                                           gCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAGGCCa          >  2:5536497/1‑150 (MQ=255)
                                                                                                                                               cGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAGGCCACGGa      <  1:4783040/150‑1 (MQ=255)
                                                                                                                                                 aTGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAGGCCACGGATc    <  1:1000912/150‑1 (MQ=255)
                                                                                                                                                   gTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAGGCCACGGATCGa  >  2:4792263/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCTGCGTCTTCAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAGGCCACGGATCGA  >  NC_002947/5876878‑5877174

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