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.

AGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGT‑TAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATTTTAGGGACTTTTTAGGGAC  >  USA300TCH1516_ALE/2085162‑2085431                                                                                                                                        |                                                                                                                                        aGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGT‑TAGTTATTCATAACGTTTATATTTTATCGACCGTTTTGTAtg                                                                                                                                    >  1:46377/1‑139 (MQ=255)   tttcgacaaTTTCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGT‑TAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCa                                                                                                                                  <  2:326233/135‑1 (MQ=255)      cagcacaTGCATTAGAAAATTTATTAATAGTTGGTGTAATAACTGACAAAAGTTACAAAGACGAATGTAAGATCAAACATAGTCGTTCCTTCC‑ATCTCATTCATAACGTTTTTATGTTCTCGACCGTTTTGTATCCAAAc                                                                                                                               <  2:70408/137‑1 (MQ=255)          acTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGT‑TAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGga                                                                                                                           <  1:15492/140‑1 (MQ=255)                 ttAGCAAAATTAATGTTAGATGGTGCACTTACTAACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGA‑TAGTTAGTCATAACGTTTATATGTTATCGACTGTATTCTATCCACACTGGCGACtgtt                                                                                                                    >  2:8736/1‑140 (MQ=255)                  tAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGT‑TGGGTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTa                                                                                                                   <  1:208257/140‑1 (MQ=255)                        aTTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGT‑TAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATaa                                                                                                             >  2:86336/1‑140 (MQ=255)                                   gTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGT‑TAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATc                                                                                                  <  1:95482/140‑1 (MQ=255)                                      ggggCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGT‑TAGTTAGTCATAACGTTTATATTTTATCGACTGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCAtt                                                                                               <  2:234495/137‑1 (MQ=255)                                         gCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGT‑TAATTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCTAAACTGCAGACATTTAACATAAAATTACTTATCAttaca                                                                                            >  2:111056/1‑137 (MQ=255)                                         gCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTT‑AGGTGGAATTAGTAATACAGTTTAAATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCaa                                                                                            <  2:319176/140‑1 (MQ=255)                                            cTTACTGACAATAGTGCAAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGT‑TAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTAC‑TATCATTCAAAAAg                                                                                        <  2:124470/140‑1 (MQ=255)                                                        gttGCCATAGCTGAATCTAAGCAAGGTGATACGTATTTGAAGT‑TAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTTTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCat                                                                             <  2:228608/138‑1 (MQ=255)                                                        aGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGT‑TAGTTAGTCATAACGTTTATATGTTATCGACCGTTTGGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCat                                                                             <  1:8548/140‑1 (MQ=255)                                                         gTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGT‑TAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCata                                                                            <  1:129108/140‑1 (MQ=255)                                                           gCCAAAGCCGAATCTAAGAAAGATCATACTTATTTGAAGT‑TAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCTAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCataat                                                                          <  1:166511/140‑1 (MQ=255)                                                                aGCCGAATCTAAGAAAGATGATACTGATTTGAAGT‑TAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTGTCATTCAAAAAGTAAAACAGCATAATATCaa                                                                     >  2:87315/1‑140 (MQ=255)                                                                aGCCGAATCTAAGAAAGATGATACTGATTTGAAGT‑TAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCTTTTAAAAAGTAAAACAGCATAATATCaa                                                                     >  1:126398/1‑140 (MQ=255)                                                                aGCCGAATCTAAGAAAGATGATACTGATTTGAAGT‑TAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCaa                                                                     >  2:269215/1‑140 (MQ=255)                                                                       tCTAAGAAAGATGATACTGATTTGAAGT‑TAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTtat                                                              <  1:216955/140‑1 (MQ=255)                                                                        cTAAGAAAGATGATACTGATTTGAAGT‑TAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTtata                                                             <  2:293809/140‑1 (MQ=255)                                                                                      aCTGATTTGAAGT‑TAGTTAGTCATAACGTTTATATGTTATCGAACGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTttatcattatc                                               <  1:2432/140‑1 (MQ=255)                                                                                       atGATTTGAAGT‑TTGTTAGTCATAACGTTTATATGAAAGGTACCGTTTTGTTTTCAAGCTTGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTttatcattatca                                              <  2:53603/139‑1 (MQ=255)                                                                                           tttGAAGT‑TAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAACAGTAAAACAGCATAATATCAATGTTTATAACTTTATCATTATCaataa                                          >  1:171488/1‑140 (MQ=255)                                                                                                            cATTACGTTTATATGTTATCGACCGTGTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCGTataaaata                          <  2:60079/140‑1 (MQ=255)                                                                                                                   tttATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAtttt                   <  1:99062/140‑1 (MQ=255)                                                                                                                     tataTGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttag                 <  1:210285/140‑1 (MQ=255)                                                                                                                                   cGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAACATTTATAACTTTAACCTTATCAATAATACCTCATATAAAATACATttttaagaactttttaggga   >  2:222696/1‑140 (MQ=255)                                                                                                                                    gTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCCATAAGTAAAACCGCCTAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAtttcaagtactttttaacgcc  >  1:218823/1‑135 (MQ=255)                                                                                                                                        |                                                                                                                                        AGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGT‑TAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATTTTAGGGACTTTTTAGGGAC  >  USA300TCH1516_ALE/2085162‑2085431

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.