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%	61.2 / NA	18	intergenic (+871/‑758)	yeeJ/shiA	inverse autotransporter adhesin/shikimate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (7/11);  total (7/11)

TGCAATCGATAATGAAAATGTACGTTCCTGCGAGCTACCTTCTAAATTGAACGCAATGTAAGATCCTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAGAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAA  >  NC_000913/2052648‑2053117                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         tGCAATCGATAATGAAAATGTACGTTCCTGCGAGCTACCTTCTAAATTGAACGCAATGTAAGATCCTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACACAGGGGAAAAGCAATTCATgg                                                                                                                                                                                                                             >  2:205111/1‑251 (MQ=255)       cGATAATGAAAATGTACGTTCCTGCGAGCTACCTTCTAAATTGAACGCAATGTAAGATCCTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATAc                                                                                                                                                                                                                       >  2:3675/1‑251 (MQ=255)                      aCGTTCCTGCGAGCTACCTTCTAAATTGAACGCAATGTAAGATCCTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATTGTTATACCAACAACCTGAAAAt                                                                                                                                                                                                        <  1:190902/251‑1 (MQ=255)                        gTTCCTGCGAGCTACCTTCTAAATTGAACGCAATGTAAGATCCTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTc                                                                                                                                                                                                      >  2:61121/1‑251 (MQ=255)                                  gCTACCTTCTAAATTGAACGCAATGTAAGATCCTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCtt                                                                                                                                                                                            <  1:186422/251‑1 (MQ=255)                                          cgaaaTTGAACGCAATGTAAGATGCTAACGATCTCATTTCCTCGCCCACCTTATCGAGTACATCTACTCTTGGTGGCAATTCATATATTGGCCGGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCGCGGAGTTTCTGATCCATACTGGATGCGTGCTGGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTTAAAATTCGTTTTTGCTTACAactga                                                                                                                                                                                    <  2:34224/249‑1 (MQ=255)                                              aTTGAACGCAATGTAAGATCCTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAActgactga                                                                                                                                                                                <  1:162969/251‑1 (MQ=255)                                                                 cctcaCGATAGCATTTCCGCGCGCGCCACAACGAGTACATCCTCTCTTAGTGGCAATTCATATATTGGCCGGACGGCTGGAAAGCCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAAGTCATGGTTATACCAACAACCTGAAAATGCGTTTTTGCTTTCAACTGACTGATAAATAACAGCGTTTTCAg                                                                                                                                                             <  2:6523/247‑1 (MQ=255)                                                                      cGATCGCATTTCCTCGCGCACCTCTTCGAGTACATCCTCACTTATTGGCAATTCATATATTGGCCTGACTGCTGGAAACCCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTTATAAATAACATCGTTTTCAGGTTCt                                                                                                                                                        <  1:112192/251‑1 (MQ=255)                                                                      cGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCt                                                                                                                                                        >  1:82972/1‑251 (MQ=255)                                                                                                cGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCaa                                                                                                                                 >  1:44788/1‑248 (MQ=255)                                                                                                      cATCCTCCCTTAGTGGCAATTCATATATGGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGc                                                                                                                        <  2:97721/251‑1 (MQ=255)                                                                                                                                                                       aaTGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTTCGTGATCCTg                                                       <  1:205111/251‑1 (MQ=255)                                                                                                                                                                                   aTAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATtgt                                           <  2:226938/251‑1 (MQ=255)                                                                                                                                                                                                             aTCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCa                  >  1:196820/1‑250 (MQ=255)                                                                                                                                                                                                               ccATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGt               >  2:16497/1‑251 (MQ=255)                                                                                                                                                                                                                   aCTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATACATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACAACGGGGTCAGGTTCTATTTACTCATCCCCTCAATGCAAAGATCCAGCGTACCGTATTGCTGTGGTATGGTCTTTGCTAAACCTTTTATTTCTTTTAATGTTGCCGGATCTTTCATAGGCATGGGGACTGATGTTAATTGAGTCTGTTCAAGTTTaa           <  1:146958/251‑1 (MQ=255)                                                                                                                                                                                                                            gTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTaa  <  1:204860/251‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         TGCAATCGATAATGAAAATGTACGTTCCTGCGAGCTACCTTCTAAATTGAACGCAATGTAAGATCCTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAGAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAA  >  NC_000913/2052648‑2053117

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