Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,045,760	C→T	K281K (AAG→AAA)	yihG ←	1‑acylglycerol‑3‑phosphate O‑acyltransferase YihG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,045,760	0	C	T	100.0%	47.9 / NA	15	K281K (AAG→AAA)	yihG	1‑acylglycerol‑3‑phosphate O‑acyltransferase YihG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (4/11);  total (4/11)

CTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGT‑CAGAAGCCGGTCTTTTTCCTGCCACAATGAGTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACATGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGCAGTTTATCGAATTGTTTACCCAGTACATTTAGCGCCATCGCAATGCCTGCCGCCTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGG  >  NC_000913/4045534‑4045998                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               atGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATGAAAAGACTCACTCGCATCGGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGAAATGAGAGACGT‑CAGAAGCCGGTCTTTTTCCGGCCACAATGAGTTCAGCCATTGCTGAAAATGGCGTTTAAAGCTGTTGTCATTGATGTAg                                                                                                                                                                                                                        <  1:39188/250‑1 (MQ=255)                    cGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGT‑CAGAAGCCGGTCTTTTTCCTGCCACAATGAGTTCAGCCATTGCTGAAAATGGCGTTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCgg                                                                                                                                                                                                     >  2:9787/1‑251 (MQ=255)                                       cgacaaGGAGATGAGCGAAGAAGAATTTATTTCTCGTGCCGGACGTTTTGTTCATTAAAATCCTCACTCGCCTCTGCGCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACAGGAGAGACGT‑CAGAAGCCGGTCTTTTTCCTGCCACAATGAGTTCAGCCATTGCTGAAAATGGCGTTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGGAAATCcac                                                                                                                                                                                  <  2:197140/246‑1 (MQ=255)                                         gactggAGATGAGCGAAGAATGATTTATGTCTCGTGCCGGACGTTTTTTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACTCTGTGATGACATGAGAGACGT‑CAGACGCCGGTCTTTTTCCGGCCACAATGAGTTCAGCCATTGCTTAAAATGTCGTTTAAAGCTTTTGTCAGTGATGTAGTAGCCATGTAACTCATCGGCAATAGGCTGTAAATCCCCAt                                                                                                                                                                                <  1:165238/247‑1 (MQ=255)                                                          gaagaaTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGT‑CAGAAGCCGGTCTTTTTCCTGCCACAATGAGTTTAGCCATTGCTGAAAATGGCGTTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACATGTACGACAATCCGCGTc                                                                                                                                                               >  2:107225/1‑251 (MQ=255)                                                                        acttgcccgacttttggtttatacaaaTCCTCACCCGTCTCGGCTCATTCAGAGGCGGGTCAAACTTACTTATTCTGTCGCGGTAATGACATTAGAGA‑GTGCCGAATCCCGTCTTTTTCCTGCCACAATGAGTTCAGCCATTGCTGAAAAAGTCGTTTAAAGCTTTTGTCATCCATGTATTCTCCATGTAACTCATCGGCAATAGGCTGTTAATCCACATGTACGACAATCCGCGTCAGTTTACCGCTTaa                                                                                                                                                 <  2:263737/227‑1 (MQ=255)                                                                                   gTTTTGTTCATTAAAATACTCCCTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGT‑CAGAAGCCGGTCTTTTTCCTGCCACAATGAGTTCAGCCATTGCTGAAAATGGCGTTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACATGTACGACTATCCGCGTCAGTTTACCGCTTAACATATCgaaga                                                                                                                                      <  1:62731/251‑1 (MQ=255)                                                                                                           cGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGT‑CAGAAGCCGGTCTTTTTCCTGCCACAATGAGTTCAGCCATTGCTGAAAATGGCGTTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACATGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCgg                                                                                                                <  2:203068/249‑1 (MQ=255)                                                                                                                     cATTCAGAGGCGAGTAAAACTTACTTATTCGGACGCTGTGATGACATGAGAGACGT‑CAGAAGCCGGTCTTTTTCCTGCCACAATGAGTTCAGCCATTGCTGAAAATGGCGTTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACATGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCg                                                                                                     <  1:107225/250‑1 (MQ=255)                                                                                                                         cAGAAGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACCTGTGTGACGT‑CAGAAACCGGTCTTTTTCCTGCCACAATGAGTTCCGCCATTGCTGAAAATGGCGTTTAAAGCTTTTGTCATTGATGTAGTCCCCATGTAACTCATCGGCAATATGCTGTAAATCCACATCTACGCCAATCCGCGTCAGTTTACCGCTTCACATATCGAAGAACGGCCGACGGTTATTGTCCGGCTAACACACCGTAAcc                                                                                                >  2:46703/1‑250 (MQ=255)                                                                                                                                                 tCTGACGCTGTGATGACATGAGAGACGT‑CAGAAGCCGGTCTTTTTCCTGCCACAATGAGTTCAGCCATTGCTGAAAATGGCGTTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACATGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGCAGTTTATCGAATTGttt                                                                         <  1:127208/250‑1 (MQ=255)                                                                                                                                                                                           ttttCCTGCCACAATGAGTTCAGCCATTGCTGAAAATGGCGTTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACATGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGCAGTTTATCGAATTGTTTACCCAGTACATTTAGCGCCATCGCAATGCCTGCCGCCttt                                 <  1:51686/249‑1 (MQ=255)                                                                                                                                                                                                               cAGCCATTGCTGAAAATGGCGTTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACATGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGCAGTTTATCGAATTGTTTACCCAGTACATTTAGCGCCATCGCAATGCCTGCCGCCTTTGGTGGCAACAAGTTTTGAAAag           <  1:124562/251‑1 (MQ=255)                                                                                                                                                                                                                 gCCATTGCTGAAAATGGCGTTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACATGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGCAGTTTATCGAATTGTTTACCCAGTACATTTAGCGCCATCGCAATGCCTGCCGCCTTTGGTGGCAACAAGTTTTGAAAagtg         <  1:222377/251‑1 (MQ=255)                                                                                                                                                                                                                        cTGAAAATGGCGTTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACATGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGCAGTTTATCGAATTGTTTACCCAGTACATTTAGCGCCATCGCAATGCCTGCCGCCTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTggg  >  1:272421/1‑251 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               CTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGT‑CAGAAGCCGGTCTTTTTCCTGCCACAATGAGTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACATGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGCAGTTTATCGAATTGTTTACCCAGTACATTTAGCGCCATCGCAATGCCTGCCGCCTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGG  >  NC_000913/4045534‑4045998

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