New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	2085296 =	18 (0.640)	6 (0.240)	6/232	NT	26.6%	coding (173/201 nt)	hlb_1	Phospholipase C
	?	USA300TCH1516_ALE	= 2128343	17 (0.680)					coding (1/825 nt)	hlb_2	Phospholipase C

AAAAAAAGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATTTTAGGGACTTTTTAGGGACTT  >  USA300TCH1516_ALE/2085156‑2085433                                                                                                                                             |                                                                                                                                          ctacatccgtgtattctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctat                                                                                                                                           <  2:174798‑M1/1‑1 (MQ=255)  aaaaaaGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTg                                                                                                                                          <  1:180213/141‑1 (MQ=255)    catccgtgtattctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTa                                                                                                                                        <  2:220606‑M1/4‑1 (MQ=255)    aaaaGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTa                                                                                                                                        <  2:189764/141‑1 (MQ=255)      aaGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATc                                                                                                                                      <  2:324544/141‑1 (MQ=255)       ccgtgtattctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttaggggcccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATcc                                                                                                                                     >  1:271076‑M1/135‑141 (MQ=255)       aGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGTTACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATcc                                                                                                                                     <  2:272539/141‑1 (MQ=255)       aGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGTTACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATcc                                                                                                                                     <  2:100631/141‑1 (MQ=255)        gTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCa                                                                                                                                    <  1:317510/141‑1 (MQ=255)         ttGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTACGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTCTATCTTATCGACCGTTTTGTATCCaa                                                                                                                                   >  2:111884/1‑141 (MQ=255)               cacTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGga                                                                                                                             <  2:9603/141‑1 (MQ=255)                          agaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTACTAACa                                                                                                                  <  2:62647‑M1/26‑1 (MQ=255)                                ttattattaGTTGGTGCACTTACTGACAACAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAAt                                                                                                            <  2:221227/141‑1 (MQ=255)                                aaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactcaaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAAt                                                                                                            <  1:197336‑M1/32‑1 (MQ=255)                                      ttaGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTa                                                                                                      >  1:33778/1‑141 (MQ=255)                                        atgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATc                                                                                                    >  2:888‑M1/102‑141 (MQ=255)                                        aGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATc                                                                                                    <  1:316534/141‑1 (MQ=255)                                                              aGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCata                                                                              >  1:243018/1‑141 (MQ=255)                                                                 aattttagaaaggatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCataata                                                                           >  1:302173‑M1/77‑141 (MQ=255)                                                                             tCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTtata                                                               <  1:26308/141‑1 (MQ=255)                                                                               tgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAAc                                                             >  1:260619‑M1/63‑141 (MQ=255)                                                                                  gaaagaTGTGACTGATTTGAAGTTAGTTAGGCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACttt                                                          <  2:4407/141‑1 (MQ=255)                                                                                   aaagaTGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACtttt                                                         >  1:40244/1‑140 (MQ=255)                                                                                     agaTGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACAGTATTGATTACAAACACGCGCTTCCGATCATAGAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTttatc                                                       <  2:333529/141‑1 (MQ=255)                                                                                        ggtgagtGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACAGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTttatcatt                                                    >  2:60908/7‑141 (MQ=255)                                                                                                     aaGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCt                                       <  1:108714/141‑1 (MQ=255)                                                                                                           gttagtCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCAtata                                 >  1:169060/1‑141 (MQ=255)                                                                                                                           atatGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttaggg                 >  1:162675/1‑141 (MQ=255)                                                                                                                               ttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttagggactt             >  2:197336‑M1/15‑141 (MQ=255)                                                                                                                                          ttttGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATtttagggactttttagggactt  <  1:440/141‑1 (MQ=255)                                                                                                                                             |                                                                                                                                          AAAAAAAGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATTTTAGGGACTTTTTAGGGACTT  >  USA300TCH1516_ALE/2085156‑2085433

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.