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%	91.4 / NA	27	D16G (GAT→GGT)	paaI	hydroxyphenylacetyl‑CoA thioesterase PaaI
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/15);  total (12/15)

TGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAAT  >  NZ_CP009273/1454688‑1454858                                                                                       |                                                                                     tGTTCATTGCTGCGCCGACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGc                                                                                   >  1:754780/1‑90 (MQ=255)       ttGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTgcgc                                                                             >  2:445762/1‑90 (MQ=255)            gcgcCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAg                                                                        >  2:290069/1‑90 (MQ=255)            gcgcCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAg                                                                        >  1:672603/1‑90 (MQ=255)                 aaCGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCtt                                                                   <  1:863765/90‑1 (MQ=255)                      gcgcTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCAt                                                              <  1:445762/90‑1 (MQ=255)                       cgcTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATc                                                             <  2:48015/90‑1 (MQ=255)                         cTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTgc                                                                               <  1:234196/70‑1 (MQ=255)                         cTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTgc                                                                               >  2:234196/1‑70 (MQ=255)                                   gCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCa                                                 >  1:958580/1‑90 (MQ=255)                                     ggTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAAt                                               <  1:290069/90‑1 (MQ=255)                                               cATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGc                                     <  1:737758/90‑1 (MQ=255)                                               cATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGc                                     <  2:629838/90‑1 (MQ=255)                                                aTAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCt                                    >  2:294028/1‑90 (MQ=255)                                                       cTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGt                             >  2:557554/1‑90 (MQ=255)                                                         ggCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAg                           >  1:582829/1‑90 (MQ=255)                                                         ggCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAg                           <  2:824071/90‑1 (MQ=255)                                                               aTGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCa                     >  1:703927/1‑90 (MQ=255)                                                                    cATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGGTGAAGGCTTTGCTGTAGTGACCATGAcc                <  2:499167/90‑1 (MQ=255)                                                                    cATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGAcc                <  1:576750/90‑1 (MQ=255)                                                                    cATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGAcc                <  2:672603/90‑1 (MQ=255)                                                                       gCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACTGTc             <  1:579944/90‑1 (MQ=255)                                                                          aTGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACt          <  1:212233/90‑1 (MQ=255)                                                                                 agaACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACaaa   <  1:543884/90‑1 (MQ=255)                                                                                  gaACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAAt                                               >  1:961329/1‑45 (MQ=255)                                                                                  gaACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAAt                                               <  2:961329/45‑1 (MQ=255)                                                                                  gaACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAAt  >  2:567050/1‑90 (MQ=255)                                                                                       |                                                                                     TGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAAT  >  NZ_CP009273/1454688‑1454858

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: ‑