New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 72783557 (0.480)12 (0.110) 10/154 3.5 18.9% coding (2797/4194 nt) rhsC rhs element protein RhsC
?NZ_CP009273 = 729729 NA (NA)pseudogene (126/1506 nt) BW25113_RS03635 RHS repeat protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACC  >  NZ_CP009273/727746‑727885
                                                                                         |                                                  
caGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGAcg                                                    >  1:1240930/1‑90 (MQ=33)
caGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGAcg                                                    <  2:2648879/90‑1 (MQ=33)
 aGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGAcgc                                                   <  2:2286460/90‑1 (MQ=33)
 aGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGAcgc                                                   >  2:1745725/1‑90 (MQ=33)
 aGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGAcgc                                                   >  2:1365917/1‑90 (MQ=33)
  gACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGAcgcg                                                  >  1:1265308/1‑90 (MQ=33)
  gACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGAcgcg                                                  >  2:2317552/1‑90 (MQ=33)
  gACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGAcgcg                                                  >  1:1404523/1‑90 (MQ=33)
  gACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGAcgcg                                                  >  2:3251059/1‑90 (MQ=33)
   aCCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACgcgc                                                 >  2:2224954/1‑90 (MQ=33)
    cccGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCa                                                <  1:3251059/90‑1 (MQ=33)
    cccGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCa                                                >  2:216156/1‑90 (MQ=33)
    cccGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCa                                                >  1:967651/1‑90 (MQ=33)
    cccGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCa                                                >  1:884094/1‑90 (MQ=33)
    cccGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCa                                                <  1:317692/90‑1 (MQ=33)
    cccGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCa                                                >  1:239138/1‑90 (MQ=33)
    cccGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCa                                                <  1:2240077/90‑1 (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:1900768/90‑1 (MQ=33)
      cGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACt                                              <  1:1757626/90‑1 (MQ=33)
      cGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACt                                              <  1:1463171/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:885417/1‑90 (MQ=33)
       ggCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTa                                             >  1:2365045/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                                           <  2:2175477/62‑1 (MQ=35)
                                     gCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATc                                           >  1:1710203/1‑62 (MQ=37)
                                     gCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATc                                           <  2:1710203/62‑1 (MQ=35)
                                     gCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATc                                           >  1:2175477/1‑62 (MQ=37)
                                     gCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCt               >  2:2758869/1‑90 (MQ=35)
                                     gCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCt               >  1:2249405/1‑90 (MQ=35)
                                     gCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCt               >  1:2143534/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)
                                                                                         |                                                  
CAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACC  >  NZ_CP009273/727746‑727885

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