New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	1803422 =	0 (0.000)	7 (0.070)	4/154	NT	100%	coding (446/807 nt)	USA300HOU_RS12765	hypothetical protein
	?	USA300TCH1516_ALE	1803448 =	NA (NA)					coding (420/807 nt)	USA300HOU_RS12765	hypothetical protein

TAATTTGGACTTTAAAGTATTCCCAATTATAATTCTTCATGATTTTCTTATTGGATTTCGAATTTGGTTTCATGCATTGATGTCTCAAAGAACATGATGAACAGTCATCACATTCATATAGTTTGAAGTCTCGTTTAAAACCATATCTATCATTACGGTATGCATATCTTTTAAAACCTA  >  USA300TCH1516_ALE/1803333‑1803512                                                                                          |                                                                                           ttatatgtcctaataacaaaagaataggttttaaaagatatgcataccgtaatgatagatatggttttaaacgagacttcaaactatatGa                                                                                           >  1:3783858‑M1/90‑91 (MQ=255) ttatatgtcctaataacaaaagaataggttttaaaagatatgcataccgtaatgatagatatggttttaaacgagacttcaaactatatGa                                                                                           >  1:3812543‑M1/90‑91 (MQ=255)  tatatgtcctaataacaaaagaataggttttaaaagatatgcataccgtaatgatagatatggttttaaacgagacttcaaactatatGaa                                                                                          <  1:528906‑M1/3‑1 (MQ=255)  tatatgtcctaataacaaaagaataggttttaaaagatatgcataccgtaatgatagatatggttttaaacgagacttcaaactatatGaa                                                                                          <  1:529523‑M1/3‑1 (MQ=255)                         taggttttaaaagatatgcataccgtaatgatagatatggttttaaacgagacttcaaactatatGAACATGATGAACAGTCATCACATTc                                                                   >  1:933747‑M1/66‑91 (MQ=255)                          aggttttaaaagatatgcataccgtaatgatagatatagttttaaacgtgacttcaagctatatGAACATGATGAACAGTCATCACATTCa                                                                  >  1:2613218‑M1/65‑91 (MQ=255)                          aggttttaaaagatatgcataccgtaatgatagatatagttttaaacgtgacttcaagctatatGAACATGATGAACAGTCATCACATTCa                                                                  >  1:664354‑M1/65‑91 (MQ=255)                           ggttttaaaagatatgcataccgtaatgatagatatggttttaaacgagacttcaaactatatGAACATGATGAACAGTCATCACATTCat                                                                 >  1:569661‑M1/64‑91 (MQ=255)                            gttttaaaagatatgcataccgtaatgatagatatggttttaaacgagacttcaaactatatGAACATGATGAACAGTCATCACATTCata                                                                >  1:3994972‑M1/63‑91 (MQ=255)                            gttttaaaagatatgcataccgtaatgatagatatggttttaaacgagacttcaaactatatGAACATGATGAACAGTCAACACATTCata                                                                >  1:2941824‑M1/63‑91 (MQ=255)                              tttaaaagatatgcataccgtaatgatagatatggttttaaacgagacttcaaactatatGAACATGATGAACAGTCATCACATTCAtata                                                              >  1:2847337‑M1/61‑91 (MQ=255)                                                                                          GAACATGATGAACAGTCATCACATTCATATAGTTTGAAGTCTCGTTTAAAACCATATCTATCATTACGGTATGCATATCTTTTAAAACCTa  >  1:753624‑M1/1‑91 (MQ=255)                                                                                          |                                                                                           TAATTTGGACTTTAAAGTATTCCCAATTATAATTCTTCATGATTTTCTTATTGGATTTCGAATTTGGTTTCATGCATTGATGTCTCAAAGAACATGATGAACAGTCATCACATTCATATAGTTTGAAGTCTCGTTTAAAACCATATCTATCATTACGGTATGCATATCTTTTAAAACCTA  >  USA300TCH1516_ALE/1803333‑1803512

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