New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	= 800928	15 (0.600)	3 (0.130)	3/266	3.0	15.7%	intergenic (+62/‑289)	hisC_1/USA300HOU_RS03915	Histidinol‑phosphate aminotransferase/Putative 5'(3')‑deoxyribonucleotidase
	?	USA300TCH1516_ALE	801948 =	18 (0.760)					intergenic (+189/+435)	USA300HOU_RS03915/USA300HOU_RS03920	Putative 5'(3')‑deoxyribonucleotidase/Putative lipid kinase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CAACTGGTGTTAGAATTACAATTGGTTTTAAAGAACAAAATGATAAAATGTTAGAAGTTTTATCAAACTTTAAATACGAATAGTAAGTGGGGAGTGGGACAGAAATGATACTTTCGCAAAATTTATTTCGTCGTCCCACCCCAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  USA300TCH1516_ALE/800784‑800928 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccccaacATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATT  >  USA300TCH1516_ALE/801948‑802009                                                                                                                                                                                                                 CAACTGGTGTTAGAATTACAATTGGTTTTAAAGAACAAAATGATAAAATGTTAGAAGTTTTATCAAACTTTAAATACGAATAGTAAGTGGGGAGTGGGACAGAAATGATACTTTCGCAAAATTTATTTCGTCGTCCCACCC                                                                    <  2:246923/141‑1             GAATTACAATTGGTTTTAAAGAACAAAATGATAAAATGTTAGAAGTTTTATCAAACTTTAAATACGAATAGTAAGTGGGGAGTGGGACAGAAATGATACTTTCGCAAAATTTATTTCGTCGTCCCACCCCAACATAAAGAA                                                        >  1:278488/1‑141                                     AAAATGATAAAATGTTAGAAGTTTTATCAAACTTTAAATACGAATAGTAAGTGGGGAGTGGGACAGAAATGATACTTTCGCAAAATTTATTTCGTCGTCCCACCCCAACATAAAGAAATACTTTTTCTTTAGAAATTACTA                                >  2:87025/1‑141                                                                   ACTTTAAATACGAATAGTAAGTGGGGAGTGGGACAGAAATGATACTTTCGCAAAATTTATTTCGTCGTCCCACCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATT  >  1:152423/1‑141                                                                                                                                                                                                                 CAACTGGTGTTAGAATTACAATTGGTTTTAAAGAACAAAATGATAAAATGTTAGAAGTTTTATCAAACTTTAAATACGAATAGTAAGTGGGGAGTGGGACAGAAATGATACTTTCGCAAAATTTATTTCGTCGTCCCACCCCAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  USA300TCH1516_ALE/800784‑800928 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccccaacATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATT  >  USA300TCH1516_ALE/801948‑802009

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