Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,337,394	A→G	S522G (AGC→GGC)	acnA →	aconitate hydratase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,337,394	0	A	G	100.0%	38.5 / NA	16	S522G (AGC→GGC)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (5/11);  total (5/11)

CTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGAC‑ACCGTATCTCGACGAACTGGG‑GTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTC  >  NC_000913/1337178‑1337633                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                cTGGCGAAAAAAGCCGTAACTCTTGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTAGTTTCTGATTATCTGGCAAAAGCGAAACTGAC‑ACCGTATCTCGACGAACTGGG‑GTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAACGGCGATTTAACCATCGGGGCGGTGCTGTCCGGCa                                                                                                                                                                                                                >  1:119313/1‑250 (MQ=255)                        gggCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGAC‑ACCGTATCTCGACGAACTGGG‑GTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatn                                                                                                                                                                                         <  1:51272/250‑2 (MQ=255)                                                                                 ttCTGATTATCTGGCAAAAGCGAAACTGAC‑GCCGTAGCTCGACGAACTGGG‑GTTTAACCTTGTGGGATACGGTGGTACCACCTGTATGGGTAACGCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATGTAGCCGTCGGTGCGGTGCGGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCt                                                                                                                                                          <  1:403740/224‑1 (MQ=255)                                                                                            tGGCAAAANCNNNANNGAC‑ACCGTATCTCGACGAACTGGG‑GTTTAACCTTGTGGGATACGNTTGTACCACCTGTATNGGNAACTNTGGNCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAAc                                                                                                                                                                                                        >  2:305540/1‑167 (MQ=255)                                                                                            tGGCAAAAGCGAAACTGAC‑TCCGTATCTCGACGAA‑TGGGTGTTTAACCTTGTGGGATACGGTTGTACCTCCTGTATGGGTAACTCTTGGCCGCTGCCCGATCCTATCGTAACGGCATTCAAAAAAGGCGAGTTAACCGTCGGTGAGGTGCTGTCCGGCAACCGTAAc                                                                                                                                                                                                        <  1:305540/167‑1 (MQ=255)                                                                                               cAAAAGCGAAACTGACCACTG‑ATCTCGACGAACTGGG‑GTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCCAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGATGTCCGGCCACCGTACCTTTGAAGGCCGGATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGACGTTGGTGGTTTCTTATGCGCGGGCGGGAAATATg                                                                                                                 <  1:340136/251‑1 (MQ=255)                                                                                                                                    g‑gTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAACGGCGATTTAAc                                                                                                                                                                                                                                       >  1:168221/1‑97 (MQ=255)                                                                                                                                       ggcgatgttGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcgccgccn                                                                                                                                                  <  1:358625/173‑2 (MQ=255)                                                                                                                                                                   ncNGTNTTGGTNACTCNGGGCCGCTGCCCGATNCTATCGAAACGGCAATCAANNANGGCGANTTAAc                                                                                                                                                                                                                                       <  2:168221/66‑1 (MQ=255)                                                                                                                                                                               aCTCTGGGCCGCTGCACGATCCTATCGAAACGGCAATCCAAACAGGCGATTTAACAGTAGGTGCGGAACTGAACGGCAACCGTAACTATGAAGGCCGTATCCATCCGCTGGTTAAACCTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATAAGAATATCACCCCGGCTTCTCAGCCTATCGGCCATGATCg                                                                           >  1:71135/1‑211 (MQ=255)                                                                                                                                                                                                      aTCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGc                                        <  1:262100/223‑1 (MQ=255)                                                                                                                                                                                                          gaactgCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGGCCGGTAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAActggctgg                                                                                                                                                             <  1:262359/97‑1 (MQ=255)                                                                                                                                                                                                          aaaCGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTGACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTAtn                                                                                                                                  <  1:80163/129‑2 (MQ=255)                                                                                                                                                                                                             cGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTTGCGGTAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATc                                                                <  1:245557/192‑1 (MQ=255)                                                                                                                                                                                                              ggCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCg                                                                                                                                                                                                            <  1:330300/51‑1 (MQ=255)                                                                                                                                                                                                                cAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGtc  >  1:371072/1‑251 (MQ=255)                                                                                                                                                                                                                      aaaaagGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGGTTATCTg                                                     >  1:369397/1‑194 (MQ=255)                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                CTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGAC‑ACCGTATCTCGACGAACTGGG‑GTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTC  >  NC_000913/1337178‑1337633

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