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

GCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGA  >  NC_000913/3625633‑3625728
                                               |                                                
ggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTAc                                                <  1:1044819‑M2/3‑1 (MQ=255)
   gtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGc                                              >  1:3756953‑M2/45‑49 (MQ=255)
      cggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCAtt                                          >  1:413214‑M2/42‑50 (MQ=255)
       ggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTg                                         <  1:1129702‑M2/10‑1 (MQ=255)
          acgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTcc                                      >  1:2250292‑M2/38‑50 (MQ=255)
          acgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTcc                                      <  1:2350027‑M2/13‑1 (MQ=255)
          acgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTcc                                      >  1:2501973‑M2/38‑50 (MQ=255)
          acgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTcc                                      >  1:295659‑M2/38‑50 (MQ=255)
          acgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTcc                                      >  1:3679195‑M2/38‑50 (MQ=255)
          acgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTcc                                      <  1:3978975‑M2/13‑1 (MQ=255)
           cgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTcc                                      >  1:1962457‑M2/37‑49 (MQ=255)
            gggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCa                                    <  1:295955‑M2/15‑1 (MQ=255)
             ggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCaa                                   >  1:1212323‑M2/35‑50 (MQ=255)
                atttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCg                                <  1:4046921‑M2/19‑1 (MQ=255)
                atttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCg                                <  1:595581‑M2/19‑1 (MQ=255)
                    gcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCt                             <  1:545485‑M2/22‑1 (MQ=255)
                    gcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCt                             <  1:3629083‑M2/22‑1 (MQ=255)
                     cagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTc                           <  1:843537‑M2/24‑1 (MQ=255)
                     cagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTc                           >  1:993005‑M2/27‑50 (MQ=255)
                     cagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTc                           <  1:1922862‑M2/24‑1 (MQ=255)
                      agtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTcc                          >  1:2288949‑M2/26‑50 (MQ=255)
                      agtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTcc                          >  1:1332712‑M2/26‑50 (MQ=255)
                       gtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCt                         >  1:649633‑M2/25‑50 (MQ=255)
                       gtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCt                         >  1:3298422‑M2/25‑50 (MQ=255)
                         gcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTg                        >  1:1644714‑M2/23‑49 (MQ=255)
                          cttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGaa                      >  1:2052341‑M2/22‑50 (MQ=255)
                             gccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATc                    <  1:4082598‑M2/31‑1 (MQ=255)
                               cgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCg                   >  1:143470‑M2/17‑48 (MQ=255)
                               cgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCg                   >  1:1257111‑M2/17‑48 (MQ=255)
                                gaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcag                <  1:4074676‑M2/35‑1 (MQ=255)
                                 aaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcag                >  1:1489873‑M2/15‑49 (MQ=255)
                                 aaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcag                >  1:1391073‑M2/15‑49 (MQ=255)
                                  aggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcagca              <  1:1195026‑M2/37‑1 (MQ=255)
                                   ggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGcagcag             >  1:3440641‑M2/13‑50 (MQ=255)
                                   ggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGcagcag             <  1:4044999‑M2/38‑1 (MQ=255)
                                    gaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCagcagc            <  1:811310‑M2/39‑1 (MQ=255)
                                    gaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCagcagc            >  1:2486470‑M2/12‑50 (MQ=255)
                                     aagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAgcg           >  1:1344794‑M2/11‑50 (MQ=255)
                                      agattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAgcg           >  1:3813010‑M2/10‑49 (MQ=255)
                                      agattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAgcg           >  1:1878004‑M2/10‑49 (MQ=255)
                                        attgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCaa        >  1:4412692‑M2/8‑50 (MQ=255)
                                        attgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCaa        >  1:655355‑M2/8‑50 (MQ=255)
                                         ttgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAAt       >  1:1770443‑M2/7‑50 (MQ=255)
                                          tgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATg      >  1:3536784‑M2/6‑50 (MQ=255)
                                          tgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATg      >  1:82816‑M2/6‑50 (MQ=255)
                                           gagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATg      <  1:208692‑M2/45‑1 (MQ=255)
                                           gagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATg      <  1:1726384‑M2/45‑1 (MQ=255)
                                           gagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATg      <  1:184895‑M2/45‑1 (MQ=255)
                                             gtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTg   <  1:2703186‑M2/48‑1 (MQ=255)
                                              tTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGa  <  1:2265916‑M2/49‑1 (MQ=255)
                                               |                                                
GCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGA  >  NC_000913/3625633‑3625728

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.