New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4607519 =93 (1.070)4 (0.050) 4/172 NT 4.2% coding (182/1032 nt) rsmC 16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
?NC_000913 4607556 = 91 (1.080)coding (145/1032 nt) rsmC 16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
Rejected: Frequency below cutoff threshold.

AGACTTACAGGATGACCTGCCCGCGCGTTTAGATACCGCGGCCAGCCGTGCTCATACCCAGCAATTCCACCACTGGCAGGTATTAAGCCGCCAGATGGGGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4607620‑4607519
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggaATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGCGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGTTCGAAATCA  >  NC_000913/4607556‑4607653
                                                                                                                                                                                                        
AGACTTACAGGATGACCTGCCCGCGCGTTTAGATCCCGCGGCCAGCCGTGCTCATACCCAGCAATTCCACCACTGGCAGGTATTAAGCCGCCAGATGGGG                                                                                                      <  1:1442877/100‑1
 GACTTACAGGATGACCTGCCCGCGCGTTTAGATACCGCGGCCAGCCGTGCTCATACCCAGCAATTCCACCACTGGCAGGTATTAAGCCGCCAGATGGGGG                                                                                                     <  1:1185600/100‑1
 GACTTACAGGATGACCTGCCCGCGCGTTTAGATACCGCGGCCAGCCGTGCTCATACCCAGCAATTCCACCACTGGCAGGTATTAAGCCGCCAGATGGGGG                                                                                                     >  1:1969041/1‑100
                                                      TACCCAGCAATTCCACCACTGGCAGGTATTAAGCCGCCAGATGGGGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGCGCGGGCAGGTC                                                >  2:927508/1‑100
                                                           AGCAATTCCACCACTGGCAGGTATTAAGCCGCCAGATGGGGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGCGCGGGCAGGTCATCCT                                           >  1:1070333/1‑100
                                                                      CACTGGCAGGTATTAAGCCGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGCGCGGGCAGGTCATCCTGTAAGTCTCCGG                               >  1:312392/1‑100
                                                                       ACTGGCAGGTATTAAGCCGCCAGATGGGGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGCGCGGGCAGGTCATCCTGTAAGTCTCCGG                               >  2:419476/1‑100
                                                                                                    GAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGCGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGTTCGAAATCA  <  2:1443214/100‑1
                                                                                                                                                                                                        
AGACTTACAGGATGACCTGCCCGCGCGTTTAGATACCGCGGCCAGCCGTGCTCATACCCAGCAATTCCACCACTGGCAGGTATTAAGCCGCCAGATGGGGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4607620‑4607519
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggaATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGCGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGTTCGAAATCA  >  NC_000913/4607556‑4607653

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