New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	1789630 =	12 (0.600)	3 (0.160)	3/238	NT	25.0%	intergenic (‑111/+59)	gapA2/coaE	Glyceraldehyde‑3‑phosphate dehydrogenase 2/Dephospho‑CoA kinase
	?	USA300TCH1516_ALE	1789673 =	7 (0.380)					intergenic (‑154/+16)	gapA2/coaE	Glyceraldehyde‑3‑phosphate dehydrogenase 2/Dephospho‑CoA kinase

AACCTTGAGAGATTGTTAGAAGAAGAAGGTTATATTGAAAAGCCGAATTACGGAGAAGAAGATTAATATTACACTATAAATAAGTCATTACTTTACGTACGCGTTGATGTATGTAAGTAATGACT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  USA300TCH1516_ALE/1789754‑1789630 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aagtaatgactTATTTATAGTGTAATATTAATCTTCTTCTCCGTAATTCGGCTTTTCAATATAACCTTCTTCTTCTAACAATCTCTCAAGGTT  >  USA300TCH1516_ALE/1789673‑1789754                                                                                                                                                                                                                 AACCTTGAGAGATTGTTAGAAGAAGAAGGTTATATTGAAAAGCCGAATTACGAAAAAGAAAATTAATATTACACTATAAATAAGTCATTACTATACGTAAGCGTTGATGTATGTAAGTAATGAATTATTTATAGTGTAATA                                                                    >  1:82008/1‑141                 aAGAAGAAGAAGGTTATATTGAAAAGAAGGATGACGGAGAAGAAGATTAATATTACACTATAAATAAGTCATTACTTTACGTACGCGTTGATGTATGTAAGTAATGACTTATTTATAGTGTAATATTAATCTTCTTCTCCG                                                    <  2:99005/140‑1                                                                   TATTACACTTATAATAAGTCATTTCTTTACGTACGCGTTGATGTATGTAAGTAATGATTTATTTATAGTGTAATATTAATCTTCTTCTCCGTAATTCGGCTTTTCAATATAACCTTCTTCTTCTAACAATCTCTCAAGGTT  <  2:82008/141‑1                                                                                                                                                                                                                 AACCTTGAGAGATTGTTAGAAGAAGAAGGTTATATTGAAAAGCCGAATTACGGAGAAGAAGATTAATATTACACTATAAATAAGTCATTACTTTACGTACGCGTTGATGTATGTAAGTAATGACT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  USA300TCH1516_ALE/1789754‑1789630 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aagtaatgactTATTTATAGTGTAATATTAATCTTCTTCTCCGTAATTCGGCTTTTCAATATAACCTTCTTCTTCTAACAATCTCTCAAGGTT  >  USA300TCH1516_ALE/1789673‑1789754

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