New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 72511422 (0.710)3 (0.110) 3/158 3.1 10.7% coding (76/4194 nt) rhsC rhs element protein RhsC
?NZ_CP009273 = 725154 30 (1.050)coding (116/4194 nt) rhsC rhs element protein RhsC
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC  >  NZ_CP009273/725073‑725239
                                                                                 |                                                                                     
tGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCcggctg                                                                              <  2:259282/91‑3 (MQ=255)
       gTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACg                                                                       >  2:820479/1‑91 (MQ=255)
            gCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCgg                                                                  >  2:865885/1‑91 (MQ=255)
            gCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACGTCCgg                                                                  >  2:422502/1‑91 (MQ=255)
               gTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCa                                                               >  1:179965/1‑91 (MQ=255)
               gTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCa                                                               >  1:837822/1‑91 (MQ=255)
                 aGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATc                                                             >  1:171249/1‑91 (MQ=255)
                 aGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATc                                                             >  1:491055/1‑91 (MQ=255)
                 aGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATc                                                             >  1:408607/1‑91 (MQ=255)
                       gTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCa                                                       >  1:779622/1‑91 (MQ=255)
                        tCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCaa                                                      <  2:491055/91‑1 (MQ=255)
                        tCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCaa                                                      >  1:383348/1‑91 (MQ=255)
                         ccAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACGTCCGGCCATCCGGTCAAt                                                     >  2:591930/1‑91 (MQ=255)
                            gggTTCAGCCGGGGTACGCATCGGTGCCCCCCCCGGCGGGGCCTGTTCGGTGTGCCCCGGCGGGGGGCCGTCCGGCCATCCGGCCAATccc                                                  >  1:447636/1‑90 (MQ=14)
                             ggTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCgc                                                 <  2:84821/91‑1 (MQ=255)
                                 cAGCCGGGGTACGCATCGGTGCCCCCCCCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGCCGTCCGGCCATCCGGTCAATCCgctgct                                             >  1:504419/1‑91 (MQ=255)
                                 cAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCgctgct                                             >  2:540910/1‑91 (MQ=255)
                                 cAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCgctgct                                             >  1:391217/1‑91 (MQ=255)
                                  aGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTc                                            <  2:171249/91‑1 (MQ=255)
                                  aGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGTGACGGCCGGCCATCCGGTCAATCCGCTGCTc                                            <  2:187396/91‑1 (MQ=255)
                                         ggtcGCCTCCGTTCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCaa                                     <  2:447636/88‑1 (MQ=14)
                                          tACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCaaa                                    <  2:179965/91‑1 (MQ=255)
                                               aTCGGTGCCCCCACCGGCGTGGCCTGGTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTcc                               <  1:422502/91‑1 (MQ=255)
                                                 cGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCtt                             >  2:492128/1‑91 (MQ=38)
                                                  ggTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTc                            <  2:408607/91‑1 (MQ=37)
                                                   gTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTcc                           >  2:474743/1‑91 (MQ=37)
                                                   gTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTcc                           >  2:324323/1‑91 (MQ=37)
                                                                  tGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACAt            <  1:686859/91‑1 (MQ=37)
                                                                   ggCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATc           >  2:675413/1‑91 (MQ=37)
                                                                        gTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCt      >  2:100053/1‑91 (MQ=37)
                                                                            ggTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgccc  <  1:100053/91‑1 (MQ=37)
                                                                                 |                                                                                     
TGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC  >  NZ_CP009273/725073‑725239

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