Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,476,605	G→A	S390N (AGC→AAC)	yjgL →	protein YjgL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,476,605	0	G	A	100.0%	26.0 / NA	12	S390N (AGC→AAC)	yjgL	protein YjgL
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (7/5);  total (7/5)

AAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCATGCAATATTTTGATGATACTAAAAAGATAAAGCATGCATGCTCCGTAATAAATCATATTAATCTTCGCAGCGATCAGTCTAAAATAGCTGAGACATTATTTTTCAATCTCGATAAAGAACCCTATAAAAATAGCCCTGAATTACAGGGGTTGATTTGGAATAAGTTGGTTGTATATGTCAATGAATTTAACTTAAGTAATCGAGAAAAAACAAATTTAATACAAAGG‑CTATTTGAT  >  NC_000913/4476365‑4476843                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                acaaaggaaaacataaagaaATAGATAGGCTCTCCGCATTAAATGATCACGGCATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCACAGATAACCTGACTCATATTCTTATTGATTCCTTCTATGATGCAAAGGTTAACCTTACACCGATCATCGATACGATTCTGGATACTTTTTAAATATCCACATATATTAATGATGTAAGAATCATGGATTGGTGTTTTATTAAAAACATGCTatat                                                                                                                                                                                                                                       >  2:112881/16‑251 (MQ=255)                     aTAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAACATGCAATATTTTGATGATACTaaaa                                                                                                                                                                                                                       >  2:22923/1‑247 (MQ=255)                                    gCATTAACTGATCACGGAATTGATAATGACAAAGTAGGTAACAATAAAATTAATGACATCAAAGTTAACATGTCTCATATTCTTATTGATTACTTAGATGATGCAAAGGTTAACCTTACACAGGTCATCGATTCGATTCTTGAGACTTTTTCAAAATCCCCATCTATTAATGCTGTAAGAATAATGGATTGGTGTATTAATAATCACATGCCATATTTTGATGATACTAAACAGATAAAgcatgcatacct                                                                                                                                                                                                    >  2:158026/1‑247 (MQ=255)                                                                           aaaaaTGAAATTACAGACATCACAGATAACCTGTCTCATATTCTTATTGATTACTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCACCATATATTAATGATGTAATAATACTGGATAGGTGTTTTATTAAAAAAATGCACTATTTTGATGATCATAAAAAGTTAAAACATGCATGCTCCGTAATAATTAATATTAATCTTCGAAGCGTTCAGTCTa                                                                                                                                                              >  2:1385/4‑250 (MQ=255)                                                                                                           gTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGCTTCGCTTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAACATGCAATATTTTGATGATACTAAAAAGATAAAGCATGCATGCTCCGTAATAAATCATATTAATCTTCGCAGCGATCAGTCTaa                                                                                                                                                             >  2:188990/1‑219 (MQ=255)                                                                                                           gTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAACATGCAATATTTTGATGATACTAAAAAGATAAAGCATGCATGCTCCGTAATAAATCATATTAATCTTCGCAGCGATCAGTCTaa                                                                                                                                                             <  1:188990/219‑1 (MQ=255)                                                                                                               cATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAACATGCAATATTTTGATGATACTAAAAAGATAAAGCATGCATGCTCCGTAATAAATCATATTAATCTTCGCAGCGATCAGTCTAAAATAGCTGAGACATTATTTTTCAATCTCGATaa                                                                                                                            <  2:3163/248‑1 (MQ=255)                                                                                                                                    gatgatGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAACATGCAATATTTTGATGATACTAAAAAGGTAAAGCATGCATGCTCCGTAATAAATCATATTAATCTTCGCAGCGATCAGTCTAAAATAGCTGAGACATTATTTTTCAATCTCGATAAAGAACCCTATAAAAATAGCCCTGa                                                                                                    <  1:82828/251‑1 (MQ=255)                                                                                                                                                                                                         atTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAACATGCAATATTTTGATGATACTAAAAAGATAAAGCATGCATGCTCCGTAATAAAGCATATTAATCTTCGCAGCGATCAGTCTAAAATAGCTGAGACATTATTTTTCAATCTCGATAAAGAACCCTATAAAAATAGCCCTGAATTACAGGGGTTGATTTGGAATAAGTTGGTTGTATATGTCAATGAATTTAACTTAAGTAATCGAGaaa                                <  1:8090/250‑1 (MQ=255)                                                                                                                                                                                                                atgTAAGAATACTTGATTGGTGTTTTAATAAAAACATGCAATATTTTGATGATACTAAAAAGATAAAGCATGCATGCTCCGTAATAAATCATATTAATCTTCGCAGCGATCAGACTAAAATAGCGAAGACATTATTTGTCAATCTCGATAAAGAACCCTATAAAAATAGCCCTGAATTACAGGGGTTGATTTGGAATAAGTTGGTTGTATATGGCAATGAATTTAACTTAAGTAATCGAGAAAAAACaaac                        <  2:123450/251‑2 (MQ=255)                                                                                                                                                                                                                              gATTGGTGTTTTAATAAAAACATGCAATATTTTGATGATACTAAAAAGATAAAGCATGCATGCTCCGTAATAAATCATATTAATCTTCGCAGCGATCAGTCTAAAATAGCTGAGACATTATTTTTCAATCTCGATAAAGAACCCTATAAAAATAGCCCTGAATTACAGTGGTTGATTTGGAATAAGTTGGTTGAATATGTCAATGAATTTAACTTAAGTAATCGAGAAAAAACAAATTTAATACAAAGG‑Ct         >  2:50147/1‑251 (MQ=255)                                                                                                                                                                                                                                     gtTTTAATAAAAACATGCAATATTTTGATGATAATAAAAAGATAAAGCATGCATGCTACATAATAAATCATATTAATCTTCGCAGCGAACAGTATGAAATAGCTGAGACATAATTTTTCAATCTCGATAAAGAACCATATAAAAATAGCCCTGAATTACAGGGGTTAATTTGTATTAAGTTGGTTGTATATGTCAATGAATTTAACTTAATAAATCGAGAACAAACGAATTTAATA‑AAAGGACTATTTGAt  >  2:276023/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                AAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCATGCAATATTTTGATGATACTAAAAAGATAAAGCATGCATGCTCCGTAATAAATCATATTAATCTTCGCAGCGATCAGTCTAAAATAGCTGAGACATTATTTTTCAATCTCGATAAAGAACCCTATAAAAATAGCCCTGAATTACAGGGGTTGATTTGGAATAAGTTGGTTGTATATGTCAATGAATTTAACTTAAGTAATCGAGAAAAAACAAATTTAATACAAAGG‑CTATTTGAT  >  NC_000913/4476365‑4476843

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