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

ATTATAATTCTTCATGATTTTCTTATTGGATTTCGAATTTGGTTTCATGCATTGATGTCTCAAAGAACATGATGAACAGTCATCACATTCATATAGTTTGAAGTCTCGTTTAAAACCATATCTATCATTACGGTATGCATATCTTTTAAAACCTATT  >  USA300TCH1516_ALE/1803358‑1803514                                                                                           |                                                                   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)                                                                  AACATGATGAACAGTCATCACATTCATATAGTTTGAAGTCTCGTTTAAAACCATATCTATCATTACGGTATGCATATCTTTTAAAACCTAt   >  1:182354‑M1/1‑91 (MQ=255)                                                                  AACATGATGAACAGTCATCACATTCATATAGTTTGAAGTCTCGTTTAAAACCATATCTATCATTACGGTATGCATATCTTTTAAAACCTAt   <  1:3143736‑M1/91‑1 (MQ=255)                                                                   ACATGATGAACAGTCATCACATTCATATAGTTTGAAGTCTCGTTTAAAACCATATCTATCATTACGGTATGCATATCTTTTAAAACCTAtt  <  1:2438413‑M1/91‑1 (MQ=255)                                                                   ACATGATGAACAGTCATCACATTCATATAGTTTGAAGTCTCGTTTAAAACCATATCTATCATTACGGTATGCATATCTTTTAAAACCTAtt  >  1:2759444‑M1/1‑91 (MQ=255)                                                                   ACATGATGAACAGTCATCACATTCATATAGTTTGAAGTCTCGTTTAAAACCATATCTATCATTACGGTATGCATATCTTTTAAAACCTAtt  <  1:297515‑M1/91‑1 (MQ=255)                                                                                           |                                                                   ATTATAATTCTTCATGATTTTCTTATTGGATTTCGAATTTGGTTTCATGCATTGATGTCTCAAAGAACATGATGAACAGTCATCACATTCATATAGTTTGAAGTCTCGTTTAAAACCATATCTATCATTACGGTATGCATATCTTTTAAAACCTATT  >  USA300TCH1516_ALE/1803358‑1803514

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.