New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	= 1912410	54 (0.690)	4 (0.050)	4/268	8.6	7.2%	coding (532/609 nt)	USA300HOU_RS15290	hypothetical protein
	?	USA300TCH1516_ALE	= 2319335	NA (NA)					coding (253/345 nt)	USA300HOU_RS15290	hypothetical protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

ATTGATGTAGAGCCTGTTTTTGGATTTATGAAGGCTATTTTGGGTTTAACTCGAACGTCAGTTCGAGGAATAAATAAAGTTAAACGAGAGCTAGGTTTTGTATTAATGGCACTTAATATAAGGAAAATAGCAGCTCAACGAGCTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  USA300TCH1516_ALE/1912266‑1912410 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gagctgCTATTTTCCTTATATTAAGTGCCATTAATACAAAACCTAGCTCTCGTTTAACTTTATTTATTCCTCGAACTGACATTCGAGTGAAACCCAAAATAGCCTTCATAAATCCAAAAACAGGCTCTACATCAATTTTTCT  <  USA300TCH1516_ALE/2319335‑2319200                                                                                                                                                                                                                                                                                           ATTGATGTAGTGCCTGTGTTTGGATTTATGAAGGCTATTTTGGGTTTAACTCGAACGTCAGTTCGAGGAATAAATAAAGTTAAACGAGAGCTAGGTTTTGTATTAATGGCACTTAATATAAGGAAAATAGCAGCTCAACGA                                                                                                                                              <  2:291661/141‑1     ATGTAGAGCCTGTTTTTGGATTTATGAAGGCTATTTTGGGTTTAACTCGAACGTCAGTTCGAGGAATAAATAAAGTTAAACGAGAGCTAGGTTTTGTATTAATGGCACTTAATATAAGGAAAATAGCAGCTCAACGAGCTG                                                                                                                                          >  1:291192/1‑141     ATGTAGAGCCTGTTTTTGGATTTATGAAGGCTATTTTGGGTTTAACTCGAACGTCAGTTCGAGGAATAAATAAAGTTAAACGAGAGCTAGGTTTTGTATTAATGGCACTTAATATAAGGAAAATAGCAGCTCAACGAGCTG                                                                                                                                          >  1:344253/1‑141                                                                                                    GTATTAATGGCACTTAATATAAGGAAAATAGCAGCTCAACGAGCTGCTATTTTCCTTATATTAAGTGCCATTAATACAAAACCTAGCTCTCGTTTAACTTTATTTATTCCTCGAACTGACATTCGAGTGAAACCCAAAATA                                           >  1:51799/1‑141                                                                                                             GCACTTAATATAAGGAAAATAGCAGCTCAACGAGCTGCTATTTTCCTTATATTAAGTGCCATTAATACAAAACCTAGCTCTCGTTTAACTTTATTTATTCCTCGAACTGACATTCGAGTGAAACCCAAAATAGCCTTCATA                                  >  2:802608/1‑141                                                                                                               ACTTAATATAAGGAAAATAGCAGCTCAACGAGCTGCTATTTTCCTTATATTAAGTGCCATTAATACAAAACCTAGCTCTCGTTTAACTTTATTTATTCCTCGAACTGACATTCGAGTGAAACCCAAAATAGCCTTCATAAA                                <  2:51799/141‑1                                                                                                                                      CTCAACGAGCTGCTATTTTCCTTATATTAAGTGCCATTAATACAAAACCTAGCTCTCCTTTAACTTTATTTATTCCTCGAACCGCCCGCGGAGTGTAACCCAACATAGCCTTCAGAAATCCAAAAACAGGCGCTACATCAA         >  1:26376/1‑141                                                                                                                                            GAGCTGCTATTTTCCTTATATTAAGTGCCATTAATACAAAACCTAGCTCTCGTTTAACTTTATTTATTCCTCGAACTGACATTCGAGTGAAACCCAAAATAGCCTTCATAAATCCAAAAACAGGCTCTACATCAATTTTTC   <  1:372367/141‑1                                                                                                                                            GAGCTGCTATTTTCCTTATATTAAGTGCCATTAATACAAAACCTAGCTCTCGTTTAACTTTATTTATTCCTCGAACTGACATTCGAGTGAAACCCAAAATAGCCTTCATAAATCCAAAAACAGGCTCTACATCAATTTTTC   >  1:583265/1‑141                                                                                                                                            GAGCTGCTATTTTCCTTATATTAAGTGCCATTAATACAAAACCTAGCTCTCGTTTAACTTTATTTATTCCTCGAACTGACATTCGAGTGAAACCCAAAATAGCCTTCATAAATCCAAAAAAAGGCTCTACATCAATTTTTt   >  1:490490/1‑140                                                                                                                                            GAGCTGCTATTTTCCTTATATTAAGTGCCATTAATACAAAACCTAGCTCTCGTTTAACTTTATTTATTCCTCGAACTGACAGTCGAGTGAAACCCAAAATAGCCTTCATAAATCCAAAAACAGGCTCTACATCAATTTTTC   >  1:438689/1‑141                                                                                                                                             AGCTGCTATTTTCCTTATATTAAGTGCCATTAATACAAAACCTAGCTCTCGTTTAACTTTATTTATTCCTCGAACTGACATTCGAGTGAAACCCAAAATAGCCTTCATAAATCCAAAAACAGGCTCTACATCAATTTTTCT  >  1:577845/1‑141                                                                                                                                                                                                                                                                                           ATTGATGTAGAGCCTGTTTTTGGATTTATGAAGGCTATTTTGGGTTTAACTCGAACGTCAGTTCGAGGAATAAATAAAGTTAAACGAGAGCTAGGTTTTGTATTAATGGCACTTAATATAAGGAAAATAGCAGCTCAACGAGCTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  USA300TCH1516_ALE/1912266‑1912410 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gagctgCTATTTTCCTTATATTAAGTGCCATTAATACAAAACCTAGCTCTCGTTTAACTTTATTTATTCCTCGAACTGACATTCGAGTGAAACCCAAAATAGCCTTCATAAATCCAAAAACAGGCTCTACATCAATTTTTCT  <  USA300TCH1516_ALE/2319335‑2319200

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