New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3235543 =28 (0.900)4 (0.140) 4/252 3.5 13.0% coding (333/1137 nt) rlmG 23S rRNA m(2)G1835 methyltransferase, SAM‑dependent
?NC_000913 3235583 = 27 (0.920)coding (293/1137 nt) rlmG 23S rRNA m(2)G1835 methyltransferase, SAM‑dependent
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CAGCACCGCCGACTACCCGCAACAGCCGGGTGTGGTGCTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3235582‑3235543
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggtgctgTCGAGAAACTTCACGCTCGACTCATCAATCCCGTTGAGGCGTAAATTCTCGCGCGTCGCCAGTTCACTGATGTATGAGTCGCCAATGCTGTACGGCTTATGTTCCGCCAGTGCGCAGCTTAACGCACCAAAGGCATCATTCAGGATC  >  NC_000913/3235583‑3235729
                                                                                                                                                                                           
CAGCACCGCCGACTACCCGCAACAGCCGGGTGTGGTGCTGTCGAGAAACTTCACGCTCGACTCATCAATCCCGTTGAG                                                                                                               >  7:7469/1‑78
CAGCACCGCCGACTACCCGCAACAGCCGGGTGTGGTGCTGTCGAGAAACTTCACGCTCGACTCATCAATCCCGTTGAG                                                                                                               <  8:7469/78‑1
   CACCGCCGACTACCCGCAACAGCCGGGTGTGGTGCTGTCGAGAAACTTCACGCTCGACTCATCAATCCCGTTGAGGCGTAAATTCTCGCGCGTCGCCAGTTCACTGATGTATGAGTCGCCAATGCTGTACGGCTTATGTTCCGC                                          <  7:105385/144‑1
   CACCGCCGACTACCCGCAACAGCCGGGTGTGGTGCTGTCGAGAAACTTCACGCTCGACTCATCAATCCCGTTGAGGCGTAAATTCTCGCGCGTCGCCAGTTCACTGATGTATGAGTCGCCAATGCTGTACGGCTTATGTTCCGC                                          >  8:105385/1‑144
                                       GTCGAGAAACTTCACGCTCGACTCATCAATCCCGTTGAGGCGTAAATTCTCGCGCGTCGCCAGTTCACTGATGTATGAGTCGCCAATGCTGTACGGCTTATGTTCCGCCAGTGCGCAGCTTAACGCACCAAAGGCATCATTCAGGATC  >  4:8579/1‑148
                                                                                                                                                                                           
CAGCACCGCCGACTACCCGCAACAGCCGGGTGTGGTGCTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3235582‑3235543
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggtgctgTCGAGAAACTTCACGCTCGACTCATCAATCCCGTTGAGGCGTAAATTCTCGCGCGTCGCCAGTTCACTGATGTATGAGTCGCCAATGCTGTACGGCTTATGTTCCGCCAGTGCGCAGCTTAACGCACCAAAGGCATCATTCAGGATC  >  NC_000913/3235583‑3235729

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