New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3334936 =137 (1.510)4 (0.050) 4/250 NT 2.9% coding (28/279 nt) sfsB malPQ operon transcriptional activator
?NC_000913 3334956 = 133 (1.520)coding (48/279 nt) sfsB malPQ operon transcriptional activator
Rejected: Frequency below/above cutoff threshold.

GCGTTGAGGAACTCAAACCATTTCTGCGAGATTCCGCCGCCATTGACGTTCCCTTCTTGCGCAAACCCGCAATGATGTCAGCGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3335020‑3334936
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgggTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGTTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGACCCCTGGGTTATCTGG  >  NC_000913/3334956‑3335103
                                                                                                                                                                                                                                         
GCGTTGAGGAACTCAAACCATTTCTGCGAGATTCCGCCGCCATTGACGTTCCCTTCTTGCGCAAACCCGCAATGATGTCAGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCG                                                                                                               >  5:115280/1‑124
GCGTTGAGGAACTCAAACCATTTCTGCGAGATTCCGCCGCCATTGACGTTCCCTTCTTGCGCAAACCCGCAATGATGTCAGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCG                                                                                                               <  6:115280/124‑1
                                                            GCAAACCCGCAATGATGTCAGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTT                                                                                                    <  7:100745/75‑1
                                                            GCAAACCCGCAATGATGTCAGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTT                                                                                                    >  8:100745/1‑75
                                                                                    GTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGTTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGACCCCTGGGTTATCTGG  <  3:293431/149‑1
                                                                                                                                                                                                                                         
GCGTTGAGGAACTCAAACCATTTCTGCGAGATTCCGCCGCCATTGACGTTCCCTTCTTGCGCAAACCCGCAATGATGTCAGCGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3335020‑3334936
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgggTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAATGGTTTGAGTTCCTCAACGCTGGCGAATGCATTATCGCGCCCATGGCCGAAAGGAGAGATGATTATTGCGAAAGCCCTGGGAACTGACCCCTGGGTTATCTGG  >  NC_000913/3334956‑3335103

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.