Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,981,030	G→A	L195L (CTG→TTG)	otsA ←	trehalose‑6‑phosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,981,030	0	G	A	100.0%	21.7 / NA	10	L195L (CTG→TTG)	otsA	trehalose‑6‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (6/4);  total (6/4)

GGCAAACCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTT‑TTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAGGAACGCCAGACGATCGTTTTCTGTCTGGAAACCCAGCAAATCATAATCACAAAGCTGTTCAAGCAAGGTGTCATATGTCGGCAGCGCGTTGAAGATTTCCGGTGTCGGGAAAGGAATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCGTAATTCATGCGCAAATGGCAACAGGTGATAATCGTGGATCCAGATAATGTCATCGTCTTGCAACAGCGGCAGTAATTTATCTG  >  NC_000913/1980801‑1981272                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   ggCAAACCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTT‑TTGCGCTACGTGTCGTGACGCGGGTCCGGTTAGAAAGACAATCCAAGAACGCCa                                                                                                                                                                                                                                            >  1:134128/1‑238 (MQ=255)                                                                                         gcAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTT‑TTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAAGAACGCCAGACGAGCGTTTTCTGTCTGGAAACCCAGCAAATCATAATCACAAAGATGTTCAAGCAAGGTGTCATATGTCGGCAGCGCGTTGAAGATTTCCGGGGTCggg                                                                                                                                       >  1:1041/1‑251 (MQ=255)                                                                                                                            ggTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTT‑TTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAAGaa                                                                                                                                                                                                                                                 >  1:8524/1‑110 (MQ=255)                                                                                                                                    gCCGATCGGGTAGACTTCTGTTCGAAATGCTTCCCATGCGGCTGTATGGCTTTCGTTC‑CTACGTGGTGTGACGCGGGTGAAGTTAGAAATAGAATAAAAGGGCGCCAGACGAGCGTTTTCTGGCTGGAAACGCAGGAACAGTGAATCACAAAGCTGTTCAGGGAAGGTGTCATATTGGGGCAGCGCGTTGAAGATGTCCGGTTACGGGAAAGGAATATGCAGAAAGAAACCAATGCGATTATTCACTcccc                                                                                            <  1:155359/251‑1 (MQ=255)                                                                                                                                        ncctacacgacgctcttccgatctGCTTTGCCCCAGGCTGTATGGCTTT‑TTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAAGAACGCCAGACGATCGTTTTCTGTCTGGAAACCCAGCAAATCATAATCACAAAGCTGTTCAAGCAAGGTGTCATATGTCGGCAGCGCGTTGAAGATTTCCGGTGTCGGGAAAGGAATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGttt                                                                                        >  1:102174/24‑251 (MQ=255)                                                                                                                                                                       ccnNNNNTGTNTGGCTTT‑TTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAAGAGCGCCAGACGATCGTTTTCTGTCTGGAAACCCAGCAAATCATAATCACAAAGCTGTTCAAGCAAGGTGTCATATGTCg                                                                                                                                                                    <  2:101325/144‑1 (MQ=255)                                                                                                                                                                                                            cgcgGGTCAGGTTAGAAAGACAATCCAAGAACGCCAGACGATCGTTTTCTGTCTGGAAACCCACCAAATCATAATCACAAAGATGTTCAAGCAAGGTGTCATATGTCGGCAGCGCGTTGAAGATTTCCGGTGTCGGGAAAGGAATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCGTAATTCATACGCAAATGGCAACAGGTGATAATCGTGGATCCAGATAATGTCATCGTCTTGCAAc                     >  1:154509/1‑251 (MQ=255)                                                                                                                                                                                                               gggTCAGGTTAGAAAGACAATCCAAGAACGCCAGACGATCGTTTTCTGTCTGGAAACCCAGCACATCATAATc                                                                                                                                                                                                    >  1:25032/1‑73 (MQ=255)                                                                                                                                                                                                                               aCAATCCAAGAACGCCAGACGATCGTTTTCTGTCTGGAAACCCAGCAAATCATAATCACAAAGCTGTTCAAGCAAGGTGTCATATGTCGGCAGCGCGTTGAAGATTTCCGGTGTCGGGAAAGGAATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCGTAATTCATGCGCAAATGGCAACAGGTGATAAGCGTGGATCCAGATAATGTCATCGTCTTGCAACAGCGGCAGTAAGTTATCTg  <  1:54068/251‑1 (MQ=255)                                                                                                                                                                                                                                   tCCAAGAACGCCAGACGATCGTTTTCTGTCTGGAAACCCAGCAAATCATAATCACAAAGCTGTTc                                                                                                                                                                                        <  1:53642/65‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   GGCAAACCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTT‑TTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAGGAACGCCAGACGATCGTTTTCTGTCTGGAAACCCAGCAAATCATAATCACAAAGCTGTTCAAGCAAGGTGTCATATGTCGGCAGCGCGTTGAAGATTTCCGGTGTCGGGAAAGGAATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCGTAATTCATGCGCAAATGGCAACAGGTGATAATCGTGGATCCAGATAATGTCATCGTCTTGCAACAGCGGCAGTAATTTATCTG  >  NC_000913/1980801‑1981272

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