New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	= 800928	25 (0.820)	4 (0.140)	4/266	3.2	16.4%	intergenic (+62/‑289)	hisC_1/USA300HOU_RS03915	Histidinol‑phosphate aminotransferase/Putative 5'(3')‑deoxyribonucleotidase
	?	USA300TCH1516_ALE	801948 =	17 (0.590)					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 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccccaacATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCATTAGTCGTTCTTTATC  >  USA300TCH1516_ALE/801948‑802085                                                                                                                                                                                                                                                                                             CAACTGGTGTTAGAATTACAATTGGTTTTAAAGAACAAAATGATAAAATGTTAGAAGTTTTATCAAACTTTAAATACGAATAGTAAGTGGGGAGTGGGACAGAAATGATACTTTCGCAAAATTTATTTCGTCGTCCCACCC                                                                                                                                                >  2:73639/1‑141     TGGTGTTAGAATTACAATTGGTTTTAAAGAACAAAATGATAAAATGTTAGAAGTTTTATCAAACTTTAAATACGAATAGTAAGTGGGGAGTGGGAAAGAAATGATACTTTCGCAAAATTTATTTTGTCGTCCCAACCCAAC                                                                                                                                            >  2:18350/1‑141         GTTAGAATTACAATTGGTTTTAAAGAACAAAATGATAAAATGTTAGCAGTTTTATCAAAAAAAAAATACGAATAGTAAGTGGGGAGTGGGACAGAAATGATACTTTCGCAAAATTTATTTCGTTGTCCCACCCCAACctgc                                                                                                                                        >  2:308024/1‑137                                     AAAATGATAAAATGTTAGAAGTGTTATCAAACTTTAAATACTAATAGTAAGTGGGGAGTGGGACAGAAATGATACTTTCGCAAAATTTATTTCGTCGTCCCACCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTA                                                                                                            <  2:272085/141‑1                                                                                   GTAAGTGGGGAGTGGGACAGAAATGATACTTTCGCAAAATTTATTTCGTCGTCCCACCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACA                                                              >  1:206555/1‑141                                                                                                                             ATTTCGTCGTCCCACCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAAT                    <  2:206555/141‑1                                                                                                                                   TCGTCCCACCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATTATACAGAATCGCTACATAATCAATCATTAa              >  2:240157/1‑140                                                                                                                                               AACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGGAGCGGCACGTAATCAGTCATTAGTCGTTCTTTATC  <  1:69033/141‑1                                                                                                                                                                                                                                                                                             CAACTGGTGTTAGAATTACAATTGGTTTTAAAGAACAAAATGATAAAATGTTAGAAGTTTTATCAAACTTTAAATACGAATAGTAAGTGGGGAGTGGGACAGAAATGATACTTTCGCAAAATTTATTTCGTCGTCCCACCCCAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  USA300TCH1516_ALE/800784‑800928 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccccaacATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCATTAGTCGTTCTTTATC  >  USA300TCH1516_ALE/801948‑802085

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