Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,419,739	C→T	25.0%	A32V (GCC→GTC)	PP_3906 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,419,739	0	C	T	25.0%	24.8 / 4.4	16	A32V (GCC→GTC)	PP_3906	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (6/6);  new base T (2/2);  total (8/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

CCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATGGTGCTCGCAGCGTTGCTTGAACGCAGCCAGGAAGCGCGCAACGAAGCCTTCGCCAGAATCGCACAGGGAGCCCCGCGCTACCAAGCATCTGGCAAGGGGAAGACCTGGGATGTGGTCGAGATCGCCACCGGCGCCGTGCAAGGTTTCGCCTACAGCTACAAGACGGCGCTGCGCTTCGTGG  >  NC_002947/4419621‑4419867                                                                                                                       |                                                                                                                                 ccTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATGGTGCTCGCAGCGTTGCTTGAACGCAGCCAGGAAGCGCGCAACGAAGCCTTCGCCAGAATCGCACAGGGAg                                                                                                                 >  1:304435/1‑136 (MQ=255)  cTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATGGTGCTCGCAGCGTTGCTTGAACGCAGCCAGGAAGCGCGCAACGAAGCCTTCGCCAGAATCGCACAGGGAGc                                                                                                                >  1:342549/1‑136 (MQ=255)    gACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATGGTGCTCGCAGCGTTGCTTGAACGCAGCCAGGAAGCGCGCAACGAAGCCTTCGCCAGAATCGCACAGGGAGccc                                                                                                              >  1:28471/1‑136 (MQ=255)     acccacccGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATGGTGCTCGCAGCGTTGCTTGAACGCAGCCAGGAAGCGCGCAACGAAGCCTTCGCCAGAATCGCACAGGGAGcccc                                                                                                             <  1:490523/136‑1 (MQ=255)          cccGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATGGTGCTCGCAGCGTTGCTTGAACGCAGCCAGGAAGCGCGCAACGAAGCCTTCGCCAGAATCGCACAGGGAGCCCCGCGCt                                                                                                        <  2:17855/136‑1 (MQ=255)             gAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATGGTGCTCGCAGCGTTGCTTGAACGCAGCCAGGAAGCGCGCAACGAAGCCTTCGTCAGAATCGCACAGGGAGCCCCGCGCTAcc                                                                                                     >  2:236384/1‑136 (MQ=255)             gAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATGGTGCTCGCAGCGTTGCTTGAACGCAGCCAGGAAGCGCGCAACGAAGCCTTCGTCAGAATCGCACAGGGAGCCCCGCGCTAcc                                                                                                     >  2:238662/1‑136 (MQ=255)                                                tttCGAGGACAAGCGAATGGTGCTCGCAGCGTTGCTTGAACGCAGCCAGGAAGCGCGCAACGAAGCCTTCGCCAGAATCGCACAGGGAGCCCCGCGCTACCAAGCATCTGGCAAGGGGAAGACCTGGGATGTGGtg                                                                  >  1:412361/1‑135 (MQ=255)                                                tttCGAGGACAAGCGAATGGTGCTCGCAGCGTTGCTTGAACGCAGCCAGGAAGCGCGCAACGAAGCCTACGCCAGAATCGCACAGGGAGCCCCGCGCTACCAAGCATCTGGCAAGGGGAAGACCTGGGATGTGGTc                                                                  >  2:392923/1‑136 (MQ=255)                                                             cGAATGGTGCTCGCAGCGTTGCTTGAACGCAGCCAGGAAGCGCGCAACGAAGCCTTCGCCAGAATCGCACAGGGAGCCCCGCGCTACCAAGCATCTGGCAAGGGGAAGACCTGGGATGTGGTCGAGATCGCCACCg                                                     <  2:437620/136‑1 (MQ=255)                                                                                                   aGCGCGCAACGAAGCCTTCGCCAGAATCGCACAGGGAGCCCCGCGCTACCAAGCATCTGGCAAGGGGAAGACCTGGGATGTGGTCGAGATCGCCACTGGCGCCGTGCAAGGTTTCGCCTACAGCTACAAGACGgcg               <  2:294475/136‑1 (MQ=255)                                                                                                   aGCGCGCAACGAAGCCTTCGCCAGAATCGCACAGGGAGCCCCGCGCTACCAAGCATCTGGCAAGGGGAAGACCTGGGATGTGGTCGAGATCGCCACCGGCGCCGTGCAAGGTTTCGCCTACAGCTACAAGACGgcg               <  1:236621/136‑1 (MQ=255)                                                                                                          aaCGAAGCCTTCGCCAGAATCGCACAGGGAGCCCCGCGCTACCAAGCATCTGGCAAGGGGAAGACCTGGGATGTGGTCGAGATCGCCACCGGCGCCGTGCAAGGTTTCGCCTACAGCTACAAGACGgcgctgcgct        >  1:53578/1‑136 (MQ=255)                                                                                                                gCCTTCGTCAGAATCGCACAGGGAGCCCCGCGCTACCAAGCATCTGGCAAGGGGAAGACCTGGGATGTGGTCGAGATCGCCACCGGCGCCGTGCAAGGTTTCGCCTACAGCTACAAGACGGCGCTGCGCTTCGTgg  <  1:236384/136‑1 (MQ=255)                                                                                                                gCCTTCGTCAGAATCGCACAGGGAGCCCCGCGCTACCAAGCATCTGGCAAGGGGAAGACCTGGGATGTGGTCGAGATCGCCACCGGCGCCGTGCAAGGTTTCGCCTACAGCTACAAGACGGCGCTGCGCTTCGTgg  <  1:238662/136‑1 (MQ=255)                                                                                                                gCCTTCGCCAGAATCGCACAGGGAGCCCCGCGCTACCAAGCATCTGGCAAGGGGAAGACCTGGGATGTGGTCGAGATCGCCACCGGCGCCGTGCAAGGTTTCGCCTACAGCTACAAGACGGCGCTGCGCTTCGTgg  <  2:108326/136‑1 (MQ=255)                                                                                                                       |                                                                                                                                 CCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATGGTGCTCGCAGCGTTGCTTGAACGCAGCCAGGAAGCGCGCAACGAAGCCTTCGCCAGAATCGCACAGGGAGCCCCGCGCTACCAAGCATCTGGCAAGGGGAAGACCTGGGATGTGGTCGAGATCGCCACCGGCGCCGTGCAAGGTTTCGCCTACAGCTACAAGACGGCGCTGCGCTTCGTGG  >  NC_002947/4419621‑4419867

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