New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3873646 =101 (1.150)4 (0.050) 3/244 NT 4.2% coding (829/879 nt) dgoK 2‑oxo‑3‑deoxygalactonate kinase
?NC_000913 3873729 = 88 (1.060)coding (746/879 nt) dgoK 2‑oxo‑3‑deoxygalactonate kinase
Rejected: Frequency below/above cutoff threshold.

TCTGGTTTGTTGATTGGCGCAGAGGTCGCCAGTATGCGCGACTATGTGGCCCATCAACACGCCATCACCCTTGTCGCCGGAACATCGCTGACCGCGCGCTACCAGCAAGCCTTTCAGGCGATGGGTTGCGACGTGACGGCGGTGGCGGGCGACA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3873799‑3873646
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgacaAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGC  >  NC_000913/3873729‑3873850
                                                                                                                                                                                                                                                                                    
TCTGGTTTGTTGATTGGCGCAGAGGTCGCCAGTATGCGCGACTATGTGGCCCCTCAACACGCCATCACCCTTGTCGCCGGAACATCGCTGACCGCGCGCTACCAGCAAGCCTTTCAGGCGATGGGTTGCGACGTGACGGCGGTGGCGGG                                                                                                                                 >  4:105236/1‑149
TCTGGTTTGTTGATTGGCGCAGAGGTCGCCAGTATGCGCGACTATGTGGCCCATCAACACGCCATCACCCTTGTCGCCGGAACATCGCTGACCGCGCGCTACCAGCAAGCCTTTCAGGCGATGGGTTGCGACGTGACGGCGGTGGCGGG                                                                                                                                 >  4:55468/1‑149
TCTGGTTTGTTGATTGGCGCAGAGGTCGCCAGTATGCGCGACTATGTGGCCCATCAACACGCCATCACCCTTGTCGCCGGAACATCGCTGACCGCGCGCTACCAGCAAGCCTTTCAGGCGATGGGTTGCGACGTGACGGCGGTGGCGGG                                                                                                                                 >  7:302922/1‑149
  TGGTTTGTTGATTGGCGCAGAGGTCGCCAGTATGCGCGACTATGTGGCCCATCAACACGCCATCACCCTTGTCGCCGGAACATCGCTGACCGCGCGCTACCAGCAAGCCTTTCAGGCGATGGGTTGCGACGTGACGGCGGTGGCGGGCG                                                                                                                               >  1:336060/1‑149
    GTTTGTTGATTGGCGCAGAGGTCGCCAGTATGCGCGACTATGTGGCCCATCAACACGCCATCACCCTTGTCGCCGGAACATCGCTGACCGCGCGCTACCAGCAAGCCTTTCAGGCGATGGGTTGCGACGTGACGGCGGTGGCGGGCGA                                                                                                                              <  7:317078/148‑1
    GTTTGTTGATTGGCGCAGAGGTCGCCAGTATGCGCGACTATGTGGCCCATCAACACGCCATCACCCTTGTCGCCGGAACATCGCTGACCGCGCGCTACCAGCAAGCCTTTCAGGCGATGGGTTGCGACGTGACGGCGGTGGCGGGCGA                                                                                                                              >  7:85575/1‑148
                                                                    CCTTGTCGCCGGAACATCGCTGACCGCGCGCTACCAGCAAGCCTTTCAGGCGATGGGTTGCGACGTGACGGCGGTGGCGGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAA                                                              >  2:360152/1‑148
                                                                    CCTTGTCGCCGGAACATCGCTGACCGCGCGCTACCAGCAAGCCTTTCAGGCGATGGGTTGCGACGTGACGGCGGTGGCGGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAA                                                              <  5:73891/148‑1
                                                                    CCTTGTCGCCGGAACATCGCTGACCGCGCGCTACCAGCAAGCCTTTCAGGCGATGGGTTGCGACGTGACGGCGGTGGCGGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAA                                                              >  6:73891/1‑148
                                                                                                                   AGGCGATGGGTTGCGACGTGACGGCGGTGGCGGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCA               <  1:360152/148‑1
                                                                                                                                                    GCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGC  <  3:177782/128‑1
                                                                                                                                                    GCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGC  >  4:177782/1‑128
                                                                                                                                                                                                                                                                                    
TCTGGTTTGTTGATTGGCGCAGAGGTCGCCAGTATGCGCGACTATGTGGCCCATCAACACGCCATCACCCTTGTCGCCGGAACATCGCTGACCGCGCGCTACCAGCAAGCCTTTCAGGCGATGGGTTGCGACGTGACGGCGGTGGCGGGCGACA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3873799‑3873646
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgacaAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGC  >  NC_000913/3873729‑3873850

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