New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 394843426 (1.030)4 (0.170) 4/238 2.4 14.0% intergenic (+10/+15) yifE/yifB UPF0438 family protein/magnesium chelatase family protein and putative transcriptional regulator
?NC_000913 = 3948439 25 (1.060)intergenic (+15/+10) yifE/yifB UPF0438 family protein/magnesium chelatase family protein and putative transcriptional regulator
Rejected: Frequency below/above cutoff threshold.

AGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3948295‑3948434
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aaaaagggcATTTCGCCCTTTTTATTAATCGTCA  <  NC_000913/3948439‑3948415
                                                                                                                                                                     
AGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAG                              >  3:5059/1‑137
AGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAG                              <  4:5059/137‑1
                                                      AAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCATTTCGCCCTTTTTAT           <  3:60841/102‑1
                                                      AAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCATTTCGCCCTTTTTAT           >  4:60841/1‑102
                                                                     CCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCATTTCGCCCTTTTTATTAATCGTCA  <  5:30694/96‑1
                                                                     CCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCATTTCGCCCTTTTTATTAATCGTCA  >  6:30694/1‑96
                                                                                                                                                                     
AGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3948295‑3948434
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aaaaagggcATTTCGCCCTTTTTATTAATCGTCA  <  NC_000913/3948439‑3948415

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