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%	79.6 / NA	24	K161E (AAA→GAA)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/11);  total (13/11)

TGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCAGGTTAACCTC  >  NC_000913/1336238‑1336385                                                                          |                                                                           tGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATggg                                                                      >  1:233699/1‑80 (MQ=255)  ggTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATgggg                                                                     <  2:233667/80‑1 (MQ=255)         atgaGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAaagcaag                                                             >  2:233707/1‑81 (MQ=255)              gCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTc                                                        <  1:233668/81‑1 (MQ=255)                   tGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTc                                                    <  1:233669/80‑1 (MQ=255)                         aaaCGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAg                                             <  2:233670/81‑1 (MQ=255)                             gTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGcgtc                                         <  2:233671/81‑1 (MQ=255)                              tACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGcgtcg                                        >  2:233710/1‑81 (MQ=255)                                cGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTg                                      >  2:233709/1‑81 (MQ=255)                                    tGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTgccgc                                  >  1:233706/1‑81 (MQ=255)                                      gAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTgccgcc                                 <  1:233673/80‑1 (MQ=255)                                      gAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCca                                <  2:233672/81‑1 (MQ=255)                                        aaTGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCcagg                              <  1:233674/81‑1 (MQ=255)                                          tGGAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCcaggca                            <  2:233675/81‑1 (MQ=255)                                           ggAGCGCAACCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCcaggca                            >  1:233703/1‑80 (MQ=255)                                                gcAACCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCcaggcacaggca                      >  2:233702/1‑81 (MQ=255)                                                  aaCCACGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCAtt                    >  1:233716/1‑81 (MQ=255)                                                      aCGAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGc                 >  1:233717/1‑80 (MQ=255)                                                        gAACGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCAt              >  1:233711/1‑81 (MQ=255)                                                          aCGTTATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCa            <  1:233676/81‑1 (MQ=255)                                                             ttATGTGTTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCAGGt         <  2:233677/81‑1 (MQ=255)                                                                 gtgtTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCAGGTTAAc     >  2:233698/1‑81 (MQ=255)                                                                  tgtTCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCAGGTTAAcc    >  2:233713/1‑81 (MQ=255)                                                                     tCCTGGAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCAGGTTAACCTc  >  2:233718/1‑80 (MQ=255)                                                                          |                                                                           TGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCAGGTTAACCTC  >  NC_000913/1336238‑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: ‑