New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? USA300TCH1516_ALE 2085296 =29 (1.040)4 (0.160) 4/250 3.0 13.8% coding (173/201 nt) hlb_1 Phospholipase C
?USA300TCH1516_ALE = 2128343 24 (0.960)coding (1/825 nt) hlb_2 Phospholipase C
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GTAATTTTAAATTTAAAGTCCCTAAAAAGTCCCTAAAATTTTATTTTATATGAGGTATTATTGATAATGATAAAGTTATAAACCTTGATATTATGCTGTTTTACTTTTTGAATGATAAGTAATTTTATGTTAAAAGTCTCCAGTTTGGATACA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  USA300TCH1516_ALE/2085448‑2085296
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccagtttggatacaTAGAAACCTTGTAACAACAGTATTTATTGGGTTTGGAGTCCCTAATGGGTCCCTAAATTACATACTTTCTAAAATTTTAGTTGTTTTTTTGTCCTCTTCATTAAATTTTTCTTCTAACAAATGAGAATA  <  USA300TCH1516_ALE/2128343‑2128215
                                                                                                                                                                                                                                                                                          
GTAATTTTAAATTTAAAGTCCCTAAAAAGTCCCTAAAATTTTATTTTATATGAGGTATTATTGATAATGATAAAGTTATAAACCTTGATATTATGCTGTTTTACTTTTTGAATGATAAGTAATTTTATGTTAAAAGTCTC                                                                                                                                                <  2:267311/140‑1
       TAAATTTAAAGTCCCTAAAAAGTCCCTAAAATTTTATTTTATATGAGGTATTATTGATAATGATAAAGTTATAAACCTTGATATTATGCTGTTTTACTTTTTGAATGATAAGTAATTTTATTTTAAAAGTCTCCAGTTTG                                                                                                                                         >  1:184796/1‑140
           TTTAAATTCCCTAAAAAGTTCCTATAATTTTATTTTATATGAGGTATTATTTATAATGATAAAGTTATTAACCTTGATATTATGCTGTTTTACTTTTTGAATGATAAGTAATTTTATGTTAAAAGTCTCCAGTTTGGATA                                                                                                                                     <  1:104571/140‑1
                                                                            TATAAACCTTGATATTATGCTGTTTTACTTTTTGAATGATAAGTAATTTTATGTTAAAAGTCTCCAGTTTGGATACATAGAAACCTTGTAACAACAGTATTTATTGGGTTTGGAGTCCCTAATGGGTCCCTAAATTACAT                                                                    >  2:124499/1‑140
                                                                                          TTATGCTGTTTTACTTTTTGAATGATAAGTAATTTTATGTTAAAAGTCTCCAGTTTGGATACATAGAAACCTTGTAACAACAGTATTTATTGGGTTTGGAGTCCCTAATGGGTCCCTAAATTACATACTTTCTAAAATTT                                                      <  1:267180/140‑1
                                                                                                                 GATAAGTAATTTTATGTTAAAAGTCTCCAGTTTGGATACATAGAAACCTTGTAACAACAGTATTTATTGGGTTTGGAGTCCCTAATGGGTCCCTAAATTACATACTTTCTAAAATTTTAGTTGTTTTTTTGTCCTCTTCA                               >  2:179635/1‑140
                                                                                                                                         CTCCAGTTTGGATACATAGAAAAGTAGTTACTACATGGTTTATTGGGTTTGGAGTCCCTAATGGGTCCCTAAATTACATACTTTCTAAAATTTTAGTTGTTTTTTTGTCCTCTTCATTAAATTTTTCTTCTAACAAATGA       >  1:312424/1‑140
                                                                                                                                              GTTTGGATACATAGAAACCTTGTAACAACAGTATTTATTGGGTTTGGAGTCCCTAATGGGTCCCTAAATTACATACTTTCTAAAATTTTATTTGTTTTTTTGTCCTCTTCATTAAATTATCATTCTAACAAATGAGAATA  >  1:43778/1‑140
                                                                                                                                                                                                                                                                                          
GTAATTTTAAATTTAAAGTCCCTAAAAAGTCCCTAAAATTTTATTTTATATGAGGTATTATTGATAATGATAAAGTTATAAACCTTGATATTATGCTGTTTTACTTTTTGAATGATAAGTAATTTTATGTTAAAAGTCTCCAGTTTGGATACA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  USA300TCH1516_ALE/2085448‑2085296
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccagtttggatacaTAGAAACCTTGTAACAACAGTATTTATTGGGTTTGGAGTCCCTAATGGGTCCCTAAATTACATACTTTCTAAAATTTTAGTTGTTTTTTTGTCCTCTTCATTAAATTTTTCTTCTAACAAATGAGAATA  <  USA300TCH1516_ALE/2128343‑2128215

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.