Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,822,536	G→T	31.2%	L168M (CTG→ATG)	PP_5526 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,822,536	0	G	T	31.2%	25.2 / 7.6	16	L168M (CTG→ATG)	PP_5526	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (3/8);  new base T (2/3);  total (5/11)
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.45e-01

ACGTAGAGTTCCTTTTCCTCGTCGGAGGCTCCGCTATGGTCGAAGAACTCCAGCAGCTCAGTGAGCTTTGTAAAGGGCTCGGCGCGCCACGTTAGCCTGGGCTTGTTTTTCTCAAGCCAACTCCATGCCAGAGCGCATTTTGTGCGTTCGTCCTTGTCTTTGAACCAGCTGAAGATCCTGTCTGACCGCGCGTGCTCGGTCCAGGCGTCTGATAGGCGCTTCAATGCACTCTCTTTCT  >  NC_002947/2822406‑2822643                                                                                                                                   |                                                                                                            aCGTAGAGTTCCTTTTCCTCGTCGGAGGCTCCGCTATGGTCGAAGAACTCCAGCAGCTCAGTGAGCTTTGTAAAGGGCTCGGCGCGCCACGTTAGCCTGGGCTTGTTTTTCTCAAGCCAACTCCATGCCAGAgcgc                                                                                                        >  1:32409/1‑136 (MQ=255)                           ggCTCCGCTATGGTCGAAGAACTCCAGCAGCTCAGTGAGCTTTGTAAAGGGCTCGGCGCGCCACGTTAGCCTGGGCTTGTTTTTCTCAAGCCAAATCCATGCCAGAGCGCATTTTGTGCGTTCGTCCTTGTCTTTg                                                                              <  1:524811/136‑1 (MQ=255)                                                cTCCAGCAGCTCAGTGAGCTTTGTATAGGGCTCGGCGCGCCACGTTAGCCTGTGCTTGTTTTTCTCAAGCCAACTCCATGCCAGAGCGCATTTTGTGCGTTCGTCCTTGTCTTTGAACCAGCTGAAGATCCTGTCt                                                         <  1:477653/136‑1 (MQ=255)                                                          tCAGTGAGCTTTGTAAAGGGCTCGGCGCGCCACGTTAGCCTGGGCTTGTTTTTCTCAAGCCAACTCCATGCCAGAGCGCATTTTGTGCGTTCGTCCTTGTCTTTGAACCAGCTGAAGATCCTGTCTGACCGCGCGt                                               <  1:45792/136‑1 (MQ=255)                                                          tCAGTGAGCTTTGTAAAGGGCTCGGCGCGCCACGTTAGCCTGGGCTTGTTTTTCTCAAGCCAACTCCATGCCAGAGCGCATTTTGTGCGTTCGTCCTTGTCTTTGAACCAGCTGAAGATCCTGTCTGACCGCGCGt                                               <  2:474086/136‑1 (MQ=255)                                                           cAGTGAGCTTTGTAAAGGGCTCGGCGCGCCACGTTAGCCTGGGCTTGTTTTTCTCAAGCCAACTCCATGCCATAGCGCATTTTGTGCGTTCGTCCTTGTCTTTGAACCAGCTGAAGATCCTGTCTGCCCGCGCGTg                                              >  2:349388/1‑136 (MQ=255)                                                           cAGTGAGCTTTGTAAAGGGCTCGGCGCGCCACGTTAGCCTGGGCTTGTTTTTCTCAAGCCAACTCCATGCCATAGCGCATTTTGTGCGTTCGTCCTTGTCTTTGAACCAGCTGAAGATCCTGTCTGACCGCGCGTg                                              >  2:352833/1‑136 (MQ=255)                                                                aGCTTTGTAAAGGGCTCGGCGCGCCACGTTAGCCTGGGCTTGTTTTTCTCAAGCCAACTCCATGCCAGAGCGCATTTTGTGCGTTCGTCCTTGTCTTTGAACCAGCTGAAGATCCTGTCTGACCGCGCGTGCTCgg                                         >  2:346538/1‑136 (MQ=255)                                                                   tttGTAAAGGGCTCGGCGCGCCACGTTAGCCTGGGCTTGTTTTTCTCAAGCCAACTCCATGCCATAGCGCATTTTGTGCGTTCGTCCTTGTCTTTGAACCAGCTGAAGATCCTGTCTGACCGCGCGTGCTCTGTcc                                      <  1:349388/136‑1 (MQ=255)                                                                   tttGTAAAGGGCTCGGCGCGCCACGTTAGCCTGGGCTTGTTTTTCTCAAGCCAACTCCATGCCATAGCGCATTTTGTGCGTTCGTCCTTGTCTTTGAACCAGCTGAAGATCCTGTCTGACCGCGCGTGCTCTGTcc                                      <  1:352833/136‑1 (MQ=255)                                                                           gggCTCGGCGCGCCACGTTAGCCTGGGCTTGTTTTTCTCAAGCCAACTCCATGCCAGAGCGCATTTTGTGCGTTCGTCCTTGTCTTTGAACCAGCTGAAGATCCTGTCTGACCGCGCGTGCTCGGTCCAGGCGTCt                              <  2:128293/136‑1 (MQ=255)                                                                              cTCGGCGCGCCACGTTAGCCTGGGCTTGTTTTTCTCAAGCCAACTCCATGCCAGAGCGCATTTTGTGCGTTCGTCCTTGTCTTTGAACCAGCTGAAGATCCTGTCTGACCGCGCGTGCTCGGTCCAGGCGTCTGAt                           <  2:32409/136‑1 (MQ=255)                                                                              cTCGGCGCGCCACGTTAGCCTGGGCTTGTTTTTCTCAAGCCAACTCCATGCCAGAGCGCATTTTGTGCGTTCGTCCTTGTCTTTGAACCAGCTGAAGATCCTGTCTGACCGCGCGTGCTCGGTCCAGGCGTCTGAt                           <  2:335185/136‑1 (MQ=255)                                                                                    gcgcCACGTTAGCCTGGGCTTGTTTTTCTCAAGCCAACTCCATGCCATAGCGCATTTTGTGCGTTCGTCCTTGTCTTTGAACCAGCTGAAGATCCTGTCTGACCGCGCGTGCTCGGTCCAGGCGTCTGATAGgcgc                     <  1:166433/136‑1 (MQ=255)                                                                                        cACGTTAGCCTGGGCTTGTTTTTCTCAAGCCAACTCCATGCCAGAGCGCATTTTGTGCGTTCGTCCTTGTCTTTGAACCAGCTGAAGATCCTGTCTGACCGCGCGTGCTCGGTCCAGGCGTCTGATAGGCGCTTCa                 <  1:343067/136‑1 (MQ=255)                                                                                                       ttGTTTTTCTCAAGCCAACTCCATGCCAGAGCGCATTTTATGCGTTCGTCCTTGTCTTTGAACCAGCTGAAGATCCTGTCTGACCGCGCGTGCTCGGTCCAGGCGTCTGATAGGCGCTTCAATGCACTCTCTTTCt  >  2:245930/1‑136 (MQ=255)                                                                                                                                   |                                                                                                            ACGTAGAGTTCCTTTTCCTCGTCGGAGGCTCCGCTATGGTCGAAGAACTCCAGCAGCTCAGTGAGCTTTGTAAAGGGCTCGGCGCGCCACGTTAGCCTGGGCTTGTTTTTCTCAAGCCAACTCCATGCCAGAGCGCATTTTGTGCGTTCGTCCTTGTCTTTGAACCAGCTGAAGATCCTGTCTGACCGCGCGTGCTCGGTCCAGGCGTCTGATAGGCGCTTCAATGCACTCTCTTTCT  >  NC_002947/2822406‑2822643

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