New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 3357043	33 (0.980)	4 (0.130)	4/264	4.1	11.1%	coding (2319/4461 nt)	gltB	glutamate synthase, large subunit
	?	NC_000913_3_bme_pgi	= 3357057	33 (1.040)					coding (2333/4461 nt)	gltB	glutamate synthase, large subunit
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCTGTGCTTCCAGGGGGCGGTCAGCCGCATTGGTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_bme_pgi/3356950‑3357043 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cctggctgATGGGCTTACGCGCCAGCCAGGCACGTTTCGACAGATTCAGCAGATCCTGCTGGAAGTCTTCAAAGCTTGCTCCACCAATGCGGCTGACCGCCCCCTGGAAGCACAGGCCCACTACATCATCGTGTAGACCGACCGCTTCAAA  <  NC_000913_3_bme_pgi/3357057‑3356915                                                                                                                                                                                                                                               CCTGTGCTTCCAGGGGGCGGTCAGCCGCATTGGTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGATGGGCTTACGCGCCAGCCAGGCA                                                                                                                         >  1:58350/1‑118 CCTGTGCTTCCAGGGGGCGGTCAGCCGCATTGGTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGATGGGCTTACGCGCCAGCCAGGCA                                                                                                                         <  2:58350/118‑1                                 GTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGATGGGCTTACGCGCCAGCCAGGCACGTTTCGACAGA                                                                                                             <  1:245638/98‑1                                 GTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGATGGGCTTACGCGCCAGCCAGGCACGTTTCGACAGA                                                                                                             >  2:245638/1‑98                                                                                         TGGCTGATGGGCTTACGCGCCAGCCAGGCACGTTTCGACAGATTCAGCAGATCCTGCTGGAAGTCTTCAAAGCTTGCTCCACCAATGCGGCTGACCGCCCCCTGGAAGCACAGGCCCACTACATCATCGTGTAGACCGACCGCTTCAAA  >  1:549625/1‑149                                                                                                                                                                                                                                               CCTGTGCTTCCAGGGGGCGGTCAGCCGCATTGGTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_bme_pgi/3356950‑3357043 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cctggctgATGGGCTTACGCGCCAGCCAGGCACGTTTCGACAGATTCAGCAGATCCTGCTGGAAGTCTTCAAAGCTTGCTCCACCAATGCGGCTGACCGCCCCCTGGAAGCACAGGCCCACTACATCATCGTGTAGACCGACCGCTTCAAA  <  NC_000913_3_bme_pgi/3357057‑3356915

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