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

CTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGG  >  NC_000913/1337160‑1337594                                                                                                                                                                                                                                           |                                                                                                                                                                                                         cTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcgg                                                                                                                                                                                          >  1:177717/1‑251 (MQ=255)                             aaGCCGAAACTCTTGGCCTCAAGCGGCAACCCTGGGTCAAAGCGACGCCGGCACAGGGATCGAAAGTCGTTTCTGCTTATCTGGCACCAGCGAAACTGACACCGTATCTCGACGAACTGGGGGTTAACCTTGTGGGATACGGATGTACCACCTGTATTGGGAACTCTGGTCCGCTGCCCGATCCTATCGAAACGGCAATCAGAAAAGGCGATTTACCCAGCGGTGCGGTGCTGTCCGGCAAcctgaactgt                                                                                                                                                              >  2:257809/1‑243 (MQ=255)                                                         aCCATGGGGCAAAGCGTCGCTGGCACCGGGGTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGAAACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGAAACTCGGGGCCGCTGCCCGAGACTAGCGAAACGGCAATCAAAAAAGGCGATTTAACCGGCGGTGCGGTGCTGTACGGAAACCGTAACTTTGAAGGCCGTATACATCCGCTGGTTaaaa                                                                                                                                  <  2:143189/251‑1 (MQ=255)                                                                    aaGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAActggctg                                                                                                                       <  1:288806/251‑1 (MQ=255)                                                                                  ccGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCt                                                                                                         >  1:320162/1‑251 (MQ=255)                                                                                          gAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTACCTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCAGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGAGGTTg                                                                                                 >  1:214014/1‑251 (MQ=255)                                                                                            aaGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTc                                                                                                                                                                                            <  1:181764/158‑1 (MQ=255)                                                                                            aaGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTc                                                                                                                                                                                            >  2:181764/1‑158 (MQ=255)                                                                                                 gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTCTgc                                                                                          >  1:309060/1‑251 (MQ=255)                                                                                                                         aaaCTGACACCGTATCTCGACGAACTGGGGTTTTACCTTGTGGGATTCGGTGTTACCACCTGAATTGGTATCTCTGGGCTGCTGTCCGATCCTATCGAACCGGCAATCAAAAAAGGCGATTTAACCTTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATACATCCGCTGGTTAAACCTTACTGGCTGGCCTCGTCGCCGCTGGTGTTTGCCTATGCGCTGGCGGGAAATATGAATATCa                                                                   <  1:257809/250‑1 (MQ=255)                                                                                                                                   cGTATCTCGACGAACTGGGGTTTAAACTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATCTGAATATCAACCTGGCTTCt                                                        >  1:356193/1‑251 (MQ=255)                                                                                                                                               aaCTGGGGTTTAACCTTGTGGGATCCGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTAGCCATCCGCTGGGTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCGGGCTTCTGAGCCTATCGGc                                            <  2:198585/251‑1 (MQ=255)                                                                                                                                                  gggggGTTATCCTTGGGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGCTCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCAt                                         <  2:178400/250‑1 (MQ=255)                                                                                                                                                           aCCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCaaa                                <  2:124378/251‑1 (MQ=255)                                                                                                                                                                gtgGGACACGGTGGTACCACCTGGATGGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCGGGCCTCGCCGCCGCGGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTTTCGGCCATGATCGCAAAGGCg                            <  2:24500/250‑1 (MQ=255)                                                                                                                                                                 tgGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGAt                          >  1:5400/1‑251 (MQ=255)                                                                                                                                                                  gggATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGACGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCATATCGGCCATGATCGCaacagcactc                         >  1:328106/1‑243 (MQ=255)                                                                                                                                                                  gggATAAGGGTGTACCACCTTTATTGGTAACTCGGGGCCGCTGCCCGATCCTATCGACACGGAAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGTCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGAt                          <  2:81123/250‑1 (MQ=255)                                                                                                                                                                          gTTGTACCACCTGGATTGGTAACTCTGGGTCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATAGCAACCTTGCTTCTGAGCCTAGCGGCCATGATCGCCAAGGCGATCCGGTTTAt                 <  2:280505/251‑1 (MQ=255)                                                                                                                                                                                        aTTGGTAACTCGTGGCCGCTGCCCGAGCCTATCGAAACGGCAAGCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCAACCGCTGGTTAAAACTAACTGGCTGGCCGCGCCGCCGCTGGTGGTTGCCTAGGCGCTGGCGGGAAATATGAATAGCAACCTGGCTTCTGAGCCGAGCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCt    <  1:117964/250‑1 (MQ=255)                                                                                                                                                                                          tGGTACCTCTGGGCAGCTGCCCGATCCTATCGAAACGGCAATTAAAAAAGGCGATTTAGCCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTgg  <  2:382768/250‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                         CTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGG  >  NC_000913/1337160‑1337594

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