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.

GCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTG  >  NZ_CP009273/4399437‑4399600
                                                                                        |                                                                           
gCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTgcggcgcg                                                                            <  2:400323/90‑1 (MQ=255)
   gACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCa                                                                         <  1:802784/90‑1 (MQ=255)
           aggctacttttttcgccgtagcgttgagctgtgcggagcgatctttcggcacaatcaccgcatgaacgcctgcggcgCGGCAGAAAGCGt                                                                 <  2:246593‑M2/13‑1 (MQ=255)
               gTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCg                                                             >  2:127925/1‑90 (MQ=255)
                 tCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCt                                                           >  2:553678/1‑90 (MQ=255)
                 tCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCt                                                           >  2:229905/1‑90 (MQ=255)
                       ggTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCg                                                     <  2:770493/90‑1 (MQ=255)
                         tGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCggg                                                   <  1:732257/90‑1 (MQ=255)
                          gATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGt                                                  >  1:303314/1‑90 (MQ=255)
                            ttGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGa                                                >  2:836638/1‑90 (MQ=255)
                              gtgCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACt                                              <  1:546091/90‑1 (MQ=255)
                              gtgCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTa                                             >  1:252692/1‑90 (MQ=255)
                                  ttgagctgtgcggagcgatctttcggcacaatcaccgcatgaacgcctgcggcgCGGCAGAAAGCGTTCCGCTGATTCGGGTGATTAAcc                                          <  2:198295‑M2/36‑1 (MQ=255)
                                        aTCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTgcgc                                    <  1:647568/90‑1 (MQ=255)
                                          cGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGGGACTAACCTTgcgcgc                                  >  1:510190/1‑90 (MQ=255)
                                           gCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCa                                 <  2:13260/90‑1 (MQ=255)
                                            cTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCAc                                <  1:534128/90‑1 (MQ=255)
                                            cTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCAc                                <  1:62582/90‑1 (MQ=255)
                                              ccGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCa                              <  1:127925/90‑1 (MQ=255)
                                                    tctttcggcaaattccccgcatgaacgcctgcggctGGGCAGAAAGCTTTCCGTTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTa                        <  1:592243‑M2/54‑1 (MQ=255)
                                                    aGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTa                        <  2:174920/90‑1 (MQ=255)
                                                     gCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTAt                       >  1:286397/1‑90 (MQ=255)
                                                          aCGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGc                  <  1:309504/90‑1 (MQ=255)
                                                               aCGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCAt                             <  2:292537/74‑1 (MQ=255)
                                                               aCGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCAt                             >  1:292537/1‑74 (MQ=255)
                                                               aCGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGaa             >  1:588944/1‑90 (MQ=255)
                                                               aCGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGaa             >  2:779761/1‑90 (MQ=255)
                                                               aCGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGaa             >  2:597771/1‑90 (MQ=255)
                                                               aCGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGaa             >  2:492989/1‑90 (MQ=255)
                                                               aCGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGaa             >  1:593172/1‑90 (MQ=255)
                                                               aCGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGaa             >  2:173973/1‑90 (MQ=255)
                                                               aCGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGaa             >  2:116777/1‑90 (MQ=255)
                                                               aCGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGaa             >  1:94115/1‑90 (MQ=255)
                                                                   cGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAGCCTTGCGCGCACCATGCGTATGTTGCAGGAAgaga         <  1:474031/90‑1 (MQ=255)
                                                                   cGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAgaga         <  1:229905/90‑1 (MQ=255)
                                                                   cGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAgaga         <  1:198295/90‑1 (MQ=255)
                                                                   cGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAgaga         <  1:359960/90‑1 (MQ=255)
                                                                          aGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTg  <  2:297877/90‑1 (MQ=255)
                                                                                        |                                                                           
GCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTG  >  NZ_CP009273/4399437‑4399600

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