Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 170,839 T→G 7.9% 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 7.9% 535.7 / 21.8 191intergenic (+20/‑550)PP_0159/PP_16SAfamily III CoA‑transferase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (78/91);  new base G (10/7);  total (88/103)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 191 total aligned reads displayed.

CGGAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGC  >  NC_002947/170692‑170987
                                                                                                                                                   |                                                                                                                                                    
cGGAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCa                                                                                                                                                    <  2:4473277/150‑1 (MQ=255)
 ggAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAg                                                                                                                                                   >  1:1224760/1‑150 (MQ=255)
  gagaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAgg                                                                                                                                                  >  1:957967/1‑150 (MQ=255)
   agaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGc                                                                                                                                                 <  2:1399943/150‑1 (MQ=255)
   agaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGc                                                                                                                                                 <  2:2214509/150‑1 (MQ=255)
   agaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGc                                                                                                                                                 <  2:5135806/150‑1 (MQ=255)
   agaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGc                                                                                                                                                 <  2:5856414/150‑1 (MQ=255)
    gaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGcc                                                                                                                                                >  2:3147624/1‑150 (MQ=255)
       gCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGCGGCCGGCAGgccgcc                                                                                                                                             >  1:1579484/1‑150 (MQ=255)
        ccAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCTGACGGCCTccc                                                                                                                                            <  2:1841498/150‑1 (MQ=255)
        ccAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTccc                                                                                                                                            <  1:2135071/150‑1 (MQ=255)
        ccAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTccc                                                                                                                                            <  2:5062460/150‑1 (MQ=255)
         cAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTcccc                                                                                                                                           <  1:3980572/150‑1 (MQ=255)
         cAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTcccc                                                                                                                                           <  1:7301511/150‑1 (MQ=255)
         cAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTcccc                                                                                                                                           <  2:1145729/150‑1 (MQ=255)
          aGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCTGACGGCCTCCCCg                                                                                                                                          >  1:1058902/1‑150 (MQ=255)
                gTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGt                                                                                                                                    >  1:1797798/1‑150 (MQ=255)
                   tCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCGTTGtttt                                                                                                                                 >  2:6614493/1‑150 (MQ=255)
                    cGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGttttt                                                                                                                                <  2:4857443/150‑1 (MQ=255)
                    cGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGttttt                                                                                                                                >  2:2963254/1‑150 (MQ=255)
                     gCCGGGCGCCACCAGTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCTGACGGCCTCCCCGCttttcatcta                                                                                                                               <  2:4818972/150‑8 (MQ=38)
                      ccGGGCGCCACCACTGCTGGGCGAGCATACCGACGGGGTGGGGGGGGACGTGTTGGGGCTGGCGGCTGGCGAGGTACAGCGTCTGCGTGTTGCTGGGGTGCTCTATATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCa                                                                                                                              <  2:5584733/150‑1 (MQ=38)
                         ggCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAgcg                                                                                                                           <  2:2406598/150‑1 (MQ=255)
                             ccaccaCTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCGACCCCCGTTGTTTTTAAGCGCTAt                                                                                                                       >  1:5727745/1‑150 (MQ=255)
                                ccacTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCGCCCCGCGTTGTTTTTCAGCCCTATCaa                                                                                                                    >  1:6855341/1‑150 (MQ=255)
                                 cacTTCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCTTGATGCTGGGGTGCTATAGATGAAAAGCGGGGTGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCaaa                                                                                                                   <  2:4569905/150‑1 (MQ=255)
                                    tgctgGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAg                                                                                                                                                   >  2:109603/1‑115 (MQ=255)
                                    tgctgGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGt                                                                                                                >  1:6533334/1‑150 (MQ=255)
                                    tgctgGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGt                                                                                                                >  2:4313416/1‑150 (MQ=255)
                                    tgctgGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCGTTGTTTTTCAGCGCTATCAAAGGt                                                                                                                >  2:559363/1‑150 (MQ=255)
                                     gctgGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTa                                                                                                               >  1:7177537/1‑150 (MQ=255)
                                      ctgGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTTCAGCGTGTGGGTGATGCTGGGGTGCTCTAGTTTAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTaa                                                                                                              <  2:3031553/150‑1 (MQ=39)
                                        gggCGAGCATACCGATGTGGTGCTGGGGGATGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAAtt                                                                                                            <  2:5437445/150‑1 (MQ=255)
                                         ggCGTGAATACGGAAGTGGTGGTGGGGGATGTGTTGGGGCTGGATGCTGGCGCGGTACAGCGTCTGCGTTATGGTGGGGTGCTAATGATCAAAAGCGGGGCGGCCGTCCGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTg                                                                                                           <  1:5294103/150‑1 (MQ=255)
                                          gCGAGCACACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGa                                                                                                          >  1:5969292/1‑150 (MQ=255)
                                            gAGCATACCGATGTGGTGGTGGGGTTTGTGTTGGGGCGGGTTGCTGGCGAGGTACTGCGTCTGCGTGAGGGTGGGGTGATCTATATCAAAAGCGGGGCGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGaaa                                                                                                        <  1:1075013/150‑1 (MQ=255)
                                            gAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGaaa                                                                                                        <  2:5786568/150‑1 (MQ=255)
                                            gAGCACCCCGACGTGGTGCTGGCGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGTTGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGaaa                                                                                                        <  2:5696493/150‑1 (MQ=39)
                                              gCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGTGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAAtt                                                                                                      >  1:2454192/1‑150 (MQ=255)
                                                 tACCGTTGTGGTGGGGGGGGTTGTTTTGGGGCTGGATGCTGGGGAGGTACAGCTTCGGGTTGTTGCGGGGGTGCTATATATTAAATGGGGGGTGGCCGTCTGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTaaa                                                                                                   <  1:2160120/150‑1 (MQ=255)
                                                 tACCGATGTGGTGCTGGGGTATGTGTTGGGGGTGGTTGGTGGCGTGGTACTGCGTCTGGGTTTTGGTGGGTTGCTCTATATCAAAAGGGGGGTGGCCGCCCGCCCCCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTaaa                                                                                                   <  1:1499005/150‑1 (MQ=255)
                                                 tACCGATGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAATGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTaaa                                                                                                   <  1:6674362/150‑1 (MQ=255)
                                                 tACCGAAGTGGTGCTGGGGGATGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTaaa                                                                                                   <  2:2317911/150‑1 (MQ=255)
                                                 tACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTaaa                                                                                                   <  2:6491772/150‑1 (MQ=255)
                                                    cGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGt                                                                                                >  2:4805575/1‑150 (MQ=255)
                                                      ttgtggtgCTGGGGTAAGTGTTGGGGCTGGATGCTGGCGTGGTACTGCGTCTGCGTGATGCTGGGGTGCTCTATATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTg                                                                                              <  1:6629160/148‑1 (MQ=255)
                                                        gtggtgCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTAAGCGCTATCAAAGGTAATTGAAATTAAATGTTTAc                                                                                            >  1:6484984/1‑150 (MQ=255)
                                                            tgCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGccc                                                                                        >  2:1082054/1‑150 (MQ=255)
                                                             gCTGGGGGATGTGTGGGGGCTGGTTGCGGGCTCGTTACCGCGTCTGGGTGTTGGGGGGGTGCTCTATTTGAAAAGGGGGGTGGCCGACCGGCCGCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCct                                                                                       <  1:3380435/150‑1 (MQ=255)
                                                             gCTGGGGGACGTGTTGGGGCTGGATGCTGGCGTGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAACGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCct                                                                                       <  2:3730505/150‑1 (MQ=255)
                                                             gCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCct                                                                                       <  2:6360844/150‑1 (MQ=255)
                                                              cTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCcta                                                                                      >  2:1932191/1‑149 (MQ=255)
                                                              cTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTACCCCCGTTGTTTTTAAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCctc                                                                                      >  2:3422911/1‑150 (MQ=255)
                                                                gggggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCTTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCtctc                                                                                    <  2:7180207/150‑1 (MQ=255)
                                                                 ggggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCCCGTTGTTTTTCAGCGCTATAAAAGGTAATTGAAATTAAAGGTTTAAGCCCTCACg                                                                                   >  1:2731736/1‑150 (MQ=255)
                                                                       gtgtTGGGGCTGTTTGCTGGCGCGGTACCGCTTCGGCTTGATGCGGGGGTGCTCTAGATGAAACGGGGGGAGGCCTTCTGGCCGCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATccc                                                                             <  1:2603718/150‑1 (MQ=255)
                                                                       gtgtTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATccc                                                                             <  2:5507397/150‑1 (MQ=255)
                                                                         gtTGGGGCTGTATGCTGGCGTGGTACAGCGTGTGCGTGTTGCTGGGGTGCTCTAGATCAAACGCGGGGAGGCCGTCTGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCt                                                                           <  1:1493139/150‑1 (MQ=255)
                                                                          ttGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCtt                                                                          <  2:5443605/150‑1 (MQ=255)
                                                                                gtGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTCGATGAAAAGCGGGGAGGCCGTCCGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATg                                                                    <  2:32824/149‑1 (MQ=255)
                                                                                gtGGATGCGGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCTGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATg                                                                    <  1:4829728/149‑1 (MQ=255)
                                                                                 tGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCGCCCCCCGTTGTTTTTCAGCGCGATCAAAGGTAATTTAAAATAAAGGTTTAAGCCCACTCGAATACCCATATAAttc                                                                   >  1:859448/1‑148 (MQ=255)
                                                                                   gATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATgcgc                                                                 <  2:3495363/150‑1 (MQ=255)
                                                                                          gCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAATGCGGGGAGGCCGTCCGGCCTCCCCGCGGTGTTGTTCAGCGCGTTCCAAGGTCATTGGAATGAGAGGTTGCAGCCCTCGGGCATCCCCTGATGCTGCGcccgccgt                                                          >  2:3907425/1‑145 (MQ=255)
                                                                                           cGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTc                                                         <  2:4719937/150‑1 (MQ=255)
                                                                                           cGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTc                                                         <  2:2474521/150‑1 (MQ=255)
                                                                                               gTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGc                                                     <  2:1026726/150‑1 (MQ=255)
                                                                                               gTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTAAGCGCTATCAAAGGTAATTTAAATTAAATGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGc                                                     >  2:2666691/1‑150 (MQ=255)
                                                                                                  cAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCGTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGAc                                                  >  1:6198924/1‑150 (MQ=255)
                                                                                                         tGCGTGATGCTGGGGGGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGaa                                           <  1:7204099/150‑1 (MQ=255)
                                                                                                           cGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTTAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTACAGCGACATAGGAACg                                         >  2:3488404/1‑150 (MQ=255)
                                                                                                           cGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACg                                         >  2:3694391/1‑150 (MQ=255)
                                                                                                                   tGGGGTGCTCTAGATGAAAAGCGGGGAGGCCCTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTTAAGCCCTCTCGAATCCCCGTATAATTCGCCCCACATCCACCGACAACCGAACGaaaaaaac                                 >  2:1258226/1‑147 (MQ=255)
                                                                                                                     gggTGCTCTCGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCt                               <  2:3225838/150‑1 (MQ=255)
                                                                                                                     gggTGCTCTAGATGAATAGGGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCt                               <  2:6964854/150‑1 (MQ=255)
                                                                                                                     gggTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGAc                                                                                            <  1:2979812/89‑1 (MQ=255)
                                                                                                                     gggTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGAc                                                                                            >  2:2979812/1‑89 (MQ=255)
                                                                                                                        tGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGTGACATCGGAACGAAAAACTCCTTga                            <  1:3662304/150‑1 (MQ=255)
                                                                                                                        tGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTga                            <  1:6692426/150‑1 (MQ=255)
                                                                                                                           tctAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGAt                         >  2:1211166/1‑150 (MQ=255)
                                                                                                                            ctAGATGAAAAGCGGGGAGGCCGTCTGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGAtt                        <  1:4679598/150‑1 (MQ=255)
                                                                                                                             tAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTc                       >  2:3945247/1‑150 (MQ=255)
                                                                                                                              aGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCCCGTTGTTTTTCAGCGCTATCAAAGGTAATTGAAAATAAAGGTTGACGCCCACTCGAATCCCCGTATAATGCGCCCCACATACAGCGACATCGGAACCAAAAACACCGTGAGAttaa                      >  2:712211/1‑148 (MQ=255)
                                                                                                                              aGATGAAAAGCGGGGAGGCCGACAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCa                      <  2:5001992/150‑1 (MQ=255)
                                                                                                                               gATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCaa                     <  2:7304725/150‑1 (MQ=255)
                                                                                                                                  gAAATGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGa                  <  1:738238/150‑1 (MQ=255)
                                                                                                                                  gAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGa                  <  1:4251396/150‑1 (MQ=255)
                                                                                                                                  gAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGa                  <  1:1149448/150‑1 (MQ=255)
                                                                                                                                   aaaaGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAg                 <  1:3056634/150‑1 (MQ=255)
                                                                                                                                       gCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTaa             <  2:957967/150‑1 (MQ=255)
                                                                                                                                       gCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTaa             >  1:5113839/1‑150 (MQ=255)
                                                                                                                                       gCGGGGAGGCCGTCAGGCCTCCCCGCTGTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTag             >  1:2807028/1‑149 (MQ=255)
                                                                                                                                          gggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAagta          >  2:600123/1‑150 (MQ=255)
                                                                                                                                          gggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAagta          <  1:4877047/150‑1 (MQ=255)
                                                                                                                                            gAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAagtagt        >  2:1564912/1‑150 (MQ=255)
                                                                                                                                               gCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGc                                                                                          <  2:1504753/65‑1 (MQ=255)
                                                                                                                                               gCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTaa     >  2:4215078/1‑150 (MQ=255)
                                                                                                                                               gCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTaa     >  1:5868198/1‑150 (MQ=255)
                                                                                                                                               gCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTaa     >  1:2651542/1‑150 (MQ=255)
                                                                                                                                                 cGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAgg   >  1:7151438/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGGAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGC  >  NC_002947/170692‑170987

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: ‑