New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 5456645 =39 (0.900)2 (0.050) 4/244 NT 5.0% coding (329/2607 nt) leuS leucine‑‑tRNA ligase
?NC_002947 5456664 = 38 (0.910)coding (348/2607 nt) leuS leucine‑‑tRNA ligase

GCCGGGTCCCAGTTCACGGTACCGTTCTTGCGGTAGATGATGCCTTTTTCGAACAGGCGGGTGAACAGCCACTGCTCCCAACGGTAATAGTCGGGTTTGCAGGTGGTGACTTCACGTGCCCAGTCGATGGCCAGGCCCAGGCTCTTGAGCTGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/5456798‑5456645
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgggCCTGGCCATCGACTGG  >  NC_002947/5456664‑5456679
                                                                                                                                                                           
GCCGGGTCCCAGTTCACGGTACCGTTCTTGCGGTAGATGATGCCTTTTTCGAACAGGCGGGTGAACAGCCACTGCTCCCAACGGTAATAGTCGGGTTTGCAGGTGGTGACTTCACGTGCCCAGTCGATGGCCAGGCCCAGGCTCTTGAGC                      >  2:270314/1‑150
                                                                ACAGCCACTGCTCCCAACGGTAATAGTCGGGTTTGCAGGTGGTGACTTCACGTGCCCAGTCGATGGCCAGGCCCAGGCTCTTGAGCTGG                   <  1:1032266/89‑1
                                                                ACAGCCACTGCTCCCAACGGTAATAGTCGGGTTTGCAGGGGGTGACTTCACGTGCCCAGTCGATGGCCAGGCCCAGGCTCTTGAGCTGG                   >  2:1032266/1‑89
                                                                        TGCTCCCAACGGTAATAGTCGGGTTTGCAGGTGGTGACTTCACGTGCCCAGTCGATGGCCAGGCCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGG  >  1:953797/1‑98
                                                                        TGCTCCCAACGGTAATAGTCGGGTTTGCAGGTGGTGACTTCACGTGCCCAGTCGATGGCCAGGCCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGG  <  2:953797/98‑1
                                                                                  GGTAATAGTCGGGTTTGCAGGTGGTGACTTCACGTGCCCAGTCGATGGCCAGGCCCAGCTCAAGAGCCTGGGCCTGG             >  1:1043954/1‑77
                                                                                  GGTAATAGTCGGGTTTGCAGGTGGTGACTTCACGTGCCCAGTCGATGGCCAGGCCCAGCTCAAGAGCCTGGGCCTGG             <  2:1043954/77‑1
                                                                                                                                                                           
GCCGGGTCCCAGTTCACGGTACCGTTCTTGCGGTAGATGATGCCTTTTTCGAACAGGCGGGTGAACAGCCACTGCTCCCAACGGTAATAGTCGGGTTTGCAGGTGGTGACTTCACGTGCCCAGTCGATGGCCAGGCCCAGGCTCTTGAGCTGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/5456798‑5456645
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgggCCTGGCCATCGACTGG  >  NC_002947/5456664‑5456679

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