New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4398991 =42 (1.920)3 (0.170) 3/140 1.5 9.4% intergenic (+79/‑101) rnr/rlmB ribonuclease R/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
?NZ_CP009273 4399052 = 23 (1.280)intergenic (+140/‑40) rnr/rlmB ribonuclease R/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Rejected: Frequency below/above cutoff threshold.

GCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACA  >  NZ_CP009273/4398904‑4399079
                                                                                       |                                                                                        
gCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                                                        >  1:126689/1‑90 (MQ=255)
  agagTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                                                      <  2:211094/90‑1 (MQ=255)
  agagTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                                                      <  2:420638/90‑1 (MQ=255)
  agagTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                                                      <  2:478914/90‑1 (MQ=255)
  agagTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                                                      <  1:149623/90‑1 (MQ=255)
   agcagggcctgcactgcgtggatgccgtaaatcatttcgctcattaatgtactcgttgtttacgttttgcctgatgcactacgtTTATcc                                                                                     >  1:409509‑M1/85‑90 (MQ=255)
   agcagggcctgcactgcgtggatgccgtaaatcatttcgctcattaatgtactcgttgtttacgttttgcctgatgcactacgtTTATcc                                                                                     >  2:62026‑M1/85‑90 (MQ=255)
     gTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                   <  1:125303/90‑1 (MQ=255)
     gTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                   <  1:363375/90‑1 (MQ=255)
     cagggcctgcactgcgtggatgccgtaaatcatttcgctcattaatgtactcgttgtttacgttttgcctgatgcactacgtTTATCCgg                                                                                   >  1:422101‑M1/83‑90 (MQ=255)
      agggcctgcactgcgtggatgccgtaaatcatttcgctcattaatgtactcgttgtttacgttttgcctgatgcactacgtTTATCCGGc                                                                                  >  1:13244‑M1/82‑90 (MQ=255)
               cccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAAt                                                                         >  1:208896/1‑90 (MQ=255)
                  tctTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCAc                                                                      >  1:293557/1‑90 (MQ=255)
                  gcgtggatgccgtaaatcatttcgctcattaatgtactcgttgtttacgttttgcctgatgcactacgtTTATcc                                                                                     <  2:577286‑M1/6‑1 (MQ=255)
                  gcgtggatgccgtaaatcatttcgctcattaatgtactcgttgtttacgttttgcctgatgcactacgtTTATcc                                                                                     >  1:577286‑M1/70‑75 (MQ=255)
                   ctTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACg                                                                     >  1:565145/1‑90 (MQ=255)
                   ctTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACg                                                                     >  1:494794/1‑90 (MQ=255)
                     ttAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCa                                                                   <  1:11366/90‑1 (MQ=255)
                     ttAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTCCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCa                                                                   <  2:297078/90‑1 (MQ=255)
                       aaaaGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCaaa                                                                 <  1:264414/90‑1 (MQ=255)
                          agaagaGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTc                                                              <  2:248021/90‑1 (MQ=255)
                            aagaGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa                                                            <  1:341902/90‑1 (MQ=255)
                             agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                  >  2:444338/1‑67 (MQ=255)
                             agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                  <  1:444338/67‑1 (MQ=255)
                             agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCa                                                                   >  1:527957/1‑82 (MQ=255)
                             agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCa                                                                   <  2:527957/82‑1 (MQ=255)
                             agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAat                                                           >  1:426619/1‑90 (MQ=255)
                             agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAat                                                           >  1:400439/1‑90 (MQ=255)
                             agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAat                                                           >  2:496648/1‑90 (MQ=255)
                             agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAat                                                           >  1:614784/1‑90 (MQ=255)
                             agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAat                                                           >  1:33747/1‑90 (MQ=255)
                                    tAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGc                                                    >  2:590248/1‑90 (MQ=255)
                                     aGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                                                   <  1:397129/90‑1 (MQ=255)
                                         tGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAgaga                                               <  2:521691/90‑1 (MQ=255)
                                         tGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAgaga                                               <  2:422566/90‑1 (MQ=255)
                                          gCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGAt                                              >  1:498317/1‑90 (MQ=255)
                                            tGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCa                                            >  1:297228/1‑90 (MQ=255)
                                                   gtactcgttgtttacgttttgcctgatgcactacgtTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGc                                     <  2:346713‑M1/54‑1 (MQ=255)
                                                   aTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGc                                     <  2:330092/90‑1 (MQ=255)
                                                    tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGcc                                    <  2:208896/90‑1 (MQ=255)
                                                       gcTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGc                                     <  1:74286/86‑1 (MQ=255)
                                                       gcTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGc                                     >  2:74286/1‑86 (MQ=255)
                                                         tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATa                               <  2:170301/90‑1 (MQ=255)
                                                         tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATa                               <  2:565145/90‑1 (MQ=255)
                                                           tgtttacgttttgcctgatgcactacgtTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATaaa                             <  1:49508‑M1/62‑1 (MQ=255)
                                                                     ggATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCAt                   <  2:498317/90‑1 (MQ=255)
                                                                      gATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATc                  <  2:33747/90‑1 (MQ=255)
                                                                           ggcgTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCa             <  1:307619/90‑1 (MQ=255)
                                                                                 ctacgtTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGGGCATCAGGCAAAACGt       <  1:102291‑M1/84‑1 (MQ=255)
                                                                                    gCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTaaa    <  1:303095/90‑1 (MQ=255)
                                                                                      cTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAaaca  <  2:604880/90‑1 (MQ=255)
                                                                                       |                                                                                        
GCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACA  >  NZ_CP009273/4398904‑4399079

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