New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1304426 =34 (1.350)3 (0.140) 3/116 2.1 9.6% coding (593/909 nt) oppC oligopeptide ABC transporter permease
?NC_000913 1304482 = 27 (1.220)coding (649/909 nt) oppC oligopeptide ABC transporter permease
Rejected: Frequency below/above cutoff threshold.

CCACACCGAGTACATTCGGCACAATGTGGCGAATAACAATGCCCGATGTCGATACACCGCCAACTTGTGCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1304497‑1304426
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttgtgccgAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATC  >  NC_000913/1304482‑1304546
                                                                                                                                         
CCACACCGAGTACATTCGGCACAATGTGGCGAATAACAATGCCCGATGTCGATACACCGCCAACTTGT                                                                       <  2:308602/68‑1
    ACCGAGTACATTCGGCACAATGTGGCGAATAACAATGCCCGATGTCGATACACCGCCAACTTGTGCCG                                                                   >  2:863575/1‑68
         tTACATTCGGCACAATGTGGCGAATAACAATGCCCGATGTCTATACACCGCCAACTTGTGCCGAATGT                                                              <  2:914166/67‑1
                                     AATGCCCGATGTCGATACACAGCCAACTTGTGCCGAATGTACAAGGTGTGGTGGTAGTCTACGCATCA                                  >  2:897729/1‑68
                                        GCCCGATGTCGATACACCGCCAACTTGTGCAGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTT                               >  2:764384/1‑68
                                                                TTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCGTCACTTCTGGTGCCCAGCAAGATCCTCTTT       >  1:649922/1‑68
                                                                TTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTT       >  1:136689/1‑68
                                                                TTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTT       >  1:24746/1‑68
                                                                TTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTT       >  1:316228/1‑68
                                                                     CCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATC  >  2:376608/1‑68
                                                                                                                                         
CCACACCGAGTACATTCGGCACAATGTGGCGAATAACAATGCCCGATGTCGATACACCGCCAACTTGTGCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1304497‑1304426
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttgtgccgAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATC  >  NC_000913/1304482‑1304546

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ 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.