Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	789,597	C→T	S128N (AGT→AAT)	galK ←	galactokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	789,597	0	C	T	100.0%	27.2 / NA	11	S128N (AGT→AAT)	galK	galactokinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (3/8);  total (3/8)

TCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCC  >  NC_000913/789411‑789779                                                                                                                                                                                           |                                                                                                                                                                                       tCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAATTTAACCCGTCACCCTGCGGCACATTGc                                                                                                                                                             <  1:373880/214‑1 (MQ=255)      cAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCCCGATTTGTGCGCCGTCCAGCGGCAGATGATAAATCTGCTCCACTACGGTTCCGTCCGCGCCGTCCAGTGGAGCGGAAGAATTTAACCCGGCGCCCTGCGGCACCTTGCCGCTGATCCCCATGTCCCCGCCGCCGAAGCTGTTGTTACGc                                                                                                                    <  1:262543/250‑1 (MQ=255)                      ttcttGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCATCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAATTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGtt                                                                                                       <  1:440454/247‑1 (MQ=255)                                                                                    ttCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAATTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTGCATGTGCGACAATGGGCGCATCGAgg                                      <  1:161978/250‑1 (MQ=255)                                                                                        gCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAATTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTCCGCAGTTGCAGATGTTTCCCCACGCCCCGAACGTAg                                                                                   >  1:97620/1‑201 (MQ=255)                                                                                               gACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAATTTAACCCGGCACCCTGCGGCACATTGCCGCTGACCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTGTCATGTGCGACAATGGGCGCATCGAGGGAAAACtcgtcg                          <  1:258425/251‑1 (MQ=255)                                                                                                ncanttANGNNCNNTTTNGGCGCCGTCCAGCGGCAGATGNNAAAGCTGCTNCAATACGGTGCCGACCGCGACTTCCAGTGAAGCGGAAGAATTTAACCCGGCACNCTGCGGCACNNTGNCGCGGNGNNNCNNNTCNACGCCGCCGAAGCTGTTGTNACNNAGTTGCNGATGTTTCACCACGCCACGAACGTAg                                                                                   <  2:97620/189‑1 (MQ=255)                                                                                                                        ntCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAATTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGcc  >  1:278670/2‑250 (MQ=255)                                                                                                                                   gatgatAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAATTTaan                                                                                                                                                                                   <  1:53026/62‑2 (MQ=255)                                                                                                                                                          gTTCCGACCGCGACTTCCAGTGAAGCGGAAGAATTTAACCCTGCTCCCTGCTGCACATTGCCGCTGATCTCCATGTCCACGCCGCCGAAGCTGGTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCAtg                                                            <  1:213719/158‑1 (MQ=255)                                                                                                                                                          gTTCCGACCGCGACTTCCAGTGAAGCGGAAGAATTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGc                                                                                                              >  1:199276/1‑108 (MQ=255)                                                                                                                                                                                           |                                                                                                                                                                                       TCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCC  >  NC_000913/789411‑789779

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑