Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,100,955	G→T	25.1%	L98F (TTG→TTT)	PP_5607 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,100,955	0	G	T	25.1%	25.8 / 4.5	16	L98F (TTG→TTT)	PP_5607	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (3/9);  new base T (2/2);  total (5/11)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

AACAACTCCCCTCTAGCCCTGCTGTTTGCCCTGTTTGCCTGGGGCCTGATGCGGCGCTGGCACAGCTCGAAGCTGGCGCTGCTGAGTGCGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCA  >  NC_002947/4100833‑4101067                                                                                                                           |                                                                                                                 aacaacTCCCCTCTAGCCCTGCTGTTTGCCCTGTTTGCCTGGGGCCTGATGCGGCGCTGGCACAGCTCGAAGCTGGCGCTGCTGAGTGCGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGt                                                                                                     >  1:72118/1‑136 (MQ=255)  actacTCCCCTCTAGCCCTGCTGTTTGCCCTGTTTGCCTGGGGCCTGTTGCGGCGCTGGCACAGCTCGAAGCTGGCGCTGCTGAGTGCGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGtt                                                                                                    <  2:163719/133‑1 (MQ=255)    aacTCCCCTCTAGCCCTGCTGTTTGCCCTGTTTGCCTGGGGCCTGATGGGGCGCTGGCACAGCTCGTAGCTGGCGCTGCTGAGTGCGGGCGTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTcc                                                                                                  <  2:250350/136‑1 (MQ=255)                 cccTGCTGTTTGCCCTGTTTGCCTGGGGCCTGATGCGGCGCTGGCACAGCTCGAAGCTGGCGCTGCTGAGTGCGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAggg                                                                                     >  2:363154/1‑136 (MQ=255)                   ctgctgTTTGCCCTGTTTGCCTGGGGCCTGATGCGGCGCTGGCACAGCTCGAAGCTGGCGCTGCTGAGTGCGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTTGGCACCGGCTGGTTCCCCTGCGACCAGGGgt                                                                                   >  1:148882/1‑136 (MQ=255)                   ctgctgTTTGCCCTGTTTGCCTGGGGCCTGATGCGGCGCTGGCACAGCTCGAAGCTGGCGCTGCTGAGTGCGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTTGGCACCGGCTGGTTCCCCTGCGACCAGGGgt                                                                                   >  1:266653/1‑136 (MQ=255)                   ctgctgTTTGCCCTGTTTGCCTGGGGCCTGATGCGGCGCTGGCACAGCTCGAAGCTGGCGCTGCTGAGTGCGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGgt                                                                                   <  2:392496/136‑1 (MQ=255)                      ctgtttgccctgtttgccTGGGGCCTGATGCGGCGCTGGCACAGCTCGAAGCTGGCGCTGCTGAGTGCGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGtgtg                                                                                >  1:217705/1‑136 (MQ=255)                      ctgtttgccctgtttgccTGGGGCCTGATGCGGCGCTGGCACAGCTCGAAGCTGGCGCTGCTGAGTGCGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGtgtg                                                                                <  1:362180/136‑1 (MQ=255)                                                              attgCTCGATGTTTGCTCTTCTGTGTGCGTGTCTTTTGTTTCTTTATGGTATGGGCTGTTTGTGCTCCGTCTGGTTCCCCTGCGATCTGGTGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCt                                        <  2:168134/133‑1 (MQ=255)                                                                   gcGAAGCTGGCGCTGCTGAGTGCGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTg                                   <  1:324125/135‑1 (MQ=255)                                                                                tgctgAGTGCGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGtt                      <  2:224762/136‑1 (MQ=255)                                                                                     aGTGCGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACGTGTTTGGTCTGCTGATGTTCctct                 <  1:247261/136‑1 (MQ=255)                                                                                          gggCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTg            <  2:153955/136‑1 (MQ=255)                                                                                                    ttGCTGCATGGCATGGGCAGTTTTGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCa  <  2:266653/136‑1 (MQ=255)                                                                                                    gtgctgcATGGCATGGGCAGTTTTGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCa  <  2:148882/135‑1 (MQ=255)                                                                                                                           |                                                                                                                 AACAACTCCCCTCTAGCCCTGCTGTTTGCCCTGTTTGCCTGGGGCCTGATGCGGCGCTGGCACAGCTCGAAGCTGGCGCTGCTGAGTGCGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCA  >  NC_002947/4100833‑4101067

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: ‑