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

GCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACT  >  NZ_CP009273/1454699‑1454850                                                                            |                                                                             gcgcCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACCGTGCCTGCGCCAAAg                                                                >  1:715544/1‑90 (MQ=255)      aaCGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCtt                                                           >  1:580624/1‑90 (MQ=255)            cgcTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACCGTGCCTGCGCCAAAGCGCTTGGCATc                                                     <  2:715544/90‑1 (MQ=255)             gcTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCg                                                    <  1:508959/90‑1 (MQ=255)             gcTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCg                                                    <  1:236178/90‑1 (MQ=255)                      gagCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTAtct                                           >  1:623333/1‑90 (MQ=255)                        gCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTAtc                                            <  2:515794/87‑1 (MQ=255)                        gCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTAtc                                            >  1:515794/1‑87 (MQ=255)                         cGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCaa                                        <  1:240980/90‑1 (MQ=255)                          ggTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAAt                                       <  1:474258/90‑1 (MQ=255)                                 aGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGaa                                <  1:82593/90‑1 (MQ=255)                                    cATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGc                             <  1:509466/90‑1 (MQ=255)                                    cATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGc                             <  1:349301/90‑1 (MQ=255)                                    cATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGc                             <  2:95724/90‑1 (MQ=255)                                     aTAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCt                            >  2:236528/1‑90 (MQ=255)                                     aTAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCt                            >  2:583752/1‑90 (MQ=255)                                         ggCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGc                        <  2:623333/90‑1 (MQ=255)                                                         cATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTCGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGAcc        <  2:367701/90‑1 (MQ=255)                                                            gCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTc     <  1:583752/90‑1 (MQ=255)                                                              aaTGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCAc   <  2:518176/90‑1 (MQ=255)                                                              aaTGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCAc   <  2:580624/90‑1 (MQ=255)                                                               aTGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACt  <  1:96240/90‑1 (MQ=255)                                                               aTGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACt  <  1:497790/90‑1 (MQ=255)                                                                         aCGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGc                             <  2:330782/53‑1 (MQ=255)                                                                            |                                                                             GCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACT  >  NZ_CP009273/1454699‑1454850

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