New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 407278478 (1.140)3 (0.040) 3/270 NT 4.0% coding (788/897 nt) yihU 3‑sulpholactaldehyde (SLA) reductase, NADH‑dependent; gamma‑hydroxybutyrate dehydrogenase, NADH‑dependent
?NC_000913 = 4072835 66 (0.990)coding (737/897 nt) yihU 3‑sulpholactaldehyde (SLA) reductase, NADH‑dependent; gamma‑hydroxybutyrate dehydrogenase, NADH‑dependent
Rejected: Frequency below/above cutoff threshold.

TGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4072636‑4072784
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cccTTGATGTCGCCAACCAGCTGCATGTGCCAATGCCGCTGGGGGCCGCCTCACGGGAGG  <  NC_000913/4072835‑4072779
                                                                                                                                                                                                              
TGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCC                                                           >  1:997359/1‑149
                                         ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGA‑GGCGGCCCCCAGCTGCATGTGCCAATGCCGCTGGGGG                 <  2:47350/149‑1
                                                        CAGGCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCTTGATGTCGCCAACCAGCTGCATGTGCCAATGCCGCTGGGGGCCGCCTCACGGGAG   <  2:323696/149‑1
                                                         cGTCATCCCCGCACCGCCACGGTCCTGTGCCAGAACGGCGGACCAGTCCAGCCGACCGCGACCCGCGGCGCGCGCCTGGCTATAAACCTCCCTTGATGTCGCCAACCAGCTGCATGTGCCAATGCCGCGGGGGGCCGCCTCACGGGAGG  <  2:789872/148‑1
                                                                                                                                                                                                              
TGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4072636‑4072784
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cccTTGATGTCGCCAACCAGCTGCATGTGCCAATGCCGCTGGGGGCCGCCTCACGGGAGG  <  NC_000913/4072835‑4072779

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.