New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? USA300TCH1516_ALE 2085296 =31 (0.810)9 (0.260) 9/252 2.5 22.0% coding (173/201 nt) hlb_1 Phospholipase C
?USA300TCH1516_ALE = 2128343 36 (1.050)coding (1/825 nt) hlb_2 Phospholipase C

AAAAAAGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATTTTAGGGACTTTTTAGGGACTTTA  >  USA300TCH1516_ALE/2085157‑2085435
                                                                                                                                           |                                                                                                                                           
tacatccgtgtattctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTg                                                                                                                                            >  1:330258‑M1/140‑141 (MQ=255)
tacatccgtgtattctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTg                                                                                                                                            <  2:102750‑M1/2‑1 (MQ=255)
tacatccgtgtattctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggcttctaTg                                                                                                                                            <  2:295990‑M1/2‑1 (MQ=255)
 acatccgtgtattctcatttgttagaagaaaaatttaatgaagcggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaatcaatactgttgttacaaggtttctaTGt                                                                                                                                           <  1:148838‑M1/3‑1 (MQ=255)
 acatccgtgtattctcatttgttagaagaaaaatttaatgaagcggacaaaaaaacaactaaaattttagaaagtatgtaatttagcgacccattagggactccacacacaataaatactgttgttaccagatttctaTGt                                                                                                                                           <  2:114940‑M1/3‑1 (MQ=255)
 acatccgtgtattctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactattgttacaatgtttctcTGt                                                                                                                                           <  1:203274‑M1/3‑1 (MQ=255)
       ttGGAACACTTGCATTAGCAAATTTATGAGTAGTTGGTGCACTTATTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCaa                                                                                                                                     <  2:248490/141‑1 (MQ=255)
        tgtatgatcaattgtaagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCaaa                                                                                                                                    >  2:363307‑M1/132‑141 (MQ=255)
          tattctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACt                                                                                                                                  <  1:317769‑M1/12‑1 (MQ=255)
          cAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACt                                                                                                                                  <  1:356868/141‑1 (MQ=255)
             cacTTGCATTAGTAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGga                                                                                                                               <  2:1462/141‑1 (MQ=255)
             aacTTGCATTAGCAAATGTATTATTAGTTTGTGCACGTACTGACAATAGTGCAAAAGCCTAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGga                                                                                                                               <  2:144020/140‑1 (MQ=255)
              acTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGgag                                                                                                                              <  1:103564/141‑1 (MQ=255)
                   cATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTCTATCCAAACTGGAGACttt                                                                                                                         >  1:245909/1‑141 (MQ=255)
                                  tttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACt                                                                                                          <  1:124801‑M1/36‑1 (MQ=255)
                                      agctGGTGCACTCACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAACCTGGAGACTTTTAACATAAAATTACTTATc                                                                                                      <  1:211644/138‑1 (MQ=255)
                                            tGCACTTACTGACAATAGTGCGAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCaa                                                                                                <  2:311466/141‑1 (MQ=255)
                                            ggCACATATTGAAAATAGTGCGAAAGCCGAATATAAAAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCaa                                                                                                <  2:83587/140‑1 (MQ=255)
                                            aggacaaaaaaacaactcacagatgagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCaa                                                                                                >  1:275904‑M1/96‑141 (MQ=255)
                                                              tGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAACAAGTAAAACAGCataat                                                                              >  2:391858/1‑141 (MQ=255)
                                                               gCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCataata                                                                             <  1:386786/141‑1 (MQ=255)
                                                                      ccGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGt                                                                      >  2:56229/1‑141 (MQ=255)
                                                                      acGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTAGCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGt                                                                      <  2:137009/140‑1 (MQ=255)
                                                                        aagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGttt                                                                    >  2:254433‑M1/68‑141 (MQ=255)
                                                                         aaTCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTta                                                                   >  1:314592/1‑141 (MQ=255)
                                                                         aaTCTAAGAAAGATGATACTGATGTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTta                                                                   <  2:283867/141‑1 (MQ=255)
                                                                          aTCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTtat                                                                  >  1:292529/1‑141 (MQ=255)
                                                                          aTCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTtat                                                                  >  1:139164/1‑141 (MQ=255)
                                                                                        atACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTttatcattat                                                    >  1:221513/1‑141 (MQ=255)
                                                                                         acccattaggtactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTttatcattatc                                                   >  2:128857‑M1/51‑141 (MQ=255)
                                                                                           ccattagggactccaaacctcataaatagtgttgttacaaggtttctaTGTATCAAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATCTCAAGGTTTATAACTTTATCATTATCaa                                                 >  2:176046‑M1/49‑141 (MQ=255)
                                                                                             gATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCaata                                               >  2:145856/1‑141 (MQ=255)
                                                                                                            agtCATAACGTTTATATGTTATCGACCGTTATGTATCCAAACTGGAGAATTTTAACATAAAATTACTTATCATTCAAAAAGGAAAACAGTATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATataaaa                                <  2:266036/141‑1 (MQ=255)
                                                                                                                  aaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTAAAAAGGTGTTGTTTTATTATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATaaaataaaat                          <  1:207666‑M1/116‑1 (MQ=255)
                                                                                                                   aCGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAtt                         >  2:166056/1‑141 (MQ=255)
                                                                                                                   aCGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAtt                         >  1:240280/1‑141 (MQ=255)
                                                                                                                    atactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttt                        <  1:61778‑M1/118‑1 (MQ=255)
                                                                                                                     cgtaatatGTTATCCACCGGTTTCTATCGAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAATAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAtttt                       <  1:93099/137‑1 (MQ=255)
                                                                                                                       ttATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttag                     >  2:65609/1‑141 (MQ=255)
                                                                                                                             ttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttagggactt               >  2:124801‑M1/15‑141 (MQ=255)
                                                                                                                                   gACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttagggactttttagg         >  2:330556/1‑141 (MQ=255)
                                                                                                                                   gACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATAGCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttagggactttttagg         <  1:1923/141‑1 (MQ=255)
                                                                                                                                   gACCGTTTCGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttagggactttttagg         <  1:117933/141‑1 (MQ=255)
                                                                                                                                     ccGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttagggactttttaggga       >  1:190138/1‑141 (MQ=255)
                                                                                                                                     ccGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttagggactttttaggga       >  1:365269/1‑141 (MQ=255)
                                                                                                                                      cGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATAACTCATATAAAATAAAAtttaaggaactttttagagaa      >  2:169986/1‑140 (MQ=255)
                                                                                                                                        ctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATtttagggactttttagggactt    <  1:163975‑M1/138‑1 (MQ=255)
                                                                                                                                          ttGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAACTCAAATTGTATGGACTTTTTAGGGACttta  >  2:70328/1‑141 (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 < 22 ≤ ATCG/ATCG < 37 ≤ 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.