Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,791,877	G→A	S225S (TCC→TCT)	kbl ←	2‑amino‑3‑ketobutyrate CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,791,877	0	G	A	100.0%	24.3 / NA	11	S225S (TCC→TCT)	kbl	2‑amino‑3‑ketobutyrate CoA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (8/3);  total (8/3)

CGCGCGTTCGCCCACAGACGGTCACGCAGTTCGCTGCCCGCTTCGACCATCTCCAGTACTTTGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAACGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGGGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCAGGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACACCATCGGTGGCGAT  >  NC_000913/3791632‑3792020                                                                                                                                                                                                                                                      |                                                                                                                                                cgcgcgTTCGCCCACAGACGGTCACGCAGTTCGCTGCCCGCTTCGACCATCTCCAGTACTTTGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAACGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGAGAAcc                                                                                                                                            <  1:9669/251‑1 (MQ=255)    ncgTTCGCCCACAGACGGTCACGCAGTTCGCTGCCCGCTTCGACCATCTCCAGTACTTTGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAACGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGAGAACCACg                                                                                                                                         >  1:260958/2‑251 (MQ=255)           nccACAGACGGTCACGCAGTTCGCTGCCCGCTTCGACCATCTCCAGTACTTTGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAACGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGAGAACCACGACCAttt                                                                                                                                  >  1:80916/2‑251 (MQ=255)             cacaGACGGTCACGCAGTTCGCTGCCCGCTTCGACCATCTCCAGTACTTTGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAACGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGAGAACCACGACCATTTTc                                                                                                                                <  1:6412/251‑1 (MQ=255)                       cACGCAGTTCGCTGCCCGCTTCGACCATCTCCAGTACTTTGATGGACGCGGCAACACTGGCCGGTGACAGCGAAGTGGAGAACAGGTACGGACGAGAACGCTGGCGCAGACAATCAACCACTTCTTTGCGCGCCGCGGTATAACAAACAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGAGAACCACGCCCATTTTCACCGCCaaa                                                                                                                       >  1:313543/1‑250 (MQ=255)                                                                                    nggTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAACGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGAGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCAGGGCATCATATTTATCTGCCAGATCgaaa                                                          >  1:197386/2‑247 (MQ=255)                                                                                                                         cGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCACCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGAGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCAggg                                                                                    >  1:106874/1‑187 (MQ=255)                                                                                                                                           nccACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGAGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCAGGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACAAACCATCGGTGGCGAt  >  1:205055/2‑251 (MQ=255)                                                                                                                                                         cgccgcGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGACCATCACATCGCAGTATTCATGAGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCAGGGCATAATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTg                        >  1:56424/1‑215 (MQ=255)                                                                                                                                                            cgcgGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGAGCACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCAGGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTg                                           >  1:163567/1‑193 (MQ=255)                                                                                                                                                                                             nntANNNNNCGTNCCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGAGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCAGGGCATCATATTTATCTGCCAGATCGCaaa                                                         <  2:197386/144‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                CGCGCGTTCGCCCACAGACGGTCACGCAGTTCGCTGCCCGCTTCGACCATCTCCAGTACTTTGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAACGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGGGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCAGGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACACCATCGGTGGCGAT  >  NC_000913/3791632‑3792020

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