Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	170,838	G→T	9.8%	intergenic (+19/‑551)	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_002947	170,838	0	G	T	9.8%	98.2 / 4.5	41	intergenic (+19/‑551)	PP_0159/PP_16SA	family III CoA‑transferase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (20/17);  new base T (2/2);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-02

GACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAG  >  NC_002947/170696‑170985                                                                                                                                               |                                                                                                                                                     gaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGc                                                                                                                                                >  3:312850/1‑149 (MQ=255) gaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCGGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGc                                                                                                                                                <  8:350798/149‑1 (MQ=255)         tgtaGTAGCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTTATGCGGGGGTGCTCTAGATGAAAAGCGGGGAGGCCTGACGGCCT‑cccc                                                                                                                                         <  8:48584/144‑1 (MQ=255)                 cGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCGGCGTGATGCGGGGGTGCTCTAGATGAAACGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGttt                                                                                                                                <  2:48401/148‑1 (MQ=255)                 cGCCGGGCCCCACCAGTGCTGGGCGTCAATACCGATGGGGTGGGGGGGGATGGGTGGGGGCTGTATGGGGGCGTGTTCCCGGGTCTGCGTTATGCGGGGGGGCTCTATATTAAACGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGttt                                                                                                                                <  6:226722/148‑1 (MQ=255)                           caccacTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTa                                                                                                                      >  5:75485/1‑148 (MQ=255)                                   ctgGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGTGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGt                                                                                                              >  2:49398/1‑148 (MQ=255)                                           gCATACCGTAGGGTTGCGGGGGTATCTGTGGGGGCGGGATGCTGGCCAGGTACAGCGTCGGCGTGATGCTGGGGTGCTCTAGATCAAACGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGaa                                                                                                       <  8:155332/147‑1 (MQ=37)                                                cggatgtggtgCGGGGGGAAGGGTTGGGGCTGGATGCTGGCGTGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTaaa                                                                                                 <  4:312850/143‑1 (MQ=255)                                                cagAAGGGGTGCTGGGGGAAGTGTTGTGGGTGTATGCTGGCTAGGTACAGCGTCTGCGTTATGCTGGGGTGCTATATATTAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTaaa                                                                                                 <  5:305932/146‑1 (MQ=38)                                                           cTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGGCGCTGTCCCGCCTGCCCCCCTTTGTTTTT‑ATCGCGTTAAAATGTTATTTGAATTTAAAGTTGACGccc                                                                                      >  1:190316/1‑148 (MQ=255)                                                              ggggAAGTGTTGGGGCTGGAGGCGGGCGAGGTCCCGCGTCTGCGTGATGCGGGGGGGCTCTAGATCAAACGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCtctc                                                                                  <  2:293620/149‑1 (MQ=38)                                                                 gatgggtTGGGGGTGGATGCGGGCGAGGTTCAGGGTGTGGTTTATGGTGGGTTTCTCTATATCAAAAGCGGGGTGGCCGTAAGGCGT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGaa                                                                               <  1:346996/144‑1 (MQ=255)                                                                      gtTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATcccc                                                                          >  6:167010/1‑149 (MQ=255)                                                                      gtTGCGGCTGGATGCTGTCGAGGTACAGCGTCAGCGTGAGGCTGGGGTGCTATAGATCAAAACCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATccc                                                                           <  7:62939/148‑1 (MQ=255)                                                                         ggggCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTta                                                                       >  5:73067/1‑149 (MQ=255)                                                                             cTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCCGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAAt                                                                   >  6:344498/1‑149 (MQ=255)                                                                                 aTGCTGGCGAGTTACAGCGTCTGCTGTATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATgcgc                                                               <  6:75485/149‑1 (MQ=255)                                                                                     tGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGcccc                                                            >  1:86318/1‑148 (MQ=255)                                                                                     tGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCCCGTTGTTTTTAAGCGCGATCAAAAGGAATTTAAATTAAATGTTGACGCCCTCTCGCATACCCTTATAATTCGCCCCa                                                           >  4:187978/1‑149 (MQ=255)                                                                                              aCAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGc                                                   >  3:99567/1‑148 (MQ=255)                                                                                                        cGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGCGGCCGGGAGGCCG‑CCCCGCGTTGTTTTTGAGAGCTATGAAAGGTAATTTGATTTTAAGGTTGCCGCGCGGGCGCAGACCCGTATAATGCGGCCCACTTCCCGCGGCATAGGTAc                                        >  6:21089/1‑149 (MQ=255)                                                                                                        cGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAAc                                        >  1:256235/1‑149 (MQ=255)                                                                                                             tGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGaaa                                    <  2:140615/148‑1 (MQ=255)                                                                                                                tGGGGTGCTCTAGATGAAAAGCGGGGCGGCCGGCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCGCTCGAGTCCCCTTATAATGGGCCCCACTTCCAGCGACATCGGAACGCAAAACt                                >  5:212266/1‑149 (MQ=25)                                                                                                                     tGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTg                           <  6:73067/149‑1 (MQ=255)                                                                                                                        tctAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCATTGTTTTTCAGCGCTATAAAAGGTAAATGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTACCGCGACATCGGAACGAAAAACACCtt                            >  1:308769/1‑145 (MQ=25)                                                                                                                         ctAGATGAAAAGCGGGGAGGCCGGCAGGCCG‑CCCCGCGGTGGTTTTCAGGGCGATCAAAGGGAATTGAAAGTAAAGGGTGCCGCCCGCGCCCATCCCCGTATAAGGCGCCCCACGTCCCGCGGCAGCGGCCCGCAAACCTCCGTgagag                       >  6:339631/1‑148 (MQ=255)                                                                                                                          tAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATcccc                                                                          >  5:373036/1‑97 (MQ=255)                                                                                                                          tAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATcccc                                                                          <  6:373036/97‑1 (MQ=255)                                                                                                                             aTGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCaa                   <  5:167010/149‑1 (MQ=255)                                                                                                                              tGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAAt                  <  5:344498/149‑1 (MQ=255)                                                                                                                              tGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAAt                  <  4:99567/149‑1 (MQ=255)                                                                                                                               gAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGGCCTCTCGAAGCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAAt                  >  7:176014/1‑148 (MQ=255)                                                                                                                                     cGGGGAGGCCTGACGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTg                           <  6:116795/133‑1 (MQ=255)                                                                                                                                     cGGGGAGGCCTGACGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTg                           >  5:116795/1‑133 (MQ=255)                                                                                                                                      ggggAGGCCTGACGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAag          >  8:69776/1‑149 (MQ=255)                                                                                                                                      ggggAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCACACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAag          >  3:313517/1‑149 (MQ=255)                                                                                                                                            gCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCCGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGtt     >  2:298353/1‑148 (MQ=255)                                                                                                                                              cGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAg  <  2:86318/149‑1 (MQ=255)                                                                                                                                              cGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAg  >  1:377705/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                     GACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAG  >  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: ‑