New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3381839 =28 (0.810)5 (0.150) 5/278 2.4 16.1% coding (1097/1368 nt) degQ serine endoprotease, periplasmic
?NC_000913 3381869 = 26 (0.800)coding (1127/1368 nt) degQ serine endoprotease, periplasmic
Rejected: Frequency below cutoff threshold.

GTTGACGCCAATGATCACATCGTCTTTTTGCAAGCCAGCCTGAGCAGCTGGGCTTCCTTTGACAACTTCATCGATTTTAATACCTTTGCCGCCATCTTTTAGCTGACCATCGCTCAACGTTGCACCTTCCAGCGCTGGCGTGATCATTTCAGCGCTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3381996‑3381839
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cagcgctggAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAA  >  NC_000913/3381869‑3382010
                                                                                                                                                                                                                                                                                                            
GTTGACGCCAATGATCACATCGTCTTTTTGCAAGCCAGCCTGAGCAGCTGGGCTTCCTTTGACAACTTCATCGATTTTAATACCTTTGGCGCCATCTTTTAGCTGACCATCGCTCAACGTTGCACCTTCCAGCGCTGGCGTGATCATTTC                                                                                                                                                        <  1:82843/150‑1
GTTGACGCCAATGATCACATCGTCTTTTTGCAAGCCAGCCTGAGCAGCTGGGCTTCCTTTGACAACTTCATCGATTTTAATACCTTTGCCGCCATCTTTTAGCTGACCATCGCTCAACGTTGCACCTTCCAGCGCTGGCGTGATCATTTC                                                                                                                                                        <  1:380617/150‑1
    ACGCCAATGATCACATCGTCTTTTTGCAAGCCAGCCTGAGCAGCTGGGCTTCCTTTGACAACTTCATCGATTTTAATACCTTTGCCGCCATCTTTTAGCTGACCATCGCTCAACGTTGCACCTTCCAGCGCTGGCGTGATCATTTCAGCGC                                                                                                                                                   <  1:451391/151‑1
     CGCCAATGATCACATCGTCTTTTTGCAAGCCAGCCTGAGCAGCTGGGCTTCCTTTGACAACTTCATCGATTTTAATACCTTTGCCGCCATCTTTTAGCTGACCATCGCTCAACGTTGCACCTTCCAGCGCTGGCGTGATCATTTCAGCG                                                                                                                                                    >  2:351374/1‑149
                   TCGTCTTTTTGCAAGCCAGCCTGAGCAGCTGGGCTTCCTTTGACAACTTCATCGATTTTAATACCTTTGCCGCCATCTTTTAGCTGACCATCGCTCAACGTTGCACCTTCCAGCGCTGGCGTGAGCATTTCAGCGCTGGAAGGTGCAACGT                                                                                                                                    <  1:437327/151‑1
                                         GAGCAGCTGGGCTTCCTTTGACAACTTCATCGATTTTAATACCTTTGCCGCCATCTTTTAGCTGACAATCGCTCAACGTTGCACCTTCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTAAGCTAAAAGA                                                                                                              >  1:516671/1‑151
                                                                         ATTTTAATACCTTTGCCGCCATCTTTTAGCTGACCATCGCTCAACGTTGCACCTTCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTG                                                                              >  2:300067/1‑151
                                                                                                             TCGCTCAACGTTGCACCTTCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAA                                          >  1:314337/1‑151
                                                                                                               GCTCAACGTTGCACCTTCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAA                                        >  2:451218/1‑151
                                                                                                                                                       GCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAA  >  1:136542/1‑149
                                                                                                                                                                                                                                                                                                            
GTTGACGCCAATGATCACATCGTCTTTTTGCAAGCCAGCCTGAGCAGCTGGGCTTCCTTTGACAACTTCATCGATTTTAATACCTTTGCCGCCATCTTTTAGCTGACCATCGCTCAACGTTGCACCTTCCAGCGCTGGCGTGATCATTTCAGCGCTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3381996‑3381839
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cagcgctggAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAA  >  NC_000913/3381869‑3382010

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