New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 863375	32 (0.730)	4 (0.100)	4/324	NT	11.4%	coding (509/1017 nt)	hemH	Ferrochelatase
	?	AE015451	= 863376	32 (0.770)					coding (510/1017 nt)	hemH	Ferrochelatase

GTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/863187‑863375 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgctggcTGCCAGCGCCTCGATGTACGCCGGGTGCTCATAGAACGGCTGCAGCACGCGCATCTGCAACGGTAACTGGTGCTCGCTGACCGTCTGTTTTGCCTGCTCGATCACAGTGGTCACCGTGCTGTCGGCGAACTGTGGATAAAGCGGTGCCAAAGTCACCTTGCGCACACCTTGGGCGGCCAGCCGCGCCAACACATCA  <  AE015451/863376‑863181                                                                                                                                                                                                                                                                                                                                                                                                   GTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCG                                                                                                                                                                                                           >  2:508665/1‑184                                                     TATCCACAGTTCGCCGACAGCACGGTGACCGCTGTGATCGAGCAGGCAAAACAGCCGGTAAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTCTGAGCACCCGGCGTACATCGAGGCG                                                                                                                                                                                                           <  2:310342/132‑1                                                                                      GTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCTGCCAGCGCCTCGATGTACGCCGGGT                                                                                                                                                                            >  2:120948/1‑130                                                                                                                                    TGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCTGCCAGCGCCTCGATGTACGCCGGGTGCTCATAGAACGGCTGCAGCACGCGCATCTGCAACGGTAACTGGTGCTCGCTGACCGTCTGTTTTGCCTGCTCGATCACAGTGGTCACCGTGCTGTCGGCGAACTGTGGATAAAGCGGTGCCAAAGTCACCTTGCGCACACCTTGGGCGGCC                    <  1:120950/236‑1                                                                                                                                              CGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCTGCCAGCGCCTCGATGTACGCCGGGTGCTCATAGAACGGCTGCAGCACGCGCATCTGCAACGGTAACTGGTGCTCGCTGACCGTCT                                                                                                                >  2:78188/1‑134                                                                                                                                                                          GCGTACATCGAGGCGCTGGCTGCCAGCGCCTCGATGTACGCCGGGTGCTCATAGAACGGCTGCAGCACGCGCATCTGCAACGGTAACTGGTGCTCGCTGACCGTCTGTTTTGCCTGCTCGATCACAGTGGTCACCGTGCTGTCGGCGAACTGTGGATAA                                                           >  2:357893/1‑159                                                                                                                                                                                        GCTGGCTGCCAGCGCCTCGATGTACGCCGGGTGCTCATAGAACGGCTGCAGCACGCGCATCTGCAACGGTAACTGGTGCTCGCTGACCGTCTGTTTTGCCTGCTCGATCACAGTGGTCACCGTGCTGTCGGCGAACTGTGGATAAAGCGGTGCCAAAGTCACCTTGCGCACACCTTGGGCGGCCAGCCGCGCCAA         >  2:728931/1‑195                                                                                                                                                                                          TGGCTGCCAGCGCCTCGATGTACGCCGGGTGCTCATAGAACGGCTGCAGCACGCGCATCTGCAACGGTAACTGGTGCTCGCTGACCGTCTGTTTTGCCTGCTCGATCACAGTGGTCACCGTGCTGTCGGCGAACTGTGGAT                                                             >  2:518742/1‑141                                                                                                                                                                                            GCTGCCAGCGCCTCGATGTACGCCGGGTGCTCATAGAACGGCTGCAGCACGCGCATCTGCAACGGTAACTGGTGCTCGCTGACCGTCTGTTTTGCCTGCTCGATCACAGTGGTCACCGTGCTGTCGGCGAACTGTGGATAAAGCGGTGCCAAAGTCACCTTGCGCACACCTTGGGCGGCCAGCCGCGCCAACACATCA  >  2:635474/1‑198                                                                                                                                                                                                                                                                                                                                                                                                   GTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/863187‑863375 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgctggcTGCCAGCGCCTCGATGTACGCCGGGTGCTCATAGAACGGCTGCAGCACGCGCATCTGCAACGGTAACTGGTGCTCGCTGACCGTCTGTTTTGCCTGCTCGATCACAGTGGTCACCGTGCTGTCGGCGAACTGTGGATAAAGCGGTGCCAAAGTCACCTTGCGCACACCTTGGGCGGCCAGCCGCGCCAACACATCA  <  AE015451/863376‑863181

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