Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 170,841 A→C 11.8% intergenic (+22/‑548) 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,8410AC11.8% 81.1 / 2.3 33intergenic (+22/‑548)PP_0159/PP_16SAfamily III CoA‑transferase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (16/12);  new base C (2/2);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAG  >  NC_002947/170696‑170985
                                                                                                                                                 |                                                                                                                                                 
gaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGATGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGc                                                                                                                                                <  1:311684/149‑1 (MQ=255)
 aCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACAGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGc                                                                                                                                                <  2:200270/148‑1 (MQ=255)
   gCCAGTGGAGTATCGCCGGGCGCCACCACTGATGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCATCAGGGCTc                                                                                                                                             >  1:294611/1‑149 (MQ=255)
      aGTGGCGTATCGCCGGGCGCCACCACTGCTGGGCGAGCTTACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCGGGCGAGGTACAGCGGCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCTGACGGCCTcccc                                                                                                                                          <  5:265100/149‑1 (MQ=255)
      aGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGGGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTTATGCTGGGGTGCTCTAGATTAAAAGCGGGGAGGCCGTCAGGCCTcccc                                                                                                                                          <  4:283447/149‑1 (MQ=255)
                 gCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGtttt                                                                                                                                >  4:240584/1‑148 (MQ=255)
                        gcCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAgcgc                                                                                                                         >  1:89865/1‑148 (MQ=255)
                              agtgctgGGCGAGAATACCCAAGTGGTGGGGGGGTACGTGTGGGGGCTGTTTGCGGGCTAGGTACCGCGTCTGGGTTTTGGTGGGGTGCTATATATCAAAAGGGGGCAGGCCGTCTGGCCTCCCCGCTTTGTTTTTCAGCGCTATCaaa                                                                                                                  <  6:17978/147‑1 (MQ=255)
                              acTGCTGGGCGAGCATACCGAAGTGATGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCACCGCTTTGTTTTTCAGAGCTATCa                                                                                                                    >  7:188738/1‑147 (MQ=255)
                                   tgGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTaa                                                                                                            <  3:240584/149‑1 (MQ=255)
                                          gattACCGATGTGTTGTTGGGGGATGTGTTGGGGTTGTTTGTTGGCTAGTTACTGCGTCTGCGTTATGGTGGGGTGCTATAGTTTAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGaaa                                                                                                      <  8:188738/145‑1 (MQ=25)
                                             tACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTaa                                                                                                  <  2:89865/149‑1 (MQ=255)
                                              aCCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTa                                                                                                   >  7:158164/1‑147 (MQ=255)
                                                  atgtggtgCTGGGGGAAGTGTGGGGCGGGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCATCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGtt                                                                                             <  2:294611/147‑1 (MQ=255)
                                                       gtgCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATGAAAGGTAATTGAAATTAAAGGTTGACGc                                                                                        >  7:25759/1‑149 (MQ=255)
                                                            gggggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCctc                                                                                    >  7:177302/1‑148 (MQ=255)
                                                              gggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCtctc                                                                                  >  1:79016/1‑148 (MQ=255)
                                                               ggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGGGGCCCTGCGGCCTGCCCCCGTTGGTTTTGACCGGAATCAAAAGTAAATTTAAATTAAAGGTGACGCCCgctcc                                                                                 >  3:20608/1‑147 (MQ=255)
                                                                   gtgtTGGGGGTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCGTTGTTTTTCAGCGCTATAAAAGGTAATTTAAAATAAAGGTTTAAGCCCTATCCAATcc                                                                            >  1:133942/1‑149 (MQ=255)
                                                                      ttGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGGCTCCCCGCTTTGTTTTTCAGGGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCt                                                                         >  6:74062/1‑149 (MQ=255)
                                                                      ttGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATcccc                                                                          <  8:161378/148‑1 (MQ=255)
                                                                                    tGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTACCCCCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAAGTTGACGCCCTCACGAATCCCCTTATAATACGCCCCa                                                           >  4:65165/1‑149 (MQ=255)
                                                                                      gCGAGCTACACCGTCTTCGTGATGCTGGTGTGCTATAGATGAAAAGCGTTGAGTCCGTCAGGCCTCCACGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGAGGCCCACt                                                         <  6:262374/149‑1 (MQ=37)
                                                                                                gCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACa                                               >  1:220195/1‑149 (MQ=255)
                                                                                                                ggggTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTc                               >  1:221504/1‑149 (MQ=255)
                                                                                                                ggggTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTc                               >  8:303467/1‑149 (MQ=255)
                                                                                                                 gggTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTc                               <  8:177302/148‑1 (MQ=255)
                                                                                                                   gTGCTGTAGATGAAATGCGGGGAGTCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTATAGGTTGACTCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCtt                            <  3:65165/149‑1 (MQ=25)
                                                                                                                      ctctAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTga                          <  7:303467/148‑1 (MQ=255)
                                                                                                                                   gCGGGGAGGCCTGACGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTg                           <  6:143451/134‑1 (MQ=255)
                                                                                                                                   gCGGGGAGGCCTGACGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTg                           >  5:143451/1‑134 (MQ=255)
                                                                                                                                           gCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGtt     >  6:153696/1‑148 (MQ=255)
                                                                                                                                             cGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAg  >  3:44613/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAG  >  NC_002947/170696‑170985

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