Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,050,159	G→T	18.2%	G72W (GGG→TGG)	PP_0905 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,050,159	0	G	T	18.2%	54.7 / 2.4	22	G72W (GGG→TGG)	PP_0905	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (8/10);  new base T (2/2);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.72e-01

GAGGCCAACCGTGTTCAGCAAGGCTGAGATGGCAGACCAGTGCGCGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAA  >  NC_002947/1050029‑1050277                                                                                                                                   |                                                                                                                       gAGGCCAACCGTGTTCAGCAAGGCTGAGATGGCAGACCAGTGCGCGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTg                                                                                                                   >  1:280268/1‑136 (MQ=255)    gacAACCGTGTTCAGCAAGGCTGAGATGGCAGACCAGTGCGCGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCgg                                                                                                                <  1:64492/134‑1 (MQ=255)      cAACCGTGTTCAGCAAGGCTGAGATGGCAGACCAGTGCGCGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTc                                                                                                              <  1:73091/136‑1 (MQ=255)          cGTGTTCAGCAAGGCTGAGATGGCAGACCAGTGCGCGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAg                                                                                                          >  1:259188/1‑136 (MQ=255)          cGTGTTCAGCAAGGCTGAGATGGCAGACCAGTGCGCGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAg                                                                                                          <  2:259188/136‑1 (MQ=255)          cGTGTTCAGCAAGGCTGAGATGGCAGACCAGTGCGCGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAg                                                                                                          <  2:138398/136‑1 (MQ=255)           gtgtTCAGCAAGGCTGAGATGGCAGACCAGTGCGCGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAgg                                                                                                         <  2:11037/136‑1 (MQ=255)                         tGAGATGGCAGACCAGTGCGCGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGATCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGt                                                                                           >  2:96695/1‑136 (MQ=255)                                    aCCAGTGCGCGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAATGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGaaaa                                                                                >  1:308716/1‑136 (MQ=255)                                            gcgTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAgg                                                                        >  2:324586/1‑136 (MQ=255)                                            gcgTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAgg                                                                        >  2:509813/1‑136 (MQ=255)                                              gTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAATGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGcc                                                                      <  2:308716/136‑1 (MQ=255)                                                                               ccGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCt                                     <  1:324586/136‑1 (MQ=255)                                                                                        gcCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCg                            <  1:96695/136‑1 (MQ=255)                                                                                        gcCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCg                            <  1:152779/136‑1 (MQ=255)                                                                                                aGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCtta                    >  1:518006/1‑135 (MQ=255)                                                                                                aGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTc                    >  1:88089/1‑136 (MQ=255)                                                                                                aGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTc                    >  1:271796/1‑136 (MQ=255)                                                                                                       tttGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCa             <  2:280268/136‑1 (MQ=255)                                                                                                                  tGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCaa  <  1:509813/136‑1 (MQ=255)                                                                                                                     ttGGCCTGCGTCAATGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTc                                      <  1:84700/97‑1 (MQ=255)                                                                                                                     ttGGCCTGCGTCAATGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTc                                      >  2:84700/1‑97 (MQ=255)                                                                                                                                   |                                                                                                                       GAGGCCAACCGTGTTCAGCAAGGCTGAGATGGCAGACCAGTGCGCGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAA  >  NC_002947/1050029‑1050277

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