Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,060,875	A→G	E143E (GAA→GAG)	yihL →	putative transcriptional regulator YihL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,060,875	0	A	G	100.0%	72.6 / NA	22	E143E (GAA→GAG)	yihL	putative transcriptional regulator YihL
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/6);  total (16/6)

TACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAAC  >  NC_000913/4060746‑4060996                                                                                                                                  |                                                                                                                          tACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAGACa                                                                                                                        >  1:582071/1‑133 (MQ=255)             ccGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAGACATATATTAATCCt                                                                                                            >  1:342256/1‑133 (MQ=255)                      cGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAGACATATATTAATCCTGAAGTTGCt                                                                                                   <  2:582071/133‑1 (MQ=255)                                      cTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAGACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGaac                                                                                   >  1:414934/1‑133 (MQ=255)                                            ccACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAGACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTg                                                                             <  2:415358/133‑1 (MQ=255)                                             cACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAGACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGa                                                                            <  1:135755/133‑1 (MQ=255)                                                       aGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAGACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCa                                                                  >  1:387641/1‑133 (MQ=255)                                                                aTCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAGACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCa                                                         >  2:317931/1‑133 (MQ=255)                                                                        aaCTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAGACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGa                                                 >  1:91745/1‑133 (MQ=255)                                                                             tATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAGACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAGACCACTCGTTTGCTACAGTCTTGGAAAAGt                                            >  1:210406/1‑133 (MQ=255)                                                                             tATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAGACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGt                                            >  2:503155/1‑133 (MQ=255)                                                                                cGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAGACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCt                                         <  1:317931/133‑1 (MQ=255)                                                                                     cACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAGACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCaa                                    <  2:91745/133‑1 (MQ=255)                                                                                     cACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAGACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCaa                                    <  2:414934/133‑1 (MQ=255)                                                                                       ccGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAGACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTTGGAAAAGTGCTACCaaaa                                  >  1:236952/1‑133 (MQ=255)                                                                                       ccGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAGACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCaaaa                                  >  2:554267/1‑133 (MQ=255)                                                                                                 ggCGCTGGAAGGACATAAAGTTTTCTATCACGAGACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGgta                        >  2:622000/1‑133 (MQ=255)                                                                                                    gcTGGAAGGACATAAAGTTTTCTATCACGAGACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGCAAATTGCTACCATAAAGGGACGgtagta                     >  2:228536/1‑133 (MQ=255)                                                                                                     cTGGAAGGACATAAAGTTTTCTATCACGAGACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTaa                    >  2:474087/1‑133 (MQ=255)                                                                                                     cTGGAAGGACATAAAGTTTTCTATCACGAGACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGAGCTACCAAGAAGAAACGGTAGTaa                    >  2:270262/1‑133 (MQ=255)                                                                                                                  aaGTTTTCTATCACGAGACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGAtttt       >  1:318885/1‑133 (MQ=255)                                                                                                                       ttCTATCACGAGACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAAc  >  1:492088/1‑133 (MQ=255)                                                                                                                                  |                                                                                                                          TACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAAC  >  NC_000913/4060746‑4060996

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