New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	855629 =	114 (1.000)	7 (0.100) +TCCCATTGGCATTTACGAAGTTTAAATGTGCAATTAGAATATATGTATAACAATATT	4/164	12.7	9.9%	intergenic (+160/‑600)	zipA_1/cggR	Cell division protein ZipA/Central glycolytic genes regulator
	?	USA300TCH1516_ALE	2830190 =	103 (0.900)					intergenic (+18/+429)	icaC/lipA_2	putative poly‑beta‑1,6‑N‑acetyl‑D‑glucosamine export protein/Lipase 1
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  USA300TCH1516_ALE/855676‑855629 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AAACACGCGGT‑GGGCGAAATCAGTTTGAATTGACTGACTTCGTTTTACCGCG  >  USA300TCH1516_ALE/2830190‑2830241                                                 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||                                                      GGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCATTGGCATTTACGAAGTTTAAATGTGCAATTAGAATATATGTATAACAATATTAAACACGCGG                                             <  1:291772/115‑1 GGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCATTGGCATTTACGAAGTTTAAATGTGCAATTAGAATATATGTATAACAATATTAAACACGCGG                                             >  2:291772/1‑115    TGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCATTGGCATTTACGAAGTTTAAATGTGCAATTAGAATATATGTATAACAATATTAAACACGCGGTAAAACGAAATCAGTTTGAATTGACTGAC                <  3:295293/141‑1           ACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCATTGGCATTTACGAAGTTTAAATGTGCAATTAGAATATATGTATAACAATATTAAACACGCGGT‑GGGCGAAATCAGTTTGAATTGACTGACTTCGTTTT        <  1:795870/141‑1            CGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCATTGGCATTTACGAAGTTTAAATGTGCAATTAGAATATATGTATAACAATATTAAACACGCGGT‑GGGCGAAATCAGTTTGAATTGACTGACTTCGTTTTA       <  2:739635/141‑1               AATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCATTGGCATTTACGAAGTTTAAATGTGCAATTAGAATATATGTATAACAATATTAAACACGCGGT‑GGGCGAAATCAGTTTGAATTGACTGACTTCGTTTTACCG    <  4:235251/141‑1                 TAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCATTGGCATTTACGAAGTTTAAATGTGCAATTAGAATATATGTATAACAATATTAAACACGCGGT‑GGGCGAAATCAGTTTGAATTGACTGACTTCGTTTTACCGCG  <  2:793136/141‑1                                                 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||                                                      GGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  USA300TCH1516_ALE/855676‑855629 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AAACACGCGGT‑GGGCGAAATCAGTTTGAATTGACTGACTTCGTTTTACCGCG  >  USA300TCH1516_ALE/2830190‑2830241

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