New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 4591278	28 (0.640)	3 (0.070)	3/324	NT	11.0%	coding (349/816 nt)	liuC	methylglutaconyl‑CoA hydratase
	?	AE015451	= 4591409	22 (0.530)					coding (480/816 nt)	liuC	methylglutaconyl‑CoA hydratase

TGGCCATCGACCAGTTGGCCGAAGACGCCAGCCTGCGCTTTGTGCTGCTGCGCGGCCGTGGCCGGCACTTCAGCGCCGGTGCCGACCTGGCCTGGATGCAGCAGTCGGCGCAACTGGACTTCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCGCGCGCCGCTACGCCCTCACCGCCGAGCGTTTCACTGGCGTGCGCGCACGCGAGCTGGGCCTGCTGGCCGAGGTGTACCCGGCCAGCGAACTGGACGACCAC  >  AE015451/4591048‑4591511                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          tGGCCATCGACCAGTTGGCCGAAGACGCCAGCCTGCGCTTTGTGCTGCTGCGCGGCCGTGGCCGGCACTTCAGCGCCGGTGCCGACCTGGCCTGGATGCAGCAGTCGGCGCAACTGGACTTCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCggc                                                                                                                                                                                                                                        <  1:456508‑M1/234‑4 (MQ=255)      aTCGACCAGTTGGCCGAAGACGCCAGCCTGCGCTTTGTGCTGCTGCGCGGCCGTGGCCGGCACTTCAGCGCCGGTGCCGACCTGGCCTGGATGCAGCAGTCGGCGCAACTGGACTTCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTg                                                                                                                                                                                                                                   <  2:685413/234‑1 (MQ=255)      aTCGACCAGTTGGCCGAAGACGCCAGCCTGCGCTTTGTGCTGCTGCGCGGCCGTGGCCGGCACTTCAGCGCCGGTGCCGACCTGGCCTGGATGCAGCAGTCGGCGCAACTGGACTTCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGAt                                                                                                                                                                                                                                 >  1:685415/1‑236 (MQ=255)                       aGACGCCAGCCTGCGCTTTGTGCTGCTGCGCGGCCGTGGCCGGCACTTCAGCGCCGGTGCCGACCTGGCCTGGATGCAGCAGTCGGCGCAACTGGACTTCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATg                                                                                                                                                                                                                 >  1:621589/1‑235 (MQ=255)                                            tctgctgcgCGGCCGTGGCCGGCACTTCAGCGCCGGTGCCGACCTGGCCTGGATGCAGCAGTCGGCGCAACTGGACTTCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCggcgcgttcgccaatggccttgaccacgaacgggctgatgaccgccgg                                                                                                                                                                                           <  2:554149‑M1/235‑49 (MQ=255)                                                ctgcGCGGCCGTGGCCGGCCCTGCAGCGCCGGGGCCGACCTGGCCTGGATGCAGCAGTCGGCGCAACTGGACTTCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTg                                                                                                                                                                                                                                         <  2:92349/186‑1 (MQ=255)                                                             gCCGGCACTTCAGCGCCGGTGCCGACCTGGCCTGGATGCAGCAGTCGGCGCAACTGGACTTCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCggcgcgttcgccaatggccttgaccacgaacgggctgatgaccgccggggccaggccgatgcgta                                                                                                                                                                          <  1:426241‑M1/236‑66 (MQ=255)                                                                                  ccGACCTGGCCTGGATGCAGCAGTCGGCGCAACTGGACTTCAACACCAACCTGGATGATGCTCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAgctgc                                                                                                                                                                                                                          >  2:303294/1‑167 (MQ=255)                                                                                  ccGACCTGGCCTGGATGCAGCAGTCGGCGCAACTGGACTTCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCCTCGGCGcc                                                                                                                                                                                                     >  1:122396/1‑188 (MQ=255)                                                                                    gACCTGGCCTGGCTGCAGCAGTCGGCGCAACTCGACTTCAACCCCACCCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGcccc                                                                                                                                                      <  1:102026/233‑1 (MQ=255)                                                                                               gATGCAGCAGTCGGCGCAACTGGACTTCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGcc                                                                                                                                                                                            <  2:702980/184‑1 (MQ=255)                                                                                                   cagcagTCGGCGCAACTGGACTTCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGt                                                                                                                                     <  1:319410/235‑1 (MQ=255)                                                                                                       agTCGGCGCAACTGGACTTCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTg                                                                                                                                >  1:676028/1‑236 (MQ=255)                                                                                                            gcgcAACTGGACTTCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGtcat                                                                                                                                             >  1:795765/1‑218 (MQ=255)                                                                                                                     gACTTCACCCCCAACCTGGATGATGCCCACGACCTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGt                                                                                                                                                                               <  2:209018/175‑1 (MQ=255)                                                                                                                                  aCCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCgc                                                                                                      <  1:356739‑M2/235‑3 (MQ=255)                                                                                                                                          gatgCCCACGAACTGGGCGAGCTGATGTACGCCCTGCGTCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCGCGCGCCGCTa                                                                                             >  1:541957/1‑236 (MQ=255)                                                                                                                                             gCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCgcgcgcc                                                                                                 >  2:531177‑M2/1‑222 (MQ=255)                                                                                                                                                cACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCCGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCGCGCGCCGCTACGCCCt                                                                                       >  1:610838/1‑236 (MQ=255)                                                                                                                                                     aCTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCGCGCGCCGCTACGCCCTCAcc                                                                                   <  1:275717/235‑1 (MQ=255)                                                                                                                                                                tGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGcct                                                                                                                                                                  >  2:222685/1‑143 (MQ=255)                                                                                                                                                                  aTGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGt                                                                                                                                     <  2:723822/172‑1 (MQ=255)                                                                                                                                                                     tccGCCCTGCATCGCCTCAAGGCCCCGACCCTGGGCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCGCGCGCCGCTACGCCCTCAccgcc                                                                                <  2:683088/220‑1 (MQ=255)                                                                                                                                                                     tACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGt                                                                                                                                                                               <  2:178752/127‑1 (MQ=255)                                                                                                                                                                      aCGCCCTGCATCGCCTCAAGCCGCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGAGCAGCTGCTGCGCCATGGCCATCg                                                                                                                                                                                                          >  2:574774/1‑99 (MQ=255)                                                                                                                                                                               aTCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCcggcggc                                                                                                                                                >  2:183617/1‑147 (MQ=255)                                                                                                                                                                               aTCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGcc                                                                                                                                        >  2:84611/1‑156 (MQ=255)                                                                                                                                                                                tCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCGCGCGCCGCTACGCCCTCACCGCCGAGCGTTTCACTGGCGTgcgcg                                                          >  1:293657/1‑233 (MQ=255)                                                                                                                                                                                         ggCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCGCGCGCCGCTACGCCCTCACCGCCGAGCGTTTCACTGGCGTGCGCGCa                                                        >  1:282802/1‑226 (MQ=255)                                                                                                                                                                                             ccGACCCTGGCCGTGGTGGAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAg                                                                                                                                                                                                 <  2:541549/85‑1 (MQ=255)                                                                                                                                                                                                                         ccTTTGGCGGCGCGCggcgcgttcgccaatggccttgaccacgaacgggctgatgaccgccggggccaggccgatgcgtacttccgacaggcacagttgggcgtcttcggcgcc                                                                                                                                        >  1:36997‑M1/1‑15 (MQ=255)                                                                                                                                                                                                                         ccTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCGCGCGCCGCTACGCCCTCACCGCCGAGc                                                                            >  2:787487/1‑174 (MQ=255)                                                                                                                                                                                                                           ttaggcggcGCGCTGGGACTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCg                                                                                                                                                                <  2:535173/85‑1 (MQ=255)                                                                                                                                                                                                                            ttGGCGGCGCGCTGGGCCTGATCCGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAgg                                                                                                                         <  2:105702/126‑1 (MQ=255)                                                                                                                                                                                                                                 ggcgcgcggcgcgttcgccaatggccttgaccacgaacgggctgatgaccgccggggccaggccgatgcgtacttccgacaggcacagttgggcgtcttcggcgccgatggccatgtcgcagcagctgatcaggcccagcgcgccgccaaaggctgcaccttgcaccacggccagggtcggggccttgaggcgatgcagggcgtacatcagctcgcccagttcgtgggc             >  2:531177‑M1/1‑7 (MQ=255)                                                                                                                                                                                                                                 ggcgCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATGGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCg                                                                                                                >  1:345770/1‑130 (MQ=255)                                                                                                                                                                                                                                 ggcgCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCgcgcgc                                                                                                  <  2:505696‑M2/144‑7 (MQ=255)                                                                                                                                                                                                                                  gcgcgcggcgcgttcgccaatggccttgaccacgaacgggctgatgaccgccggggccaggccgatgcgtacttccgacaggcacagttgggcgtcttcggcgccgatggccatgtcgcagcagctgatcaggcccagcgcgcc                                                                                                 <  2:505696‑M1/144‑139 (MQ=255)                                                                                                                                                                                                                                      gcggcgcgttcgccaatggccttgaccacgaacgggctgatgaccgccggggccaggccgatgcgtacttccgacaggcacagttgggcgtcttcggcgccgatggccatgtcgcagcagctgatcaggcccagcgcgccgccaaaggctgcaccttgcaccacggccagggtcggggccttgaggcgatgcagggcgtacatcagctcgcccagttcgtgggcatcatccaggt  <  1:356739‑M1/235‑234 (MQ=255)                                                                                                                                                                                                                                      gcTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCGCGCGCCGCTACGCCCTCACCGCCGAGCGTTTCACTGGCGTGCGCGCACGCGAGCTGGGCCTGCTGGCCGAGGTGTACCCGGCCAGCGAACTGGACGACCAc  >  1:801780/1‑235 (MQ=255)                                                                                                                                                                                                                                      gcTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGCCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCGCGCGCCGCTACGCCCt                                                                                       >  2:725011/1‑150 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          TGGCCATCGACCAGTTGGCCGAAGACGCCAGCCTGCGCTTTGTGCTGCTGCGCGGCCGTGGCCGGCACTTCAGCGCCGGTGCCGACCTGGCCTGGATGCAGCAGTCGGCGCAACTGGACTTCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCGCGCGCCGCTACGCCCTCACCGCCGAGCGTTTCACTGGCGTGCGCGCACGCGAGCTGGGCCTGCTGGCCGAGGTGTACCCGGCCAGCGAACTGGACGACCAC  >  AE015451/4591048‑4591511

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ 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.