Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	789,358	C→T	V208I (GTC→ATC)	galK ←	galactokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	789,358	0	C	T	100.0%	61.2 / NA	19	V208I (GTC→ATC)	galK	galactokinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/11);  total (8/11)

ACATCACGCAGGGCTGGCTGCTGGAAGAAACGCGCACCGGTTTCGCACTGTTCACGACGGGTGTTGTATTCGCTGCCAACCAGGGTACGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAAC  >  NC_000913/789246‑789488                                                                                                                 |                                                                                                                                   aCATCACGCAGGGCTGGCTGCTGGAAGAAACGCGCACCGGTTTCGCACTGTTCACGACGGGTGTTGTATTCGCTGCCAACCAGGGTACGTTTGAAGTTACTGTTGATGATGATGACAGCCACACCTTTGGGCATGGAAAc                                                                                                         >  1:43763/1‑140 (MQ=255)    tCACGCAGGGCTGGCTGCTGGAAGAAACGCGCACCGGTTTCGCACTGTTCACGACGGGTGTTGTATTCGCTGCCAACCAGGGTACGTTTGAAGTTACTGTTGATGATGATGACAGCCACACCTTTg                                                                                                                    <  2:634842/126‑1 (MQ=255)    tCACGCAGGGCTGGCTGCTGGAAGAAACGCGCACCGGTTTCGCACTGTTCACGACGGGTGTTGTATTCGCTGCCAACCAGGGTACGTTTGAAGTTACTGTTGATGATGATGACAGCCACACCTTTg                                                                                                                    >  1:634842/1‑126 (MQ=255)     cACGCAGGGCTGGCTGCTGGAAGAAACGCGCACCGGTTTCGCACTGTTCACGACGGGTGTTGTATTCGCTGCCAACCAGGGTACGTTTGAAGTTACTGTTGATGATGATGACAGCCACACCTTTGGGCATGGAAACTGCt                                                                                                     <  1:610030/140‑1 (MQ=255)       cGCAGGGCTGGCTGCTGGAAGAGACGCGCCCCGGTTTGGCACTGTTCACGACGGGTGTTGTATTCGCTGCCAACCAGGGTACGTTTGAAGTTACTGTTGATGATGATGACAGCCACACCTTTGGGCATGGAAACTGCttt                                                                                                   <  1:240075/140‑1 (MQ=255)                                             gCACTGTTCACGACGGGTGTTGTATTCGCTGCCAACCAGGGTACGTTTGAAGTTACTGTTGATGATGATGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGAt                                                             <  2:558280/140‑1 (MQ=255)                                                   ttCACGACGGGTGTTGTATTCGCTGCCAACCAGGGTACGTTTGAAGTTACTGTTGATGATGATGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGAtctttct                                                       >  2:67141/1‑140 (MQ=255)                                                       cgacgGGTGTTGTATTCTCTGCCAACCAGGGAACGTTTGCAGTTACTGTTGATGATGATGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGcc                                                   <  2:200784/140‑1 (MQ=255)                                                            ggTGTTGTATTCGCTGCCAACCAGGGTACGTTTGAAGTTACTGTTGATGATGATGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAgcg                                              <  1:412294/140‑1 (MQ=255)                                                             gtgtTGTATTCGCTGCCAACCAGGGTACGTTTGAAGTTACTGTTGATGATGATGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGAt                                                             >  2:146013/1‑124 (MQ=255)                                                             gtgtTGTATTCGCTGCCAACCAGGGTACGTTTGAAGTTACTGTTGATGATGATGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGAt                                                             <  1:146013/124‑1 (MQ=255)                                                             gtgtTGTATTCGCTGCCAACCAGGGTACGTTTGAAGTTACTGTTGATGATGATGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAgcgc                                             >  2:344913/1‑140 (MQ=255)                                                                       cGCTGCCAACCAGGGTACGTTTGAAGTTACTGTTGATGATGATGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGc                                   <  2:550143/140‑1 (MQ=255)                                                                       cGCTGCCAACCAGGGTACGTTTGAAGTTACTCTTGCTGATGATGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGc                                   <  2:170269/140‑1 (MQ=255)                                                                                         gTTTGAAGTTACTGTTGATGATGATGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAg                 >  2:533536/1‑140 (MQ=255)                                                                                             gAAGTTACTGTTGATGATGATGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTAc             <  1:67141/140‑1 (MQ=255)                                                                                                   aCTGTTGATGATGATGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTa       >  2:615659/1‑140 (MQ=255)                                                                                                     tGTTGATGATGATGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTCCCGAGCGCCGAAATTAGCTGATCCATGATCGCGCAGTTACAGCCTACa     >  2:271656/1‑140 (MQ=255)                                                                                                       ttGATGATGATGACAGCCACACCTTGGGGGATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACaaa   <  1:264223/140‑1 (MQ=255)                                                                                                        tgatgatgatgaCAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAAc  <  1:615659/140‑1 (MQ=255)                                                                                                                 |                                                                                                                                   ACATCACGCAGGGCTGGCTGCTGGAAGAAACGCGCACCGGTTTCGCACTGTTCACGACGGGTGTTGTATTCGCTGCCAACCAGGGTACGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAAC  >  NC_000913/789246‑789488

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