Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,666,371 T→A D322E (GAT→GAA hcaD → phenylpropionate dioxygenase ferredoxin reductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,666,3710TA100.0% 39.5 / NA 11D322E (GAT→GAAhcaDphenylpropionate dioxygenase ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (4/7);  total (4/7)

TGGCAATCACTCGTCTTGATAATGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGT  >  NZ_CP009273/2666232‑2666496
                                                                                                                                           |                                                                                                                             
tgGCAATCACTCGTCTTGATAATGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGG‑‑‑‑‑‑GCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGAACAGTa                                                                                                                          <  1:141316/139‑1 (MQ=255)
      tCACTCGTCTTGATAATGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGAACAGTa                                                                                                                          <  1:377317/139‑1 (MQ=255)
            gTCTTGATAATGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGAACAGTACAGTGa                                                                                                                    <  1:250637/139‑1 (MQ=255)
                  aTAATGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGAACAGTACAGTGATAACt                                                                                                               >  2:208131/1‑138 (MQ=255)
                                                                    gcgcAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGAACAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAAcc                                                            >  2:182258/1‑139 (MQ=255)
                                                                                               gggCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGAACAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTa                                 <  1:232973/139‑1 (MQ=255)
                                                                                                    accgctaccgctACTGCCGCCGCCGTGGTTCTGGAGCGAACAGTACAGTGATAACTTACAGTTTATTGGc                                                                                                 >  2:417429/1‑70 (MQ=255)
                                                                                                    accgctaccgctACTGCCGCCGCCGTGATTCTGGAGCGAACAGTACAGTGATAACTTACAGTTTATTGGc                                                                                                 <  1:417429/70‑1 (MQ=255)
                                                                                                     ccgctaccgctaCTGCCGCCGCCGTGGTTCTGGAGCGAACAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTg                            >  2:134526/1‑138 (MQ=255)
                                                                                                             gctacTGCCGCCGCCGTGGTTCTGGAGCGAACAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGc                   <  2:152975/139‑1 (MQ=255)
                                                                                                                               gTTCTGGAGCGAACAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGt  <  1:182258/138‑1 (MQ=255)
                                                                                                                                           |                                                                                                                             
TGGCAATCACTCGTCTTGATAATGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGT  >  NZ_CP009273/2666232‑2666496

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