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
*	AM260479	2,955,195	0	C	.	100.0%	115.5 / NA	26	coding (1032/1071 nt)	h16_A2735	phenylpropionate dioxygenase or related ring‑hydroxylating dioxygenase, large terminal subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (11/15);  total (11/15)

CGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGC  >  AM260479/2955075‑2955284                                                                                                                         |                                                                                          cgtcgAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATg                                                                              >  1:326123/1‑133 (MQ=255)        gCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACgg                                                                      <  2:88782/134‑1 (MQ=255)         ccGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGc                                                                     <  2:103690/134‑1 (MQ=255)               gcgcCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGCATgcag                                                               <  2:612497/134‑1 (MQ=255)                cgccgcCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGCATg                                                                  >  1:715649/1‑130 (MQ=255)                cgccgcCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGCATg                                                                  <  2:715649/130‑1 (MQ=255)                                aCCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTa                                              <  1:847667/134‑1 (MQ=255)                                   tGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATcg                                           >  2:431740/1‑134 (MQ=255)                                         gAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCgcgc                                     >  2:992243/1‑134 (MQ=255)                                            gacgacGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGCATgc                                                                 >  2:984480/1‑103 (MQ=255)                                            gacgacGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGCATgc                                                                 <  1:984480/103‑1 (MQ=255)                                                               gCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTcc               <  2:286851/134‑1 (MQ=255)                                                                     aTGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGa         <  2:402655/134‑1 (MQ=255)                                                                          ggccggcTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCcgg             <  1:466317/125‑1 (MQ=255)                                                                          ggccggcTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCcgg             >  2:466317/1‑125 (MQ=255)                                                                            ccggcTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGc  <  1:431740/134‑1 (MQ=255)                                                                               gcTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGt                 >  2:741986/1‑116 (MQ=255)                                                                               gcTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGt                 <  1:741986/116‑1 (MQ=255)                                                                                                     aCCAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAg                  >  1:993527/1‑93 (MQ=255)                                                                                                     aCCAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAg                  <  2:993527/93‑1 (MQ=255)                                                                                                       cAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTa                                              <  1:896716/63‑1 (MQ=255)                                                                                                       cAGCGAAGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTa                                              >  2:896716/1‑63 (MQ=255)                                                                                                             aGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTcc               >  2:485657/1‑88 (MQ=255)                                                                                                             aGTCGGGCCGTA‑CCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTcc               <  1:485657/88‑1 (MQ=255)                                                                                                                  ggCCGTA‑CCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCcgg             >  2:857649/1‑85 (MQ=255)                                                                                                                  ggCCGTA‑CCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCcgg             <  1:857649/85‑1 (MQ=255)                                                                                                                         |                                                                                          CGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGC  >  AM260479/2955075‑2955284

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: ‑