New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? USA300TCH1516_ALE 2085296 =19 (0.930)4 (0.220) 4/252 2.0 17.7% coding (173/201 nt) hlb_1 Phospholipase C
?USA300TCH1516_ALE = 2128343 20 (1.090)coding (1/825 nt) hlb_2 Phospholipase C
Rejected: Frequency below/above cutoff threshold.

AAAAAGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATTTTAGGGACTTTTTAG  >  USA300TCH1516_ALE/2085158‑2085427
                                                                                                                                          |                                                                                                                                   
aaaaaGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGt                                                                                                                                   >  2:214297/1‑141 (MQ=255)
 catccgtgtattctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTa                                                                                                                                  <  1:41347‑M1/4‑1 (MQ=255)
 aaaccgtgaattctaatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTa                                                                                                                                  >  2:213308‑M1/138‑141 (MQ=255)
      ttGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCaa                                                                                                                             <  1:191487/141‑1 (MQ=255)
            tctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGga                                                                                                                       >  1:197072‑M1/127‑141 (MQ=255)
                   aTTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTCGTCATAACGTTTATATGTTATTGAACGTTTTGTAGTCAAACTGGAGTCgttt                                                                                                                >  1:13085/1‑141 (MQ=255)
                   aTTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACtttt                                                                                                                >  1:82377/1‑141 (MQ=255)
                   aTTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACtttt                                                                                                                >  1:5229/1‑141 (MQ=255)
                     tAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTaa                                                                                                              <  1:143812/141‑1 (MQ=255)
                           aTTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATaaa                                                                                                        <  2:226955/141‑1 (MQ=255)
                              tattattaGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAAtt                                                                                                     <  1:219303/141‑1 (MQ=255)
                                     aGTGGGGGCACTGACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATc                                                                                              <  2:5302/141‑1 (MQ=255)
                                                ttACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGt                                                                                   >  2:35409/1‑141 (MQ=255)
                                                  aaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAATaa                                                                                 <  2:192341‑M1/53‑1 (MQ=255)
                                                        aaTAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGc                                                                           >  1:210274/1‑141 (MQ=255)
                                                             tGCCAACGCTGAATCTAATAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCataat                                                                      <  2:140411/141‑1 (MQ=255)
                                                             ggCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCataat                                                                      <  1:50991/140‑1 (MQ=255)
                                                              gCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCataata                                                                     >  1:121132/1‑141 (MQ=255)
                                                                              agaaagaTGATACTGATTTGAAGTGAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACtt                                                     <  1:197216/141‑1 (MQ=255)
                                                                                          ccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCaa                                         >  2:53063‑M1/49‑141 (MQ=255)
                                                                                           tGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCaat                                        >  1:29742/1‑141 (MQ=255)
                                                                                           tGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCaat                                        <  1:195075/141‑1 (MQ=255)
                                                                                                                           tGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttagggact        <  2:210274/141‑1 (MQ=255)
                                                                                                                               caaggtttctaTGTATCCAGACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttagggacttttt    <  2:151752‑M1/130‑1 (MQ=255)
                                                                                                                                 cGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttagtgactttttag  >  1:204109/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                   
AAAAAGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATTTTAGGGACTTTTTAG  >  USA300TCH1516_ALE/2085158‑2085427

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.