New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? USA300TCH1516_ALE 2085296 =17 (0.750)9 (0.440) 8/250 1.1 38.9% coding (173/201 nt) hlb_1 Phospholipase C
?USA300TCH1516_ALE = 2128343 13 (0.640)coding (1/825 nt) hlb_2 Phospholipase C

ACATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATT  >  USA300TCH1516_ALE/2128206‑2128464
                                                                                                                                         |                                                                                                                         
aCATCCGTGTATTCTCATTTGTTAGAAGAAAACTTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtg                                                                                                                         <  1:38754‑M2/140‑3 (MQ=255)
   tcggtgtATTCTTATTTGTTAGAAGAAAAATTTAATGAAGAGTACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTAt                                                                                                                      <  2:203532/137‑1 (MQ=255)
      gtgtATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAATATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatcca                                                                                                                   >  2:108555‑M2/1‑132 (MQ=255)
      gtgtATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatcca                                                                                                                   <  1:240823‑M2/140‑9 (MQ=255)
        gtATTCTCATTTGTTAGAAGAAAAATTGAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaa                                                                                                                 >  2:196271‑M2/1‑130 (MQ=255)
          aTTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaact                                                                                                               >  1:106493‑M2/1‑128 (MQ=255)
                       agaagaAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaac                                                                                                  <  1:6264‑M2/140‑26 (MQ=255)
                              aaaTTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaat                                                                                           <  1:121850‑M2/140‑33 (MQ=255)
                                          gagGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGa                                                                               >  1:45744/1‑140 (MQ=255)
                                           agGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAATCCAATTTTTTCTGTTTTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCTGAc                                                                              <  1:97582/140‑1 (MQ=255)
                                           agGACAAAAAAACAACTAAAAGTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAATCCCAAGAAATCCTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGCCCATCTATATGTGCTGATTTAATTGGAc                                                                              <  2:46312/140‑1 (MQ=255)
                                              aCAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGCACAAt                                                                           >  1:60402/1‑140 (MQ=255)
                                                      acaacTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGAGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATAGAACAATCTCCTtat                                                                   >  2:242831/1‑140 (MQ=255)
                                                               aTTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataata                                                          >  1:250346‑M2/1‑75 (MQ=255)
                                                                          tATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttat                                               >  2:238544‑M2/1‑64 (MQ=255)
                                                                          tATGTAAGTTAGGGACCCATTAGGGACTACACAAACAGTAAGTACTTTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttat                                               >  1:113159‑M2/1‑64 (MQ=255)
                                                                                       gaaccATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGaa                                  <  2:258722/137‑1 (MQ=255)
                                                                                       gACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTCATATTTAATGaa                                  >  1:40180/1‑140 (MQ=255)
                                                                                             ttAGTGACTCCAAAGCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCAttt                            <  2:256486/140‑1 (MQ=255)
                                                                                              tAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattataaataa                           <  2:100655‑M2/140‑97 (MQ=255)
                                                                                               aGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattatcaataat                          <  1:55394‑M2/140‑98 (MQ=255)
                                                                                                         aaCCCAATAAATACTGTTGTTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaaagtaaaacagcataatatcaaggtttataactttatcattatcaataatacgtcatata                <  1:249128‑M2/140‑108 (MQ=255)
                                                                                                          aCCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACAGGGGTGGATATACACGCGCTTATGTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCCTTTGATACTGGTGAAt               >  2:232208/1‑140 (MQ=255)
                                                                                                            ctaataagtaCTGTTGTTACAATGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCa             <  2:232815/138‑1 (MQ=255)
                                                                                                            ccAATAAATACTGTTATTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACCATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCa             <  2:171071/140‑1 (MQ=255)
                                                                                                                 aaagaCTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACaa        <  1:26155/140‑1 (MQ=255)
                                                                                                                      cTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAattat   <  2:221440/140‑1 (MQ=255)
                                                                                                                       tgttgtTACAAGGTTTCTAtgtatccaaactggagacttttaacataaaattacttatcattcaaaatgttcaacagcataatatcaagggttataactttatcattatcaataatacctcatataaaataaatttttag  <  1:100959‑M2/140‑122 (MQ=255)
                                                                                                                                         |                                                                                                                         
ACATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGCATCAGACAAATTATT  >  USA300TCH1516_ALE/2128206‑2128464

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.