Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,508,936	G→A	27.3%	A164V (GCG→GTG)	glk ←	glucokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,508,936	0	G	A	27.3%	34.3 / 12.4	22	A164V (GCG→GTG)	glk	glucokinase
Reads supporting (aligned to +/- strand):  ref base G (10/6);  new base A (3/3);  total (13/9)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

TCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCAGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCAACCCCAAGCCCCGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAA  >  NC_000913/2508795‑2509072                                                                                                                                              |                                                                                                                                         tCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCATGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCGCAAAATCa                                                                                                                                  <  2:117442/150‑1 (MQ=255) tCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCAGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCGCAAAATCa                                                                                                                                  <  1:85912/150‑1 (MQ=255)         gCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCAGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGAcc                                                                                                                          >  1:38908/1‑150 (MQ=255)                      ttCACAATTGCGCGATACAAATTCACCAGCCCAGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGc                                                                                                             <  2:87110/150‑1 (MQ=255)                      ttCACAATTGCGCGATACAAATTCACCAGCCCAGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCACAAAATCAACGTGACCGCCTTCGCCTGGc                                                                                                             >  2:85024/1‑150 (MQ=255)                      ttCACAATTGCGCGATACAAATTCACCAGCCCAGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCACAAAATCAACGTGACCGCCTTCGCCTGGc                                                                                                             >  2:209961/1‑150 (MQ=255)                      ttCACAATTGCGCGATACAAATTCACCAGCCCAGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCACAAAATCAACGTGACCGCCTTCGCCTGGc                                                                                                             >  2:209960/1‑150 (MQ=255)                                 tccaTACAAAATCACCAGCCCAGGGCCAGAAACCAAGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGTAATCATCTTCACTATACGGCGCAAAATCAACGTGACAGCCTTCGCCAGGCAAACTAACCCa                                                                                                  >  2:42477/4‑150 (MQ=255)                                          aTTCACCAGCCCAGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATc                                                                                         >  2:59275/1‑150 (MQ=255)                                                cAGCCCAGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATg                                                                                   >  1:158424/1‑150 (MQ=255)                                                cAGCCCAGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATg                                                                                   >  1:58322/1‑150 (MQ=255)                                                              aaaGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTCTTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGAGGCGCac                                                                     >  2:140564/1‑149 (MQ=255)                                                                    cgcgCTCCGCCGAAACATGACCATTTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCACAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCAACCCCaa                                                               <  1:209960/150‑1 (MQ=255)                                                                    cgcgCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCACAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCAACCCCaa                                                               <  1:85024/150‑1 (MQ=255)                                                                    cgcgCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCACAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCAACCCCaa                                                               <  1:209961/150‑1 (MQ=255)                                                                              cgAAACATGACCAATTTCCGCACTCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCAACCCCAAGCCCCGTTcc                                                     <  2:223057/150‑1 (MQ=255)                                                                              cgAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCAACCCCAAGCCCCGTTcc                                                     <  1:83425/150‑1 (MQ=255)                                                                                       aCCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCAACCCCAAGCCCCGTTCCGGCACCGTa                                            >  2:19387/1‑150 (MQ=255)                                                                                       aCCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACACAACGCTAATCGACAAGGACCAGATGCGCAACCCCAAGCCCCGTTCCGGCACCGTa                                            >  1:30783/1‑150 (MQ=255)                                                                                            tttCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCAACCCCAAGCCCCGTTCCGGCACCGTAAACCg                                       <  2:14107/150‑1 (MQ=255)                                                                                                                            tcttctTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCAACCCCAAGCCCCGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGcca       >  2:156242/1‑150 (MQ=255)                                                                                                                                 ttcACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCAACCCCAAGCCCCGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCaaa  >  2:202498/1‑150 (MQ=255)                                                                                                                                              |                                                                                                                                         TCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCAGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCAACCCCAAGCCCCGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAA  >  NC_000913/2508795‑2509072

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