New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3823258 =246 (3.720)3 (0.050) 3/264 NT 1.2% coding (859/2109 nt) spoT bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
?NC_000913 3823343 = 238 (3.680)coding (944/2109 nt) spoT bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Rejected: Frequency below/above cutoff threshold.

CGGCAACACCCATCTCCGCCATCTGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTGGAATGGCGATATAGTCTTTCACGCGGCCCGGACG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3823408‑3823258
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑acgGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTACCGCACAAATCCGCGCCCAGCGCTGGATGCAAAGCCTGCTG  >  NC_000913/3823343‑3823488
                                                                                                                                                                                                                                                                                                         
CGGCAACACCCATCTCCGCCATCTGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTGGAATGGCGATATAGTCTTTCACGCGGCCCGGA                                                                                                                                                      <  2:387670/149‑1
CGGCAACACCCATCTCCGCCATCTGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTGGAATGGCGATATAGTCTTTCACGCGGCCCGGA                                                                                                                                                      <  2:463193/149‑1
CGGCAACACCCATCTCCGCCATCTGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTGGAATGGCGATATAGTCTTTCACGCGGCCCGGA                                                                                                                                                      <  2:491628/149‑1
 GGCAACACCCATCTCCGCCATCTGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTGGAATGGCGATATAGTCTTTCACGCGGCCCGGAC                                                                                                                                                     <  4:307198/149‑1
                                                           GGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTGGAATGGCGATATAGTCTTTCACGCGGCCCGGACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGG                                                                                           >  3:315670/1‑149
                                                                 CCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTGGAATGGCGATATAGTCTTTCACGCGGCCCGGACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTTG                                                                                     >  3:538410/1‑149
                                                                                                            TCGCTTTTGGAATGGCGATATAGTCTTTCACGCGGCCCGGACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTAC                                          <  4:315670/149‑1
                                                                                                                                                    ACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTACCGCACAAATCCGCGCCCAGCGCTGGATGCAAAGCCTGCTG  >  3:171113/1‑149
                                                                                                                                                                                                                                                                                                         
CGGCAACACCCATCTCCGCCATCTGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTGGAATGGCGATATAGTCTTTCACGCGGCCCGGACG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3823408‑3823258
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑acgGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTACCGCACAAATCCGCGCCCAGCGCTGGATGCAAAGCCTGCTG  >  NC_000913/3823343‑3823488

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