New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 3554817 =49 (0.790)3 (0.050) 3/254 NT 6.4% coding (433/1203 nt) PP_3139 group 1 family glycosyl transferase
?NC_002947 3554849 = 41 (0.690)coding (465/1203 nt) PP_3139 group 1 family glycosyl transferase

CTTCAGGGCCGTGCACGGTAAAGCTGTATGCCGGCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/3554853‑3554817
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggcccTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTGGGTGGCGCCCGAGCAGTGGGGCAAGGTGAAGGTGGTGCATTG  >  NC_002947/3554849‑3554995
                                                                                                                                                                                        
CTTCAGGGCCGTGCACGGTAAAGCTGTATGCCGGCCCTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCGAGTTGTACCGTTGGGTGGCG                                     >  7:318910/1‑149
     GGGCCGTGCACGGTAAAGCTGTATGCCGGCCCTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTC                                                                       <  3:178377/110‑1
     GGGCCGTGCACGGTAAAGCTGTATGCCGGCCCTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTC                                                                       >  4:178377/1‑110
                                GGCCCTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGAGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTGGGTTGCGCGCGAGCAGTGGGGCAAGGTGAAGGTCGAGCA     >  3:139098/1‑149
                                  CCCTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTGGGTGGCGCCCGAGCAGTGGGGCAAGGTGAAGGTGGTGCATT   <  2:339643/149‑1
                                   CCTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTGGGTGGCGCCCGAGCAGTGGGGCAAGGTGAAGGTGGTGCATTG  <  6:15690/149‑1
                                                                                                                                                                                        
CTTCAGGGCCGTGCACGGTAAAGCTGTATGCCGGCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/3554853‑3554817
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggcccTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTGGGTGGCGCCCGAGCAGTGGGGCAAGGTGAAGGTGGTGCATTG  >  NC_002947/3554849‑3554995

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