New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4606703103 (0.910)9 (0.080) 6/168 NT 8.1% coding (998/1032 nt) rsmC 16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
?NC_000913 = 4606716 108 (1.010)coding (985/1032 nt) rsmC 16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
Rejected: Frequency below cutoff threshold.

AACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4606602‑4606703
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgcTTCAAGGTGTATCGCGCCATTATGACCCGCCAGGCGAAGAAAGGTTAATTATCTCGCCGAATACCGTGTCGGATGCGGCGTGAACGCCTTATCCGACCT  <  NC_000913/4606716‑4606618
                                                                                                                                                                                                         
AACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGC                                                                                                       >  1:2245480/1‑100
AACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGC                                                                                                       >  2:2027203/1‑100
 ACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCG                                                                                                      <  2:2786604/100‑1
  CGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGC                                                                                                     >  2:843369/1‑100
  CGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGC                                                                                                     <  2:2100368/100‑1
                     GGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCTTCAAGGTGTATCGCGCCA                                                                                  <  1:2057670/100‑1
                      GATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCTTCAAGGTGTATCGCGCCAT                                                                                 <  1:802525/100‑1
                      GATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCTTCAAGGTGTATCGCGCCAT                                                                                 <  2:2950438/100‑1
                       ATAAGGCGTTCACGCCGCATCCGACACGGTATTCTGCGAGATAATTAACCTTTCTGCGCCTGGCGGGTCATAATGGCGCTTCAAGGTGTATCGCGCCAT                                                                                 <  1:155250/99‑1
                        TAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCTTCAAGGTGTATCGCGCCATTA                                                                               <  1:2882418/100‑1
                         AAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCTTCAAGGTGTATCGCGCCATTAT                                                                              <  2:2401951/100‑1
                          AGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCTTCAAGGTGTATCGCGCCATTATG                                                                             <  2:198528/100‑1
                          AGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCTTCAAGGTGTATCGCGCCATTATG                                                                             <  1:2844266/100‑1
                              GTTCTCGCCGAACACGACACGGTATTCGGCGAGATAATTCACCTTTCTTCGCCTGGCGGGTCATAATGGCGCTTCAAGGTGTATCGCGCCATTATGACC                                                                          <  2:606556/99‑1
                                                                                                  GCGCTTCAAGGTGTATCGCGCCATTATGACCCGCCAGGCGAAGAAAGGTTAATTATCTCGCCGAATACCGTGTCGGATGCGGCGTGAACGCCTTATCCGA     >  2:2652575/1‑100
                                                                                                    GCTTCAAGGTGTATCGCGCCATTATGACCCGCCAGGCGAAGAAAGGTTAATTATCTCGCCGAATACCGTGTCGGATGCGGCGTGAACGCCTTATCCGACC   <  2:687215/100‑1
                                                                                                    GCTTCAAGGTGTAGCGCGCCATTATGAACCGCCAGGACGATAAAGGCTAATTATCTCGCCGAATACCGTGTCGGATGCGGCGTGAACGCCTTA          >  1:1887969/1‑93
                                                                                                     CTTCAAGGTGTATCGCGCCATTATGACCCGCCAGGCGAAGAAAGGTTAATTATCTCGCCGAATACCGTGTCGGATGCGGCGTGAACGCCTTATCCGACCT  >  2:1278060/1‑100
                                                                                                     CTTCAAGGTGTATCGCGCCATTATGACCCGCCAGGCGAAGAAAGGTTAATTATCTCGCCGAATACCGTGTCGGATGCGGCGTGAACGCCTTATCCGACCT  >  2:2916881/1‑100
                                                                                                                                                                                                         
AACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4606602‑4606703
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgcTTCAAGGTGTATCGCGCCATTATGACCCGCCAGGCGAAGAAAGGTTAATTATCTCGCCGAATACCGTGTCGGATGCGGCGTGAACGCCTTATCCGACCT  <  NC_000913/4606716‑4606618

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