Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,602,469	G→T	25.0%	L779M (CTG→ATG)	PP_4071 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,602,469	0	G	T	25.0%	32.2 / 3.8	16	L779M (CTG→ATG)	PP_4071	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (4/8);  new base T (2/2);  total (6/10)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

CGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGTCTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGTTGAGCACGTTCCCCCAGTGCCCACTGGTCGATCATGGCCATTTCCGTAAAGCGCTGACTACGCGGTTCAAAGTACTCCTT  >  NC_002947/4602350‑4602601                                                                                                                        |                                                                                                                                     cgccgcTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGTCTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACATGTTGCGCAGGCGCTCg                                                                                                                      >  2:280191/1‑136 (MQ=255) cgccgcTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGTCTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACATGTTGCGCAGGCGCTCg                                                                                                                      >  2:34324/1‑136 (MQ=255)      cTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGTCTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAg                                                                                                                 >  1:139291/1‑136 (MQ=255)                         aaGGATCGTGACGCCATGATCCAGGTCTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCg                                                                                              <  2:445061/136‑1 (MQ=255)                                   aCGCCATGATCCAGGTCTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCa                                                                                    >  2:182244/1‑136 (MQ=255)                                                        aaaTCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGTTGAGCACGTTCCCCCAGTg                                                               >  1:56571/1‑136 (MQ=255)                                                              gCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGTTGAGCACGTTCCCCCAGTGCCCACt                                                         <  1:182244/136‑1 (MQ=255)                                                                   ggAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGTTGAGCACGTTCCCCCAGTGCCCACTGGTCg                                                    <  1:271057/136‑1 (MQ=255)                                                                                         cgcCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGTTGAGCACGTTCCCCCAGTGCCCACTGGTCGATCATGGCCATTTCCGTAAAgc                              >  2:325916/1‑136 (MQ=255)                                                                                                tGTCGATGTAGTCGGCGCTGTACATGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGTTGAGCACGTTCCCCCAGTGCCCACTGGTCGATCATTGCCATTTCCGTAAAGCGCTGACt                       <  1:280191/136‑1 (MQ=255)                                                                                                tGTCGATGTAGTCGGCGCTGTACATGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGTTGAGCACGTTCCCCCAGTGCCCACTGGTCGATCATTGCCATTTCCGTAAAGCGCTGACt                       <  1:34324/136‑1 (MQ=255)                                                                                                tGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGTTGAGCACGTTCCCCCAGTGCCCACTGGTCGATCATGGCCATTTCCGTAAAGCGCTGACt                       <  1:492761/136‑1 (MQ=255)                                                                                                tGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCATTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGTTGAGCACGTTCCCCCAGTGCCCACTGGTCGATCATGGCCATTTCCGTAAAGCGCTGACt                       <  2:230457/136‑1 (MQ=255)                                                                                                            cGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGTGTAATCCAGTTGAGCACGTTCCCCCAGTGCCCACTGGTCGATCATGGCCATTTCCGTAAAGCGCTGACTACGCGGTTCaaa           <  2:474982/136‑1 (MQ=255)                                                                                                            cGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGTTGAGCACGTTCCCCCAGTGCCCACTGGTCGATCATGGCCATTTCCGTAAAGCGCTGACTACGCGGTTCaaa           <  2:52002/136‑1 (MQ=255)                                                                                                                     aCAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGTTGAGCACGTTCCCCCAGTGCCCACTGGTCGATCATGGCCATTTCCGTAAAGCGCTGACTACGCGGTTCAAAGTACTCCtt  <  2:56571/136‑1 (MQ=255)                                                                                                                        |                                                                                                                                     CGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGTCTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGTTGAGCACGTTCCCCCAGTGCCCACTGGTCGATCATGGCCATTTCCGTAAAGCGCTGACTACGCGGTTCAAAGTACTCCTT  >  NC_002947/4602350‑4602601

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