New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	2085296 =	15 (0.810)	9 (0.550)	9/240	0.7	37.2%	coding (173/201 nt)	hlb_1	Phospholipase C
	?	USA300TCH1516_ALE	= 2128343	17 (1.030)					coding (1/825 nt)	hlb_2	Phospholipase C

GTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATTTTAGGGACTTTTTAGG  >  USA300TCH1516_ALE/2085163‑2085428                                                                                                                                      |                                                                                                                                     cgtgtattctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTg                                                                                                                                     <  1:32896‑M1/2‑1 (MQ=255)      attctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATcc                                                                                                                                <  1:145983‑M1/7‑1 (MQ=255)      attctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATcc                                                                                                                                >  2:173007‑M1/129‑135 (MQ=255)       ttctcatttgttagaagaaaaatttaatgaagaggacaaaaaaaaaactaaaaatttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCa                                                                                                                               <  1:15939‑M1/8‑1 (MQ=255)         ctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaactcaataaatactgttgttacaaggtttctaTGTATCCaaa                                                                                                                             >  1:134118‑M1/126‑135 (MQ=255)          cTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGCTGATACTGATTTGAAGTTAGTTAGTCATCACGTTTATATGTTATCGACCATTTTCTATCCCAAc                                                                                                                            >  2:47539/1‑135 (MQ=255)          cTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAAc                                                                                                                            >  2:91027/1‑135 (MQ=255)          cTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATATAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAAc                                                                                                                            >  2:132421/1‑135 (MQ=255)           catttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACt                                                                                                                           <  2:146410‑M1/12‑1 (MQ=255)           catttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacacattagggactccaaacccaataaatactgttgttacaaggtttctaAGTAtccaagtt                                                                                                                           <  2:141945‑M1/12‑4 (MQ=255)             tttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagtgtctccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTgg                                                                                                                         >  1:210412‑M1/122‑135 (MQ=255)             tgtgttagaagagaaatgtaatgaagaagagaaaaagacacctatagttttagaaagtatgtaatttagggacccattatggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTgg                                                                                                                         >  2:185382‑M1/122‑135 (MQ=255)                gttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattaggcactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGgaga                                                                                                                      >  2:145592‑M1/119‑135 (MQ=255)                     aagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaatttatactgttgttacaaggtttctaTGTATCCAAACTGGAGACtttt                                                                                                                 >  1:187217‑M1/114‑135 (MQ=255)                     aaaTTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACtttt                                                                                                                 <  1:65331/135‑1 (MQ=255)                          tattattaGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACAt                                                                                                            <  1:7714/135‑1 (MQ=255)                                           gttACTGAAAATAGTTACAAAGCTGAAGTTAAGAAAGATGATACGGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTc                                                                                           <  2:42408/134‑1 (MQ=255)                                            ttACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCa                                                                                          <  2:109457/135‑1 (MQ=255)                                                  acaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGt                                                                                    >  2:145983‑M1/85‑135 (MQ=255)                                                      tAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCGAAAAGTaaaa                                                                                <  1:136850/135‑1 (MQ=255)                                                       taaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAAc                                                                               >  1:177722‑M1/80‑135 (MQ=255)                                                                gCCGAATCTAAGAAAGATGATACTGATTTGAAGTTCGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATACACAATGGAGACTTTTAACATAACATTACTTATCATTCAAAAAATAAAACAACataaaa                                                                      >  2:86366/1‑133 (MQ=255)                                                                   agattgaagatagaTGATAGTGTTTTGTAGTTCGTTTTTTATTTCGTCTATTTCGCATCGTCCGTTTTGTATTCATACTTGTGACTTTTATTATAAAATTACTTATCATTCAAACAGTAAAACAGCCTAATATCa                                                                   <  2:74065/129‑1 (MQ=255)                                                                          agaaagaTGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATAGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTta                                                            >  1:89634/1‑135 (MQ=255)                                                                                         aTTTGAAGTTAGTTAGTCATAACGTTTATATGTTATAGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTttatcattat                                             >  1:174517/1‑135 (MQ=255)                                                                                          tttGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTttatcattatc                                            >  2:152910/1‑135 (MQ=255)                                                                                              gactccaaacccaataaataatgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCaata                                        <  2:102036‑M1/95‑1 (MQ=255)                                                                                                 tccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAATTAAAACAGCATAATATCAAGGTTTATACCTTTATCATTATCAataata                                     <  2:68409‑M1/98‑1 (MQ=255)                                                                                                      tagtCATAACGTTTATATTTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTttatcatta                                              <  1:149132/121‑1 (MQ=255)                                                                                                      tagtCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTttatcatta                                              >  2:149132/1‑121 (MQ=255)                                                                                                             aaatactgttgttacaagtttcctaTGTAGCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATataaaa                         >  1:79003‑M1/26‑135 (MQ=255)                                                                                                                   cgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATaaaataaaat                   >  2:185475‑M1/20‑135 (MQ=255)                                                                                                                                    taTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttagggactttttagg  <  1:124459‑M1/133‑1 (MQ=255)                                                                                                                                      |                                                                                                                                     GTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATTTTAGGGACTTTTTAGG  >  USA300TCH1516_ALE/2085163‑2085428

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.