New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	2085296 =	24 (0.980)	13 (0.590)	12/252	0.7	40.9%	coding (173/201 nt)	hlb_1	Phospholipase C
	?	USA300TCH1516_ALE	= 2128343	16 (0.730)					coding (1/825 nt)	hlb_2	Phospholipase C

TTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATTTTAGGGACTTTTTAGGGACTTTA  >  USA300TCH1516_ALE/2085164‑2085435                                                                                                                                     |                                                                                                                                            ttGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCaa                                                                                                                                     >  1:250522/1‑141 (MQ=255)   gCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTCAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAAc                                                                                                                                   <  1:100898/141‑1 (MQ=255)    cAACACTTGCATTAGCAAATTTATTATTAGTTGGGGCACTGACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACt                                                                                                                                  <  1:106661/141‑1 (MQ=255)       cacTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCTTTTTGTATCCAAACTGga                                                                                                                               >  2:29458/1‑141 (MQ=255)                  agaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACa                                                                                                                    <  1:210890‑M1/26‑1 (MQ=255)                            tttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacacattagggattccaaaaccaataaatactgttgttaaaaagtttctaTGTATCCAAAATTAAAACTTTTAACATAAAATTACt                                                                                                          <  2:42009‑M1/36‑1 (MQ=255)                                       gCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTGTTGTAAACAAACTGGAGACTTTTAACATAAAATTACTTATCATTCaaa                                                                                               >  2:244991/1‑141 (MQ=255)                                               tGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTaaaa                                                                                       <  1:34736/141‑1 (MQ=255)                                                aacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAAc                                                                                      <  2:223414‑M1/56‑1 (MQ=255)                                                    actaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAACCTGGAGACTTTTAACATAAAATTACTAATCATTCAAAAAGTAAAACAGCa                                                                                  <  1:126825‑M1/60‑1 (MQ=255)                                                       gTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCataa                                                                               <  1:174990/141‑1 (MQ=255)                                                            aaaGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCa                                                                          >  1:242051/1‑141 (MQ=255)                                                                ccGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGt                                                                      >  2:10925/1‑141 (MQ=255)                                                                    ttatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTtat                                                                  >  2:3022‑M1/66‑141 (MQ=255)                                                                    aTCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTtat                                                                  >  1:238018/1‑141 (MQ=255)                                                                        gtaatttagggacacattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCAAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGAATAATATCAAGGTTTATAACt                                                              <  2:45858‑M1/80‑1 (MQ=255)                                                                                   tACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTATTTATCATTCAAAAAGTAATACAGCATAATATCAAGGTTTATAACTttatcattatc                                                   >  2:2619/1‑141 (MQ=255)                                                                                     ccattagggactccaaacccaataagtactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACCCCCCTGTTTCAAGGTTTATAACTTTATCATTATCaa                                                 >  2:129327‑M1/49‑141 (MQ=255)                                                                                     ccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCaa                                                 >  1:96356‑M1/49‑141 (MQ=255)                                                                                         tttGAAGTTAGTTAGTCATAACGTTTATCTGTTATCGACCGTTTTGTATCCACACTGGAAACTTTTAACATCAAATTAATTATCTTTCAAAAAGTAAAACAGCATAATATAAATGTTTATAACTTTGTTATCATCagtaat                                             >  2:124683/1‑141 (MQ=255)                                                                                               ctccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCa                                       >  1:272293‑M1/39‑141 (MQ=255)                                                                                                 tagttagtCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCata                                     >  1:245352/1‑141 (MQ=255)                                                                                                  agttagtCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATCCCTCctat                                    >  2:154494/1‑141 (MQ=255)                                                                                                  agttagtCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCatat                                    >  2:44549/1‑141 (MQ=255)                                                                                                  agttagtCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATAACTCatat                                    >  2:264391/1‑141 (MQ=255)                                                                                                    aacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATataa                                  <  2:195815‑M1/108‑1 (MQ=255)                                                                                                     tagtCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATataaa                                 >  2:222626/1‑141 (MQ=255)                                                                                                         cATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATataaaataa                             <  1:248536/141‑1 (MQ=255)                                                                                                         aataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATataaaataa                             <  2:16957‑M1/113‑1 (MQ=255)                                                                                                            aaCCTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATaaaataaaat                          <  1:256629/141‑1 (MQ=255)                                                                                                                actgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAtttta                      <  2:119387‑M1/120‑1 (MQ=255)                                                                                                                   atatGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttaggg                   >  1:128923/1‑141 (MQ=255)                                                                                                                            cGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttagggactttttag          <  1:122052/141‑1 (MQ=255)                                                                                                                              aCCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttagggactttttaggg        >  2:66203/1‑141 (MQ=255)                                                                                                                               ccGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTCTCAATAATACCTCATATAAAATAAAAttttagcgactttttaggga       >  1:165423/1‑141 (MQ=255)                                                                                                                               ccGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttagggactttttaggga       >  1:125054/1‑141 (MQ=255)                                                                                                                                    aTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATTTTAGGGACTTTTTAGGGACttta  >  1:223414‑M1/2‑141 (MQ=255)                                                                                                                                     |                                                                                                                                            TTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATTTTAGGGACTTTTTAGGGACTTTA  >  USA300TCH1516_ALE/2085164‑2085435

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.