New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AM260480 = 288777837 (1.120)6 (0.190) 6/208 2.0 14.0% coding (1513/1683 nt) h16_B2557 3‑(3‑hydroxy‑phenyl) propionate hydroxylase
?AM260480 2887792 = 38 (1.190)coding (1527/1683 nt) h16_B2557 3‑(3‑hydroxy‑phenyl) propionate hydroxylase
Rejected: Frequency below/above cutoff threshold.

GACACACGTGGTGGTCACCGGCGGCAACGTGGCCGGCGCCACGGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AM260480/2887732‑2887778
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgcCGGCCAGGCCCACGCCCGCTACGGCGCACGACCCGGCACGCTCTACCTGATCCGCCCCGATGGCTATGTCCTGGCACGCTGGCGCGAACCGTCATGGGACGCCGTGCGCGCCACCCTTGCCCCCCTGATGC  >  AM260480/2887792‑2887922
                                                                                                                                                                                  
GACACACGTGGTGGTCACCGGCGGCAACGTGGCCAGCGCCGTGGCGCCGGCCAGGC                                                                                                                            >  1:1023172/1‑56
GACACACGTGGTGGTCACCGGCGGCAACGTGGCCAGCGCCGTGGCGCCGGCCAGGC                                                                                                                            <  2:1023172/56‑1
          GTGGTCACCGGCGGCAACGTGGCCAGCGCCGTGGCGCCGGCCAGGCCCACGCCCGCTACGGCGCACGACCCGGCACGCTCTACCTGATCCGCCCCGATGGCTATGTCCTGGCACGC                                                      <  1:1108671/116‑1
          GTGGTCACCGGCGGCAACGTGGCCAGCGCCGTGGCGCCGGCCAGGCCCACGCCCGCTACGGCGCACGACCCGGCACGCTCTACCTGATCCGCCCCGATGGCTATGTCCTGGCACGC                                                      >  2:1108671/1‑116
           TGGTCACCGGCGGCAACGTGGCCAGCGCCGTGGCGCCGGCCAGGCCCACGCCCGCTACGGCGCACGACCCGGCACGCTCTACCTGATCCGCCCCGAT                                                                        <  1:186885/97‑1
           TGGTCACCGGCGGCAACGTGGCCAGCGCCGTGGCGCCGGCCAGGCCCACGCCCGCTACGGCGCACGACCCGGCACGCTCTACCTGATCCGCCCCGAT                                                                        >  2:186885/1‑97
                                             GCCGGCCAGGCCCACGCCCGCTACGGCGCACGACCCGGCACGCTCTACCTGATCCGCCCCGATGGCTATGTCCTGGCACGCTGGCGCGAACCGTCATGGGACGCCGTGCGCGCCACCCTTGCCCCCCTGATGC  >  2:901478/1‑133
                                             GCCGGCCAGGCCCACGCCCGCTACGGCGCACGACCC                                                                                                   >  1:601889/1‑36
                                             GCCGGCCAGGCCCACGCCCGCTACGGCGCACGACCC                                                                                                   >  1:878104/1‑36
                                             GCCGGCCAGGCCCACGCCCGCTACGGCGCACGACCC                                                                                                   <  2:601889/36‑1
                                             GCCGGCCAGGCCCACGCCCGCTACGGCGCACGACCC                                                                                                   <  2:878104/36‑1
                                                                                                                                                                                  
GACACACGTGGTGGTCACCGGCGGCAACGTGGCCGGCGCCACGGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AM260480/2887732‑2887778
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgcCGGCCAGGCCCACGCCCGCTACGGCGCACGACCCGGCACGCTCTACCTGATCCGCCCCGATGGCTATGTCCTGGCACGCTGGCGCGAACCGTCATGGGACGCCGTGCGCGCCACCCTTGCCCCCCTGATGC  >  AM260480/2887792‑2887922

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ 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.