New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 440959535 (0.860)6 (0.150) 6/284 2.6 14.6% coding (321/732 nt) rlmB 23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent
?NC_000913 = 4409617 37 (0.950)coding (343/732 nt) rlmB 23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent
Rejected: Frequency below cutoff threshold.

CCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTT  >  NC_000913/4409468‑4409765
                                                                                                                                                     |                                                                                                                                                    
ccGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAg                                                                                                                                                     >  2:604493‑M2/1‑150 (MQ=255)
ccGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAg                                                                                                                                                     >  2:553474‑M2/1‑150 (MQ=255)
          gggCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAgtgcggatcgg                                                                                                                                           <  1:438920‑M2/151‑12 (MQ=255)
          gggCATTAGCCCCCGCGTGACGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGAcgc                                                                                                                                            >  2:19475/1‑150 (MQ=255)
             cATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAgtgcggatcggtta                                                                                                                                        >  2:620496‑M2/1‑137 (MQ=255)
                      ccGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCAt                                                                                                                               <  2:203015/151‑1 (MQ=255)
                           tGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGt                                                                                                                          >  1:173900/1‑151 (MQ=255)
                             aaGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTg                                                                                                                         <  2:524669/150‑1 (MQ=255)
                                ccAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGAt                                                                                                                       <  1:421793/149‑1 (MQ=255)
                                    gACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATtgtg                                                                                                                   <  2:363713/149‑1 (MQ=255)
                                      cGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGa                                                                                                               >  2:78493/1‑151 (MQ=255)
                                         cAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTACCGATCCGCACTACGCAGGCACGCGCCGAGGTTgtgcggatcggttacgccgtcgaggatcagcaggaacggttga                                                                                                           <  2:438911‑M2/151‑44 (MQ=255)
                                                 ggAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTc                                                                                                     >  2:605558/1‑150 (MQ=255)
                                                  gAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTc                                                                                                     <  1:460872/149‑1 (MQ=255)
                                                     aaCGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGcac                                                                                                <  1:477492/151‑1 (MQ=255)
                                                     aaCGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGcac                                                                                                <  2:345837/151‑1 (MQ=255)
                                                     aaCGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGcac                                                                                                <  2:259530/151‑1 (MQ=255)
                                                       cGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAg                                                                                              >  2:366335/1‑151 (MQ=255)
                                                          tCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCt                                                                                            >  1:19410/1‑150 (MQ=255)
                                                          tCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTc                                                                                           <  1:54134/151‑1 (MQ=255)
                                                             gCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCaa                                                                                         >  1:368/1‑150 (MQ=255)
                                                              ccGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACTACGCAGGCACGCGCCGAGGTTgtgcggatcggttacgccgtcgaggatcagcaggaacggttgatcgagcgaagcgatcagatc                                                                                       >  1:333227‑M2/1‑88 (MQ=255)
                                                              ccGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAAc                                                                                        <  2:532141/150‑1 (MQ=255)
                                                              ccGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAAc                                                                                        <  2:529662/150‑1 (MQ=255)
                                                               cGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGc                                                                                      <  2:166753/151‑1 (MQ=255)
                                                               cGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGc                                                                                      <  1:466053/151‑1 (MQ=255)
                                                                   tCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACg                                                                                  >  1:506893/1‑151 (MQ=255)
                                                                              cgctCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGt                                                                       <  2:536703/151‑1 (MQ=255)
                                                                                                  ctgATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTgcggcgtccgcact                                                    <  1:454131/150‑9 (MQ=255)
                                                                                                     aTCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCg                                                >  1:153392/1‑151 (MQ=255)
                                                                                                     aTCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCg                                                >  2:283446/1‑151 (MQ=255)
                                                                                                             cGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCg                                           >  1:108864/1‑148 (MQ=255)
                                                                                                                  tAACCGATCCGCACAACCTCGGCGCGTGCCTGCGCAGTGCGGACCCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGAGAGAAAAATAGCCTGCGGCGCGGCAGAAAGCGATCCGCTGATTCg                                   >  2:533692/1‑151 (MQ=255)
                                                                                                                    aCCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGAGGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCggg                                 <  1:177253/151‑1 (MQ=255)
                                                                                                                             cacaACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAAc                         <  1:90354/150‑1 (MQ=255)
                                                                                                                                 aCCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTgcggcg                                                            >  2:528723/1‑111 (MQ=255)
                                                                                                                                 aCCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTgcggcg                                                            <  1:528733/111‑1 (MQ=255)
                                                                                                                                 aCCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTg                     >  1:567980/1‑150 (MQ=255)
                                                                                                                                 aCCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTg                     >  2:191680/1‑150 (MQ=255)
                                                                                                                                       gcgcGTGCCTGCGTAgtgcggatcggttacgccgtcgaggatcagcaggaacggttgatcgagcgaagcgatcagatccggcagatcgttttcctgatactgacgtcctggcttcacgcgggcgataatgccctgatgcacggcaccgtcg              <  1:439202‑M2/151‑137 (MQ=255)
                                                                                                                                            tGCCTGCGTAgtgcggatcggttacgccgtcgaggatcagcaggaacggttgatcgagcgaagcgatcagatccggcagatcgttttcctgatactgacgtcctggcttcacgcgggcgataatgccctgatgcacggcaccgtcgctttt         <  1:97162‑M2/151‑142 (MQ=255)
                                                                                                                                            tGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATg          >  2:588948/1‑150 (MQ=255)
                                                                                                                                              ccTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGt       >  2:410117/1‑151 (MQ=255)
                                                                                                                                              ccTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGt       >  2:76386/1‑151 (MQ=255)
                                                                                                                                                   gTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGtt  >  2:425933/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTT  >  NC_000913/4409468‑4409765

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