New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1881698 =11 (0.430)3 (0.130) 3/122 2.3 20.6% coding (112/924 nt) dmlR transcriptional activator of dmlA
?NC_000913 1881770 = 13 (0.540)coding (40/924 nt) dmlR transcriptional activator of dmlA

ACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACC  >  NC_000913/1881634‑1881765
                                                                |                                                                   
ctggtgaaacataaaacacgaatcatgaataatttaccgctgctgaatgatttgcgcgtctttaTGCt                                                                  <  2:380941‑M1/4‑1 (MQ=255)
ctggtgaaacataaaacacgaatcatgaataatttaccgctgctgaatgatttgcgcgtctttaTGCt                                                                  <  2:309621‑M1/4‑1 (MQ=255)
 tggtgaaacataaaacacgaatcatgaataatttaccgctgctgaatgatttgcgcgacttttTGCTg                                                                 <  1:355842‑M1/5‑1 (MQ=255)
        ggtggtGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACg                                                          <  1:335264/68‑1 (MQ=255)
              gCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTg                                                    >  2:703858/1‑68 (MQ=255)
                   gaatcatgaataatttaccgctgccgaatgatttgcgcggctttaTGCTGACGAACGCCGGTGAAACg                                               >  2:518783‑M1/46‑68 (MQ=255)
                        aacaCCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGCCGAACGCCGGTGCCACGCCCAg                                          >  2:145499/1‑68 (MQ=255)
                        aacaCCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAg                                          >  1:801041/1‑68 (MQ=255)
                          caccacGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGtt                                        >  2:903087/1‑68 (MQ=255)
                          caccacGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGtt                                        >  2:812867/1‑68 (MQ=255)
                          caccacGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGtt                                        >  1:309621/1‑68 (MQ=255)
                          caccacGTTTAGCGTTAGCTCCAGCAAGGAGATGCGCTTGCTGACGAACGCCGGTGAAGAGCCCAGta                                        >  1:620457/1‑67 (MQ=255)
                                      cGTTTGCTCCAGCAAGGCGGTGCGCTTGCTGACGACCGCCGGTGAAACGCCCAGTTCTTCCGCCAcgg                            >  2:693264/1‑68 (MQ=255)
                                      cGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGACCGCCGGTGAAACGCCCAGTTCTTCCGCCAcgg                            >  1:11012/1‑68 (MQ=255)
                                      cGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCAcgg                            >  2:256381/1‑68 (MQ=255)
                                                gatttgcgcgtctttaTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCaaaa                    >  1:722805‑M1/17‑66 (MQ=255)
                                                gatttgcgcgtctttaTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCaaaa                    <  2:722803‑M1/50‑1 (MQ=255)
                                                                TGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGccacc  >  2:883059‑M1/1‑68 (MQ=255)
                                                                |                                                                   
ACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACC  >  NC_000913/1881634‑1881765

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