New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? USA300TCH1516_ALE 2085296 =24 (0.980)13 (0.590) 12/252 0.7 40.9% coding (173/201 nt) hlb_1 Phospholipase C
?USA300TCH1516_ALE = 2128343 16 (0.730)coding (1/825 nt) hlb_2 Phospholipase C

ACTACATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAATGTGAAAAAAGAA  >  USA300TCH1516_ALE/2128203‑2128483
                                                                                                                                            |                                                                                                                                            
aCTACATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTa                                                                                                                                              >  2:136653/1‑141 (MQ=255)
                      tGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAata                                                                                                                        <  1:34393/141‑1 (MQ=255)
                         tAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTTCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATaa                                                                                                                     <  2:35933/141‑1 (MQ=255)
                          agaagaAAAATTTAATGAAGAGGACAACACAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAACTACTGTTGGCACAAGGTTTCTATCTATCCAAACTGGGGGACAtataca                                                                                                                    >  2:9817/1‑139 (MQ=255)
                          agaagaAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaaca                                                                                                                    <  1:210890‑M2/141‑27 (MQ=255)
                          agaagaAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGAACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAACGTTTCTATGTATCCAAACTGGCGGCAATATaaa                                                                                                                    >  2:240419/1‑141 (MQ=255)
                                    tttAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAGTATAAACGCGCTGAtt                                                                                                          >  2:240710/1‑141 (MQ=255)
                                    tttAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACACATTAGGGATTCCAAAACCAATAAATACTGTTGTTAAAAAGTTTCTAtgtatccaaaattaaaacttttaacataaaattact                                                                                                          <  2:42009‑M2/141‑37 (MQ=255)
                                    tttAAGGAAGAGGACAAAAAGACAACGAGTGTTTTGGGTAGGAGATAAGGTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGAtt                                                                                                          <  1:120236/141‑1 (MQ=255)
                                               ggACAAAAAAAGAACTAAAGTTTTAGAAAGGATGTAATTGAGGGATCCACTAGGGACTCCAAACCCAATAAATACTGGTGTTACAAGGTGTCTATGTATCCAAACGGGGGGCAATATAAACGCGCTGATTTAATCGGACaa                                                                                               <  2:48447/141‑1 (MQ=255)
                                                        aacaacTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaac                                                                                      <  2:223414‑M2/141‑57 (MQ=255)
                                                            acTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaacctggagacttttaacataaaattactaatcattcaaaaagtaaaacagca                                                                                  <  1:126825‑M2/141‑61 (MQ=255)
                                                                 aaTGTTAGAAAGTATGTACTTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAaat                                                                             <  1:126561/141‑1 (MQ=255)
                                                                       aGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAatgat                                                                       >  2:9013/1‑141 (MQ=255)
                                                                            ttaTGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttat                                                                  >  2:3022‑M2/2‑65 (MQ=255)
                                                                                gTAATTTAGGGACACATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatcaaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagaataatatcaaggtttataact                                                              <  2:45858‑M2/141‑81 (MQ=255)
                                                                                             ccATTAGGGACTCCAAACCCAATAAGTACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacccccctgtttcaaggtttataactttatcattatcaa                                                 >  2:129327‑M2/1‑48 (MQ=255)
                                                                                             ccATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattatcaa                                                 >  1:96356‑M2/1‑48 (MQ=255)
                                                                                               aTTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCAttt                                               <  1:83496/141‑1 (MQ=255)
                                                                                               aTTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCAttt                                               <  2:13033/141‑1 (MQ=255)
                                                                                                 tAGGGACTCCAAACCCAATTAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGa                                             <  1:158063/141‑1 (MQ=255)
                                                                                                       cTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattatcaataatacctca                                       >  1:272293‑M2/1‑38 (MQ=255)
                                                                                                            aaCCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattatcaataatacctcatataa                                  <  2:195815‑M2/141‑109 (MQ=255)
                                                                                                                 aataaataCTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattatcaataatacctcatataaaataa                             <  2:16957‑M2/141‑114 (MQ=255)
                                                                                                                        aCTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattatcaataatacctcatataaaataaaatttta                      <  2:119387‑M2/141‑121 (MQ=255)
                                                                                                                        aCTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAattat                      <  1:239594/141‑1 (MQ=255)
                                                                                                                                    aGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAATGTGa          >  1:161905/1‑141 (MQ=255)
                                                                                                                                        ttCTATTTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAATGTGaaaaa      <  2:256692/141‑1 (MQ=255)
                                                                                                                                            atgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattatcaataatacctcatataaaataaaattttagggactttttagggacttta  >  1:223414‑M2/1‑1 (MQ=255)
                                                                                                                                            aTGTATTCAAATTGAGTGTAATATAAACGTGTTGATTTAATCTGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTTAATGAAGCATTTTATAATGGTGCATCAGACAAATTATTAAGTAATGTGAAAAAAGaa  <  2:167457/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                            
ACTACATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAATGTGAAAAAAGAA  >  USA300TCH1516_ALE/2128203‑2128483

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.