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 28 [0] [0] 28 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)27 (0.920) 23/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

GCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACAC  >  NC_000913/3618943‑3619040
                                                |                                                 
gcggcgTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTt                                                  <  1:1150076‑M1/50‑2 (MQ=255)
gcggcgTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTt                                                  <  1:1877054‑M1/50‑2 (MQ=255)
    cgTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgc                                              >  1:1462609‑M1/1‑45 (MQ=255)
     gTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgcc                                             <  1:601878‑M1/50‑7 (MQ=255)
      tGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgcca                                            >  1:1210379‑M1/1‑43 (MQ=255)
          gCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgt                                        <  1:1816175‑M1/50‑12 (MQ=255)
            aCGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcc                                      <  1:2200500‑M1/50‑14 (MQ=255)
              gggCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtccc                                     >  1:1947790‑M1/1‑35 (MQ=255)
                 cATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatc                                 <  1:1440415‑M1/50‑19 (MQ=255)
                  aTTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcg                                >  1:304420‑M1/1‑31 (MQ=255)
                  aTTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcg                                >  1:1509537‑M1/1‑31 (MQ=255)
                      gCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtct                             <  1:189284‑M1/49‑23 (MQ=255)
                        aGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcc                          >  1:841355‑M1/1‑25 (MQ=255)
                            cTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaa                      >  1:940247‑M1/1‑21 (MQ=255)
                               gCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatcg                   >  1:1962732‑M1/1‑18 (MQ=255)
                               gCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatc                    <  1:661633‑M1/49‑32 (MQ=255)
                                 cGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatcgca                 <  1:1437606‑M1/50‑35 (MQ=255)
                                  gaaggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcag                >  1:1232369‑M1/1‑15 (MQ=255)
                                   aaggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagc               <  1:1782049‑M1/50‑37 (MQ=255)
                                   aaggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagc               <  1:1063158‑M1/50‑37 (MQ=255)
                                   aaggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcag                >  1:23333‑M1/1‑14 (MQ=255)
                                           ttGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaat       <  1:1148258‑M1/50‑45 (MQ=255)
                                           ttGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaat       <  1:908478‑M1/50‑45 (MQ=255)
                                            tGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatg      >  1:1061099‑M1/1‑5 (MQ=255)
                                              aGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatggt    <  1:547524‑M1/50‑48 (MQ=255)
                                              aGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatggt    >  1:1133965‑M1/1‑3 (MQ=255)
                                               gTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatggtg   <  1:926918‑M1/50‑49 (MQ=255)
                                                ttacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatggtga  <  1:1731731‑M1/50‑50 (MQ=255)
                                                |                                                 
GCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACAC  >  NC_000913/3618943‑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.