New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? USA300TCH1516_ALE 2085296 =20 (1.080)15 (0.910) 13/228 0.3 52.8% coding (173/201 nt) hlb_1 Phospholipase C
?USA300TCH1516_ALE = 2128343 9 (0.550)coding (1/825 nt) hlb_2 Phospholipase C

ATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAAT  >  USA300TCH1516_ALE/2128216‑2128471
                                                                                                                               |                                                                                                                                
aTTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTTTGTAATTTAGGGACCCATTAGGGACGGCAAACCCAAGAAATACTGTTGTTACAAGGTTTCTAt                                                                                                                                 >  1:137935‑M2/1‑128 (MQ=255)
  tctcATTTGTTAGAAGAAAAATTTTAGGATGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgt                                                                                                                               >  1:133034‑M2/1‑126 (MQ=255)
      aTTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTAAGAAAGTATGTAATTTAGGGACCCATTAGGGACACTAAACCCAATAAATACTGGTGTTACAAGGTGTGTAGGTactc                                                                                                                           >  2:198609/1‑126 (MQ=255)
       tttGTTAGAAGAAAAAATTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAAGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGATTCAAtatatcca                                                                                                                          <  2:99724‑M2/129‑9 (MQ=255)
            tAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactg                                                                                                                     <  1:176238‑M2/129‑14 (MQ=255)
                        ttAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGAGACACACTAGGGACTCCAAACCCAATAAATACTGTTATTACAAGATTTCTAtgtatccaaaatgaagacttttaac                                                                                                         <  2:101977‑M2/129‑26 (MQ=255)
                               agagGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACACAATAAATACTGTTGTTACAACGTTTCTAtgtatccaaactggagacttttaacatacaat                                                                                                  <  1:228195‑M2/129‑33 (MQ=255)
                                           aacaacTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattc                                                                                      <  1:228735‑M2/129‑45 (MQ=255)
                                              aacTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGAACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACaa                                                                                   >  1:195526/1‑129 (MQ=255)
                                                                      aTTTAGTGACCCATTTGTGACTCCAAACCCAATAAAGACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAatgat                                                           <  2:198101/129‑1 (MQ=255)
                                                                         tAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggttt                                                        <  2:118172‑M2/129‑75 (MQ=255)
                                                                                  aTTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGTGGCAATATATACGCGCTTATTTAATTTGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTc                                               <  1:20348/129‑1 (MQ=255)
                                                                                  aTTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGTGCTGATTTAATTTGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTc                                               <  1:20345/129‑1 (MQ=255)
                                                                                    tAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatc                                             <  2:235264‑M2/129‑86 (MQ=255)
                                                                                     aGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattaaaaaagtaaaacagcataatatcaaggtttataactttatca                                            <  1:91027‑M2/129‑87 (MQ=255)
                                                                                       ggACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGa                                          <  2:108628/129‑1 (MQ=255)
                                                                                         aCTCCAAACCCAATAAATACTGTTGTTACAAGGTTTAAAAATAACCAAAATAGAAGCAAAATAAACGCACAGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAg                                        >  1:189771/1‑129 (MQ=255)
                                                                                            ccAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactgaagacttttaacataaaattacttataattcaaaaagtaaaacagcataatatcaaggtttataactttatcattatcaa                                     <  1:94068‑M2/129‑94 (MQ=255)
                                                                                               aaCCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttatcataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattatcaataa                                  >  1:134678‑M2/1‑33 (MQ=255)
                                                                                               aaCCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGAGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTg                                  >  1:158124/1‑129 (MQ=255)
                                                                                                aCCCAATAAATACTGTTGTTACAAGGTCTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGa                                 >  2:233178/1‑129 (MQ=255)
                                                                                                   cAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAGACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATaa                              >  1:187064/1‑129 (MQ=255)
                                                                                                    aataaataCTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattatcaataatacct                             <  1:189658‑M2/129‑102 (MQ=255)
                                                                                                    aataaataCTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattatcaataatacct                             <  2:199464‑M2/129‑102 (MQ=255)
                                                                                                    aataaataCTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaaatacaaagaactaagatcaaggtttataactttatcattattaataatacct                             <  2:166500‑M2/129‑102 (MQ=255)
                                                                                                                  ttACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattatcaataatacctcatataaaataaaa               >  2:91027‑M2/1‑14 (MQ=255)
                                                                                                                     cAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattccaaaagtaaaacaccataatatcaagctttataactttatcattctcaataatacctcatataaaataaaatct            <  2:9080‑M2/129‑119 (MQ=255)
                                                                                                                           ttCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtctataactttatcattatcaataatacctcatataaaataaaattttaggga      >  1:118172‑M2/1‑5 (MQ=255)
                                                                                                                               atgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattatcaataatacctcatataaaataaaattttagggacttt  <  2:211888‑M2/129‑129 (MQ=255)
                                                                                                                               |                                                                                                                                
ATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATTAAGTAAT  >  USA300TCH1516_ALE/2128216‑2128471

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