New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi 2968371 =14 (0.970)3 (0.220) 3/248 1.6 18.0% coding (64/2247 nt) ptsP PEP‑protein phosphotransferase enzyme I; GAF domain containing protein
?NC_000913_3_bme_pgi 2968388 = 14 (1.020)coding (47/2247 nt) ptsP PEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Rejected: Frequency below/above cutoff threshold.

ACGCATCTGCTTATCGACGTAAAAGAGGTTAAGTCACGCCAATTATGCTCACTCGCCTGCGCGAAATAGTCGAAAAGGTAGCCAGCGCACCACGCCTGAATGAGGCGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913_3_bme_pgi/2968478‑2968371
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aggcgtGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATCACCACACTCGCGAACTCTT  >  NC_000913_3_bme_pgi/2968388‑2968532
                                                                                                                                                                                                                                                             
ccgacTCTGCTTATCGACGTAAAAGAGGTTAAGTCACGCAAATTATGCTCACTCGCCTGCGCGAAATAGTCGAAAAGGTAGCCAGCGCACCACGCCTGAATGAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTG                                                                                                         <  4:64299/145‑1
                                                             CGAAATAGTCGAAAAGGTAGCCAGCGCACCACGCCTGAATGAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCG                                                                      <  3:39160/124‑1
                                                             CGAAATAGTCGAAAAGGTAGCCAGCGCACCACGCCTGAATGAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCG                                                                      >  4:39160/1‑124
                                                                                                       GGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATCACCACACTCGCGAACTCTT  <  3:64299/150‑1
                                                                                                                                                                                                                                                             
ACGCATCTGCTTATCGACGTAAAAGAGGTTAAGTCACGCCAATTATGCTCACTCGCCTGCGCGAAATAGTCGAAAAGGTAGCCAGCGCACCACGCCTGAATGAGGCGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913_3_bme_pgi/2968478‑2968371
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aggcgtGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATCACCACACTCGCGAACTCTT  >  NC_000913_3_bme_pgi/2968388‑2968532

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