Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,666,371	T→A	D322E (GAT→GAA)	hcaD →	phenylpropionate dioxygenase ferredoxin reductase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,666,371	0	T	A	100.0%	62.4 / NA	16	D322E (GAT→GAA)	hcaD	phenylpropionate dioxygenase ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (6/10);  total (6/10)

GTCTTGATAATGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCAC  >  NZ_CP009273/2666244‑2666499                                                                                                                                |                                                                                                                                 gTCTTGATAATGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGAACAGTACAGTg                                                                                                                        <  1:615359/138‑1 (MQ=255)  tCTTGATAATGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGAACAGTACAGTGAt                                                                                                                      <  1:729056/139‑1 (MQ=255)                     cacCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGAACAGTACAGTGAt                                                                                                                      >  1:376514/1‑120 (MQ=255)                     cacCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGAACAGTACAGTGAt                                                                                                                      <  2:376514/120‑1 (MQ=255)                             cGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGAACAGTACAGTGATAACTTACAGTTTATTGGCGATATGc                                                                                             <  1:629960/137‑1 (MQ=255)                             cGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGAACAGTACAGTGATAACTTACAGTTTATTGGCGATATGc                                                                                             >  2:629960/1‑137 (MQ=255)                                  gCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGAACAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGa                                                                                      >  2:47691/1‑139 (MQ=255)                                            cGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGAACAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCtt                                                                            >  2:849567/1‑139 (MQ=255)                                                                  gccgctgccgcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGAACAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACt                                                       <  2:690884/138‑1 (MQ=255)                                                                                 ttAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGAACAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGt                                       <  1:849567/139‑1 (MQ=255)                                                                                         accgctaccgctACTGCCGCCGCCGTGGTTCTGGAGCGAACAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTg                               >  1:263263/1‑139 (MQ=255)                                                                                                ccgctacTGCCGCCGCCGTGGTTCTGGAGCGAACAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACg                        <  2:263263/139‑1 (MQ=255)                                                                                                         ccgccgccgTGGTTCTGGAGCGAACAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCtt                <  1:583459/138‑1 (MQ=255)                                                                                                            ccgccgTGGTTCTGGAGCGAACAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCg            >  2:700402/1‑139 (MQ=255)                                                                                                               ccgTGGTTCTGGAGCGAACAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTg         <  1:700402/139‑1 (MQ=255)                                                                                                                       cTGGAGCGAACAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCAc  <  1:707033/138‑1 (MQ=255)                                                                                                                                |                                                                                                                                 GTCTTGATAATGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCAC  >  NZ_CP009273/2666244‑2666499

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