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

TGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4606944‑4607047
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcggc‑T‑GTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGCC  <  NC_000913/4607088‑4606992
                                                                                                                                                                                                           
TGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTT                                                                                                         >  1:2868009/1‑100
 GCTGGCAAAGACTTCAGCTTCAACACCGTTGGCCGGAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTC                                                                                                        <  1:2929128/100‑1
  CTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCG                                                                                                       >  2:3655528/1‑100
  CTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCG                                                                                                       <  1:1881442/100‑1
  CTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCG                                                                                                       <  1:3646885/100‑1
  CTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCG                                                                                                       <  1:3414212/100‑1
  CTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCG                                                                                                       >  1:3325994/1‑100
   TGGCAAAGACTTCACCTTCAACACCGTTGGCTGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGG                                                                                                      <  1:2025203/100‑1
   TGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGG                                                                                                      <  2:2388943/100‑1
   TGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCCCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGTCGAATTTTCGG                                                                                                      <  1:2820768/100‑1
    GGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC                                                                                                     <  1:4279299/100‑1
    GGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC                                                                                                     >  1:2410059/1‑100
    GGCAAAGACTTCACCTTCAACACCGTTGGCCACAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC                                                                                                     >  2:2101344/1‑100
                                       GTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑T‑GTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCG                                                                <  2:1276692/100‑1
                                          TTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGA‑ATGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCC                                                              <  2:2201378/100‑1
                                              GCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGA‑ATGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTC                                                          <  1:32412/100‑1
                                                GGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGA‑ATGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGn                                                        <  1:4290210/100‑2
                                                GGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑T‑GTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCC                                                       <  1:2671575/100‑1
                                                GGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑T‑GTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCC                                                       <  1:2054278/100‑1
                                                GGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑T‑GTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCC                                                       <  1:1424098/100‑1
                                                GGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑T‑GTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCC                                                       <  2:4027011/100‑1
                                                 GCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGA‑ATGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGac                                                       <  2:2990232/100‑3
                                                 GCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGA‑ATGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCC                                                       <  1:617448/100‑1
                                                 GCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑T‑GTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCG                                                      <  1:2200207/100‑1
                                                  CTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑GAATGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCG                                                      <  1:2745623/100‑1
                                                      CTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑T‑GTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAAT                                                 <  1:1104563/100‑1
                                                      CTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑T‑GTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAAT                                                 <  1:218503/100‑1
                                                                                                   TCGGC‑T‑GTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAG    >  1:3054982/1‑100
                                                                                                    CGGC‑T‑GTGGCGCGGGCGTGCTTTCCGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGC   <  2:4378823/100‑1
                                                                                                    CGGC‑T‑GTGGCGCGGGCGTGCTTTCAGTTGCCGTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAG    <  1:554001/99‑1
                                                                                                     GGC‑T‑GTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGCC  >  1:4359274/1‑100
                                                                                                     GGC‑T‑GTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGCC  >  1:1186876/1‑100
                                                                                                                                                                                                           
TGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4606944‑4607047
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcggc‑T‑GTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGCC  <  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.