New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	2085296 =	20 (0.890)	5 (0.250)	5/250	2.0	24.4%	coding (173/201 nt)	hlb_1	Phospholipase C
	?	USA300TCH1516_ALE	= 2128343	13 (0.640)					coding (1/825 nt)	hlb_2	Phospholipase C

CTACATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAATGTGAAAAAAG  >  USA300TCH1516_ALE/2128204‑2128481                                                                                                                                            |                                                                                                                                           cTACATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTa                                                                                                                                            >  2:177651/1‑140 (MQ=255)   aCATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAATGTTTCTAtg                                                                                                                                          <  1:35516‑M2/140‑3 (MQ=255)                          agaagaAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCATACTGGGGGCAatatcc                                                                                                                   >  1:213531/1‑138 (MQ=255)                            aagaaAAATTTAATGAAGAGGAAAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacat                                                                                                                 >  1:71099‑M2/1‑113 (MQ=255)                             agaaAAATTTAATGAAGAGGACAAAAAAACAACTAAAATAATAGAAAGTATGTAATTTAGGAACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGAAGGCAATATAAAcg                                                                                                                >  2:148559/1‑140 (MQ=255)                              gaaAAATTTAATGAAGAGTACAATAAGACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataa                                                                                                               <  2:161159‑M2/140‑30 (MQ=255)                                aaaaTTTAATGAAGAGAACAAAAAACCAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACGGGGGGCAATATAAACgcgc                                                                                                             <  1:110712/140‑1 (MQ=255)                                     tgaaagAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGAtt                                                                                                        <  2:121333/135‑1 (MQ=255)                                               ggACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaa                                                                                              >  1:2499‑M2/1‑94 (MQ=255)                                               ggACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCACTATAAACGCGCTGATTTAATCGGACa                                                                                              >  2:47125/1‑140 (MQ=255)                                                   aaaaaaaCAACTAAAATTTTAGAAAGTATGTAATATAGGGCCCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaag                                                                                          >  2:197251‑M2/1‑90 (MQ=255)                                                                             tATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGt                                                                <  2:143070/140‑1 (MQ=255)                                                                                 tAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCATATCTTAAAAATAATGATGTCGTAata                                                            <  1:63405/140‑1 (MQ=255)                                                                                                      aCTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACAGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAAt                                       <  1:104612/140‑1 (MQ=255)                                                                                                             aCCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCAt                                <  1:76356/140‑1 (MQ=255)                                                                                                             aCCCAATAAATACTGTTGTTACAAAGTTTCTATGTATCCAAACTGAGGGAAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCAt                                <  1:4864/140‑1 (MQ=255)                                                                                                                        aCTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAatta                     >  1:11449/1‑140 (MQ=255)                                                                                                                                        ttCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattatcaataatacctcatataaaataaaattttagggactttttaggga     <  1:196593‑M2/140‑136 (MQ=255)                                                                                                                                          cTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAATGTGaaaaaa   <  1:208600/140‑1 (MQ=255)                                                                                                                                           tATGGATCCAAAGTGGGGGGAATATAAACGAGCTTATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAATGTGAAAAAAg  <  2:37203/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                           CTACATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAATGTGAAAAAAG  >  USA300TCH1516_ALE/2128204‑2128481

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