Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	AE015451	3,718,088	Δ28 bp	6.2%	coding (478‑505/792 nt)	paaG ←	1,2‑epoxyphenylacetyl‑CoA isomerase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 3718087	44 (0.620)	3 (0.040)	3/350	NT	6.2%	coding (506/792 nt)	paaG	1,2‑epoxyphenylacetyl‑CoA isomerase
	?	AE015451	3718116 =	47 (0.660)					coding (477/792 nt)	paaG	1,2‑epoxyphenylacetyl‑CoA isomerase

GGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAATGTTGGCACCGGCACCGGCCGCCACGCCGTTGACCGCACATATCACCGGCAAAGGCAGGTCGCGCAGGGTGCGCA  >  AE015451/3717857‑3718315                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                     ggTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGcc                                                                                                                                                                                                                                  >  1:642248/1‑235 (MQ=255)                  ccTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGcc                                                                                                                                                                                                                        >  2:1160753/1‑228 (MQ=255)                                    tCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGAc                                                                                                                                                                                              >  1:276892/1‑236 (MQ=255)                                                             cAAAACCGTTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCaa                                                                                                                                                                                           <  1:1114387/214‑1 (MQ=255)                                                              aaaaCCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTAccg                                                                                                                                                                     <  2:766246/235‑1 (MQ=255)                                                                        tCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGGGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGc                                                                                                                                                                                                             >  1:433017/1‑185 (MQ=255)                                                                                   agtGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATg                                                                                                                                                                                                                           <  2:279866/158‑1 (MQ=255)                                                                                      gAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAgc                                                                                                                                                                                                                     <  1:1237633/163‑1 (MQ=255)                                                                                                       gtgtCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGc                                                                                                                                                                                                                         <  2:418769/142‑1 (MQ=255)                                                                                                       gtgtCAGGCCGTAGGTGGGCCGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGa                                                                                                                                               <  1:513957/216‑1 (MQ=255)                                                                                                                cGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGc                                                                                                                   >  1:1096004/1‑235 (MQ=255)                                                                                                                    ggTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAg                                                                                                               <  1:1160758/235‑1 (MQ=255)                                                                                                                     gTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCCCCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAgcg                                                                                                             >  1:333906/1‑236 (MQ=255)                                                                                                                     gTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGAt                                                                                                                                              >  2:787964/1‑203 (MQ=255)                                                                                                                     gTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGACCTTGCAGAAGGCCTGGATGAAACTGGcc                                                                                                                  >  2:564335/1‑231 (MQ=255)                                                                                                                               cggcgAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAg                                                                                                    >  1:687441/1‑235 (MQ=255)                                                                                                                                ggcgAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTAccgccg                                                                                                                                                                  >  2:941319/1‑172 (MQ=255)                                                                                                                                ggcgAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGCTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCa                                                                                                  >  1:705360/1‑236 (MQ=255)                                                                                                                                       tGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGc                                                                                                                                                                                                                              >  2:632258/1‑105 (MQ=255)                                                                                                                                                aGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCt                                                                                                                                                                                                                                      <  2:634376/88‑1 (MQ=255)                                                                                                                                                      aGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAAt                                                                            <  1:629862/236‑1 (MQ=255)                                                                                                                                                               tcaCGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGAGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCgg                                                                                                                                                                 >  2:693688/1‑142 (MQ=255)                                                                                                                                                                caCGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAgg                                                                                             >  2:1000254/1‑209 (MQ=255)                                                                                                                                                                         cTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAACGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCa                                                                                                                                                     <  2:717059/144‑1 (MQ=255)                                                                                                                                                                         cTGCATCACCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAgcg                                                                                                                                                                                                                    <  2:953305/81‑1 (MQ=255)                                                                                                                                                                          tGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCg                                                                                          >  2:1245892/1‑202 (MQ=255)                                                                                                                                                                          tGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAATGTTGGCACCGGCACCGgccg                                                        >  1:568135/1‑236 (MQ=255)                                                                                                                                                                            catcatCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAATGTTGGCACCGGCACCGgccgc                                                       >  1:266028/1‑235 (MQ=255)                                                                                                                                                                                    accacGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTggccatgccgaccaagcgcggcagcagccaagtaccgccggagtccggcaccaggccgatcttgcagaaggcctggatgaaactg                                                                                                                                                 >  2:811510‑M1/1‑52 (MQ=255)                                                                                                                                                                                    accacGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGc                                                                                                                                         >  2:1027094/1‑145 (MQ=255)                                                                                                                                                                                    accacGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAATGTTGGCACCGGCACCGGCCGCCACCCCGt                                               >  1:852914/1‑235 (MQ=255)                                                                                                                                                                                    accacGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGAACCGCCGGa                                                                                                                                                                >  2:242059/1‑122 (MQ=255)                                                                                                                                                                                       acGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAATGTTGGCACCGGCACCGGCCGCCACGCCGTTGAc                                           >  1:292660/1‑236 (MQ=255)                                                                                                                                                                                           ggTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGCCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAATGTTGGCACCGGCACCGGCCGCCACGCCGTTGACCGca                                       >  1:385512/1‑236 (MQ=255)                                                                                                                                                                                              aGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTggccatgccgaccaagcgcggcagcagccaagtaccgccggagtccggcaccaggccgatcttgcagaaggcctggatgaaactggccgagcgggcggccagcaccagg                                                                                                                         <  2:977203‑M1/151‑110 (MQ=255)                                                                                                                                                                                              aGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTggccatgccgaccaagcgcggcagcagccaagtaccgccggagtccggcaccaggccgatcttgcagaaggcctggatgaaactggccgagcgggcggccagcaccagg                                                                                                                         >  1:977208‑M1/1‑42 (MQ=255)                                                                                                                                                                                                   aGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAgc                                                                                                              <  2:226493/157‑1 (MQ=255)                                                                                                                                                                                                        ccATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGTCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAATGTTGGCACCGGCACCGGCCGCCACGCCGTTGACCGCACATATCACCGGCa                          >  1:629380/1‑236 (MQ=255)                                                                                                                                                                                                         cATGGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGa                                                                                                                                                                <  2:40552/101‑1 (MQ=255)                                                                                                                                                                                                          aTTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGaa                                                                                                                         >  1:1235857/1‑139 (MQ=255)                                                                                                                                                                                                          aTTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAATGTTGGCACCGGCACCGGCCGCCACGCCGTTGACCGCACATATCACCGGCaa                         >  1:592492/1‑235 (MQ=255)                                                                                                                                                                                                          aTTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAATGTTGGCACCGGCACCGGCCGCCACGCCGTGGACCg                                         >  1:316743/1‑219 (MQ=255)                                                                                                                                                                                                              cTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAATGTTGGCACCGGCACCGGCCGCCACGCCGTTGACCGCACATATCACCGGCa                          >  2:802469/1‑230 (MQ=255)                                                                                                                                                                                                                 ggCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAgc                                                                                                              <  2:276888/143‑1 (MQ=255)                                                                                                                                                                                                                    cTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAATGTTGGCACCGGCACCGGCCGCCACGCCGTTGACCGCACATATCACCGGCAAAGGCAGGTcgc              <  1:787969/236‑1 (MQ=255)                                                                                                                                                                                                                      gTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCa                                                                                                     >  2:490558/1‑147 (MQ=255)                                                                                                                                                                                                                         cGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAAt                                                                            >  2:205938/1‑169 (MQ=255)                                                                                                                                                                                                                                aaGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAATGTTGGCACCGGCACCGGCCGCCACGCCGTTGACCGCACATATCACCGGCAAAGGCAGGTCGCGCAGGGTGCGCa  >  1:1204003/1‑236 (MQ=255)                                                                                                                                                                                                                                aaGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGAAGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAATGTTGGCACCGGCACCGGCCGCCAc                                                    >  1:733401/1‑186 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                     GGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAATGTTGGCACCGGCACCGGCCGCCACGCCGTTGACCGCACATATCACCGGCAAAGGCAGGTCGCGCAGGGTGCGCA  >  AE015451/3717857‑3718315

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.