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%	42.8 / NA	15	S390N (AGC→AAC)	yjgL	protein YjgL
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (5/10);  total (5/10)

ACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCATGCAATATTTTGATGATACTAAAAAGATAAAGCATGCATGCTCCGTAATAAATCATATTAATCTTCGCAGCGATCAGTCTAAAATAGCTGAGACATTATTTTTCAATCTCGATAAAGAACCCTATAAAAATAGCCCTGAATTACAGGGGTTGATTTGGAATAAGTTGGTTGTATATGTCAATGAATTTAACTTAAGTAATCGAGAAAAAACAAATTTAATACAA  >  NC_000913/4476375‑4476831                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                   aCTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAACATGCAATATTTTgatgat                                                                                                                                                                                                                 <  2:84757/250‑1 (MQ=255)     gaagaaATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAACATGCAATATTTTGATGATACTaa                                                                                                                                                                                                            <  1:246628/251‑1 (MQ=255)                 aaGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAACATGCAATATTTTGATGATACTAAAAAGATAAAgca                                                                                                                                                                                                <  1:146306/251‑1 (MQ=255)                        aCGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAACATGCAATATTTTGATGATACTAAAAAGATAAAgcatgcat                                                                                                                                                                                           <  2:3719/249‑1 (MQ=255)                                                aTAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAACATGCAATATTTTGATGATACTAAAAAGATAAAGCATGCATGCTCCGTAATAAATCATATTAATCtt                                                                                                                                                                 >  2:67016/1‑251 (MQ=255)                                                   aTGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAACATGCAATATTTTGATGATACTAAAAAGATAAAGCATGCATGCTCCGTAATAAATCATATTAATCTTCGc                                                                                                                                                              <  2:74497/251‑1 (MQ=255)                                                                       gAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAACATGCAATATTTTGATGATACTAAAAAGATAAAGCATGCATGCTCCGTAATAAATCATATTAATCTTCGCAGCGATCAGTCTAAAATAg                                                                                                                                           <  1:67016/250‑1 (MQ=255)                                                                                                                  aTTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAACATGCAATATTTTGATGATACTAAAAAGATAAAGCATGCATGCTCCGTAATAAATCATATTAATCTTCGCAGCGATCAGTCTAAAATAGCTGAGACATTATTTTTCAATCTCGATAAAGAACCCTATAAAAAt                                                                                               <  1:44400/251‑1 (MQ=255)                                                                                                                                gCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAACATGCAATATTTTGATGATACTAAAAAGATAAAGCATGCATGCTCCGTAATAAATCATATTAATCTTCGCAGCGATCAGTCTAAAATAGCTGAGACATTATTTTTCAATCTCGATAAAGAACCCTATAAAAATAGCCCTGAATTAc                                                                                  >  1:186986/1‑250 (MQ=255)                                                                                                                                                      gcaacgattcgattGTGGAGACTTTTTCAAAATCCCCATATATTTATGATGTAAGAATGATGGATTGGTGTCTTAATATAAACATGCAATATTTTGATGATACTAAAAAGATAAATTATGCATGCTGCGTAATATAACATATTAATCGTCGCTGCGATCAGTCTAAAATAGTTGAGTCATTATTGTTCAATGTCGATAATGAGCCCTAGAAAAAGAGCTCAGAATTACAgggtttcatgttgattaagtt                                                            <  2:107583/246‑19 (MQ=255)                                                                                                                                                                           atttttCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAACATGCAATATTTTGATGATACTAAAAAGATAAAGCATGCATGCTCCGTAATAAATCATATTAATCTTCGCAGCGATCAGTCTAAAATAGCTGAGACATTATTTTTCAATCTCGATAAAGAAACCTATAAAAATAGCCCTGAATTACAGGGGTTGATTTGGAATAAGTTGGTTGTATATGTCAAGGAAtgt                                      <  2:232024/250‑3 (MQ=255)                                                                                                                                                                            tttttCAAAATCCCCATATATTCATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAACATGCAATATTTTGATGATACTAAAAAGATAAAGCATGCATGCTCCGTAATAAATCATATTAATCTTCGCAGCGATCAGTCTAAAATAGCTGAGACATTATTTTTCAATCTCGATAAAGAACCCTATAAAAATAGCCCTGAATTACAGGGGTTGATTTGGAATAAGTTGGTTGTATATGTCAATGAAttt                                      >  2:112346/1‑250 (MQ=255)                                                                                                                                                                                                   atgatgTAAGAATACTGGATTGGTGTTTTAATAAAAACATGCAATATTTTGATGATACTAAAAAGATAAAGCATGCATGCTCCGTAATAAATCATATTAATCTTCGCAGCGATCAGTCTAAAATAGCTGAGACATTATTTTTCAATCTCGATAAAGAACCCTATAAAAATAGCCCTGAATTACAGGGGTTGATTTGGAATAAGTTGGTTGTATATGTCAATGAATTTAACTTAAGTAATCGAGAAAAAAc               >  2:198845/1‑250 (MQ=255)                                                                                                                                                                                                            gAATACTGGATTGGTGTTTTAATAAAAACATGCAATATTTTGATGATACTAAAAAGATAAAGCATGCATGCTCCGTAATAAATCATATTAATCTTCGCAGCGATCAGTCTAAAATAGCTGAGACATTATTTTTCAATCTCGATAAAGAACCCTATAAAAATAGCCCTGAATTACAGGGGTTGATTTGGAATAAGTTGGTTGTATATGTCAATGAATTTAACTTAAGTAATCGAGAAAAAACAAATTTAAt      >  1:33982/1‑250 (MQ=255)                                                                                                                                                                                                                aCTGGATTGGTGTTTTAATAAAAACATGCAATATTTTGATGATACTAAAAAGATAAAGCATGCATGCTCCGTAATAAATCATATTAATCTTCGCAGCGATCAGTCTAAAATAGCTGAGACATTATTTTTCAATCTCGATAAAGAACCCTATAAAAATAGCCCTGAATTACAGGGGTTGATTTGGAATAAGTTGGTTGTATATGTCAATGAATTTAACTTAAGTAATCGAGAAAAAACAAATTTAATACaa  <  1:51822/250‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                   ACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCATGCAATATTTTGATGATACTAAAAAGATAAAGCATGCATGCTCCGTAATAAATCATATTAATCTTCGCAGCGATCAGTCTAAAATAGCTGAGACATTATTTTTCAATCTCGATAAAGAACCCTATAAAAATAGCCCTGAATTACAGGGGTTGATTTGGAATAAGTTGGTTGTATATGTCAATGAATTTAACTTAAGTAATCGAGAAAAAACAAATTTAATACAA  >  NC_000913/4476375‑4476831

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