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

TTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAA  >  NZ_CP009273/1454694‑1454857
                                                                                |                                                                                   
ttGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTgcgc                                                                            >  2:340131/1‑90 (MQ=255)
   ctgcGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCaa                                                                         <  2:217298/90‑1 (MQ=255)
                       ggAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTa                                                     <  2:239332/90‑1 (MQ=255)
                         agagCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTAtc                                                   <  2:224436/90‑1 (MQ=255)
                            gCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCAt                                                             <  1:229673/77‑1 (MQ=255)
                            gCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCAt                                                             >  2:229673/1‑77 (MQ=255)
                            gCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCg                                                           <  1:275132/79‑1 (MQ=255)
                            gCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCg                                                           >  2:275132/1‑79 (MQ=255)
                            gCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCa                                                >  1:371702/1‑90 (MQ=255)
                            gCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCa                                                >  1:31284/1‑90 (MQ=255)
                              ggTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAAt                                              >  2:157728/1‑90 (MQ=255)
                                         aTAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCt                                   >  1:344943/1‑90 (MQ=255)
                                             ggCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGc                               <  2:152704/90‑1 (MQ=255)
                                                  ggCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAg                          <  2:201340/90‑1 (MQ=255)
                                                                          agaACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACaaa  <  2:313140/90‑1 (MQ=255)
                                                                             aCGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCa                    <  2:418423/69‑1 (MQ=255)
                                                                             aCGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATg                  >  1:399124/1‑71 (MQ=255)
                                                                             aCGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATg                  <  2:399124/71‑1 (MQ=255)
                                                                                |                                                                                   
TTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAA  >  NZ_CP009273/1454694‑1454857

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