New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4631582 =150 (1.140)4 (0.030) 4/168 NT 2.8% coding (850/1938 nt) slt lytic murein transglycosylase, soluble
?NC_000913 4631612 = 133 (1.070)coding (880/1938 nt) slt lytic murein transglycosylase, soluble
Rejected: Frequency below cutoff threshold.

GGCGTTCAATAAGCGAAGTAGATTGCGAGCGCATAATGGCGTCATCGCGCCATTTCGCCTGCTCGTCGGTGACATCGTTGCCCATCAAACGCCAGGCGACGAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4631684‑4631582
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑acgatGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCAC  >  NC_000913/4631612‑4631706
                                                                                                                                                                                                      
GGCGTTCAATAAGCGAAGTAGATTGCGAGCGCATAATGGCGTCATCGCGCCATTTCGCCTGCTCGTCGGTGACATCGTTGCCCATCAAACGCCAGGCGAC                                                                                                    <  1:561996/100‑1
 ttccgatcTAAGCGAAGTAGATTGCGAGCGCATAATGGCGTCATCGCGCCATTTCGCCTGCTCGTCGGTGACATCGTTGCCCATCAAACGCCAGGCGACG                                                                                                   <  1:3252628/92‑1
 GCGTTCAATAAGCGAAGTAGATTGCGAGCGCATAATGGCGTCATCGCGCCATTTCGCCTGCTCGTCGGTGACATCGTTGCCCATCAAACGCCAGGCGACG                                                                                                   <  1:1892712/100‑1
   GTTCAATAAGCGAAGTAGATTGCGAGCGCATAATGGCGTCATCGCGCCATTTCGCCTGCTCGTCGGTGACATCGTTGCCCATCAAACGCCAGGCGACGAT                                                                                                 >  1:2417079/1‑100
   GTTCAATAAGCGAAGTAGATTGCGAGCGCATAATGGCGTCATCGCGCCATTTCGCCTGCTCGTCGGTGACATCGTTGCCCATCAAACGCCAGGCGACGAT                                                                                                 >  2:1667226/1‑100
                                             CGCGCCATTTCGCCTGCTCGTCGGTGACATCGTTGCCCATCAAACGCCAGGCGACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTAT                                                        >  1:972248/1‑99
                                                  CATTTCGCCTGCTCGTCGGTGACATCGTTGCCCATCAAACGCCAGGCGACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTC                                                  >  1:1753065/1‑100
                                                                       ACATCGTTGCCCATCAAACGCCAGGCGACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGA                             >  1:1728077/1‑100
                                                                                          GCCAGGCGACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCG          <  2:1728084/100‑1
                                                                                                  ACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCAC  <  1:3501402/100‑1
                                                                                                                                                                                                      
GGCGTTCAATAAGCGAAGTAGATTGCGAGCGCATAATGGCGTCATCGCGCCATTTCGCCTGCTCGTCGGTGACATCGTTGCCCATCAAACGCCAGGCGACGAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4631684‑4631582
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑acgatGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCAC  >  NC_000913/4631612‑4631706

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