New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4607578109 (0.880)7 (0.060) 3/172 NT 6.1% coding (123/1032 nt) rsmC 16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
?NC_000913 = 4607601 110 (0.920)coding (100/1032 nt) rsmC 16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
Rejected: Frequency below cutoff threshold.

GACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4607479‑4607578
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgCCCGCGCGTTTAGATACCGCGGCCAGCCGTGCTCATACCCAGCAATTCCACCACTGGCAGGTATTAAGCCGCCAGATGGGGGATAACGCCCGTTTCAGT  <  NC_000913/4607601‑4607503
                                                                                                                                                                                                       
GACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTG                                                                                                     <  2:1231252/100‑1
GACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTG                                                                                                     >  2:2125356/1‑100
GACGTCATCCACCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTG                                                                                                     <  1:1844488/100‑1
                           GAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGCCCGCGCGTTTAGATACCGCGGCCAGC                                                                          <  1:1063319/100‑1
                           GAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGCCCGCGCGTTTAGATACCGCGGCCAGC                                                                          <  2:1389137/100‑1
                           GAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGCCCGCGCGTTTAGATACCGCGGCCAGC                                                                          <  2:1595421/100‑1
                           GAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGCCCGCGCGTTTAGATACCGCGGCCAGC                                                                          <  2:1723948/100‑1
                           GAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGCCCGCGCGTTTAGATACCGCGGCCAGC                                                                          <  2:576620/100‑1
                               CGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGCCCGCGCGTTTAGATACCGCGGCCAGCCGT                                                                       <  2:107374/99‑1
                                   CGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGCCCGCGCGTTTAGATACCGCGGCCAGCCGTGCTCA                                                                  <  1:1501320/100‑1
                                                                                                  TGCCCGCGCGTTTAGATACCGCGGCCAGCCGTGCTCATACCCAGCAATTCCACCACTGGCAGGTATTAAGCCGCCAGATGGGGGATAACGCCCGTTTCAG   >  1:2379838/1‑100
                                                                                                  TGCCCGCGCGTTTAGATACCGCGGCCAGCCGTGCTCATACCCAGCAATTCCACCACTGGCAGGTATTAAGCCGCCAGATGGGGGATAACGCCCGTTTCAG   >  1:265152/1‑100
                                                                                                  TGCCCGCGCGTTTAGATACCGCGGCCAGCCGTGCTCATACCCAGCAATTCCACCACTGGCAGGTATTAAGCCGCCAGATGGGGGATAACGCCCGTTTCAG   <  2:1872127/100‑1
                                                                                                  TGCCCGCGCGTTTAGATACCGCGGCCAGCCGTGCTCATACCCAGCAATTCCACCACTGGCAGGTATTAAGCCGCCAGATGGGGGATAACGCCCGTTTCAG   >  2:504692/1‑100
                                                                                                   GCCCGCGCGTTTAGATACCGCGGCCAGCCGTGCTCATACCCAGCAATTCCACCACTGGCAGGTATTAAGCCGCCAGATGGGGGATAACGCCCGTTTCAGT  >  2:2613702/1‑100
                                                                                                                                                                                                       
GACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4607479‑4607578
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgCCCGCGCGTTTAGATACCGCGGCCAGCCGTGCTCATACCCAGCAATTCCACCACTGGCAGGTATTAAGCCGCCAGATGGGGGATAACGCCCGTTTCAGT  <  NC_000913/4607601‑4607503

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.