New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 460704794 (1.090)7 (0.090) 6/168 NT 6.9% coding (654/1032 nt) rsmC 16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
?NC_000913 = 4607088 99 (1.210)coding (613/1032 nt) rsmC 16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
Rejected: Frequency below cutoff threshold.

GCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4606945‑4607047
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcggc‑TGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGC  <  NC_000913/4607088‑4606993
                                                                                                                                                                                                        
GCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTC                                                                                                      >  1:395870/1‑100
 CTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCG                                                                                                     >  2:2018577/1‑100
  TGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGG                                                                                                    >  1:840478/1‑100
  TGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGG                                                                                                    >  1:918284/1‑100
   GGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC                                                                                                   <  1:452969/100‑1
   GGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC                                                                                                   >  1:758515/1‑100
   GGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC                                                                                                   >  2:1624479/1‑100
                          CGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑TGTGGCGCGGGCGTGCTTTCAGT                                                                           <  1:1798060/100‑1
                                             GCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTC                                                         <  1:2041988/100‑1
                                               GGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑TGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCC                                                      <  2:582282/100‑1
                                               GGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑TGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCC                                                      <  2:264730/100‑1
                                                  TGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑TGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAA                                                   <  1:196871/100‑1
                                                     CTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑TGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAAT                                                <  1:820983/100‑1
                                                        CTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑TGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCG                                             <  1:1677925/100‑1
                                                                                                  TCGGC‑TGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAG   <  1:884743/100‑1
                                                                                                  TCGGC‑TGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGACGGTAGAAGCCA    >  1:8957/1‑99
                                                                                                   CGGC‑TGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGC  >  2:2201590/1‑100
                                                                                                   CGGC‑TGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCA    >  2:413700/1‑98
                                                                                                                                                                                                        
GCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4606945‑4607047
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcggc‑TGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGC  <  NC_000913/4607088‑4606993

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