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_CP0092731,454,7740AG100.0% 42.0 / NA 13D16G (GAT→GGT) paaIhydroxyphenylacetyl‑CoA thioesterase PaaI
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (8/5);  total (8/5)

CTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAA  >  NZ_CP009273/1454697‑1454857
                                                                             |                                                                                   
ctgcGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCaa                                                                         <  1:304784/90‑1 (MQ=255)
  gcgcCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTAGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAg                                                                       >  2:216862/1‑90 (MQ=255)
       aaCGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCtt                                                                  >  2:173608/1‑90 (MQ=255)
                         gCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCa                                                >  2:59955/1‑90 (MQ=255)
                                     cATAAGGCCTGGCAAAATGCCCATGCAATGTTTGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGc                                    <  1:188451/90‑1 (MQ=255)
                                     cATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGc                                    <  1:154824/90‑1 (MQ=255)
                                      aTAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCt                                   >  1:274122/1‑90 (MQ=255)
                                      aTAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCt                                   >  2:180381/1‑90 (MQ=255)
                                             cTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGt                            >  1:236246/1‑90 (MQ=255)
                                                   aaaTGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGAc                      >  2:21649/1‑90 (MQ=255)
                                                   aaaTGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTCGCTGTAGTGAc                      <  2:351277/90‑1 (MQ=255)
                                                               aaTGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCAc          <  2:115638/90‑1 (MQ=255)
                                                                       agaACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACaaa  >  2:174788/1‑90 (MQ=255)
                                                                             |                                                                                   
CTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAA  >  NZ_CP009273/1454697‑1454857

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