New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4399442 =39 (1.340)3 (0.120) 3/150 2.0 8.0% coding (351/732 nt) rlmB 23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
?NZ_CP009273 4399525 = 35 (1.360)coding (434/732 nt) rlmB 23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Rejected: Frequency below/above cutoff threshold.

CTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCA  >  NZ_CP009273/4399353‑4399529
                                                                                         |                                                                                       
cTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGAGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGAc                                                                                         <  2:13949/90‑1 (MQ=255)
         tcgctCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGc                                                                                >  2:604746/1‑90 (MQ=255)
           gctCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGt                                                                              <  1:265690/90‑1 (MQ=255)
               gATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCAt                                                                          >  2:359960/1‑90 (MQ=255)
                  cAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCg                                                                       <  2:431954/90‑1 (MQ=255)
                   aaCCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCgg                                                                      >  1:384592/1‑90 (MQ=255)
                         tCCTGCTGATCCTCGGCGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATACGGTGATtg                                                                >  1:378152/1‑90 (MQ=255)
                         tCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATtg                                                                >  2:284071/1‑90 (MQ=255)
                         tCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATtg                                                                >  1:700749/1‑90 (MQ=255)
                          ccTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCCCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATtgt                                                               <  1:250432/90‑1 (MQ=255)
                            tgctgATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGc                                                             >  1:60071/1‑90 (MQ=255)
                             gctgTTCCTCGATGGCGTAACCGATCCGCACAACCTCGGCGCTTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGcc                                                            <  1:589863/90‑1 (MQ=255)
                               tgATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGa                                                          <  1:468727/90‑1 (MQ=255)
                                  tCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGCGATTGTGCCGAAAg                                                       >  1:200359/1‑90 (MQ=255)
                                   tacgcatggtgcgcgcaaggttagtcacccgaatcaacggaaagctttctgcccAGCCGCAGGCGTTCATGCGGGGAATTTGCCGAAAGa                                                      <  1:592243‑M1/36‑1 (MQ=255)
                                    cTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGAt                                                     <  1:83278/90‑1 (MQ=255)
                                       gACGGCGTAACCGATCCGCACAACCTCGGGGCGGGCCTGCGTATTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGAgcgc                                                  >  1:26079/1‑90 (MQ=255)
                                       gACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCg                                                                       <  1:404448/69‑1 (MQ=255)
                                       gACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCg                                                                       >  2:404448/1‑69 (MQ=255)
                                       gACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCAGGCGGTGATTGTGCCGAAAGATCGc                                                  >  1:185704/1‑90 (MQ=255)
                                         cGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTc                                                <  1:604746/90‑1 (MQ=255)
                                           gCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCg                                              >  2:572965/1‑90 (MQ=255)
                                            cGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGc                                             <  2:641454/90‑1 (MQ=255)
                                              tAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGcac                                           <  2:80299/90‑1 (MQ=255)
                                              tAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGcac                                           <  2:60071/90‑1 (MQ=255)
                                               aaCCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGcaca                                          >  2:575512/1‑90 (MQ=255)
                                                     ggttaatcacccgaatcagcggaacgctttctgccgCGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCaa                                    <  2:198295‑M1/54‑1 (MQ=255)
                                                          acaACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTa                               >  1:153719/1‑90 (MQ=255)
                                                            aaCCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACg                             <  2:427993/90‑1 (MQ=255)
                                                             aCCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACgg                            <  2:678239/90‑1 (MQ=255)
                                                             aCCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACgg                            <  2:700749/90‑1 (MQ=255)
                                                                tCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGa                         >  1:750099/1‑90 (MQ=255)
                                                                 cGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGaa                        <  2:33352/90‑1 (MQ=255)
                                                                       gTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGt                  >  2:546091/1‑90 (MQ=255)
                                                                       gTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGt                  >  2:283886/1‑90 (MQ=255)
                                                                       gTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGt                  >  1:29163/1‑90 (MQ=255)
                                                                            acgctttctgccgCGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCt             <  2:246593‑M1/77‑1 (MQ=255)
                                                                                  gTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTgcggcg       >  1:753369/1‑90 (MQ=255)
                                                                                  gTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTgcggcg       >  2:647568/1‑90 (MQ=255)
                                                                                  gTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTgcggcg       >  2:802784/1‑90 (MQ=255)
                                                                                    gCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTgcggcgcg     <  2:400323/90‑1 (MQ=255)
                                                                                       gACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCa  <  1:802784/90‑1 (MQ=255)
                                                                                         |                                                                                       
CTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCA  >  NZ_CP009273/4399353‑4399529

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