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

TGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACC  >  NZ_CP009273/1454688‑1454844
                                                                                      |                                                                      
tGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGc                                                                     >  2:322416/1‑90 (MQ=255)
               cAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCt                                                      <  2:59019/90‑1 (MQ=255)
                                    ggTTATGGGTCATAAGGCCTGGGAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAAt                                 <  1:470387/90‑1 (MQ=255)
                                    ggTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAAt                                 <  1:60559/90‑1 (MQ=255)
                                    ggTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCCCCAAAGCGCTTGGCATCGACATTATCTCAAt                                 >  2:212237/1‑90 (MQ=255)
                                           aGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGaa                          <  1:100341/90‑1 (MQ=255)
                                             tCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAgg                        >  2:56955/1‑90 (MQ=255)
                                               aTAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCtt                                                     >  1:462903/1‑59 (MQ=255)
                                               aTAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCtt                                                     <  2:462903/59‑1 (MQ=255)
                                                      cTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGt               >  2:143389/1‑90 (MQ=255)
                                                      cTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGt               >  2:500755/1‑90 (MQ=255)
                                                                   cATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGAcc  <  1:212237/90‑1 (MQ=255)
                                                                                   aCGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGc                  <  2:39301/58‑1 (MQ=255)
                                                                                      |                                                                      
TGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACC  >  NZ_CP009273/1454688‑1454844

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: