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 34 [0] [0] 33 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)30 (0.820) 24/96 0.2 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:1323473‑M1/50‑1 (MQ=255)
  cggcgTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTta                                                 >  1:861136‑M1/1‑48 (MQ=255)
  cggcgTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTta                                                 <  1:805909‑M1/50‑3 (MQ=255)
     cgTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgc                                              >  1:963133‑M1/1‑45 (MQ=255)
      gggcgGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGCtacgc                                              >  1:1693972‑M1/3‑44 (MQ=255)
      gTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgc                                              >  1:1105742‑M1/1‑44 (MQ=255)
          ggCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattg                                         >  1:2794323‑M1/1‑40 (MQ=255)
          gacacGGGCATTTGCAGTGCTTGCCGAAGGTAGATTGAGTtacgccattg                                         >  1:2024917‑M1/3‑40 (MQ=255)
           gCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgt                                        >  1:2496496‑M1/1‑39 (MQ=255)
            cACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtc                                       >  1:3959‑M1/1‑38 (MQ=255)
               gggCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtccc                                     >  1:688811‑M1/1‑35 (MQ=255)
                 gCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaat                                  <  1:1275359‑M1/50‑18 (MQ=255)
                  cATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatc                                 <  1:1953056‑M1/50‑19 (MQ=255)
                   aTTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcg                                <  1:1142511‑M1/50‑20 (MQ=255)
                      tGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtct                             >  1:265934‑M1/1‑28 (MQ=255)
                               tGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatc                    >  1:2254465‑M1/1‑19 (MQ=255)
                                  cGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatcgca                 <  1:479657‑M1/50‑35 (MQ=255)
                                    aaggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagc               <  1:1865683‑M1/50‑37 (MQ=255)
                                    aaggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcag                >  1:1109947‑M1/1‑14 (MQ=255)
                                    aaggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcag                >  1:2636948‑M1/1‑14 (MQ=255)
                                     aggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagca              <  1:1336067‑M1/50‑38 (MQ=255)
                                      ggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcag             >  1:1753871‑M1/1‑12 (MQ=255)
                                      ggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcag             >  1:1510161‑M1/1‑12 (MQ=255)
                                        aagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcg           <  1:2302409‑M1/50‑41 (MQ=255)
                                        aagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcg           <  1:949759‑M1/50‑41 (MQ=255)
                                           aTTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaa        <  1:581286‑M1/50‑44 (MQ=255)
                                           aTTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaa        >  1:1012160‑M1/1‑7 (MQ=255)
                                           aTTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaa        >  1:1757826‑M1/1‑7 (MQ=255)
                                             tGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatg      >  1:2443489‑M1/1‑5 (MQ=255)
                                              gAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatgg     >  1:1566729‑M1/1‑4 (MQ=255)
                                                gTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatggtg   <  1:1124781‑M1/50‑49 (MQ=255)
                                                 ttacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatggtga  <  1:2644140‑M1/50‑50 (MQ=255)
                                                 ttacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatggtg   >  1:2178998‑M1/1‑1 (MQ=255)
                                                 ttacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatggtg   >  1:183981‑M1/1‑1 (MQ=255)
                                                 |                                                 
CGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACAC  >  NC_000913/3618942‑3619040

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