New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? USA300TCH1516_ALE 2085296 =19 (0.650)4 (0.150) 4/232 NT 17.4% coding (173/201 nt) hlb_1 Phospholipase C
?USA300TCH1516_ALE = 2128343 21 (0.800)coding (1/825 nt) hlb_2 Phospholipase C

AATTTAAAGTCCCTAAAAAGTCCCTAAAATTTTATTTTATATGAGGTATTATTGATAATGATAAAGTTATAAACCTTGATATTATGCTGTTTTACTTTTTGAATGATAAGTAATTTTATGTTAAAAGTCTCCAGTTTGGATACA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  USA300TCH1516_ALE/2085439‑2085296
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccagtttggatacaTAGAAACCTTGTAACAACAGTATTTATTGGGTTTGGAGTCCCTAATGGGTCCCTAAATTACATACTTTCTAAAATTTTAGTTGTTTTTTTGTCCTCTTCATTAAATTTTTCTTCTAACAAATGAGAATACACGGATGT  <  USA300TCH1516_ALE/2128343‑2128206
                                                                                                                                                                                                                                                                                          
AATTTAAAGTCCCTAAAAAGTCCCTAAAATTTTATTTTATATGAGGTATTATTGATAATGATAAAGTTATAAACCTTGATATTATGCTGTTTTACTATTTTAATGGTAAGTAATTTTATTTTAAAAGTCTCCAGTTTGGAT                                                                                                                                               >  1:154400/1‑141
                                          GAGGTATTATTGATAATGATAAAGTTATAAACCTTGATATTATGCTGTTTTACTTTTTGAATGATAAGTAATTTTATGTTAAAAGTCTCCAGTTTGGCTACATAGAAACCTTGTAACAACAGGATTTATTGGGTATGGAtt                                                                                                     >  2:341024/1‑139
                                                 TATTGATAATGATAAAGTTATAAACCTTGATATTATGCTGTTTTACTTTTTGAATGATAAGTAATTTTATGTTAAAAGTCTCCAGTTTGGATACATAGAAACCTTGTAACAACATTAGTTATTGGTTTTGGAGGCCCTAAT                                                                                              >  1:162160/1‑141
                                                                    ATAAACCTTGATATTATGCTGTTTTACTTTTTGAATGATAAGTAATTTTATGTTAAAAGTCTCCAGTTTGGATACATAGAAACCTTGTAACAACAGTATTTATTGGGTTTGGAGTCCCTAATGGGTCCCTAAATTACATAC                                                                           >  1:181835/1‑141
                                                                                                           AAGTAATTTTATGTTAAAAGTCTCCAGTTTGGATACATAGAAACCTTGTAACAACAGTATTTATTGGGTTTGGAGTCCCTAATGGGTCCCTAAATTACATACTTTCTAAAATTTTAGTTGTTTTTTTGTCCTCTTCATTAA                                    <  2:353944/141‑1
                                                                                                                             AGTCTCCAGTTTGGATACATAGAAACCTTGTCACAACAGTAGTTAGTGTGTTTGGAGTCCCTAATGGGTCCCTAAATTACATACTTTCTAAAATTTTAGTTGTTTTTTTGTCCTCTTCATTAAATTTTTCTTCTAACAAAT                  <  2:162160/141‑1
                                                                                                                                     GTTTGGATACATAGAAACCTTGTAACAACAGTATTTATTGGGTTTGGAGTCCCTAATGGGTCCCTAAATTACATACTTTCTAAAATTTTAGTTGTTTTTTTGTCCTCTTCATTAAATTTTTCTTCTAACAAATGAGAATAC          <  1:19937/141‑1
                                                                                                                                          GATACATAGAAACCTTGTAACAACAGTATTTATTGGGTTTGGAGTCCCTAATGGGTCCCTAAATTACATACTTTCTAAAATTTTAGTTGTTTTTTTGTCCTCTTCATTAAATTTTTCTTCTAACAAATGAGAATACACGGA     >  2:316219/1‑141
                                                                                                                                            TACATAGAAACCTTGTAACAACAGTATTTATTGGGTTTGGAGTCCCTAATGGGTCCCTAAATTACATACTTTCTAAAATTTTAGTTGTTTTTTTGTCCTCTTCATTAAATTTTTCTTCTAACAAATGAGAATACACGGATG   <  1:233885/141‑1
                                                                                                                                             ACATAGAAACCTTGTAACAACAGTATTTATTGGGTTTGGAGTCCCTAATGGGTCCCTAAATTACATACTTTCTAAAATTTTAGTTGTTTTTTTGTCCTCTTCATTAAATTTTTCTTCTAACAAATGAGAATACACGGATGT  <  2:181835/141‑1
                                                                                                                                                                                                                                                                                          
AATTTAAAGTCCCTAAAAAGTCCCTAAAATTTTATTTTATATGAGGTATTATTGATAATGATAAAGTTATAAACCTTGATATTATGCTGTTTTACTTTTTGAATGATAAGTAATTTTATGTTAAAAGTCTCCAGTTTGGATACA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  USA300TCH1516_ALE/2085439‑2085296
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccagtttggatacaTAGAAACCTTGTAACAACAGTATTTATTGGGTTTGGAGTCCCTAATGGGTCCCTAAATTACATACTTTCTAAAATTTTAGTTGTTTTTTTGTCCTCTTCATTAAATTTTTCTTCTAACAAATGAGAATACACGGATGT  <  USA300TCH1516_ALE/2128343‑2128206

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.