New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	2085296 =	26 (1.030)	6 (0.270)	6/232	NT	23.0%	coding (173/201 nt)	hlb_1	Phospholipase C
	?	USA300TCH1516_ALE	= 2128343	17 (0.750)					coding (1/825 nt)	hlb_2	Phospholipase C

CATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAA  >  USA300TCH1516_ALE/2128207‑2128458                                                                                                                                         |                                                                                                                    cATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgta                                                                                                                 >  1:89454‑M2/1‑137 (MQ=255)      gtgtATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaa                                                                                                            <  2:96457‑M2/141‑10 (MQ=255)      gtgtATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaa                                                                                                            <  2:157337‑M2/141‑10 (MQ=255)         tATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaact                                                                                                         <  1:268685‑M2/141‑13 (MQ=255)            tctcATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTggg                                                                                                      >  1:272709/1‑141 (MQ=255)                    gTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAatat                                                                                              <  1:45699/141‑1 (MQ=255)                       agaagaAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAAACTTATAGTTATCCAAACTGGGGGCAATATaaa                                                                                           >  2:151064/1‑141 (MQ=255)                                 tttAATGAGGAGGACAAAAAAACAACTAAAATTTTAGAAAGTAGGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGAtt                                                                                 <  1:125407/141‑1 (MQ=255)                                          gagTACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGAc                                                                        <  2:115791/141‑1 (MQ=255)                                          gagGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGAc                                                                        <  2:82646/141‑1 (MQ=255)                                          gagGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGAc                                                                        >  2:72090/1‑141 (MQ=255)                                               cAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAAtct                                                                   >  1:250977/1‑141 (MQ=255)                                                     aacaacTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTtaa                                                             >  2:264823/1‑140 (MQ=255)                                                       taacTAAAATTTTAGAGAGGACTGAACTTTGGGTGCCATGAGTGACTCCAAACTCATTGAGCAGTGACTTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTatat                                                           <  2:299311/140‑1 (MQ=255)                                                         aataaAATTATAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagca                                                         >  2:43210‑M2/3‑81 (MQ=255)                                                                         gTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGt                                         <  1:65723/141‑1 (MQ=255)                                                                               aaTTTAGAGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataacttt                                   <  2:215538‑M2/141‑83 (MQ=255)                                                                                 tttAGGGACCAATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCa                                 <  2:295118/141‑1 (MQ=255)                                                                                     gggACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGa                             >  2:178087/1‑141 (MQ=255)                                                                                       gccccatTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattat                           >  1:139569‑M2/3‑51 (MQ=255)                                                                                         cccATTAGGGACTTCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattatca                         >  1:58702‑M2/1‑49 (MQ=255)                                                                                          ccATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCAt                        <  1:151064/141‑1 (MQ=255)                                                                                           cATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattatcaat                       >  2:123716‑M2/1‑47 (MQ=255)                                                                                                  gACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattatcaataatacct                >  2:67416‑M2/1‑40 (MQ=255)                                                                                                     tCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGt             >  1:140495/1‑141 (MQ=255)                                                                                                     tCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGt             >  1:150968/1‑141 (MQ=255)                                                                                                             aaataaataATGTTGTTACAAGGTTTCTATGTATTCAAACTGGCGGCAATATAAACGCGGTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGa     <  1:3408/140‑1 (MQ=255)                                                                                                                taaataCTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTATTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGATGCATCAGACaa  >  1:217546/1‑141 (MQ=255)                                                                                                                                         |                                                                                                                    CATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAA  >  USA300TCH1516_ALE/2128207‑2128458

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