New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 561977475 (1.440)3 (0.060) 3/252 NT 4.0% coding (727/1140 nt) PP_4938 glycosyl transferase
?NC_002947 = 5619802 73 (1.480)coding (699/1140 nt) PP_4938 glycosyl transferase
Rejected: Frequency below/above cutoff threshold.

GGAACGCGAAGGTACCAATAGCCAGTCGAATGCTCGATACAGCTGCTGCAGGTCTTCGCGGTGACCGCGCAGATGAACGCGCTCGGCCAGGCCCAACGCGTCGATACGCTGCTGCAGCACGGCATGCTGCGGCCCCTCGCCAATGATCGCCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5619622‑5619774
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcgccagGGCCAGTGCACGGCAACCGGGCTTGCAGCTGGCGATCATTGGCGAGGGGCCGCAGCATGCCGTGCTGCAGCAGCGTATCGACGCGTTGGGCCTGGCCGAGCGCGTTCATCTGCGCGGTCACCGCGAAGACCTGCAGCAGCTGTATCGA  <  NC_002947/5619802‑5619655
                                                                                                                                                                                                                                                                                                             
GGAACGCGAAGGTACCAATAGCCAGTCGAATGCTCGATACAGCTGCTGCAGGTCTTCGCGGTGACCGCGCAGATGAACGCGCTCGGCCAGGCCCAACGCGTCGATACGCTGCTGCAGCACGGCATGCTGCGGCCCCTCGCCAATGATC                                                                                                                                                           <  2:274882/148‑1
 GAACGCGAAGGTACCAATAGCCAGTCGAATGCTCGATACAGCTGCTGCAGGTCTTCGCGGTGACCGCGCAGATGAACGCGCTCGGCCAGGCCCAACGCGTCGATACGCTGCTGCAGCACGGCATGCTGCGGCCCCTCGCCAATGATCGC                                                                                                                                                         >  2:121733/1‑149
    CGCGAAGGTACCAATAGCCAGTCGAATGCTCGATACAGCTGCTGCAGGTCTTCGCGGTGACCGCGCAGATGAACGCGCTCGGCCAGGCCCAACGCGTCGATACGCTGCTGCAGCACGGCATGCTGCGGCCCCTCGCCAATGATCGCCAG                                                                                                                                                      >  1:291461/1‑149
                                    ATACAGCTGCTGCAGGTCTTCGCGGTGACCGCGCAGATGAACGCGCTCGGCCAGGCCCAACGCGTCGATACGCTGCTGCAGCACGGCATGCTGCGGCCCCTCGCCAATGATCGCCAGGGCCAGTGCACGGCAACCGGGCTTGCAGCTGG                                                                                                                      <  8:85708/149‑1
                                             CTGCAGGTCTTCGCGGTGACCGCGCAGATGAACGCGCTCGGCCAGGCCCAACGCGTCGATACGCTGCTGCAGCACGGCATGCTGCGGCCCCTCGCCAATGATCGCCAGGGCCAGTGCACGGCAACCGGGCTTGCAGCTGGCGATCATT                                                                                                              <  8:262242/148‑1
                                                CAGGTCTTCGCGGTGACCGCGCAGATGAACGCGCTCGGCCAGGCCCAACGCGTCGATACGCTGCTGCAGCACGGCATGCTGCGGCCCCTCGCCAATGATCGCCAGGGCCAGTGCACGGCAACCGGGCTTGCAGCTGGCGATCATTGG                                                                                                            <  4:165453/147‑1
                                                                                                                                                        GGGCCAGTGCACGGCAACCGGGCTTGCAGCTGGCGATCATTGGCGAGGGGCCGCAGCATGCCGTGCTGCAGCAGCGTATCGACGCGTTGGGCCTGGCCGAGCGCGTTCATCTGCGCGGTCACCGCGAAGACCTGCAGCAGCTGTATCGA  <  7:161876/149‑1
                                                                                                                                                                                                                                                                                                             
GGAACGCGAAGGTACCAATAGCCAGTCGAATGCTCGATACAGCTGCTGCAGGTCTTCGCGGTGACCGCGCAGATGAACGCGCTCGGCCAGGCCCAACGCGTCGATACGCTGCTGCAGCACGGCATGCTGCGGCCCCTCGCCAATGATCGCCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5619622‑5619774
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcgccagGGCCAGTGCACGGCAACCGGGCTTGCAGCTGGCGATCATTGGCGAGGGGCCGCAGCATGCCGTGCTGCAGCAGCGTATCGACGCGTTGGGCCTGGCCGAGCGCGTTCATCTGCGCGGTCACCGCGAAGACCTGCAGCAGCTGTATCGA  <  NC_002947/5619802‑5619655

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