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 24 [0] [0] 24 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)24 (0.890) 19/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

CGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACA  >  NC_000913/3618944‑3619039
                                               |                                                
cggcgTGCGGCACGGGCATTTGCAGTGCTTGCCTAAGGAAGATTGAGTta                                                >  1:239681‑M1/1‑48 (MQ=255)
    gTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgc                                             >  1:2243370‑M1/1‑44 (MQ=255)
    gTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgc                                             >  1:1090256‑M1/1‑44 (MQ=255)
     tGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgcca                                           <  1:1391503‑M1/50‑8 (MQ=255)
        ggCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattg                                        >  1:1396425‑M1/1‑40 (MQ=255)
         gCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgt                                       >  1:545516‑M1/1‑39 (MQ=255)
                cATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatc                                >  1:704017‑M1/1‑32 (MQ=255)
                   ttGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtc                             >  1:2091180‑M1/1‑29 (MQ=255)
                    tGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtc                             <  1:1772764‑M1/49‑22 (MQ=255)
                    tGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtc                             <  1:1497324‑M1/49‑22 (MQ=255)
                       aGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcc                         >  1:1694547‑M1/1‑25 (MQ=255)
                         tGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcct                        >  1:938151‑M1/1‑23 (MQ=255)
                                cGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatcgca                <  1:856088‑M1/50‑35 (MQ=255)
                                cGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatcg                  >  1:1476397‑M1/1‑16 (MQ=255)
                                 gaaggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcag               >  1:1979675‑M1/1‑15 (MQ=255)
                                 gaaggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcag               >  1:487720‑M1/1‑15 (MQ=255)
                                  aaggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagc              <  1:1618171‑M1/50‑37 (MQ=255)
                                   aggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagca             <  1:1679009‑M1/50‑38 (MQ=255)
                                   aggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagca             <  1:1881643‑M1/50‑38 (MQ=255)
                                   aggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagca             <  1:1518610‑M1/50‑38 (MQ=255)
                                     gaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagc           >  1:1644980‑M1/1‑11 (MQ=255)
                                       agATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcg          >  1:14934‑M1/1‑9 (MQ=255)
                                             aGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatggt   >  1:1600303‑M1/1‑3 (MQ=255)
                                              gTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatggtg  <  1:242491‑M1/50‑49 (MQ=255)
                                               |                                                
CGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACA  >  NC_000913/3618944‑3619039

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: