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

TTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGA  >  NC_000913/3625631‑3625728
                                                 |                                                
gcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtt                                                  <  1:1150076‑M2/1‑1 (MQ=255)
gcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtt                                                  <  1:1877054‑M2/1‑1 (MQ=255)
    cgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGc                                              >  1:1462609‑M2/46‑50 (MQ=255)
     gtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGcc                                             <  1:601878‑M2/6‑1 (MQ=255)
      tgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCa                                            >  1:1210379‑M2/44‑50 (MQ=255)
          gcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGt                                        <  1:1816175‑M2/11‑1 (MQ=255)
            acgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTcc                                      <  1:2200500‑M2/13‑1 (MQ=255)
              gggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTccc                                     >  1:1947790‑M2/36‑49 (MQ=255)
                 catttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATc                                 <  1:1440415‑M2/18‑1 (MQ=255)
                  atttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCg                                >  1:304420‑M2/32‑50 (MQ=255)
                  atttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCg                                >  1:1509537‑M2/32‑50 (MQ=255)
                      gcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCt                             <  1:189284‑M2/22‑1 (MQ=255)
                        agtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTcc                          >  1:841355‑M2/26‑50 (MQ=255)
                            cttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGaa                      >  1:940247‑M2/22‑50 (MQ=255)
                               gccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATc                    <  1:661633‑M2/31‑1 (MQ=255)
                               gccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCg                   >  1:1962732‑M2/19‑50 (MQ=255)
                                 cgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgca                 <  1:1437606‑M2/34‑1 (MQ=255)
                                  gaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcag                >  1:1232369‑M2/16‑50 (MQ=255)
                                   aaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcagc               <  1:1782049‑M2/36‑1 (MQ=255)
                                   aaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcagc               <  1:1063158‑M2/36‑1 (MQ=255)
                                   aaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcag                >  1:23333‑M2/15‑49 (MQ=255)
                                           ttgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAAt       <  1:1148258‑M2/44‑1 (MQ=255)
                                           ttgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAAt       <  1:908478‑M2/44‑1 (MQ=255)
                                            tgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATg      >  1:1061099‑M2/6‑50 (MQ=255)
                                              agtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGt    <  1:547524‑M2/47‑1 (MQ=255)
                                              agtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGt    >  1:1133965‑M2/4‑50 (MQ=255)
                                               gtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTg   <  1:926918‑M2/48‑1 (MQ=255)
                                                tTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGa  <  1:1731731‑M2/49‑1 (MQ=255)
                                                 |                                                
TTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGA  >  NC_000913/3625631‑3625728

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: