New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4607047127 (1.170)9 (0.090) 8/168 NT 7.4% coding (654/1032 nt) rsmC 16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
?NC_000913 = 4607088 106 (1.030)coding (613/1032 nt) rsmC 16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
Rejected: Frequency below cutoff threshold.

GCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4606945‑4607047
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcggc‑TGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGCC  <  NC_000913/4607088‑4606992
                                                                                                                                                                                                         
GCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTC                                                                                                       >  1:379767/1‑100
 CTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCG                                                                                                      <  1:1505271/100‑1
 CTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCG                                                                                                      >  2:1572576/1‑100
   GGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC                                                                                                    <  2:755901/100‑1
   GGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC                                                                                                    >  1:2144279/1‑100
         GACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCTGTGGC                                                                                             >  2:1172674/1‑100
                                               GGCTGGCTTCTACCGCCGGGGCAGAGACATCGCATAAGGTGAGACGAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCC                                                       <  2:1639558/100‑1
                                               GGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGC                                                        <  1:2978638/100‑1
                                               GGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCC                                                       <  2:2402918/100‑1
                                                     CTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAAT                                                 <  1:1582814/100‑1
                                                      TTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATT                                                <  1:20609/100‑1
                                                       TCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTC                                               <  1:1492695/100‑1
                                                                           acCGCACAAGGTGAGACGACTTTTCGGCCGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATCCCCCGAAAATTCGTCTCCCCTTGTGCGATG                             <  2:2005860/96‑1
                                                                             CGCACAAGGTGAGACGAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTC                         <  1:1172681/100‑1
                                                                                                  TCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAG    >  2:1306505/1‑100
                                                                                                   CGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGC   >  2:2102576/1‑100
                                                                                                   CGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGC   >  1:2149035/1‑100
                                                                                                    GGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGCC  >  1:556816/1‑100
                                                                                                                                                                                                         
GCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4606945‑4607047
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcggc‑TGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGCC  <  NC_000913/4607088‑4606992

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