New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4607047158 (1.140)13 (0.100) 7/168 NT 7.7% coding (654/1032 nt) rsmC 16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
?NC_000913 = 4607088 164 (1.260)coding (613/1032 nt) rsmC 16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
Rejected: Frequency below cutoff threshold.

TGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4606944‑4607047
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcggc‑TGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGCCG  <  NC_000913/4607088‑4606991
                                                                                                                                                                                                            
TGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTT                                                                                                          <  2:2418674/100‑1
 GCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTC                                                                                                         <  1:3698063/100‑1
 GCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTC                                                                                                         <  2:1651749/100‑1
 GCTGGCAAAGACGTCACCTTCATCACCGTGGGCCGCGAGTGTTGCGCGGCTGTCTTCTACCGCCGGGTCAGAGACATCGCACAAGGTGAGACGAATTTTC                                                                                                         <  1:1711422/100‑1
  CTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACCTCGCACAAGGTGAGACGAATTTTCG                                                                                                        >  2:947074/1‑100
  CTGGCAAAGACTTCACCTTCAACACCGGTGGCCGCAAGTGTTGCGCGGCTGGCTCCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTC                                                                                                         >  2:1921663/1‑99
   TGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGG                                                                                                       <  1:139541/100‑1
   TGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGG                                                                                                       >  1:562523/1‑100
   TGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGG                                                                                                       >  1:562658/1‑100
   TGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCGGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGG                                                                                                       <  1:3492915/100‑1
    GGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC                                                                                                      >  2:3517900/1‑100
    GGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC                                                                                                      >  2:783921/1‑100
    GGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC                                                                                                      >  2:1559909/1‑100
      CAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC                                                                                                      >  1:3735752/1‑98
                                              GCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTC                                                           <  2:31592/100‑1
                                              GCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTC                                                           <  1:2870966/100‑1
                                                GGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCC                                                        <  2:2286200/100‑1
                                                GGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCC                                                        <  1:1232133/100‑1
                                                 GCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCC                                                        <  2:2083395/100‑1
                                                 GCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCC                                                        <  1:721973/100‑1
                                                  CTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGA                                                      <  1:1009946/100‑1
                                                   TGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGA                                                      <  1:1269459/100‑1
                                                      CTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAAT                                                  <  1:845644/100‑1
                                                         tctCCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCG                                               <  1:2430115/97‑1
                                                         CTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCG                                               <  2:3582523/100‑1
                                                          TACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGT                                              <  1:2790514/100‑1
                                                            CCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGaga                                            >  2:2430133/1‑97
                                                                                                   TCGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAG     <  1:2002352/100‑1
                                                                                                    CGGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGC    >  1:1245177/1‑100
                                                                                                     GGCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCAATGTCTCTGCCCCGGCGGTAGAAGCCAGC    >  1:788963/1‑99
                                                                                                      GCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAACATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGCCG  >  1:182729/1‑100
                                                                                                      GCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGCCG  <  2:1083262/100‑1
                                                                                                      GCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGCCG  >  2:1305610/1‑100
                                                                                                      GCTGTGGCGCGGGCGTGCTTTCAGTTGCCTTTGCGCGCCATTCGCCGAAAATTCGTCTCACCTTGTGCGATGTCTCTGCCCCGGCGGTAGAAGCCAGCCG  >  1:149241/1‑100
                                                                                                                                                                                                            
TGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4606944‑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 < 32 ≤ 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.