New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? USA300TCH1516_ALE 2085296 =20 (0.890)5 (0.250) 5/250 2.0 24.4% coding (173/201 nt) hlb_1 Phospholipase C
?USA300TCH1516_ALE = 2128343 13 (0.640)coding (1/825 nt) hlb_2 Phospholipase C

AAAAGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATTTTAGGGACTTTTTAGGGACTT  >  USA300TCH1516_ALE/2085159‑2085433
                                                                                                                                          |                                                                                                                                         
acatccgtgtattctcatttgttagaagaaaaatttaatgaagaggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaatgtttctaTg                                                                                                                                          <  1:35516‑M1/2‑1 (MQ=255)
 aaaGTTGCAACACTTGCATTAGCAAATTTATCATCCAGGCAGACTCTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATCACGTTTATATGTTCTCGACCGTTTTGTa                                                                                                                                        >  2:188076/1‑140 (MQ=255)
   aGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATc                                                                                                                                      >  1:116792/1‑140 (MQ=255)
    gTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGGCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATcc                                                                                                                                     <  2:71673/140‑1 (MQ=255)
             cTTGTATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGgag                                                                                                                            <  2:20140/140‑1 (MQ=255)
              ttccATTAGCAAATTTATTACGAGTTTGTGCACTTACTCACAATACTGCCAAAGCCGAATCTAAGAAAGCTGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGgaga                                                                                                                           <  1:170752/137‑1 (MQ=255)
                         aagaaaaatttaatgaagaggaaaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACAt                                                                                                                 >  1:71099‑M1/114‑140 (MQ=255)
                          aTTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAAATGGAGACGTTTCACATaa                                                                                                               >  2:54017/1‑140 (MQ=255)
                           gaaaaatttaatgaagagtacaataagacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATaa                                                                                                               <  2:161159‑M1/29‑1 (MQ=255)
                                            ggacaaaaaaacaactaaaattttagaaagtatgtaatttagggacccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCaa                                                                                              >  1:2499‑M1/95‑140 (MQ=255)
                                                aaaaaaacaactaaaattttagaaagtatgtaatatagggccccattagggactccaaacccaataaatactgttgttacaaggtttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAg                                                                                          >  2:197251‑M1/91‑140 (MQ=255)
                                                      agaaTAGTGCTAGAGGCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAAc                                                                                    <  2:251649/138‑1 (MQ=255)
                                                           aGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAAAATAAATTTACTTATCATTCAAAAAGGAAAACAGCat                                                                               >  2:117211/1‑140 (MQ=255)
                                                            gTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACAGTTTTGTATCCAAACTGGAGACTTTTAAAATAAAATTACTTATCATTAAACAAGTaaaaacccatc                                                                              >  2:89896/1‑133 (MQ=255)
                                                            gTGCCAAAGCCGAATATAAGAAAGACGATACTGATTTGAAGTTAGTTAGTCATAAAGATTATATGTTATCGACCGTTCTGTATCCAAACTGTAGACTATTCACATAACATTACTTATCATTTGAAAAGTAACACAGCata                                                                              >  2:252450/1‑140 (MQ=255)
                                                               ccAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCataata                                                                           <  1:24693/140‑1 (MQ=255)
                                                                 aaaGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATTATTCAAAAAGTAAAATAGCataataac                                                                         >  2:258144/1‑138 (MQ=255)
                                                                           cTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTtata                                                               <  2:220889/140‑1 (MQ=255)
                                                                                       atACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTttatcatta                                                   <  2:95390/140‑1 (MQ=255)
                                                                                        tACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTttatcattat                                                  <  1:54017/140‑1 (MQ=255)
                                                                                         aCTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTACAACAGCATAATATCAAGGTTTATAACTttatcattatc                                                 >  2:250283/1‑140 (MQ=255)
                                                                                                                        atatGTTATCGACCGTTTTGTATCCAAACTGTAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAATAGCATAATATCAATGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATTAAAttttagg                  <  2:7330/140‑1 (MQ=255)
                                                                                                                                     ttctaTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAAttttagggactttttaggga     <  1:196593‑M1/135‑1 (MQ=255)
                                                                                                                                     cGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACATCATATAAAATAAAAttttagggactttttagcga     >  2:163959/1‑140 (MQ=255)
                                                                                                                                        tttGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATtttagggactttttagggactt  >  2:258797/1‑140 (MQ=255)
                                                                                                                                        tttGTATCCAAACTGGAGACGTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATtttagggactttttagggactt  <  1:145365/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                         
AAAAGTTGCAACACTTGCATTAGCAAATTTATTATTAGTTGGTGCACTTACTGACAATAGTGCCAAAGCCGAATCTAAGAAAGATGATACTGATTTGAAGTTAGTTAGTCATAACGTTTATATGTTATCGACCGTTTTGTATCCAAACTGGAGACTTTTAACATAAAATTACTTATCATTCAAAAAGTAAAACAGCATAATATCAAGGTTTATAACTTTATCATTATCAATAATACCTCATATAAAATAAAATTTTAGGGACTTTTTAGGGACTT  >  USA300TCH1516_ALE/2085159‑2085433

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.