New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4398991 =38 (1.830)3 (0.180) 3/136 1.6 10.3% intergenic (+79/‑101) rnr/rlmB ribonuclease R/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
?NZ_CP009273 4399052 = 21 (1.240)intergenic (+140/‑40) rnr/rlmB ribonuclease R/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Rejected: Frequency below/above cutoff threshold.

TGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAG  >  NZ_CP009273/4398902‑4399057
                                                                                         |                                                                  
tGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGATAGATTGCTGACAAAAGGCGCTGTGGTCATTCCGGATGCGGCGTGCACGCCt                                                                    >  2:509850/1‑90 (MQ=37)
    cagcagggcctgcactgcgtggatgccgtaaatcatttcgctcattaatgtactcgttgtttacgttttgcctgatgcactacgtTTATc                                                                >  2:166524‑M1/86‑90 (MQ=255)
    agagTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGTCTTATc                                                                <  1:302221/90‑1 (MQ=255)
     agcagggcctgcactgcgtggatgccgtaaatcatttcgctcattaatgtactcgttgtttacgttttgcctgatgcactacgtTTATcc                                                               >  1:621864‑M1/85‑90 (MQ=255)
       cagggcctgcactgcgtggatgccgtaaatcatttcgctcattaatgtactcgttgtttacgttttgcctgatgcactacgtTTATCCgg                                                             >  2:97887‑M1/83‑90 (MQ=255)
         gggcctgcactgcgtggatgccgtaaatcatttcgctcattaatgtactcgttgtttacgttttgcctgatgcactacgtTTATCCGGcc                                                           <  2:533419‑M1/10‑1 (MQ=255)
         gggcctgcactgcgtggatgccgtaaatcatttcgctcattaatatactcgttgtttacgttttgcctgatgcactacgtTTATCCGGcc                                                           >  2:166809‑M1/81‑90 (MQ=255)
          ggcctgcactgcgtggatgccgtaaatcatttcgctcattaatgtactcgttgtttacgttttgcctgatgcactacgtTTATCCGGCCt                                                          >  1:34063‑M1/80‑90 (MQ=255)
           gcctgcactgcgtggatgccgtaaatcatttcgctcattaatgtactcgttgtttacgttttgcctgatgcactacgtTTATCCgg                                                             >  2:368458‑M1/79‑86 (MQ=255)
           gcctgcactgcgtggatgccgtaaatcatttcgctcattaatgtactcgttgtttacgttttgcctgatgcactacgtTTATCCgg                                                             <  1:368458‑M1/8‑1 (MQ=255)
           gcctgcactgcgtggatgccgtaaatcatttcgctcattaatgtactcgttgtttacgttttgcctgatgcactacgtTTATCCGGc                                                            <  1:367191‑M1/9‑1 (MQ=255)
           gcctgcactgcgtggatgccgtaaatcatttcgctcattaatgtactcgttgtttacgttttgcctgatgcactacgtTTATCCGGc                                                            >  2:367191‑M1/79‑87 (MQ=255)
           gcctgcactgcgtggatgccgtaaatcatttcgctcattaatgtactcgttgtttacgttttgcctgatgcactacgtTTATCCGGCCTa                                                         <  1:453997‑M1/12‑1 (MQ=255)
                 cccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAAt                                                   >  1:459660/1‑90 (MQ=255)
                   ctctTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCa                                                 >  1:96794/1‑90 (MQ=255)
                   ctctTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCa                                                 >  2:482712/1‑90 (MQ=255)
                    gcgtggatgccgtaaatcatttcgctcattaatgtactcgttgtttacgttttgcctgatgcactacgtTTATCCGGc                                                            <  1:587149‑M1/9‑1 (MQ=255)
                    gcgtggatgccgtaaatcatttcgctcattaatgtactcgttgtttacgttttgcctgatgcactacgtTTATCCGGc                                                            >  2:587149‑M1/70‑78 (MQ=255)
                     ctTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACg                                               >  1:535597/1‑90 (MQ=255)
                     ctTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACg                                               >  2:375894/1‑90 (MQ=255)
                      tttAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGc                                              <  2:260831/90‑1 (MQ=255)
                      tttAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGc                                              <  1:482712/90‑1 (MQ=255)
                       ttAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCa                                             <  1:603693/90‑1 (MQ=255)
                        tAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCaa                                            >  2:159634/1‑90 (MQ=255)
                            agaagaGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTc                                        <  1:607921/90‑1 (MQ=255)
                              aagaGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa                                      <  2:210494/90‑1 (MQ=255)
                              aagaGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa                                      <  2:459660/90‑1 (MQ=255)
                              aagaGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa                                      <  1:375894/90‑1 (MQ=255)
                               agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                             <  1:452917/66‑1 (MQ=255)
                               agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                             >  2:452917/1‑66 (MQ=255)
                               agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                        <  2:501569/71‑1 (MQ=255)
                               agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                        >  1:501569/1‑71 (MQ=255)
                               agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAat                                     >  2:175197/1‑90 (MQ=255)
                               agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAat                                     >  2:93188/1‑90 (MQ=255)
                               agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAat                                     >  1:393256/1‑90 (MQ=255)
                                gagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAata                                    <  1:449106/90‑1 (MQ=255)
                                   ggTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatatat                                 <  2:146684/90‑1 (MQ=255)
                                         aTTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCaga                           <  2:349294/90‑1 (MQ=255)
                                            gCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCaa                                            >  1:579903/1‑70 (MQ=255)
                                            gCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCaa                                            <  2:579903/70‑1 (MQ=255)
                                            gCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCa                                             >  2:59689/1‑69 (MQ=255)
                                            gCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCa                                             <  1:59689/69‑1 (MQ=255)
                                            gCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                             <  1:563058/85‑1 (MQ=255)
                                            gCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                             >  2:563058/1‑85 (MQ=255)
                                            gCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGAt                        >  2:249150/1‑90 (MQ=255)
                                                      tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGcc              <  2:96794/90‑1 (MQ=255)
                                                      tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGcc              <  1:267760/90‑1 (MQ=255)
                                                          cgttgtttacgttttgcctgatgcactacgtTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGAt          <  1:365816‑M1/59‑1 (MQ=255)
                                                          cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGAt          <  2:393256/90‑1 (MQ=255)
                                                            ttgtttacgttttgcctgatgcactacgtTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATaa        <  2:275004‑M1/61‑1 (MQ=255)
                                                                  acgttttgcctgatgcactacgtTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAg  <  2:233205‑M1/67‑1 (MQ=255)
                                                                  aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAg  <  2:488654/90‑1 (MQ=255)
                                                                   tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa                                      <  1:337165/53‑1 (MQ=255)
                                                                   tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa                                      >  2:337165/1‑53 (MQ=255)
                                                                                         |                                                                  
TGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAG  >  NZ_CP009273/4398902‑4399057

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