New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	1965737 =	103 (1.220)	3 (0.040)	3/246	NT	2.9%	coding (71/924 nt)	hemH	Ferrochelatase
	?	USA300TCH1516_ALE	1965782 =	101 (1.260)					coding (26/924 nt)	hemH	Ferrochelatase
Rejected: Frequency below/above cutoff threshold.

CACATATACACAAAGATAATCATAGTTTTACATTGAAGCACATGAAAGATAATTCATTTAAAGGAAGGTATTATCAATGACTAAAAAAATGGGATTATTAGTTATGGCTTATGGCACACCTTATAAAGAAAGTGACATAGAGCCATA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  USA300TCH1516_ALE/1965883‑1965737 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agccataACTAATAATCCCATTTTTTTAGTCATTGATAATACCTTCCTTTAAATGAATTATCTTTCATGTGCTTCAATGTAAAACTATGATTATCTTTGTGTATATGTGTGTACGAATTCGCTTACTTTACGTAACGTCTCTGGTTG  >  USA300TCH1516_ALE/1965782‑1965921                                                                                                                                                                                                                                                                                                 CACATATACACAAAGATAATCATAGTTTTACATTGAAGCACATGAAAGATAATTCATTTAAAGGAAGGTATTATCAATGACTAAAAAAATGGGATTATTAGTTATGGCTTATGGCACACCTTATAAAGAAAGTGACATAGA                                                                                                                                                    <  1:391186/141‑1  ACATATACACAAAGATAATCATAGTTTTACATTGAAGCACATGAAAGATAATTCATTTAAAGGAAGGTATTATCAATGACTAAAAAAATGGGATTCTTCGTCAAGGCTTATGGCAAACCATATAAAGAAAGTGACATAGAG                                                                                                                                                   >  2:727375/1‑141  ACATATACACAAAGATAATAATAGTTATACATTGAAGCAAATGAAAGATAATTCATTTAAAAGAAGGTATTATCAATGACTAAAAAAATGGGATTATTAGTTATGGCTTATGGCACACCTTATAAAGAAAGTGACATAGAG                                                                                                                                                   <  1:547674/141‑1   CATATAGACAGATATAATCATAGTTTTACATTGAAGCACATGAAAGATAATTCATTTAAAGGAAGGTATTATCAATGACTAAAAAAATGGGATTATTAGTTATGGCTTATGGCACACCTTATAAAGAAAGTGACATAGAGC                                                                                                                                                  <  1:381603/141‑1    ATATACAAAAATTAAATCATAGTTTTAAATTGAAGCAAATGAAAGATAATTCATTTAAAGGAAGGTATTATCAATGACTAAAAAAATGGGATTATTAGTTATGGCTTATGGCACACCTTATAAAGAAAGTGACATAGAGCC                                                                                                                                                 <  2:601117/141‑1       TACACAAAGATAATCATAGTTTTACATTGAAGCACATGAAAGATAATTCATTTAAAGGAAGGTATTATCAATGACTAAAAAAATAGGATTATTAGTTATGGCTTATGGCACACCTTATAAAAAACGTAACATAGAGCCATA                                                                                                                                              >  2:664709/1‑141                   ATCATAGTTTTACATTGAAGCACATGAAAGATAATTCATTTAAAGGAAGGTATTATCAATGACTAAAAAAATGGGATTATTAGTTATGGCTTATGGCACACCTTATAAAGAAAGTGACATAGAGCCATAACTAATAATCCC                                                                                                                                  >  1:251614/1‑141                                                                                           GGGATTATTAGTTATGGCTTATGGCACACCTTATAAAGAAAGTGACATAGAGCCATAACTAATAATCCCATTTTTTTAGTCATTGATAATACCTTCATTTAAATGAATTATCTTTCATGTGCTTCAATGTAAAACTATGAT                                                          <  2:251614/141‑1                                                                                                               ATGGCACACCTTATAAAGAAAGTGACATAGAGCCATAACTAATAATCCCATTTTTTTAGTCATTGATAATACCTTCCTTTAAATGAATTATCTTTCATGTGCTTCAATGTAAAACTATGATTATCTTTGTGTATATGTGTG                                      >  2:554521/1‑141                                                                                                                                                 ATAACTAATAATCCCATTTTTTTAGTCATTGATAATACCTTCCTTTAAATGAATTATCTTTCATGTGCTTCAATGTAAAACTATGATTATCTTTGTGTTTGTTTGTGTACGAATTCGCTTACTTTACGTAACGTCTCTGGT    >  1:34342/1‑141                                                                                                                                                  TAACTAATAATCCCATTTTTTTAGTCATTGATAATACCTTCCTTTAAATGAATTATCTTTCATGTGCTTCAATGTAAAACTATGATTATCTTTGTGTATATTTGTGTTCAAATTTGCTTACTTTACGTACTGTCTCTGTTT   >  2:960375/1‑141                                                                                                                                                  TAACTAATAATCCCATTTTTTTAGTCATTGATAATACCTTCCTTTAAATGAATTATCTTTCATGTACTTCAATGTACAACTATGATTATCTTTGTGTATATGTGTTTACGTCTTAGTTTACTTTATGTAACGTCTGTGGTT   >  2:639489/1‑141                                                                                                                                                   AACTAATAATCCCATTTTTTTAGTCATTGATAATACCTTCCTTTAAATGAATTATCTTTCATGTGCTTCAATGTAAAACTATCATTATCTTTGTGTATATGTGTGTACGAATTCGCTTACTTTACGTAACGTCTCTGGTTG  <  1:522118/141‑1                                                                                                                                                   AACGAATAATCACATTTTTTTAGTCATTGATAAGACCTTCAGTTAAATTAATTATCTTTCATGTGCTTCAATGTAAAACGATGATTATCTTTGTGTATATGTGTGTACGAATTCGCTTACTTTACGTAACGCCGCAGGTTG  <  2:25389/141‑1                                                                                                                                                                                                                                                                                                 CACATATACACAAAGATAATCATAGTTTTACATTGAAGCACATGAAAGATAATTCATTTAAAGGAAGGTATTATCAATGACTAAAAAAATGGGATTATTAGTTATGGCTTATGGCACACCTTATAAAGAAAGTGACATAGAGCCATA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  USA300TCH1516_ALE/1965883‑1965737 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agccataACTAATAATCCCATTTTTTTAGTCATTGATAATACCTTCCTTTAAATGAATTATCTTTCATGTGCTTCAATGTAAAACTATGATTATCTTTGTGTATATGTGTGTACGAATTCGCTTACTTTACGTAACGTCTCTGGTTG  >  USA300TCH1516_ALE/1965782‑1965921

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.