Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,454,774	A→G	D16G (GAT→GGT)	paaI →	hydroxyphenylacetyl‑CoA thioesterase PaaI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,454,774	0	A	G	100.0%	51.8 / NA	16	D16G (GAT→GGT)	paaI	hydroxyphenylacetyl‑CoA thioesterase PaaI
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/7);  total (9/7)

TGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAA  >  NZ_CP009273/1454688‑1454857                                                                                       |                                                                                    tGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGc                                                                                  >  2:90237/1‑90 (MQ=255)          ctgcGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCaa                                                                         >  2:126495/1‑90 (MQ=255)            gcgcCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAg                                                                       >  1:287956/1‑90 (MQ=255)            gcgcCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAg                                                                       >  2:231078/1‑90 (MQ=255)                                 gagCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTAtct                                                  >  2:340115/1‑90 (MQ=255)                                     ggTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAAt                                              >  2:126981/1‑90 (MQ=255)                                     ggTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAAt                                              >  2:179666/1‑90 (MQ=255)                                     ggTTATGAGTCATAAGGCCTGGCAAAATGCCCAAGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAAt                                              <  2:294405/90‑1 (MQ=255)                                               cATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGc                                    <  1:69551/90‑1 (MQ=255)                                                aTAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCt                                   >  1:315154/1‑90 (MQ=255)                                                 tAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCtt                                  <  1:126495/90‑1 (MQ=255)                                                           cAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTg                        <  1:340115/90‑1 (MQ=255)                                                               aTGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCa                    <  1:64991/90‑1 (MQ=255)                                                                          aTGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACt         >  2:340764/1‑90 (MQ=255)                                                                                 agaACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACaaa  <  1:90237/90‑1 (MQ=255)                                                                                 agaACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACaaa  <  2:315906/90‑1 (MQ=255)                                                                                       |                                                                                    TGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAA  >  NZ_CP009273/1454688‑1454857

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