Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	213,169	T→C	26.7%	V6225V (GTT→GTC)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	213,169	0	T	C	26.7%	17.2 / 6.2	15	V6225V (GTT→GTC)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (8/3);  new base C (2/2);  total (10/5)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

GACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAAC  >  NC_002947/213035‑213304                                                                                                                                       |                                                                                                                                        gACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATTTCGATCGCAACCGTTGAAGGGGAGGACGTTg                                                                                                                                        >  1:229723/1‑136 (MQ=255) gACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTg                                                                                                                                        >  2:337600/1‑136 (MQ=255)                    ccGTGCCGGCTAATGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGAc                                                                                                                     <  1:217695/136‑1 (MQ=255)                        gCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTg                                                                                                                 >  1:22366/1‑136 (MQ=255)                           ggCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGCTCGCAACCGTTGAAGGGGCGGACGTTGGCAAGTTCGAGCACCTGACGCTGGAc                                                                                                              >  1:288923/1‑136 (MQ=255)                               aaCGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGa                                                                                                          <  2:259558/136‑1 (MQ=255)                                            gCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCAc                                                                                             >  1:109789/1‑136 (MQ=255)                                                    aCCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGa                                                                                     <  1:337600/136‑1 (MQ=255)                                                             gCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGc                                                                            >  1:343257/1‑136 (MQ=255)                                                             gCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGc                                                                            >  1:449594/1‑136 (MQ=255)                                                                    aCAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTAc                                                                     <  2:343257/136‑1 (MQ=255)                                                                    aCAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTAc                                                                     <  2:449594/136‑1 (MQ=255)                                                                           cTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTTAAGGCGCGGACGTTGGCAAGTTCGAGCAACGGACGCTGGACAAGACTCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGc                                                              >  2:105617/1‑136 (MQ=255)                                                                                            aCCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGt                                             >  1:345584/1‑136 (MQ=255)                                                                                                               cGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGcc                          >  2:669190/1‑136 (MQ=255)                                                                                                                                       tGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAAc  >  1:290205/1‑136 (MQ=255)                                                                                                                                       |                                                                                                                                        GACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAAC  >  NC_002947/213035‑213304

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