New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	2085296 =	29 (1.040)	4 (0.160)	4/250	3.0	13.8%	coding (173/201 nt)	hlb_1	Phospholipase C
	?	USA300TCH1516_ALE	= 2128343	24 (0.960)					coding (1/825 nt)	hlb_2	Phospholipase C
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CTACATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACA‑ATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTA  >  USA300TCH1516_ALE/2128204‑2128469                                                                                                                                            |                                                                                                                                cTACATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAAGCTTCTa                                                                                                                                 >  2:59003/1‑140 (MQ=255)            tataCTCATTTGTTAGAAGAAAAATTGAAGGAAAACGACAAGTTGGGTACCAAAAAATTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAAc                                                                                                                      <  1:74668/140‑1 (MQ=255)            tATTCTCATTTGTTAGAATGATAATTTAATGAAGAGGACAAAAAAACAAATAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAAc                                                                                                                      <  1:43778/140‑1 (MQ=255)             aTTCTCATTTGCTAGAAGAAAAATATAATGCTGAGGACAAAAAAACAACTAAAATCTTAGAAAGTATGTAATTTAGGGTCCCATTAAGGACGCCAAACCCAGTAAATCCTGTGGTTACCAGGGTCCTATGTATCCCAACt                                                                                                                     >  2:183214/1‑140 (MQ=255)                 tcATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAACCATGTAGTAACTACTTTTCTATGTATCCAAACTggag                                                                                                                 <  1:312424/140‑3 (MQ=255)                          agaagaAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATGTTAGAAAGTATGTAATTTAGGGCCCCATTAGGGACTCAAAACCCAATAAATACTGTTGTTACAAGGTATCTATGTATCCAAACTGGGGGCAATATaa                                                                                                        >  1:118218/1‑140 (MQ=255)                            aagaaAAATTTAATGAAGAGGACAAAACAACAACTAAAATTTTAGAAAGTATGAAATTTAGAGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAAc                                                                                                      >  1:284946/1‑140 (MQ=255)                                       aaTGAAGAGGACAAAAAAACAAATAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTa                                                                                           <  1:36340/140‑1 (MQ=255)                                            agagGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCgg                                                                                      <  1:29306/140‑1 (MQ=255)                                            agagGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCgg                                                                                      <  1:29929/140‑1 (MQ=255)                                             gagGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTATGAGCAATATAAAAGCGCTGATTTAATCGGa                                                                                     >  1:3739/1‑140 (MQ=255)                                                        aacaacTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCACACTGAGGGCAATATAAACGCGCTGATTTAATCGGACA‑ATCTTCTta                                                                         >  2:103320/1‑140 (MQ=255)                                                        aacaacTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACA‑ATCTTCTta                                                                         >  1:12897/1‑140 (MQ=255)                                                          caacTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTAAAAGGTATATATGTATACAAAATAGGGGCAATATAAACGCGCTGATTTAATCGGACA‑ATATTCTtata                                                                       >  1:20645/1‑140 (MQ=255)                                                            acTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACA‑ATCTTCTTATAtt                                                                     >  1:193010/1‑140 (MQ=255)                                                            acTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACACAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGAAA‑ATCTTCTTATAtt                                                                     >  2:308670/1‑140 (MQ=255)                                                                  aTTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATCTAAACGCGCTGATTTAATCGGACA‑ATCTTCTTATATTAAAaat                                                               <  1:172168/140‑1 (MQ=255)                                                                                    tttAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACA‑ATCTTCTTATATTAAAAATAATGATGTCGTAATATTc                                             >  2:215244/1‑140 (MQ=255)                                                                                      tAGGGACCCATGAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGTAATATAAACGCGCTGATTTAATCGGACA‑ATCTTCTTATATTAAAAATAATGATGTCGTAATATTCaa                                           <  1:51402/140‑1 (MQ=255)                                                                                          gACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACA‑ATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGaa                                       <  1:206993/140‑1 (MQ=255)                                                                                                  aGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACA‑ATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGa                               >  2:71851/1‑140 (MQ=255)                                                                                                       cTCCAACCCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATAGGACA‑ATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATg                          >  1:313857/1‑140 (MQ=255)                                                                                                        tCCAAATCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACA‑ATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATgg                         <  2:13238/140‑1 (MQ=255)                                                                                                            aaCCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACA‑ATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCAGTAGATAATGGTACa                     >  1:191383/1‑140 (MQ=255)                                                                                                             aCCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACA‑ATCTTCTTATATTAAAAATAATGATGTCTTAATATTCAATGAAGCATTTGATAATGGTGCAt                    >  1:190777/1‑140 (MQ=255)                                                                                                                 aataaataCTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACA‑ATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGa                >  1:258892/1‑140 (MQ=255)                                                                                                                   taaataCTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACA‑ATCTTCTTATATTAAAACTAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCacaaa              >  2:166664/1‑136 (MQ=255)                                                                                                                   taaataCTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACA‑ATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACa              >  1:131214/1‑140 (MQ=255)                                                                                                                               ttACAAGGTTTCTATGTATCCAAACTGGGGGCCGACTACACGCGGCGTGGTAATAGGAGATAACGTAATA‑ATGTAAAATAATGATGTAGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGAAAAATTATTAAGTa  >  1:327093/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                                CTACATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACA‑ATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTA  >  USA300TCH1516_ALE/2128204‑2128469

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.