Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	795,434	G→A	A116T (GCC→ACC)	modA →	molybdate ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	795,434	0	G	A	100.0%	29.7 / NA	12	A116T (GCC→ACC)	modA	molybdate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (3/9);  total (3/9)

GCAGTTTAAAAAAGAGAAAGGCGTGGATGTGGTTTCTTCTTTCGCTTCGTCATCTACTCTCGCCCGTCAGATTGAAGCGGGTGCGCCTGCGGATCTGTTTATTTCTGCCGATCAGAAATGGATGGATTATGCGGTTGATAAAAAAGCGATCGATACAGCTACGCGTCAGACACTGCTCGGCAATAGCCTGGTC‑GTTGTAGCACCGAAAGCCAGCGTGCAGAAAGATTTCACCATCGACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCAAAAACTGGGCGCATGGGATACGCTCTCTCCGAAACTGGCCCCAGCGGAAGATGTTCGTGGGGCGCTGGCGCTGGTC  >  NC_000913/795226‑795637                                                                                                                                                                                                                  |                                                                                                                                                                                                            gCAGTTTAAAAAAGAGAAAGGCGTGGATGTGGTTTCTTCTTTCGCTTCGTCATCTACTCTCGCCCGGCAGATTGAAGCGGGTGCGACTGCGGATCTGTTTATTTCTGCCGATCATAAATGGAGGGATTATGCGTTTTATAAAAAAGCGGGCGATACAGCTACGCGTCAGACACTGCTCGGCAATAGCCTGGTCTGTT‑TAGCACCGAAAACCAGCGTGCAGAAAGATTTCACCATCGACAGCAAAACCAACt                                                                                                                                                                   <  1:81593/251‑1 (MQ=255)          aaaaGAGAAAGGCGTGGATGTGGTTTCTTCTTTCGCTTCGTCATCTACGCTCGCCCGTCAGATTGAAGCGGGTGCGCCTGCGGATCTGTTTATTTCTGCCGATCAGAAATGGATGGATTATGCGGTTGATAAAAAAGCGATCGATACAGCTACGCGTCAGACACTGCTCGGCAATAGCCTGGTC‑GTTGTAGCACCGAAAACCAGCGTGCAGAAAGATTTCACCATCGACAGCAAAACCAACTGGACTTCAc                                                                                                                                                          <  1:129065/251‑1 (MQ=255)                                                                              nnnnnnnnnntGNNNNNCTGNTTATGTCTGCCGATCAGAAATGGATGGATTATGCGGTTGATAAAAAAGCGATCGATACAGCTACGCGTCAGACACTGCTCGGCAATAGCCTGGTC‑GTTGTAGCACCGAAAACCAGCGTGTAGAAAGATTTCAc                                                                                                                                                                                       <  2:107463/144‑1 (MQ=255)                                                                                 gTGCGCCTGCGGATCTGTTTATTTCTGCCGATCAGAAATGGATGGATTATGCGGTTGATAAAAAAGCGATCGATACAGCTACGCGTCAGACACTGCTCGGCAATAGCCTGGTC‑GTTGTAGCACCGAAAACCAGCGTGCAGAAAGATTTCACCATCGACAGCAAAACCAACTGGACTTCACTGCTGAATGGc                                                                                                                                               <  1:100343/191‑1 (MQ=255)                                                                                                              gATCAGAAATGGATGGATTATGCGGTTGATAAAAAAGCGATCGAAACAGCTACGCGTCAGAGACTGCTCGGCAATAGCCTGGTC‑GTTGTAGCGCCGAAAACCAGCGTGCAGAAAGAGTTCACCATCGTCAGCAAATCCAACGGGGCTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATg                                                                                                                <  1:171803/193‑1 (MQ=255)                                                                                                                                                  gCGATCGATACAGCTACGCGTCAGACACTGCTCGGCAATAGCCTGGTC‑GTTGTAGCACCGAAAACCAGCGTGCAGAAAGATTTCACCATCGACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCAAAAACTGGGCGCATGGGATACGCTCTCTCCGAAACTGGCCCCAGCGGAAGATGTTCGTgn                  <  1:35231/251‑2 (MQ=255)                                                                                                                                                    gatcgatACAGCTACGCGTCAGACACTGCTCGGCAATAGCCTGGTC‑GTTGTAGCACAGAAAACCAGCGTGCCGAAAGATTTCACCATCGACATCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATc                                                                                                                         >  1:111136/1‑146 (MQ=255)                                                                                                                                                                              ngNTCNGCAATNGCCTNGTC‑GTTGTAGCACCGNAAACCAGCGTGGAGAAGGANNTNAGCATNGACAg                                                                                                                                                                              <  2:65510/66‑1 (MQ=255)                                                                                                                                                                                                 c‑GTTGTAGCACCGAAAACCAGCGTGCAGAAAGATTTCACCATCGACAGCAAAACCAACTGGACTTCACTGCTGACTGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCAAAAACTGGGCGCATGGGATACGCTCTCTCCGAAACTGGCCCCAGCGGAAGATGTTCGTGGGGCGCTGGCGCTGGTc  >  1:63190/1‑220 (MQ=255)                                                                                                                                                                                                        aGCACCGAAAACCAGCGTGCAGAAAGATTTCACCATCGACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCAAAAACTGGGCGCATGGGATACGCTCTCTCCGAAACTGGCCCCAGCGGAAGATGTTCGTgggg                <  1:126106/200‑1 (MQ=255)                                                                                                                                                                                                         gCACCGAAAACCAGCGTGCAGAAAGATTTCACCATCGACAGCAAAACCAACTGGACTTc                                                                                                                                                            >  1:8880/1‑59 (MQ=255)                                                                                                                                                                                                         gCACCGAAAACCAGCGTGCAGAAAGATTTCACCATCGACAGCAAAACCAACTGGACTTc                                                                                                                                                            <  2:8880/59‑1 (MQ=255)                                                                                                                                                                                                                  |                                                                                                                                                                                                            GCAGTTTAAAAAAGAGAAAGGCGTGGATGTGGTTTCTTCTTTCGCTTCGTCATCTACTCTCGCCCGTCAGATTGAAGCGGGTGCGCCTGCGGATCTGTTTATTTCTGCCGATCAGAAATGGATGGATTATGCGGTTGATAAAAAAGCGATCGATACAGCTACGCGTCAGACACTGCTCGGCAATAGCCTGGTC‑GTTGTAGCACCGAAAGCCAGCGTGCAGAAAGATTTCACCATCGACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCAAAAACTGGGCGCATGGGATACGCTCTCTCCGAAACTGGCCCCAGCGGAAGATGTTCGTGGGGCGCTGGCGCTGGTC  >  NC_000913/795226‑795637

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