New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 23256454 (1.450)8 (0.240) 8/226 2.9 13.8% intergenic (+15/+33) yafE/mltD putative S‑adenosyl‑L‑methionine‑dependent methyltransferase/putative membrane‑bound lytic murein transglycosylase D
?NC_000913 = 232571 51 (1.510)intergenic (+22/+26) yafE/mltD putative S‑adenosyl‑L‑methionine‑dependent methyltransferase/putative membrane‑bound lytic murein transglycosylase D
Rejected: Frequency below/above cutoff threshold.

CGCCAGAAGCGTTAGTAGACGCTATTCGCATTTACCAACAGAGCGCATCGACAGAGGTGAGAACGTATTTTGCCTTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/232413‑232564
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aaaaaggcaccgATTCCCCCGGTGCCTTTTTTATTTATGCCGCTTTATGTGCATCTACCATGATGATATCACTGGTGAAAAAGCCATCATTCTGCAAGGCAAAATACGTTCTCACCTCTGTCGATGCGCTC  <  NC_000913/232571‑232453
                                                                                                                                                                                                                                                                               
CGCCAGAAGCGTTAGTAGACGCTATTCGCATTTACCAACAGAGCGCATCGACAGAGGTGAGAACGTATTTTGCCTTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGC                                                                                                                             <  8:15371/148‑1
                          CGCATTTACCAACAGAGCGCATCGACAGAGGTGAGAACGTATTTTGCCTTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCA                                                                                                                            >  1:1779/1‑123
                          CGCATTTACCAACAGAGCGCATCGACAGAGGTGAGAACGTATTTTGCCTTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCA                                                                                                                            <  2:1779/123‑1
                                                                          TTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGATTCCCCCGGTGCCTTTTTTA                                                                                                    >  3:12568/1‑99
                                                                          TTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGATTCCCCCGGTGCCTTTTTTA                                                                                                    <  4:12568/99‑1
                                                                            GCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGATTCCCCCGGTGCCTTTTTTATTTAT                                                                                               <  5:99998/102‑1
                                                                            GCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGATTCCCCCGGTGCCTTTTTTATTTAT                                                                                               >  6:99998/1‑102
                                                                                   GATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGATTCCCCCGGTGCCTTTTTTATTT                                                                                                 >  3:161006/1‑93
                                                                                   GATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGATTCCCCCGGTGCCTTTTTTATTT                                                                                                 <  4:161006/93‑1
                                                                                                               TGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTTATTTATGCCGCTTTATGTGCATCTACCATGATGATATCACTGGTGAAAAAGCCATCATTCTGCAAGGCAAAATACGTTCTCACCTCTG             >  2:165180/1‑149
                                                                                                                          CATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTTATTTATGCCGCTTTATGTGCATCTACCATGATGATATCACTGGTGAAAAAGCCATCATTCTGCAAGGCAAAATACGTTCTCACCTCTGTCGATGCGCTC  <  1:165180/149‑1
                                                                                                                                                                                                                                                                               
CGCCAGAAGCGTTAGTAGACGCTATTCGCATTTACCAACAGAGCGCATCGACAGAGGTGAGAACGTATTTTGCCTTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/232413‑232564
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aaaaaggcaccgATTCCCCCGGTGCCTTTTTTATTTATGCCGCTTTATGTGCATCTACCATGATGATATCACTGGTGAAAAAGCCATCATTCTGCAAGGCAAAATACGTTCTCACCTCTGTCGATGCGCTC  <  NC_000913/232571‑232453

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