Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	170,841	A→C	9.4%	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_002947	170,841	0	A	C	9.4%	106.1 / 2.8	41	intergenic (+22/‑548)	PP_0159/PP_16SA	family III CoA‑transferase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (16/20);  new base C (2/2);  total (18/23)
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

GGAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAG  >  NC_002947/170693‑170985                                                                                                                                                     |                                                                                                                                                 ggAGACGCTAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCa                                                                                                                                                  <  2:213309/149‑1 (MQ=255)  gagaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGAGGGGGGGCCGTCAg                                                                                                                                                 >  8:271452/1‑149 (MQ=255)    gaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAgga                                                                                                                                               >  5:340346/1‑148 (MQ=255)         caggtgAGTAGCGCCGGGCGCCACCACTGCTGGGCGAGAATACCGACGTGGTGGGGGGGTATGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCTGACGGCCTccc                                                                                                                                          <  1:276021/144‑1 (MQ=255)                 tATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCTGACGGCCTccc                                                                                                                                          <  4:65658/141‑1 (MQ=255)                 tATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCTGACGGCCTccc                                                                                                                                          >  3:65658/1‑141 (MQ=255)                       cGGGCGCCACCACTGCTGGGCGAGCATACCTAAGTGGTGCGGGTGGAAGTGTTGGGGCTGGATGCTGGCTAGGTACAGCTTCTGCGTGATGCTGGGGGGCTCTATATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCa                                                                                                                            <  6:268544/149‑1 (MQ=38)                         ggggcaaccacTGGTGGGCGAGAATAGCGAAGTGGTGCTGGGGGATGTTTTGGGGCTGGATGCTGGCTAGTTACAGCGTCTGCGTGATGCTGGGGTGCTCTTGATTAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAgc                                                                                                                          <  6:340346/143‑1 (MQ=37)                              caccacTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTTTTTTTAAGCGCTa                                                                                                                      >  8:51520/1‑148 (MQ=255)                              caccacTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTa                                                                                                                      >  1:26661/1‑148 (MQ=255)                              caccacTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTAAGCGCTa                                                                                                                      >  1:221975/1‑148 (MQ=255)                               accacTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCCGGCCTCCCCGCGTTGTTTTTGAgggggtt                                                                                                                     >  5:345112/1‑142 (MQ=255)                                ccacTGCTGGGCGAGCAACCCGAAGTGGGGCGGGGGGAAGTGTGGGGGCGGGAGGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGGGCTCAACATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCa                                                                                                                   <  5:148266/149‑1 (MQ=38)                                      ctgGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAgg                                                                                                               <  3:92485/147‑1 (MQ=255)                                      ctgGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTa                                                                                                             <  5:194380/149‑1 (MQ=255)                                            gAGTATACGGATTTGGTGGTTGTGGATGGGGTGGGGCTGGATGATGGGGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTTGTTCAAAATGGGGGTGGCCGTCTGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGaa                                                                                                       <  2:221975/149‑1 (MQ=255)                                            gAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAAtt                                                                                                          <  8:289437/146‑1 (MQ=255)                                            gAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAAtt                                                                                                          >  7:289437/1‑146 (MQ=255)                                             atcATACCGAAGTTGTTGTTGGGTAATGGTTGTGGCTTGATGCTGGCGAGGTACAGCGTCTGCTTGATGCTGGTGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGaaa                                                                                                      <  7:280869/147‑1 (MQ=37)                                                   ccGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAg                                                                                                <  6:92161/149‑1 (MQ=255)                                                     gAAGTGGTGCTGGGGGAAGTGTTGGGGCTGTATGCTGGGGAGGTACTGCGTCTGCTTGATGCTGGGTTGCTCTAGATTAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAg                                                                                                <  7:350304/147‑1 (MQ=39)                                                     gAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTAAGGCGTCCCCGCTTTGTTTTTCAGGGCTATCAAAGGTAAATGAAATTAAAgg                                                                                               >  8:350304/1‑148 (MQ=255)                                                           gtggtgGGGGATGTGTTGGGGCTGGTTGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAATGGGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGc                                                                                        <  2:26661/149‑1 (MQ=255)                                                               tGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGGCCTCAGGGCGACCCCCATTTTTTTTAAGCGCTATCAAAAGTAATTGAAATTAAAGGTTGACGCCctc                                                                                    >  8:133585/1‑149 (MQ=255)                                                                         gtTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTTATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATccc                                                                           <  1:246379/148‑1 (MQ=255)                                                                          ttGGGGCTGTATGGTGGGCAGGTACAGCGTCTGGGTGATGCTTGGGTGCTCTTGATTAAAAGCGTGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATcccc                                                                          <  1:52364/148‑1 (MQ=255)                                                                          ttGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCtctc                                                                                  <  5:66062/140‑1 (MQ=255)                                                                          ttGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTAAGCGCTATCAAAGGTAATTTAAATTAAAGGTTGACGCCCTCACGagatag                                                                           >  6:66062/1‑142 (MQ=255)                                                                              ggCTGGATGCTGGCGAGGTACTGGGTCTGCGTGATGGTGGGGTGCTCTATATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAAtt                                                                                                          <  2:132307/112‑1 (MQ=38)                                                                              ggCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTccc                                                                                                                                          >  3:248886/1‑80 (MQ=255)                                                                              ggCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTccc                                                                                                                                          <  4:248886/80‑1 (MQ=255)                                                                              ggCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAAtt                                                                                                          >  1:132307/1‑112 (MQ=255)                                                                                                    gCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACa                                               >  1:254450/1‑149 (MQ=255)                                                                                                           cgttaggCTGGGGTGCGATCGATCAAACGCGGGCCGGCCGTCAGGCCGCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAAc                                        <  6:345112/143‑1 (MQ=21)                                                                                                                     gggTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTAAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCGCTCGAATCCCCTTATAAAGCGCCCCACTTCCAGCGCAATCGGAACGAAAAACTcc                              >  8:325167/1‑149 (MQ=25)                                                                                                                               taTGAAAAGCGGGTAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCa                    <  7:51520/148‑1 (MQ=255)                                                                                                                                 tGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAAt                  <  2:254450/149‑1 (MQ=255)                                                                                                                                           ggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTaa           <  7:152993/146‑1 (MQ=255)                                                                                                                                           ggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTaa           >  8:152993/1‑146 (MQ=255)                                                                                                                                               gCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGtt     >  2:174690/1‑148 (MQ=255)                                                                                                                                               gCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGtt     <  1:174690/148‑1 (MQ=255)                                                                                                                                                 cGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAg  >  1:141394/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GGAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAG  >  NC_002947/170693‑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: ‑