New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	= 800928	25 (0.820)	4 (0.140)	4/266	3.2	16.4%	intergenic (+62/‑289)	hisC_1/USA300HOU_RS03915	Histidinol‑phosphate aminotransferase/Putative 5'(3')‑deoxyribonucleotidase
	?	USA300TCH1516_ALE	801948 =	17 (0.590)					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.

TTGTGACTGAGATGAACTTTTATGTCTTAGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCATTAGTCGTTCTTTATC  >  USA300TCH1516_ALE/801816‑802085                                                                                                                                     |                                                                                                                                          ttGTGACTGAGATGAACTTTTATGTCTTAGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGaaa                                                                                                                                   <  2:103755/141‑1 (MQ=255)   gtgACTGAGATGAACTTTTAGGTCTTAGATACTACATCACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATa                                                                                                                                 <  2:142524/141‑1 (MQ=255)   gcgACTGCGATGAACTTTTATGTCTTAGACACCACAACACTATATTGGCAGTAGGTGACTGCGTGGCCCCAACATAGAGAAATTTGAACTCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATa                                                                                                                                 <  2:324972/139‑1 (MQ=25)    gtaCTGATATGAACTTTTAGTTCTTAGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATAc                                                                                                                                <  2:237533/139‑1 (MQ=255)        tGAGATGAACTTTTATGTCTTAGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCGACAATATAGGAATTGCAACGCCAATCTTTACAGACAATGCAAGTTGGGGTGGGCCCCAACCTAAAGAAAGACtttt                                                                                                                            >  1:152213/1‑141 (MQ=25)           gaTGAACTTTTATGTCTTAGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCt                                                                                                                         <  2:245179/141‑1 (MQ=255)                     tATGTCTTAGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTa                                                                                                               >  1:189343/1‑141 (MQ=255)                     tATGTCTTAGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGCGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGCACTTa                                                                                                               >  1:250411/1‑141 (MQ=255)                           ttAGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTAttt                                                                                                         >  2:69033/1‑141 (MQ=255)                                  cTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTCTTTCTTATgt                                                                                                  >  2:300655/1‑140 (MQ=255)                                   tACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAACTACTATTTCTTTAGAAACTAGTCTCTCTTATGCa                                                                                                 >  2:334942/1‑141 (MQ=18)                                            atatTGGCAGTAGTTGACTGCAGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTATAAATTAGTATTTCTTATGCATGAGTGTaa                                                                                        >  2:6141/1‑141 (MQ=255)                                                               gCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATAttttt                                                                     >  1:3903/1‑141 (MQ=255)                                                                 ggggCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCGTGAGTGTAACTCCTGCATTCATCTTTTTat                                                                   >  1:106426/1‑140 (MQ=255)                                                                                 gtaaTTGGAAAGCCAATTTCTACAGACACTGCAAGTTGGGGTGGGCCCCAACATACAGAACTACTGTTTCTTTAGAACTTCGTATTTCTTATGCATGCGTGTAACTCCTGCGTTGATATCTTGAAGTACACTTTAGCtgtg                                                   >  1:56305/3‑141 (MQ=255)                                                                                                   tatACAGGCAATGTAAGGTGGGGTGGGGCGCAAGATAAAGAAATACTTTTTCTTTATAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAACCGCACGGGAGGTGTGGCTGATGATAAAGAATCg                                 <  2:328532/139‑1 (MQ=14)                                                                                                          aCAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATa                          <  1:160138/141‑1 (MQ=255)                                                                                                                                  aaCATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGGAGCGGCACGTAATCAGTCATTAGTCGTTCTTTATc  <  1:69033/141‑1 (MQ=18)                                                                                                                                     |                                                                                                                                          TTGTGACTGAGATGAACTTTTATGTCTTAGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCATTAGTCGTTCTTTATC  >  USA300TCH1516_ALE/801816‑802085

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