New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	= 800928	158 (1.210)	5 (0.040)	4/280	9.3	3.4%	intergenic (+62/‑289)	hisC_1/USA300HOU_RS03915	Histidinol‑phosphate aminotransferase/Putative 5'(3')‑deoxyribonucleotidase
	?	USA300TCH1516_ALE	801948 =	134 (1.080)					intergenic (+189/+435)	USA300HOU_RS03915/USA300HOU_RS03920	Putative 5'(3')‑deoxyribonucleotidase/Putative lipid kinase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

ACCATTTCCAACTGGTGTTAGAATTACAATTGGTTTTAAAGAACAAAATGATAAAATGTTAGAAGTTTTATCAAACTTTAAATACGAATAGTAAGTGGGGAGTGGGACAGAAATGATACTTTCGCAAAATTTATTTCGTCGTCCCACCCCAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  USA300TCH1516_ALE/800776‑800928 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccccaacATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCATTAGTCGTTCTTTATCATTT  >  USA300TCH1516_ALE/801948‑802089                                                                                                                                                                                                                                                                                                         ACCATTTCCAACTGGTGTTAGAATTACAATTGGTTTTAAAGAACAAAATGATAAAATGTTAGAAGTGTTATCAAACTTTAAATACGAATAGTAAGTGGGGAGTGGGACAGAAATGATACTTTCGCAAAATTTATTTCGTCGTCCCACC                                                                                                                                                     >  1:1238049/1‑148           ACTGGTGTTAGAATTACAATTGGTTTTAAAGAACAAAATGATAAAATGTTAGAAGTTTTATCAAACTTTAAATACGAATAGTAAGTGGGGAGTGGGACAGAAATGATACTTTCGCAAAATTTATTTCGTCGTCCCACCCCAACATAAA                                                                                                                                           <  1:160891/148‑1                                                  GATAAAATGTTAGAAGTTTTATCAAACTTTAAATACGAATAGTAAGTGGGGAGTGGGACAGAAATGATACTTTCGCAAAATTTATTTCGTCGTCCCACCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCAT                                                                                                    >  1:1309309/1‑148                                                                                                 gggggtTGGGACACAAATGATACTTTCGCAAAATTTATTTCGTCGTCCCACCCCAACATAAAGAAATACCTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGAC                                                     <  1:586360/142‑1                                                                                                                AATGATACTTTCGCAAAATTTATTTCGTCGTCCCACCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATC                                      <  2:1309309/148‑1                                                                                                                                              gCCCACCCCAACATAAAGAAATACTTTTTCTTCAGAAATTAGTATTTCTTATGCATGAGTGTAAATCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCATTAGTCGTTCTTTA        >  1:1029219/2‑148                                                                                                                                                    cctAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCGTTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCATTAGTCGTTCTTTATCATTT  <  2:1049278/145‑1                                                                                                                                                    aCCAACAGTAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAGTCGCTACATAATCAATCATTAGTCGTTCTTTATCATTT  <  2:900955/147‑1                                                                                                                                                    CCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCATTAGTCGTTCTTTATCATTT  <  2:543688/148‑1                                                                                                                                                    CCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCATTAGTCGTTCTTTATCATTT  <  2:801045/148‑1                                                                                                                                                                                                                                                                                                         ACCATTTCCAACTGGTGTTAGAATTACAATTGGTTTTAAAGAACAAAATGATAAAATGTTAGAAGTTTTATCAAACTTTAAATACGAATAGTAAGTGGGGAGTGGGACAGAAATGATACTTTCGCAAAATTTATTTCGTCGTCCCACCCCAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  USA300TCH1516_ALE/800776‑800928 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccccaacATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCATTAGTCGTTCTTTATCATTT  >  USA300TCH1516_ALE/801948‑802089

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