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

CTACATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAATGTGAAAAAAGAA  >  USA300TCH1516_ALE/2128204‑2128483
                                                                                                                                           |                                                                                                                                             
cTACATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAt                                                                                                                                              <  2:174798‑M2/141‑2 (MQ=255)
   cATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgta                                                                                                                                           <  2:220606‑M2/141‑5 (MQ=255)
      ccGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGGCCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatcc                                                                                                                                        >  1:271076‑M2/1‑134 (MQ=255)
                tcATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTggggg                                                                                                                              <  1:224964/141‑1 (MQ=255)
                    ttGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAat                                                                                                                          >  2:257856/1‑141 (MQ=255)
                      gTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAatat                                                                                                                        >  1:38719/1‑141 (MQ=255)
                         agaagaAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagactactaaca                                                                                                                     <  2:62647‑M2/141‑27 (MQ=255)
                         agaagaAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATaaa                                                                                                                     >  1:242884/1‑141 (MQ=255)
                               aaaaTTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCAAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaat                                                                                                               <  1:197336‑M2/141‑33 (MQ=255)
                                       aTGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatc                                                                                                       >  2:888‑M2/1‑101 (MQ=255)
                                                        acaacTAAAATTTTAGAAAGTATGTCATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCAGACCATATTCTtata                                                                                     >  2:198747/1‑141 (MQ=255)
                                                        acaacTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTtata                                                                                     >  1:228952/1‑141 (MQ=255)
                                                         aaaacaaaaTGTTAGAAAGTATGTAATTTAGGTACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTTGGGGCAATATAATCGCGCTGATTTAATCGGACAATCTTCTTatat                                                                                    <  1:20418/136‑1 (MQ=255)
                                                                aaTTTTAGAAAGGATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataata                                                                              >  1:302173‑M2/1‑76 (MQ=255)
                                                                     tAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATCTTAAAAATAatga                                                                        >  1:77790/1‑141 (MQ=255)
                                                                      aGAAAGGATGTGATTGAGGGACACATTGGGGACTCCAAACCCAATAAATACTGTTATTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAatgat                                                                       <  2:224382/141‑1 (MQ=255)
                                                                           gTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGt                                                                  <  1:187684/141‑1 (MQ=255)
                                                                           atATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGTGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGt                                                                  <  1:220606/140‑1 (MQ=255)
                                                                              tGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataac                                                                >  1:260619‑M2/1‑62 (MQ=255)
                                                                               gTAATTTAGGGCCCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAata                                                              <  1:257856/141‑1 (MQ=255)
                                                                                                       tCCAAACCCAATAAATACTGTTGTTACAAAGTTTCTATGTATCCAAACTGGGGGCAATATAAAAGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGAATTTAATAAAGGt                                      >  2:253663/1‑141 (MQ=255)
                                                                                                          aaaCCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAAACGGTCATGGGAAGCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCa                                   <  1:161932/141‑1 (MQ=255)
                                                                                                           aaCCCAATAAATACTGTTGTTACAAGGTCTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCAt                                  >  2:186245/1‑141 (MQ=255)
                                                                                                                             gttACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTa                <  2:77790/141‑1 (MQ=255)
                                                                                                                              ttACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattatcaataatacctcatataaaataaaattttagggactt                >  2:197336‑M2/1‑14 (MQ=255)
                                                                                                                               tACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAAt              >  2:183240/1‑141 (MQ=255)
                                                                                                                                   aGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAAGATTCAATCAAGCATTTGATAATGGTGCATCAGCCAAATTATTAATTCATGTGa          >  2:207154/1‑141 (MQ=255)
                                                                                                                                           aTGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAATGTGAAAAAAGaa  <  1:186245/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                             
CTACATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAATGTGAAAAAAGAA  >  USA300TCH1516_ALE/2128204‑2128483

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