New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 727836 =NA (NA)9 (0.080) 9/154 3.8 16.3% coding (2798/4194 nt) rhsC rhs element protein RhsC
?NZ_CP009273 727859 = 51 (0.430)coding (2821/4194 nt) rhsC rhs element protein RhsC
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGA  >  NZ_CP009273/727772‑727947
                                                                                       |                                                                                        
gACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGt                                                                                        >  1:199643/1‑90 (MQ=35)
      gAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGc                                                                                  <  2:1265308/90‑1 (MQ=33)
          tGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGc                                                                              <  1:1564838/90‑1 (MQ=33)
           gCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCt                                                                             >  1:2143534/1‑90 (MQ=35)
           gCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCt                                                                             >  1:2249405/1‑90 (MQ=35)
           gCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCt                                                                             >  2:2758869/1‑90 (MQ=35)
               ccGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACa                                                                         <  1:333515/90‑1 (MQ=33)
                cGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACag                                                                        >  2:2634209/1‑90 (MQ=35)
                cGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCGGACag                                                                        >  2:2745674/1‑90 (MQ=18)
                  gACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCACACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACagag                                                                      >  2:2814546/1‑90 (MQ=18)
                      gcgcCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGaaaa                                                                  >  1:704120/1‑90 (MQ=33)
                      gcgcCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGaaaa                                                                  >  1:741963/1‑90 (MQ=33)
                      gcgcCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGaaaa                                                                  >  2:1924406/1‑90 (MQ=33)
                      gcgcCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGaaaa                                                                  >  1:1416515/1‑90 (MQ=33)
                       cgcCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAAc                                                                 <  2:483535/90‑1 (MQ=32)
                        gcCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAAcc                                                                >  1:2128471/1‑90 (MQ=31)
                        gcCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAAcc                                                                >  2:645997/1‑90 (MQ=31)
                                                                  cACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAgcgc                      <  2:3168891/90‑1 (MQ=31)
                                                                  cACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAgcgc                      >  2:400787/1‑90 (MQ=31)
                                                                  cACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAgcgc                      >  2:1513243/1‑90 (MQ=31)
                                                                  cACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAgcgc                      <  1:934864/90‑1 (MQ=31)
                                                                    cTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcac                    <  1:1365917/90‑1 (MQ=33)
                                                                    cTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcac                    <  1:3345291/90‑1 (MQ=33)
                                                                    cTATCTTTACCGGTAGGACCGTCACGGCAGGCTGACAGAGAAAACCGCCCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcac                    <  2:659219/90‑1 (MQ=255)
                                                                     tATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcacc                   <  1:2656536/90‑1 (MQ=33)
                                                                        cTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcacccac                <  2:1404523/90‑1 (MQ=33)
                                                                             ccGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTa           <  2:775503/90‑1 (MQ=33)
                                                                              cGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGt            >  1:527716/1‑88 (MQ=35)
                                                                              cGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGt            <  2:527716/88‑1 (MQ=33)
                                                                              cGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTAc          >  1:54945/1‑90 (MQ=35)
                                                                              cGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTAc          >  1:2490705/1‑90 (MQ=35)
                                                                               ggTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTAcc         <  1:2285625/90‑1 (MQ=33)
                                                                                gTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcaccc                  <  2:2206916/80‑1 (MQ=34)
                                                                                gTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcaccc                  >  1:2206916/1‑80 (MQ=34)
                                                                                gTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCa        >  2:1937640/1‑90 (MQ=33)
                                                                                gTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCa        >  1:1025267/1‑90 (MQ=35)
                                                                                gTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCa        >  1:2048907/1‑90 (MQ=35)
                                                                                gTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCa        >  1:2172523/1‑90 (MQ=35)
                                                                                gTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCa        >  2:2920094/1‑90 (MQ=35)
                                                                                gTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCa        >  2:1624565/1‑90 (MQ=33)
                                                                                gTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGAAAGGGGTTATCCGCACGGATGATAAGCGCACCCACCGGTACCa        >  2:1437736/1‑90 (MQ=255)
                                                                                 tATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcac                    <  2:2725941/77‑1 (MQ=34)
                                                                                 tATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcac                    >  1:2725941/1‑77 (MQ=34)
                                                                                 tATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCAt       <  1:3346285/90‑1 (MQ=33)
                                                                                  aTGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCAtt      >  2:3193474/1‑90 (MQ=33)
                                                                                    gACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTAc    <  2:1345788/90‑1 (MQ=33)
                                                                                    gACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTAc    <  1:216156/90‑1 (MQ=33)
                                                                                    gACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTAc    >  1:793832/1‑90 (MQ=35)
                                                                                    gACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTAc    >  1:3283087/1‑90 (MQ=35)
                                                                                    gACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAACGCACCCACCGGTACCATTAc    >  2:1968455/1‑90 (MQ=18)
                                                                                      ccGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGa  >  2:2409038/1‑90 (MQ=33)
                                                                                       |                                                                                        
GACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGA  >  NZ_CP009273/727772‑727947

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