New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	4913709 =	22 (0.550)	4 (0.110)	4/316	NT	17.5%	coding (990/1989 nt)	ccmF	holocytochrome c synthetase
	?	AE015451	4913769 =	17 (0.450)					coding (930/1989 nt)	ccmF	holocytochrome c synthetase

GCTGGCGGTGACCGAGAAACGCGGGGTGTTCAAGAGCTGGACCGTGCTGCTGGCAATTGCCGCATTCTCCCTGAGCCTGCTGGGTACCTTCCTGGTGCGTTCCGGAGTGCTGACCTCGGTGCACGCGTTTGCCGCCGACCCGTCGCGGGGTATTTTCATCCTGTTTTTCCTGCTGTTCGTGGTCGGCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/4913898‑4913709 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggtcggcggcAAACGCGTGCACCGAGGTCAGCACTCCGGAACGCACCAGGAAGGTACCCAGC  >  AE015451/4913769‑4913820                                                                                                                                                                                                                                                    GCTGGCGGTGACCGAGAAACGCGGGGTGTTCAAGAGCTGGACCGTGCTGCTGGCAATTGCCGCATTCTCCCTGAGCCTGCTGGGTACCTTCCTGGTGCGTTCCGGAGTGCTGACCTCGGTGCACGCGTTTGCCGCCGACCCGTCGCGGGGTATTTTCATCCTGTTTTTCCTGCTGTTCGTGGTCGG                                                          <  2:435365/186‑1                     GCGGGGTGTGCGAGAGCTGGCCCGTGCTGCTGGCAATTGCCGCATTCTCCCTGAGCCTGCTGGGTACCTTCCTGGTGCGTTCCGGAGTGCTGACCTCGGTGCACGCGTTTGCCGCCGACCCGTCGCGGGGTATTTTCATCCTGTTTTTCCTGCTGTTCGTGGTCGGCGGCAAACGCGTGCACCGAGGTCAGCACTCCGGAACGCACCAGGAAGGTACC      <  1:101173/218‑1                         GGTGTTCAAGAGCTGGACCGTGCTGCTGGCAATTGCCGCATTCTCCCTGAGCCTGCTGGGTACCTTCCTGGTGCGTTCCGGAGTGCTGACCTCGGTGCACGCGTTTGCCGCCGACCCGTCGCGGGGTATTTTCATCCTGTTTTTCCTGCTGTTCGTGGTCGG                                                          <  2:511887/162‑1                                                                          cGCCTGCTGGGTACCTTCCTGGTGCGTTCCGGAGTGCTGACCTCGGTGCACGCGTTTGCCGCCGACCCGTCGCGGGGTATTTTCATCCTGTTTTTCCTGCTGTTCGTGGTCGGCGGCAAACGCGTGC                                            <  1:3333/126‑1                                                                                                    TTCCGGAGTGCTGACCTCGGTGCACGCGTTTGCCGCCGACCCGTCGCGGGGTATTTTCATCCTGTTTTTCCTGCTGTTCGTGGTCGGCGGCAAACGCGTGCACCGAGGTCAGCACTCCGGAACGCACCAGGAAGGTACCCAGC  <  2:522043/143‑1                                                                                                                              CGTTTGCCGCCGACCCGTCGCGGGGTCTGTTCATCCTGTTTTTCCTGCTGTTCGTGGTCGGCGGCAAACGCGTGCACCGAGGTCAGCACTCCGGAACGCACCAGGAAGGTACCCA    <  2:265691/115‑1                                                                                                                                                                                                                                                    GCTGGCGGTGACCGAGAAACGCGGGGTGTTCAAGAGCTGGACCGTGCTGCTGGCAATTGCCGCATTCTCCCTGAGCCTGCTGGGTACCTTCCTGGTGCGTTCCGGAGTGCTGACCTCGGTGCACGCGTTTGCCGCCGACCCGTCGCGGGGTATTTTCATCCTGTTTTTCCTGCTGTTCGTGGTCGGCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/4913898‑4913709 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggtcggcggcAAACGCGTGCACCGAGGTCAGCACTCCGGAACGCACCAGGAAGGTACCCAGC  >  AE015451/4913769‑4913820

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