New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3873274120 (1.180)3 (0.030) 3/270 NT 2.5% coding (339/618 nt) dgoA 2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase
?NC_000913 = 3873314 120 (1.210)coding (299/618 nt) dgoA 2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase
Rejected: Frequency below/above cutoff threshold.

TGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3873122‑3873274
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggctACGGCATGACCGTCTGCCCCGGCTGCGCGACGGCGACCGAAGCCTTTACCGCGCTCGAAGCGGGCGCGCAGGCGCTGAAAATATTTCCGTCATCGGCTTTTGGTCCGCAATACATCAAAGCGTTAAAAGCGGTATTGCCATCGGACATC  <  NC_000913/3873314‑3873166
                                                                                                                                                                                                                                                                                                              
TGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAG                                                                                                                                                          >  1:1147005/1‑150
TGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAG                                                                                                                                                          >  1:1147021/1‑150
 GCGCCAGGTTTTCTGGCGTCACGCCGCCATCGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTTTTGCGGACCAAAAGCCGATGACTGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGG                                                                                                                                                         <  2:901034/150‑1
 GCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTCCCGCTTTGATGTAGTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGG                                                                                                                                                         <  1:738336/150‑1
 GCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGG                                                                                                                                                         <  1:738341/150‑1
 GCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGG                                                                                                                                                         <  2:901040/150‑1
                                         TGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTACGGCATGACCGTCTGCCCCGGCTGCGCGACGGCGACC                                                                                                                 >  1:862380/1‑150
                                                   GATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTACGGCATGACCGTCTGCCCCGGCTGCGCGACGGCGACCGAAGCCTTTA                                                                                                       <  2:1703811/150‑1
                                                    ATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTACGGCATGACCGTCTGCCCCGGCTGCGCGACGGCGACCGAAGCCTTTAC                                                                                                      <  2:862380/150‑1
                                                                                                                                                     GGCTACGGCATGACCGTCTGCCCCGGCTGCGCGACGGCGACCGAAGCCTTTACCGCGCTCGAAGCGGGCGCGCAGGCGCTGAAAATATTTCCGTCATCGGCTTTTGGTCCGCAATACATCAAAGCGTTAAAAGCGGTATTGCCATCGGAC     <  1:1576987/150‑1
                                                                                                                                                      GCTACGGCATGACCGTCTGCCCCGGCTGCGCGACGGCGACCGAAGCCTTTACCGCGCTCGAAGCGGGCGCGCAGGCGCTGAAAATATTTCCGTCATCGGCTTTTGGTCCGCAATACATCAAAGCGTTAAAAGCGGTATTGCCATCGGACA    >  2:90690/1‑150
                                                                                                                                                        TACGGCATGACCGTCTGCCCCGGCTGCGCGACGGCGACCGAAGCCTTTACCGCGCTCGAAGCGGGCGCGCAGGCGCTGAAAATATTTCCGTCATCGGCTTTTGGTCCGCAATACATCAAAGCGTTAAAAGCGGTATTGCCATCGGACATC  <  1:1680425/150‑1
                                                                                                                                                        TACGGCATGACCGTCTGCCCCGGCTGCGCGACGGCGACCGAAGCCTTTACCGCGCTCGAAGCGGGCGCGCAGGCGCTGAAAATATTTCCGTCATCGGCTTTTGGTCCGCAATACATCAAAGCGTTAAAAGCGGTATTGCCATCGGACATC  >  1:863500/1‑150
                                                                                                                                                        TACGGCATGACCGTCTGCCCCGGCTGCGCGACGGCGACCGAAGCCTTTACCGCGCTCGAAGCGGGCGCGCAGGCGCTGAAAATATTTCCGTCATCGGCTTTTGGTCCGCAATACATCAAAGCGTTAAAAGCGGTATTGCCATCGGACATC  >  1:863503/1‑150
                                                                                                                                                                                                                                                                                                              
TGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3873122‑3873274
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggctACGGCATGACCGTCTGCCCCGGCTGCGCGACGGCGACCGAAGCCTTTACCGCGCTCGAAGCGGGCGCGCAGGCGCTGAAAATATTTCCGTCATCGGCTTTTGGTCCGCAATACATCAAAGCGTTAAAAGCGGTATTGCCATCGGACATC  <  NC_000913/3873314‑3873166

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