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%	59.9 / NA	23	S522G (AGC→GGC)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/14);  total (9/14)

CCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTT‑GTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCCAC  >  NC_000913/1337150‑1337638                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     cccAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAACAAGCCGTAACTCTGGGCCTCAAGCGGCAACCCTGGGTCAAAGCGTCGCTGGCCCCGGGTTCGAAAGTCGTTTCTGATTATCTGGAAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTT‑GTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAACCGGCCATCAAAAAAGGCGAtt                                                                                                                                                                                                                                                >  2:167875/1‑251 (MQ=255)             gATGGCCGCAGGCTTGCTGGCGCAAAAAGCCGTCACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTGCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTT‑GTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggtg                                                                                                                                                                                                                                    >  1:115358/1‑251 (MQ=255)                 gCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTT‑GTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggtgcggt                                                                                                                                                                                                                                <  1:242640/251‑1 (MQ=255)                                                             gCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGTTTCGAAAGTCGTTGCTGATGACCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGTGTTTAACCTTGTTGGATACGGTT‑GTACCACCTGTATTGGTAACTCGGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTTTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTgg                                                                                                                                                                                    <  1:144633/251‑1 (MQ=255)                                                                 gaaCCATGGGACAAAACGTCGCGGGCCCCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGGGAAACGGACACCGTATCTCGACGATCTGGGGTTTAACCTTGTGGGATACGGTT‑GTACCTCCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTaa                                                                                                                                                                                <  2:229973/250‑1 (MQ=255)                                                                                    cGCTGGCACCGGGTTCTAATGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGAGTCCGTCTCTCGATGGCCTGGGGTTTACCCCTGTGTGATACGGTTTGT‑CCACCTGTATTGGAAACTCTGGGCCGCGGCCCGATCCTATCCAACCGGAAATCAAAAAAGGCGATTTAACCGTCGGTGCTGTGCTGTTCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGTTGAACACTAACTGGCTGGCCTc                                                                                                                                                              <  2:8340/250‑1 (MQ=255)                                                                                    cGCTGGCACCGGGTTCGAAAGTCGGTTCTGATGATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGTTTTAACCTTGTGGGATACGGTT‑GTACCGCCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGACCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTc                                                                                                                                                              <  2:90767/250‑1 (MQ=255)                                                                                                gTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTT‑GTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtggtg                                                                                                                                                 <  2:104567/251‑1 (MQ=255)                                                                                                       aGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTT‑GTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCt                                                                                                                                          <  1:368638/251‑1 (MQ=255)                                                                                                             ttCTGATTATCTGGCACAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACATTGTGGGATACGGTT‑GTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCACAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATgcgc                                                                                                                                    >  2:367992/1‑251 (MQ=255)                                                                                                                                    aaCTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTT‑GTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATccac                                                                                                             >  2:23683/1‑248 (MQ=255)                                                                                                                                      cTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTT‑GTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGGGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCTGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCaaact                                                                                                           >  1:115035/1‑248 (MQ=255)                                                                                                                                            ccGTATCGCGACGAACTGGGGTTTAACCTTGGGGGATACGGTT‑GTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGCTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTc                                                                                                     <  2:47970/251‑1 (MQ=255)                                                                                                                                                       cgaACTGGGGTTTAACCTTGTGGGATACGGTT‑GTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTTAGCCTATCg                                                                                          <  2:105723/251‑1 (MQ=255)                                                                                                                                                         aacgggggTTTAACCTTGTGGGATACGGTT‑GTACCACCTGTATTTGTAACGCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGTCTTCTGAGCCTATCGGc                                                                                        <  2:398747/247‑1 (MQ=255)                                                                                                                                                                gTTTAACCTTGTGGGATGCGGTT‑GTCCCACCTGGATTGTGAACGCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGGCCGGCAACCGTAACGTTGAAGGCCGTATCCTTCCGCTGGTTAAAACTAGCGTTCTGGCCTCGGCGCCGTTGTTGGTGGCCTATTCGCTTGCGGGAAATAGGAATATTAGCCTGGCTGCGCAGCCGATCGGCCATGagc                                                                                 <  2:60448/251‑3 (MQ=255)                                                                                                                                                                    aaCCTTGTGGGATACGGTT‑GTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCaa                                                                             <  1:167875/251‑1 (MQ=255)                                                                                                                                                                             ggATACGGTT‑GTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATcc                                                                    <  1:23683/251‑1 (MQ=255)                                                                                                                                                                                                                           gATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGACGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGCAATATGACTATCAACCCGGCTTCTGAGCCTATCGGCCCTGATCGCAAAGTCGATCCGGTTTATCTGATAGATTTCCCGCCATCGGCACAAGaaaattcccg                       >  2:343475/1‑244 (MQ=255)                                                                                                                                                                                                                                tATCGAAACGGCAATCACAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAat                                                                                                                   >  2:355578/1‑154 (MQ=255)                                                                                                                                                                                                                                tATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAat                                                                                                                   <  1:355578/154‑1 (MQ=255)                                                                                                                                                                                                                                   cGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCCCAAGAAATTTCCCGTGCGGt                 >  2:379508/1‑249 (MQ=255)                                                                                                                                                                                                                                              tCAAAAAAGGCGATTTACCCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGATCCTATCGGCCCTGATCGAAAAGGCGATCCGGGTTATCTGAAAGGTATCTGGCCATAGGCACAAGAAATTGCCCGTGCGGTCGAACAAGtcacc    >  2:381735/1‑248 (MQ=255)                                                                                                                                                                                                                                                aaaaaagGCGATTTAACCGTCGGTGTGGTGCTGTCCGGCAACTGTCACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCcac  <  2:8797/251‑1 (MQ=255)                                                                                                                                                                                                                                                aaaaaagGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCcac  <  1:86413/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     CCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTT‑GTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCCAC  >  NC_000913/1337150‑1337638

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