Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,052,885	G→A	intergenic (+871/‑758)	yeeJ → / → shiA	inverse autotransporter adhesin/shikimate:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,052,885	0	G	A	100.0%	75.7 / NA	23	intergenic (+871/‑758)	yeeJ/shiA	inverse autotransporter adhesin/shikimate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/13);  total (10/13)

GCAATCGATAATGAAAATGTACGTTCCTGCGAGCTACCTTCTAAATTGAACGCAATGTAAGATCCTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAGAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCT‑TTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTC  >  NC_000913/2052649‑2053124                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 gCAATCGATAATGAAAATGTACGTTCCTGCGAGCTACCTTCTAAATTGAACGCAATGTAAGATCCTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGt                                                                                                                                                                                                                                    >  1:198856/1‑251 (MQ=255)             gAAAATGTACGTTCCTGCGAGCTACCTTCTAAATTGAACGCAATGTAAGATCCTAACGATCGCATTTCCTCGCGCTCCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACcaacaa                                                                                                                                                                                                                         <  2:22367/250‑1 (MQ=255)                                aGCTACCTTCTAAATTGAACGCAATGTAAGATCCTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCAAGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCt                                                                                                                                                                                                     >  1:332864/1‑251 (MQ=255)                                           aaaTTGAACGCAATGTAAGATCCTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAactg                                                                                                                                                                                             <  1:137218/248‑1 (MQ=255)                                                gAACGCAATGTAAGATCCTAACGAGCGCATGTACTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCGTAAACGAAGGGGAAAAGCAATTCATGGTTATACcaaca                                                                                                                                                                                                                          >  2:146333/1‑214 (MQ=255)                                                 aaCGCAATGTAAGATCCTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGataa                                                                                                                                                                                     >  1:430220/1‑250 (MQ=255)                                                        tGTAAGATCCTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACCAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATc                                                                                                                                                                             >  2:139833/1‑251 (MQ=255)                                                         gTAAGATCCTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCg                                                                                                                                                                            >  1:189834/1‑251 (MQ=255)                                                                ccTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAg                                                                                                                                                                     <  1:139833/251‑1 (MQ=255)                                                                    aCGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTc                                                                                                                                                                 >  1:111811/1‑251 (MQ=255)                                                                     ggATCGCATTTCCTCGCGCCCCTCTTCGAGTACATCCTCACTTAGTGGCAATGCATATATTGGCCTGACTGGTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTGTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGTAAAAGCAAATCAATGTTATACCAACGACCTGTAAAGTCGTTTTTGCTTTCAAGTGACGGATAAATAACATCGTTTTCTGGTTCt                                                                                                                                                                <  2:206415/250‑1 (MQ=255)                                                                                     cgcACCTCAGCGAGTACATCCTCACTTAGTGGCAATTCATATATTGTCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCT‑TTTTACGCATCCCCTc                                                                                                                                               <  1:529735/251‑1 (MQ=255)                                                                                             aTCGAGTACATCCTCGCGTAGTGGAAATTCATATATTGGCCTGCCGGCTGGAAACCCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCT‑TTTTACGCATCCCCTCAATGCaaa                                                                                                                                       <  2:459203/251‑1 (MQ=255)                                                                                                    aCATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCT‑TTTTACGCATCCCCTCAATGCAAAGATCCg                                                                                                                                 <  1:500827/250‑1 (MQ=255)                                                                                                        ccTCACTTAGTGGCAATTCATAAATGTGCCTGACTGCTGGAAAACCCTCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCT‑TTTTACGCATCCCCTCAATGCAAAGATCCGGCGTa                                                                                                                            <  2:332864/251‑1 (MQ=255)                                                                                                                     caatggacatatTGGCCGGACGGCAGGAAAACCCGCCTCTCGCATCATAGATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGAGCCGGACGGGACGCGTGCTTTCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCGTAAAATTCGTTTCTGCTTTGAACGGACGGATAAATAACAGCGGTTTCATGTTCTGTGTTGC‑CAGCCCCGCAGTGCAAAGCTCCGGCGTTCCGTGTTGCtgtg                                                                                                               <  2:218834/243‑1 (MQ=255)                                                                                                                          cATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAACTGCCCCTGTCATAGGTACGGGAGCCAGGAGTTTCTGATCAATAATGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGCCTGATAACTAACATCGTTTTCAGGTTCT‑TTTTACGCATCCCATCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATg                                                                                                          >  1:29713/1‑251 (MQ=255)                                                                                                                                           gcgggAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACATCCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCT‑TTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCttt                                                                                         <  1:539302/248‑1 (MQ=255)                                                                                                                                                                atcatAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCT‑TTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTg                                                                    >  1:521187/1‑251 (MQ=255)                                                                                                                                                                   atAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCT‑TTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATc                                                                 <  1:312960/251‑1 (MQ=255)                                                                                                                                                                         gCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCT‑TTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCa                                                            <  1:93432/250‑1 (MQ=255)                                                                                                                                                                                                    agtctCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCT‑TTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCGTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATGGTGACTAATGTTAATTTAGt                                <  2:548924/247‑1 (MQ=255)                                                                                                                                                                                                                                  cacaAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCT‑TTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTc  >  2:287541/1‑251 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 GCAATCGATAATGAAAATGTACGTTCCTGCGAGCTACCTTCTAAATTGAACGCAATGTAAGATCCTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAGAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCT‑TTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTC  >  NC_000913/2052649‑2053124

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