New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4607047182 (1.160)11 (0.070) 9/168 NT 6.0% coding (654/1032 nt) rsmC 16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
?NC_000913 = 4607088 175 (1.180)coding (613/1032 nt) rsmC 16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
Rejected: Frequency below cutoff threshold.

GCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4606945‑4607047
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcggc‑TGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGCCG  <  NC_000913/4607088‑4606991
                                                                                                                                                                                                           
GCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTC                                                                                                         >  2:2184426/1‑100
GCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTC                                                                                                         >  1:3867241/1‑100
 CTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCG                                                                                                        <  1:1358626/100‑1
 CTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCG                                                                                                        <  2:907190/100‑1
 CTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCG                                                                                                        >  2:2717106/1‑100
   GGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC                                                                                                      >  1:3966246/1‑100
    GCAACGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGCGAGACGAATTTTCGGC                                                                                                      <  1:3803877/99‑1
                                          TGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCA                                                              <  2:207937/100‑1
                                             GCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTC                                                           <  2:4147711/100‑1
                                               GGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCC                                                        <  1:2923517/100‑1
                                                GCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCG                                                       <  2:369043/100‑1
                                                GCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCC                                                        <  1:991442/100‑1
                                                 CTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCC                                                        <  1:1794588/98‑1
                                                   GGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAA                                                    <  1:2760413/100‑1
                                                     CTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAAT                                                  <  2:3537900/100‑1
                                                       TCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTC                                                <  1:474409/100‑1
                                                        CTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCG                                               <  1:626263/100‑1
                                                          ACCGCCGGGGCAGAGACATCGCACAAGGTGAGACTAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTC                                             <  2:4335251/100‑1
                                                                                                  TCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAG     >  2:3672832/1‑100
                                                                                                   CGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCCGCGGTAAAAGCCAGC    >  1:4016373/1‑100
                                                                                                   CGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTATGCGATGTATCTTCCCCGGCGGTAGAAGCCAac    >  1:745633/1‑98
                                                                                                    GGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGCC   >  1:496343/1‑100
                                                                                                    GGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGCC   >  2:2378040/1‑100
                                                                                                    GGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGC    >  1:1082236/1‑99
                                                                                                     GCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGCCG  >  1:3990016/1‑100
                                                                                                                                                                                                           
GCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4606945‑4607047
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcggc‑TGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGCCG  <  NC_000913/4607088‑4606991

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.