New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? USA300TCH1516_ALE 2085296 =29 (0.900)9 (0.310) 8/232 NT 22.8% coding (173/201 nt) hlb_1 Phospholipase C
?USA300TCH1516_ALE = 2128343 35 (1.210)coding (1/825 nt) hlb_2 Phospholipase C

AAAGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATTTTAGGGACTTTTTAGGGACTTT  >  USA300TCH1516_ALE/2085160‑2085434
                                                                                                                                        |                                                                                                                                          
atccgtgtattctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTAt                                                                                                                                        >  1:301375‑M1/137‑141 (MQ=255)
   cgtgtattctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttataaggctcccaTGTATCCa                                                                                                                                     <  2:95643‑M1/8‑1 (MQ=255)
   cgtgtattctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCa                                                                                                                                     <  2:297842‑M1/8‑1 (MQ=255)
    gtgtattctcattcgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaaatttagaaagtatgtaatatagggacccattagggactccaaaccaaataaatgctgttactattatgtttcaaTGTATCCaa                                                                                                                                    <  2:257517‑M1/9‑1 (MQ=255)
     tgtattctcatttattagtagaaaaatttaatgaaaaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCaaa                                                                                                                                   <  2:347644‑M1/10‑1 (MQ=255)
      gtattctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaagctttctaTGTATccaacc                                                                                                                                  <  2:314803‑M1/11‑3 (MQ=255)
                     gCAAATTTATTGTTAGATTTGATACTTACTGACAATAGTGCAAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACa                                                                                                                   <  2:360585/141‑1 (MQ=255)
                     gCAAATTTATTATTAGTTGGTGAAATTACTGACAATAGTGCAAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACa                                                                                                                   <  2:168747/141‑1 (MQ=255)
                     agaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACa                                                                                                                   <  1:162945‑M1/26‑1 (MQ=255)
                       aaaTTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATa                                                                                                                 <  2:293546/141‑1 (MQ=255)
                           ttactattaTTATGTGCAGTTAGTGACAATAGTGCGAAAGCCGAAGCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAAt                                                                                                             <  2:175964/141‑1 (MQ=255)
                             attattaGTTGTTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTa                                                                                                           <  2:235534/141‑1 (MQ=255)
                                            aCTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCaaaaa                                                                                            >  1:306374/1‑141 (MQ=255)
                                             cTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCAAAACTGGAGACTTTTAACATAAAATTACTTATCCTTaaaaaaa                                                                                           >  2:215002/1‑140 (MQ=255)
                                                  tGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCACAAAGTcaaa                                                                                      >  2:3116/1‑141 (MQ=255)
                                                  tGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTaaaa                                                                                      <  2:115197/141‑1 (MQ=255)
                                                     cAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAg                                                                                   <  1:109648/141‑1 (MQ=255)
                                                       aTAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCa                                                                                 <  1:271498/141‑1 (MQ=255)
                                                          gTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCataa                                                                              >  2:356125/1‑141 (MQ=255)
                                                             ccAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAatat                                                                           <  1:49769/141‑1 (MQ=255)
                                                                ttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCaa                                                                        >  2:120384‑M1/73‑141 (MQ=255)
                                                                    cGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGtt                                                                    <  1:211282/141‑1 (MQ=255)
                                                                     aagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGttt                                                                   <  2:188892‑M1/74‑1 (MQ=255)
                                                                        tatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTtata                                                                >  2:251658‑M1/65‑141 (MQ=255)
                                                                        tatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTtata                                                                >  1:121894‑M1/65‑141 (MQ=255)
                                                                          tAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAAc                                                              >  2:333912/1‑141 (MQ=255)
                                                                           gtaatgtagggacccagtagggacttcaaaaccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGATTAAGATCAAGGTTTATAACt                                                             >  2:144515‑M1/62‑141 (MQ=255)
                                                                               aagaTGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTttat                                                         <  1:175983/141‑1 (MQ=255)
                                                                                ttagggatccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTttatc                                                        >  2:235771‑M1/57‑141 (MQ=255)
                                                                                agaTGATACTGATTTGAAGTTAGTTAGTCATATCGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTttatc                                                        <  2:111581/141‑1 (MQ=255)
                                                                                agaTGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTttatc                                                        <  1:317836/141‑1 (MQ=255)
                                                                                 gatgatACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACATCATAATATCAAGGTTTATAAATttataa                                                       >  1:18993/1‑139 (MQ=255)
                                                                                    gatACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTGTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATAAAGGTTTATAACTttatcatta                                                    >  2:163252/1‑141 (MQ=255)
                                                                                       cccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAACAGTAAAACAGCATAATATCAAGGTTTATAACTttatcattatca                                                 >  1:86469‑M1/50‑141 (MQ=255)
                                                                                        cTGATTTGAATGTAGTTAGTCATAGCTTTTATATGTTATCGACCGTTTTATATCCAAACTTGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCaa                                                <  2:288869/141‑1 (MQ=255)
                                                                                         tGATTTGAAGTTAGTTAGTCAAAACGTTTATATGTTATAGAACGTTTAATATCCAAACTGGAGACTTTTAACATAAAATTAAAAAAAAAAAAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCaat                                               <  2:216443/141‑1 (MQ=255)
                                                                                             ttGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTCGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATAATTACCAataata                                           >  2:225295/1‑141 (MQ=255)
                                                                                              gggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATAc                                          <  1:311024‑M1/99‑1 (MQ=255)
                                                                                                            cATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACATTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTGATataaaataa                            >  2:330583/1‑141 (MQ=255)
                                                                                                                     tataTGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttagg                   <  1:82324/141‑1 (MQ=255)
                                                                                                                      atatGTTATCGACCGTTTGGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCGCATATAAAATAAAAttttaggg                  <  2:137991/141‑1 (MQ=255)
                                                                                                                          gTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttagggactt              <  1:356125/141‑1 (MQ=255)
                                                                                                                               cGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttagtcattttttag         >  2:273359/1‑141 (MQ=255)
                                                                                                                               cGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttagggactttttag         >  2:48335/1‑141 (MQ=255)
                                                                                                                                    gTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttagggactttttagtgcct    >  1:195028/1‑136 (MQ=255)
                                                                                                                                      taTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATTttagggactttttagggacttt  >  1:347644‑M1/3‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                          
AAAGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATTTTAGGGACTTTTTAGGGACTTT  >  USA300TCH1516_ALE/2085160‑2085434

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