Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,672,123	G→A	G26E (GGG→GAG)	hcaD →	phenylpropionate dioxygenase, ferredoxin reductase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,672,123	0	G	A	100.0%	52.7 / NA	18	G26E (GGG→GAG)	hcaD	phenylpropionate dioxygenase, ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/9);  total (9/9)

TTCGCGGCATTCGCCACGTAGCGGCTGGGCTGGATCTTTAAGGAAGCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGT  >  NC_000913/2671997‑2672221                                                                                                                               |                                                                                                   ttCGCGGCATTCGCCACGTAGCGGCTGGGCTGGATCTTTAAGGAAGCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGa                                                                                        <  2:281211/139‑1 (MQ=255)                                              gCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAg                                                                  <  2:422503/116‑1 (MQ=255)                                              gCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAg                                                                  >  1:422503/1‑116 (MQ=255)                                              gCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGa                                                 <  1:321132/133‑1 (MQ=255)                                              gCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGa                                                 >  2:321132/1‑133 (MQ=255)                                               cACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTAt                                                   <  1:301437/130‑1 (MQ=255)                                               cACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTAt                                                   >  2:301437/1‑130 (MQ=255)                                               cACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGaa                                                <  1:480180/133‑1 (MQ=255)                                               cACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGaa                                                >  2:480180/1‑133 (MQ=255)                                               cACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCt                                          >  2:131481/1‑139 (MQ=255)                                               cACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCt                                          >  2:178301/1‑139 (MQ=255)                                                       aagaaaAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACtctc                                  <  2:287419/139‑1 (MQ=255)                                                            aaaaCGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAAt                             >  2:301316/1‑139 (MQ=255)                                                                           gTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATGTTGCTGGAAg              >  1:652034/1‑139 (MQ=255)                                                                            tCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCAt                                                                                 <  1:485166/71‑1 (MQ=255)                                                                            tCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCAt                                                                                 >  2:485166/1‑71 (MQ=255)                                                                                       cAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGt  <  2:652034/139‑1 (MQ=255)                                                                                       cAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGCTTTCCGATGCGCGACATCTTCCTTATGAACGACCTCCTCTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGt  <  1:286849/139‑1 (MQ=255)                                                                                                                               |                                                                                                   TTCGCGGCATTCGCCACGTAGCGGCTGGGCTGGATCTTTAAGGAAGCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGT  >  NC_000913/2671997‑2672221

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