Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 2,955,197 Δ1 bp coding (1034/1071 nt) h16_A2735 → phenylpropionate dioxygenase or related ring‑hydroxylating dioxygenase, large terminal subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604792,955,1950C.100.0% 68.4 / NA 16coding (1032/1071 nt)h16_A2735phenylpropionate dioxygenase or related ring‑hydroxylating dioxygenase, large terminal subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (9/7);  total (9/7)

AGCGCGAGTTCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACGCACGCAAGAACGGACCGC  >  AM260479/2955065‑2955317
                                                                                                                                  |                                                                                                                          
aGCGCGAGTTCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTAccag                                                                                                                         <  1:877166/134‑3 (MQ=255)
          cgtcgAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATg                                                                                                               >  1:13619/1‑133 (MQ=255)
                     aGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATg                                                                                                   <  2:13619/134‑1 (MQ=255)
                                                tCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATcgccgc                                                                        <  2:527508/134‑1 (MQ=255)
                                                tCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATcgccgc                                                                        >  2:1128241/1‑134 (MQ=255)
                                                                             gATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCcggcg                                            >  2:730836/1‑133 (MQ=255)
                                                                                            gcgcTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATcgccg                                                                         <  1:87879/89‑1 (MQ=255)
                                                                                            gcgcTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATcgccg                                                                         >  2:87879/1‑89 (MQ=255)
                                                                                            gcgcTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGc                                                          >  2:25812/1‑104 (MQ=255)
                                                                                            gcgcTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGc                                                          <  1:25812/104‑1 (MQ=255)
                                                                                            gcgcTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAAcc                            >  1:516781/1‑134 (MQ=255)
                                                                                                  ctCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGcac                        <  1:730836/132‑1 (MQ=255)
                                                                                                      aGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTa                                                                               >  2:360388/1‑73 (MQ=255)
                                                                                                      aGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTa                                                                               <  1:360388/73‑1 (MQ=255)
                                                                                                      aGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAacgca                  >  1:788356/1‑134 (MQ=255)
                                                                                                              aCCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACGCACGCAAGaa          <  2:516781/134‑1 (MQ=255)
                                                                                                                      aGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACGCACGCAAGAACGGACCGc  >  1:166647/1‑134 (MQ=255)
                                                                                                                                  |                                                                                                                          
AGCGCGAGTTCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACGCACGCAAGAACGGACCGC  >  AM260479/2955065‑2955317

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