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

AAAAAAGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATTTTAGGGACTTTTTAGGGACTTTA  >  USA300TCH1516_ALE/2085157‑2085435
                                                                                                                                           |                                                                                                                                           
aaaaaaGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGtttt                                                                                                                                             >  1:220582/1‑140 (MQ=255)
 acatccgtgtattctcatttgttagaagaaaactttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTg                                                                                                                                            <  1:38754‑M1/2‑1 (MQ=255)
   aaaGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTa                                                                                                                                          >  1:253336/1‑140 (MQ=255)
       gtgtattctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaatattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCa                                                                                                                                      >  2:108555‑M1/133‑140 (MQ=255)
       gtgtattctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCa                                                                                                                                      <  1:240823‑M1/8‑1 (MQ=255)
         gtattctcatttgttagaagaaaaattgaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCaaa                                                                                                                                    >  2:196271‑M1/131‑140 (MQ=255)
           attctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACt                                                                                                                                  >  1:106493‑M1/129‑140 (MQ=255)
                   cATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACtt                                                                                                                          <  1:242993/140‑1 (MQ=255)
                    aTTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCGGAGAAAGATGATAATGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACttt                                                                                                                         >  2:3792/1‑140 (MQ=255)
                        agaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAAc                                                                                                                     <  1:6264‑M1/25‑1 (MQ=255)
                               aaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAAt                                                                                                              <  1:121850‑M1/32‑1 (MQ=255)
                                            tGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATCACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCa                                                                                                 >  1:67190/1‑140 (MQ=255)
                                                  tACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGt                                                                                           <  1:203731/140‑1 (MQ=255)
                                                    cTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATAATTCAAAAAGTaa                                                                                         >  1:13171/1‑140 (MQ=255)
                                                               gCCAAAGCCGAATCTAAGAAAGATGATACTTATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGAATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCataat                                                                              <  1:69891/140‑1 (MQ=255)
                                                                attttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCataata                                                                             >  1:250346‑M1/76‑140 (MQ=255)
                                                                   aaGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCa                                                                          >  2:99523/1‑140 (MQ=255)
                                                                     gctgAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTTATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAg                                                                        <  2:8097/137‑1 (MQ=255)
                                                                     gCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCAGAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTAGAAACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACATCATAATATCaat                                                                        >  2:184997/1‑139 (MQ=255)
                                                                           tatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTtat                                                                  >  2:238544‑M1/65‑140 (MQ=255)
                                                                           tatgtaagttagggacccattagggactacacaaacagtaagtacttttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTtat                                                                  >  1:113159‑M1/65‑140 (MQ=255)
                                                                              aagaaagaTGATACTGATTTGAAGGTAGTTAGTCATAACGTTTATATGGTATCGACAGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATaaa                                                               >  2:259988/1‑139 (MQ=255)
                                                                                gaaagaTGATACTGATTTGAAGTTAGTTAGTCATAATGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACtt                                                             <  2:65253/140‑1 (MQ=255)
                                                                                          aCTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTttatcattatc                                                   >  1:262674/1‑140 (MQ=255)
                                                                                          aCTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTttatcattatc                                                   >  1:262626/1‑140 (MQ=255)
                                                                                               tagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTataaataa                                              <  2:100655‑M1/96‑1 (MQ=255)
                                                                                                agggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCaataat                                             <  1:55394‑M1/97‑1 (MQ=255)
                                                                                                          aacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACGTCAtata                                   <  1:249128‑M1/107‑1 (MQ=255)
                                                                                                                        tgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAATGTTCAACAGCATAATATCAAGGGTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAATttttag                     <  1:100959‑M1/121‑1 (MQ=255)
                                                                                                                                      cGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGCTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttagggactttttaggga       >  1:120084/1‑140 (MQ=255)
                                                                                                                                         tttGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTTTCAGTGTCAATAATACCTCATATAAACTAACATtttagggactttttatggcctt    >  1:50713/1‑140 (MQ=255)
                                                                                                                                           TGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATTTTAGGAACTTATTAGGGAAttta  <  2:9988‑M1/140‑1 (MQ=255)
                                                                                                                                           TGTATCCAAACTGCAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATTTTAGGGACTTTTTAGGGACttta  >  1:124626‑M1/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                           
AAAAAAGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATTTTAGGGACTTTTTAGGGACTTTA  >  USA300TCH1516_ALE/2085157‑2085435

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.