New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4528081345 (1.000)7 (0.020) 5/278 NT 1.9% intergenic (+15/+30) ryjB/sgcQ novel sRNA, function unknown/putative nucleoside triphosphatase
?NC_000913 = 4528108 370 (1.080)intergenic (+42/+3) ryjB/sgcQ novel sRNA, function unknown/putative nucleoside triphosphatase
Rejected: Frequency below/above cutoff threshold.

GCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTAGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4527952‑4528081
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TGGGGATAAGCAGAGCGCTTATCCCCGGCAAAGA  <  NC_000913/4528108‑4528075
                                                                                                                                                                    
GCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTAGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAAGCAGAGCGCTT                <  1:2163128/150‑1
GCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTAGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAAGCAGAGCGCTT                <  1:2590361/150‑1
GCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTAGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAAGCAGAGCGCTT                <  1:2590467/150‑1
     TACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTAGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAAGCAGAGCGCTTATCCC           <  2:2318132/150‑1
              AATATCTTTAATCATCCGTCTTTGACTCCATGCCGATTCGGTTTAATCTGTTAGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCTGGGTATAAGCAGAGCGCTTATCCCCGGCAAAGA  <  1:245888/150‑1
                                                                                             AGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCTGGGGATAAGCAGAGCGCTTATCC            >  1:4341141/1‑61
                                                                                             AGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCTGGGGATAAGCAGAGCGCTTATCC            <  2:4341141/61‑1
                                                                                                                                                                    
GCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTAGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4527952‑4528081
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TGGGGATAAGCAGAGCGCTTATCCCCGGCAAAGA  <  NC_000913/4528108‑4528075

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