Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 170,839 T→G 8.8% intergenic (+20/‑550) PP_0159 → / → PP_16SA family III CoA‑transferase/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947170,8390TG 8.8% 406.1 / 20.8 157intergenic (+20/‑550)PP_0159/PP_16SAfamily III CoA‑transferase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (58/83);  new base G (10/5);  total (68/89)
Fisher's exact test for biased strand distribution p-value = 9.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

TCGGAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGCA  >  NC_002947/170691‑170988
                                                                                                                                                    |                                                                                                                                                     
tCGGAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTc                                                                                                                                                      <  1:3176356/150‑1 (MQ=255)
 cGGAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCa                                                                                                                                                     <  2:2958786/150‑1 (MQ=255)
  ggAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAg                                                                                                                                                    >  2:3899816/1‑150 (MQ=255)
   gagaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGGGGCCGTCAgg                                                                                                                                                   >  1:1514474/1‑150 (MQ=255)
   gagaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGCGGCCGTCAgg                                                                                                                                                   >  1:1514270/1‑150 (MQ=255)
    agaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGc                                                                                                                                                  <  2:3620554/150‑1 (MQ=255)
    agaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGc                                                                                                                                                  <  2:3363666/150‑1 (MQ=255)
    agaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGc                                                                                                                                                  <  2:725784/150‑1 (MQ=255)
     gaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGTGGCCGTCAGGcg                                                                                                                                                 >  2:3920332/1‑149 (MQ=255)
     gaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGcc                                                                                                                                                 >  2:3920333/1‑150 (MQ=255)
        gCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTcc                                                                                                                                              <  2:4314076/150‑1 (MQ=255)
        gCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTcc                                                                                                                                              <  2:2388894/150‑1 (MQ=255)
          cAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTcccc                                                                                                                                            >  2:3788898/1‑150 (MQ=255)
          cAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTcccc                                                                                                                                            >  2:3788856/1‑150 (MQ=255)
            gTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCTGACGGCCTCCCCGc                                                                                                                                          <  2:3475195/150‑1 (MQ=255)
            gTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGc                                                                                                                                          <  2:1786538/150‑1 (MQ=255)
              ggAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGTGGGGAGGCCTGACGGCCTCCCCGCtt                                                                                                                                        <  2:1810195/150‑1 (MQ=255)
              ggAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGTGGGGAGGCCTGACGGCCTCCCCGCtt                                                                                                                                        <  2:1709196/150‑1 (MQ=255)
                  tATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGtt                                                                                                                                    >  2:1085032/1‑150 (MQ=255)
                         gggCGCCACCACTGCTGGGCGAGCATCCCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAgc                                                                                                                             <  2:3238817/150‑1 (MQ=255)
                         gggCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAgc                                                                                                                             <  2:987057/150‑1 (MQ=255)
                             gcCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGGGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTa                                                                                                                         >  1:2663150/1‑150 (MQ=255)
                             gcCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTAAGCGCTa                                                                                                                         >  1:3797605/1‑150 (MQ=255)
                             gcCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGAGGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTa                                                                                                                         <  2:3175455/150‑1 (MQ=255)
                                accacTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCa                                                                                                                      >  2:3185578/1‑150 (MQ=255)
                                accacTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCGTTGTTTTTCCGCGCTATCa                                                                                                                      >  2:4149871/1‑150 (MQ=255)
                                 ccacTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCCGCGCTATCaa                                                                                                                     >  2:578986/1‑150 (MQ=255)
                                  cacgggtgGGCGCGCATACCTCAGTGTGTCTGGGGGTTGTGTTGGGGCTGGTTGGTGGCGTGGTACCGGGTCGGCGTGATGGTGGGGTGCTCTATTTGAAACGGGTGGTGGCCTTCTGGCCTCCCCGCTTTGTTTTTCAGCGCTATCaaa                                                                                                                    <  1:3638144/144‑1 (MQ=255)
                                  cacTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCaaa                                                                                                                    >  1:850794/1‑150 (MQ=255)
                                    ggggggggCGCGCATACGGTGGGGGTGGGGGGGGAGGTGTGGGGGGGGGTTGCGGGGGTGGTACCGGGTCTGGGTGATGGGGGGGTGCTCTAGATGAAAGGGGGGGCGGCCGTCCGGCCCCCCCGCTTTGTTTTTCAGCGCTATCAAAgg                                                                                                                  <  1:2133487/145‑1 (MQ=255)
                                       ctgGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCGGTGTTTTTCAGCGCTATCAAAGGTaa                                                                                                               >  2:4200935/1‑150 (MQ=255)
                                       ctgGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCGCTAGATGAAAAGCGGGGGGGCCGGCCGGCCGCCCCGCGTTGGTTTTCAGCGCTATCAAAGGtca                                                                                                               >  1:2405410/1‑148 (MQ=255)
                                         gggCGAGCATACCGATGTGGTGCTGGGGGATGTGTTGGGGGTGGTTGCTGGCGAGTTACTGCGTCTGCGTGATGGTGGGGTGCTCTAGATTAAAAGGGGGGAGGCGGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAAtt                                                                                                             <  1:2874220/150‑1 (MQ=255)
                                         gggCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCGTTGTTTTTCAGCGCTATCAAAGGTAAtt                                                                                                             >  2:1828846/1‑150 (MQ=255)
                                             gAGCATACCGAAGTGGTGCTGGGGGAATTGTTGGGGCTGGATGCTGGCGAGGTAGAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGaaa                                                                                                         <  2:4119585/150‑1 (MQ=255)
                                             gAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGaaa                                                                                                         <  2:3740431/150‑1 (MQ=255)
                                             gAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTGTGTTTTTCAGCGCTATCAAAGGTAATTGaaa                                                                                                         >  2:2426192/1‑150 (MQ=255)
                                              aGCATACCGACGGGGTGCGGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGAGGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAAt                                                                                                        <  2:2466908/150‑1 (MQ=255)
                                                  tACCGATGTGGTGCTGGGGGATGTGTGGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTCGATGAAACGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTaaa                                                                                                    <  2:1618603/150‑1 (MQ=255)
                                                  tACCGAAGTGGTGGTGGGGGATGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTaaa                                                                                                    <  2:62946/150‑1 (MQ=255)
                                                  tACCGAAGTGGTGCTGGGGGATGTGTTGGGGGTGGATGCTGGCGAGTTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGCCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTaaa                                                                                                    <  2:3781231/150‑1 (MQ=255)
                                                  tACCGAAGTGGTGCTGGGGGATGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTaaa                                                                                                    <  2:3216016/150‑1 (MQ=255)
                                                  tACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGGGTGATGCTGGGGTGCTCTTGATGAAAAGCGGGGCGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTaaa                                                                                                    <  1:4425479/150‑1 (MQ=255)
                                                   agcGAAGTGGTGCTGGGGGTTGTGTTGGGGGTGGATGCTGGCGAGGTACTGCGTCTGCGTGTTGCTGGGGTGCTCTTGATGAAAAGCGGGGAGGCCGTCAGGGCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAg                                                                                                   <  1:3346747/148‑1 (MQ=255)
                                                   aCCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAg                                                                                                   <  2:2945720/150‑1 (MQ=255)
                                                   aCCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCGCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAg                                                                                                   >  1:2895570/1‑150 (MQ=255)
                                                       aaGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAATGTAATTTAAATTAAAGGttt                                                                                               >  2:894417/1‑149 (MQ=255)
                                                         gtggtgCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGAc                                                                                             <  2:3557161/150‑1 (MQ=255)
                                                                 gggggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCtctc                                                                                     <  2:1428024/150‑1 (MQ=255)
                                                                  ggggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCg                                                                                    >  2:2039743/1‑150 (MQ=255)
                                                                     gAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAAt                                                                                 >  2:832612/1‑150 (MQ=255)
                                                                        gtgtTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATccc                                                                              <  2:2560861/150‑1 (MQ=255)
                                                                         tgtTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATcccc                                                                             <  2:3294293/150‑1 (MQ=255)
                                                                          gtTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCGATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCt                                                                            >  1:2607393/1‑150 (MQ=255)
                                                                           ttGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGGGGCCGTCACGCCTCCCCGCGTTGTTTTTAAGCGCTATCAAATGTAAGTGTAATTAAATGTTGTCGCCCACTCGAATacccgg                                                                           >  1:122850/1‑148 (MQ=255)
                                                                           ttGGGGCTGGATGCCGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCtt                                                                           <  2:2580320/150‑1 (MQ=255)
                                                                                 cTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCCGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCGTATAATg                                                                     >  1:2212300/1‑150 (MQ=255)
                                                                                  tGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGGGGCCGTAAGGCCACCCCGCGTTGTTTTTCAGCCCGTAAAAAAGTAAATGAAATTAAATGTTTAAGCCCACACCAAAACCCATATaaata                                                                    >  1:848713/1‑147 (MQ=255)
                                                                                    gATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATgcgc                                                                  <  1:4516090/150‑1 (MQ=255)
                                                                                       gCTGGCGTGGTACAGGGTCTGCGTGATGGTGGGGTGCTCTAGATGAAACGCGGGGAGGCCGTCTGGCCCCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGcccc                                                               <  1:2107468/150‑1 (MQ=38)
                                                                                         tGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGCGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATAAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATACCCGTATAATGCGCCCCAc                                                             >  1:3583204/1‑150 (MQ=255)
                                                                                         tGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTGATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCAc                                                             >  2:4496229/1‑150 (MQ=255)
                                                                                         tGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCAc                                                             >  2:2705079/1‑150 (MQ=255)
                                                                                         gggCGTGGTACCGGGTGTGGGTGTTGGGGGGGTGCTCTAGTTGAACCGGGGGGTGGCGGTCCGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCAc                                                             <  1:2054718/149‑1 (MQ=21)
                                                                                              tggTACTGCGTCTGCGTTTTGCGGGGGTGCTCTCGTTGAAAAGGGGGGCGGCCGTCAGGCGTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCa                                                        <  1:2125244/149‑1 (MQ=25)
                                                                                                gTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGc                                                      >  2:1460743/1‑150 (MQ=255)
                                                                                                gTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGc                                                      <  2:2677657/150‑1 (MQ=255)
                                                                                                      cGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATc                                                <  2:208221/150‑1 (MQ=255)
                                                                                                      cGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATc                                                <  2:2465980/150‑1 (MQ=255)
                                                                                                          tGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGaa                                            <  1:3223825/150‑1 (MQ=255)
                                                                                                           gCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTGAGGCCTCCCCCCGTTGTTTTTCATCGCTATCAAAAGTAAATTAAATTAAAGGTTTAAGCCCAAAAGAATACCCGTAAAATGCGCCCCACATACAGCGAAAAAGGAAc                                           >  1:3233691/1‑150 (MQ=255)
                                                                                                                 tGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGGCTCCCCGGTTTGTTTTTGAGGGCTATCAAAGGTAATTGAAATTAAATGTTGTCGCCCTCTCGAATACCCTTATAATGCGCCCCACTTCCAGCGAAATAGGAAAGaaaaa                                     >  1:4519190/1‑150 (MQ=25)
                                                                                                                 tGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGaaaaa                                     <  1:3711342/150‑1 (MQ=255)
                                                                                                                 tGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAACTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGaaaaa                                     <  1:3711330/150‑1 (MQ=255)
                                                                                                                     ggggTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTcc                                 <  1:4450522/150‑1 (MQ=255)
                                                                                                                     ggggTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTcc                                 <  1:1460743/150‑1 (MQ=255)
                                                                                                                      gggTGCTCTAGTTGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCt                                <  2:1431289/150‑1 (MQ=255)
                                                                                                                       gttGCTATATATTAAATGCGGGGAGGCGGTCTGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCtt                               <  1:1457367/148‑1 (MQ=21)
                                                                                                                           ctctAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTgaga                           <  2:2234415/150‑1 (MQ=255)
                                                                                                                            tctAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAATGTAAATTAAATTAAATGTTGAAGCCCTCTCGAATACCCATATAATTCGCCCCACATCCAGCGACATCGGAACGAAAAACTCCATGAGAt                          >  1:2585134/1‑150 (MQ=25)
                                                                                                                            tctAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTAAGCGCTATAAAAGGTAATTGAAATTAAAGGTTGAAGCCCAAACGAATCCCCTTATAATGCGCCCCACTTACAGCGACAACGGAACGAAAAACTCCATGAGAt                          >  1:1861683/1‑150 (MQ=25)
                                                                                                                              ttgATTAAACGGGGGGAGGCCGTCAGGCCGCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTc                        <  1:3788898/148‑1 (MQ=255)
                                                                                                                              tAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATTCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTc                        >  2:4411571/1‑150 (MQ=255)
                                                                                                                                gATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCaa                      <  2:1241922/150‑1 (MQ=255)
                                                                                                                                gATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCaa                      <  1:2469026/150‑1 (MQ=255)
                                                                                                                                aaaacaaaGCGGGGAGGGCTGACGGCCTCCCCGCTTTGTTTTTCAGGGCTATGAAATGTAATTTTAATTAAAGGTTTACGCCCTCTCGAATACCCTTATAATGGGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCaa                      >  2:864066/6‑150 (MQ=17)
                                                                                                                                   gAAAAGCGGGGCGGCCGTCAGGCCGCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGa                   <  1:4233259/150‑1 (MQ=255)
                                                                                                                                   gAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGa                   <  2:3660709/150‑1 (MQ=255)
                                                                                                                                     aaaGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGt                 <  1:1085032/150‑1 (MQ=255)
                                                                                                                                        gCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCACTGAGTTaa              >  2:3841852/1‑150 (MQ=255)
                                                                                                                                         cGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAag             <  2:2732999/150‑1 (MQ=255)
                                                                                                                                         cGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCAATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAag             >  1:571441/1‑150 (MQ=255)
                                                                                                                                          ggggAGGCCTGACGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATCATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAagt            >  2:4505319/1‑150 (MQ=18)
                                                                                                                                          ggggAGGCCTGACGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAagt            >  2:3887767/1‑150 (MQ=255)
                                                                                                                                            ggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAagtag          >  1:2485066/1‑150 (MQ=255)
                                                                                                                                              aGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGtt        <  2:95115/150‑1 (MQ=255)
                                                                                                                                                gCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTaa      <  1:3185578/150‑1 (MQ=255)
                                                                                                                                                  cGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAgg    >  1:2390138/1‑150 (MQ=255)
                                                                                                                                                  cGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAgg    <  2:1514474/150‑1 (MQ=255)
                                                                                                                                                  cGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAgg    <  2:1514270/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCGGAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGCA  >  NC_002947/170691‑170988

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