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.

AGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACC  >  NZ_CP009273/727747‑727885
                                                                                         |                                                 
aGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGAcgc                                                   <  2:2286460/90‑1 (MQ=33)
aGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGAcgc                                                   >  2:1745725/1‑90 (MQ=33)
aGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGAcgc                                                   >  2:1365917/1‑90 (MQ=33)
 gACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGAcgcg                                                  >  1:1404523/1‑90 (MQ=33)
 gACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGAcgcg                                                  >  2:3251059/1‑90 (MQ=33)
 gACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGAcgcg                                                  >  1:1265308/1‑90 (MQ=33)
 gACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGAcgcg                                                  >  2:2317552/1‑90 (MQ=33)
  aCCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACgcgc                                                 >  2:2224954/1‑90 (MQ=33)
   cccGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCa                                                >  1:239138/1‑90 (MQ=33)
   cccGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCa                                                <  1:317692/90‑1 (MQ=33)
   cccGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCa                                                <  1:2240077/90‑1 (MQ=33)
   cccGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCa                                                <  1:3251059/90‑1 (MQ=33)
   cccGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCa                                                >  1:884094/1‑90 (MQ=33)
   cccGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCa                                                >  1:967651/1‑90 (MQ=33)
   cccGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCa                                                >  2:216156/1‑90 (MQ=33)
    ccGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATACCCGTATCGCCCGTGACGCGCAc                                               <  2:16700/90‑1 (MQ=18)
    ccGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCAc                                               <  1:3346218/90‑1 (MQ=33)
    ccGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCAc                                               <  1:3230047/90‑1 (MQ=33)
    ccGGCAGGAAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCAc                                               <  2:2302841/90‑1 (MQ=18)
     cGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACt                                              <  1:1757626/90‑1 (MQ=33)
     cGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACt                                              <  1:1463171/90‑1 (MQ=33)
     cGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACt                                              <  1:1900768/90‑1 (MQ=33)
      ggCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGCCAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTa                                             >  1:2351409/1‑90 (MQ=18)
      ggCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCAc                                               >  1:144537/1‑88 (MQ=33)
      ggCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCAc                                               <  2:144537/88‑1 (MQ=33)
      ggCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTa                                             >  1:2365045/1‑90 (MQ=33)
      ggCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTa                                             >  1:885417/1‑90 (MQ=33)
           gTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCttt                                        >  2:2453350/1‑90 (MQ=33)
             aaCCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTAc                                      >  1:2626051/1‑90 (MQ=33)
              gccGCCTGCCCGACCCGGAGCTGCACCCGGATAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTAcc                                     <  2:3239751/89‑1 (MQ=18)
              aCCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTAcc                                     <  2:1408701/90‑1 (MQ=33)
                    tGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGCCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATg                               >  1:378025/1‑90 (MQ=18)
                     gCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGa                              <  1:215665/90‑1 (MQ=33)
                     gCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGa                              <  1:2490537/90‑1 (MQ=33)
                         gACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGt                          >  1:199643/1‑90 (MQ=35)
                               gAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGc                    <  2:1265308/90‑1 (MQ=33)
                                   tGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGc                <  1:1564838/90‑1 (MQ=33)
                                    gCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCa                                                >  2:1343258/1‑57 (MQ=38)
                                    gCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCa                                                <  1:1343258/57‑1 (MQ=35)
                                    gCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATc                                           >  1:1710203/1‑62 (MQ=37)
                                    gCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATc                                           <  2:1710203/62‑1 (MQ=35)
                                    gCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATc                                           >  1:2175477/1‑62 (MQ=37)
                                    gCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATc                                           <  2:2175477/62‑1 (MQ=35)
                                    gCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCt               >  1:2143534/1‑90 (MQ=35)
                                    gCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCt               >  2:2758869/1‑90 (MQ=35)
                                    gCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCt               >  1:2249405/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    >  2:1924406/1‑90 (MQ=33)
                                               gcgcCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGaaaa    >  1:741963/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)
                                                                                         |                                                 
AGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACC  >  NZ_CP009273/727747‑727885

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