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

ATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAATGTGAAAAAAGA  >  USA300TCH1516_ALE/2128208‑2128482
                                                                                                                                       |                                                                                                                                           
aTCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtat                                                                                                                                        >  1:301375‑M2/1‑136 (MQ=255)
   cGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTATAAGGCTCCCAtgtatcca                                                                                                                                     <  2:95643‑M2/141‑9 (MQ=255)
   cGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatcca                                                                                                                                     <  2:297842‑M2/141‑9 (MQ=255)
    gtgtATTCTCATTCGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAAATTTAGAAAGTATGTAATATAGGGACCCATTAGGGACTCCAAACCAAATAAATGCTGTTACTATTATGTTTCAAtgtatccaa                                                                                                                                    <  2:257517‑M2/141‑10 (MQ=255)
     tgtATTCTCATTTATTAGTAGAAAAATTTAATGAAAAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaa                                                                                                                                   <  2:347644‑M2/141‑11 (MQ=255)
      gtATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGCTTTCTAtgtatccaacc                                                                                                                                  <  2:314803‑M2/141‑12 (MQ=255)
             cagtggTTAGAAGAAAAAGGTAAGGAAGAGGACAAAAAAACAACTAAAAGTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGc                                                                                                                           <  1:372288/136‑1 (MQ=255)
              atatgatAGCAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCa                                                                                                                          <  1:202205/135‑1 (MQ=255)
                     agaagaAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaaca                                                                                                                   <  1:162945‑M2/141‑27 (MQ=255)
                     agaagaAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATCAGAGACTCAAAACCAAATAAATACTATTGGTCAAGGTTTCTATGTATCCAAACTGGGGGCAGTATaaaa                                                                                                                  >  2:383860/1‑140 (MQ=255)
                          aaaaaTTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACgcgc                                                                                                              <  2:206519/141‑1 (MQ=255)
                              aTTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGAGGGCAATATAAACGTGCTggt                                                                                                          >  2:13374/1‑139 (MQ=255)
                                ttAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACCGGGGGCAATATACACGCGCTGAttt                                                                                                        >  1:139728/1‑141 (MQ=255)
                                ttAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTAAGGGACCAAGTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGAAATATCAACGCGCTGAttt                                                                                                        >  2:373999/1‑141 (MQ=255)
                                      aagaagACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCgg                                                                                                  <  2:106040/141‑1 (MQ=255)
                                        gagGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGAc                                                                                                <  1:22755/141‑1 (MQ=255)
                                        aagGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGAc                                                                                                <  2:152911/140‑1 (MQ=255)
                                               aaaaaaCAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAAtcttc                                                                                         >  2:365844/1‑141 (MQ=255)
                                                  aaaCAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTta                                                                                      <  1:348760/141‑1 (MQ=255)
                                                     gaacGAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGTACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTatat                                                                                   <  2:54283/140‑1 (MQ=255)
                                                     caacTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTatat                                                                                   <  1:157576/141‑1 (MQ=255)
                                                                ttAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaa                                                                        >  2:120384‑M2/1‑72 (MQ=255)
                                                                     aaGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggttt                                                                   <  2:188892‑M2/141‑75 (MQ=255)
                                                                        tATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttata                                                                >  1:121894‑M2/1‑64 (MQ=255)
                                                                        tATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttata                                                                >  2:251658‑M2/1‑64 (MQ=255)
                                                                        tATGAAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTa                                                                <  2:342795/141‑1 (MQ=255)
                                                                         aTGAAATTTAGGGACCCATTAGAGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGAAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTaa                                                               <  1:279642/141‑1 (MQ=255)
                                                                           gTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAata                                                             <  1:13245/141‑1 (MQ=255)
                                                                           gTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAata                                                             <  1:25903/141‑1 (MQ=255)
                                                                           gTAATGTAGGGACCCAGTAGGGACTTCAAAACCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagattaagatcaaggtttataact                                                             >  2:144515‑M2/1‑61 (MQ=255)
                                                                            tAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGCCAATCTTCTAATATTAAAAATAATGATGTCGTAatat                                                            >  2:75069/1‑141 (MQ=255)
                                                                                ttAGGGATCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatc                                                        >  2:235771‑M2/1‑56 (MQ=255)
                                                                                     gACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAg                                                   >  2:21611/1‑141 (MQ=255)
                                                                                       cccATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaacagtaaaacagcataatatcaaggtttataactttatcattatca                                                 >  1:86469‑M2/1‑49 (MQ=255)
                                                                                          aTTAGAGACTGCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAGAAATAATGATGTCGTAATATTCAATGAAGCAttt                                              <  2:135788/141‑1 (MQ=255)
                                                                                           ttAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTg                                             >  1:335904/1‑141 (MQ=255)
                                                                                            tAGAGACTCCAAACCCAATAAAGACTGGTGTTAAAAGGTTGCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGa                                            <  2:164421/141‑1 (MQ=255)
                                                                                             aGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCATTATAAACGCGCTGATTTAATCGGACAATCTTCTTATCTTAAAAATACTGATGTCGTCATATTCAATGCAGCATTTCAt                                           >  1:151672/1‑141 (MQ=255)
                                                                                             aGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGAt                                           >  2:323941/1‑141 (MQ=255)
                                                                                              gggACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattatcaataatac                                          <  1:311024‑M2/141‑100 (MQ=255)
                                                                                                      aaaCCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCa                                  >  1:193083/1‑141 (MQ=255)
                                                                                                          ccAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAg                              <  2:202411/141‑1 (MQ=255)
                                                                                                          ccAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAg                              <  1:368787/141‑1 (MQ=255)
                                                                                                               aaataCTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACaaa                         >  1:88231/1‑141 (MQ=255)
                                                                                                               aaataCTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACaaa                         <  2:181681/141‑1 (MQ=255)
                                                                                                                  taCTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGACGATATAAACGGGCTGATTTAAGCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAatta                      <  1:73549/141‑1 (MQ=255)
                                                                                                                               aGGTTTCTATGTATGCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAATGTGa         <  2:18139/141‑1 (MQ=255)
                                                                                                                               aGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAATGTGa         >  1:261890/1‑141 (MQ=255)
                                                                                                                                ggTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAATGTGaa        <  1:297842/141‑1 (MQ=255)
                                                                                                                                ggTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAATGTGaa        >  1:102763/1‑141 (MQ=255)
                                                                                                                                ggTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAATGTGaa        <  1:95643/141‑1 (MQ=255)
                                                                                                                                  tttCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAATGTGaaaa      >  2:120414/1‑141 (MQ=255)
                                                                                                                                      tAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattatcaataatacctcatataaaataaaattttagggactttttagggacttt  >  1:347644‑M2/1‑2 (MQ=255)
                                                                                                                                       |                                                                                                                                           
ATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAATGTGAAAAAAGA  >  USA300TCH1516_ALE/2128208‑2128482

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.