New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 5625659 =61 (1.020)4 (0.070) 4/232 NT 6.6% coding (153/1131 nt) PP_4943 glycosyl transferase
?NC_002947 5625684 = 55 (0.950)coding (128/1131 nt) PP_4943 glycosyl transferase

CCTCACCGGCGCCGCCGATGCGCAGGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/5625686‑5625659
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggtgAGGAATACCGTGGTCACCTGGTACCCACGCCCCTGAAACAGGCT  >  NC_002947/5625684‑5625727
                                                                        
CCTCACCGGCGCCGCCGATGCGCAGGTGAGGAATACCGTGGTCACCTGGTACCCACGCCCCTGAAA        >  1:1439516/1‑66
CCTCACCGGCGCCGCCGATGCGCAGGTGAGGAATACCGTGGTCACCTGGTACCCACGCCCCTGAAA        <  2:1439516/66‑1
  TCACCGGCGCCGCCGATGCGCAGGTGAGGAATACCGTGGTCACCTGGTACCCACGCCCCTGAAACAGGCT  >  1:55152/1‑70
  TCACCGGCGCCGCCGATGCGCAGGTGAGGAATACCGTGGTCACCTGGTACCCACGCCCCTGAAACAGGCT  <  2:55152/70‑1
                                                                        
CCTCACCGGCGCCGCCGATGCGCAGGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/5625686‑5625659
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggtgAGGAATACCGTGGTCACCTGGTACCCACGCCCCTGAAACAGGCT  >  NC_002947/5625684‑5625727

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