Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,618,992 Δ6,688 bp rhsB[yhhJ] rhsB, yhhH, yhhI, [yhhJ]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3618992 3625679 6688 51 [0] [0] 50 rhsB–[yhhJ] rhsB,yhhH,yhhI,[yhhJ]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 36189910 (0.000)49 (0.860) 33/96 0.1 100% intergenic (+2/‑201) nikR/rhsB DNA‑binding transcriptional repressor NikR/rhs element protein RhsB
?NC_000913 3625680 = NA (NA)coding (1124/1125 nt) yhhJ ABC transporter family protein YhhJ

CGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACAC  >  NC_000913/3618942‑3619040
                                                 |                                                 
cgcgGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGT                                                   >  1:362820‑M1/1‑50 (MQ=255)
   ggcgTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtac                                                <  1:1044819‑M1/50‑4 (MQ=255)
      gTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgc                                              >  1:3756953‑M1/1‑44 (MQ=255)
         cGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccatt                                          >  1:413214‑M1/1‑41 (MQ=255)
          ggCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattg                                         <  1:1129702‑M1/50‑11 (MQ=255)
             aCGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcc                                      <  1:3978975‑M1/50‑14 (MQ=255)
             aCGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcc                                      >  1:3679195‑M1/1‑37 (MQ=255)
             aCGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcc                                      >  1:295659‑M1/1‑37 (MQ=255)
             aCGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcc                                      >  1:2501973‑M1/1‑37 (MQ=255)
             aCGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcc                                      <  1:2350027‑M1/50‑14 (MQ=255)
             aCGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcc                                      >  1:2250292‑M1/1‑37 (MQ=255)
              cGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcc                                      >  1:1962457‑M1/1‑36 (MQ=255)
               gggCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtccca                                    <  1:295955‑M1/50‑16 (MQ=255)
                ggCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaa                                   >  1:1212323‑M1/1‑34 (MQ=255)
                   aTTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcg                                <  1:595581‑M1/50‑20 (MQ=255)
                   aTTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcg                                <  1:4046921‑M1/50‑20 (MQ=255)
                       gCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtct                             <  1:545485‑M1/49‑23 (MQ=255)
                       gCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtct                             <  1:3629083‑M1/49‑23 (MQ=255)
                        cAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttc                           <  1:843537‑M1/50‑25 (MQ=255)
                        cAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttc                           >  1:993005‑M1/1‑26 (MQ=255)
                        cAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttc                           <  1:1922862‑M1/50‑25 (MQ=255)
                         aGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcc                          >  1:2288949‑M1/1‑25 (MQ=255)
                         aGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcc                          >  1:1332712‑M1/1‑25 (MQ=255)
                          gTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcct                         >  1:649633‑M1/1‑24 (MQ=255)
                          gTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcct                         >  1:3298422‑M1/1‑24 (MQ=255)
                            gCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctg                        >  1:1644714‑M1/1‑22 (MQ=255)
                             cTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaa                      >  1:2052341‑M1/1‑21 (MQ=255)
                                gCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatc                    <  1:4082598‑M1/49‑32 (MQ=255)
                                  cGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatcg                   >  1:1257111‑M1/1‑16 (MQ=255)
                                  cGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatcg                   >  1:143470‑M1/1‑16 (MQ=255)
                                   gaaggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcag                <  1:4074676‑M1/50‑36 (MQ=255)
                                    aaggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcag                >  1:1391073‑M1/1‑14 (MQ=255)
                                    aaggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcag                >  1:1489873‑M1/1‑14 (MQ=255)
                                     aggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagca              <  1:1195026‑M1/50‑38 (MQ=255)
                                      ggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcag             <  1:4044999‑M1/50‑39 (MQ=255)
                                      ggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcag             >  1:3440641‑M1/1‑12 (MQ=255)
                                       gaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagc            <  1:811310‑M1/50‑40 (MQ=255)
                                       gaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagc            >  1:2486470‑M1/1‑11 (MQ=255)
                                        aagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcg           >  1:1344794‑M1/1‑10 (MQ=255)
                                         agATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcg           >  1:3813010‑M1/1‑9 (MQ=255)
                                         agATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcg           >  1:1878004‑M1/1‑9 (MQ=255)
                                           aTTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaa        >  1:4412692‑M1/1‑7 (MQ=255)
                                           aTTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaa        >  1:655355‑M1/1‑7 (MQ=255)
                                            ttGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaat       >  1:1770443‑M1/1‑6 (MQ=255)
                                             tGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatg      >  1:3536784‑M1/1‑5 (MQ=255)
                                             tGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatg      >  1:82816‑M1/1‑5 (MQ=255)
                                              gAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatg      <  1:1726384‑M1/49‑46 (MQ=255)
                                              gAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatg      <  1:184895‑M1/49‑46 (MQ=255)
                                              gAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatg      <  1:208692‑M1/49‑46 (MQ=255)
                                                gTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatggtg   <  1:2703186‑M1/50‑49 (MQ=255)
                                                 ttacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatggtga  <  1:2265916‑M1/50‑50 (MQ=255)
                                                 |                                                 
CGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACAC  >  NC_000913/3618942‑3619040

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