Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,336,311	A→G	K161E (AAA→GAA)	acnA →	aconitate hydratase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,336,311	0	A	G	100.0%	70.4 / NA	21	K161E (AAA→GAA)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/9);  total (12/9)

TTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCAGGTTAACCTC  >  NC_000913/1336237‑1336385                                                                           |                                                                           ttGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATggg                                                                      >  2:258378/1‑81 (MQ=255)    gTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGaa                                                                   <  1:258348/81‑1 (MQ=255)    gTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGa                                                                    <  1:258347/80‑1 (MQ=255)          atgaGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAaagcaag                                                             >  1:258374/1‑81 (MQ=255)          atgaGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAaagcaag                                                             >  2:258377/1‑81 (MQ=255)           tgaGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAaagcaagc                                                            <  2:258349/81‑1 (MQ=255)               gCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGtt                                                         >  1:258381/1‑80 (MQ=255)                 aTTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAg                                                      <  2:258350/81‑1 (MQ=255)                     gaagaaAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCgg                                                  >  1:258380/1‑81 (MQ=255)                              gTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGcgtc                                         >  2:258382/1‑81 (MQ=255)                                       gAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTgccgcc                                 <  2:258351/80‑1 (MQ=255)                                         aaTGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCcagg                              >  2:258390/1‑81 (MQ=255)                                          aTGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCcaggc                             <  2:258352/81‑1 (MQ=255)                                             gAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCcaggca                            >  1:258386/1‑79 (MQ=255)                                                      cACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGc                 >  2:258414/1‑81 (MQ=255)                                                          aaCGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCAt              >  1:258409/1‑80 (MQ=255)                                                             gTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCAgg          <  1:258353/81‑1 (MQ=255)                                                              ttATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCAGGt         <  1:258354/81‑1 (MQ=255)                                                                aTGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCAGGtt        >  2:258391/1‑80 (MQ=255)                                                                    gtTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCAGGTTAACCt   >  1:258388/1‑81 (MQ=255)                                                                     ttCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCAGGTTAACCTc  <  1:258355/81‑1 (MQ=255)                                                                           |                                                                           TTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCAGGTTAACCTC  >  NC_000913/1336237‑1336385

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