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 57 [0] [0] 57 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)55 (0.770) 40/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

TGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGA  >  NC_000913/3625632‑3625728
                                                |                                                
cggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTa                                                 <  1:1087644‑M2/2‑1 (MQ=255)
cggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTa                                                 >  1:379223‑M2/49‑50 (MQ=255)
 ggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTAc                                                <  1:5358629‑M2/3‑1 (MQ=255)
   cgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGc                                              <  1:1409919‑M2/5‑1 (MQ=255)
   cgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGc                                              >  1:3632848‑M2/46‑50 (MQ=255)
    gtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGc                                              >  1:1092105‑M2/45‑49 (MQ=255)
    gtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGc                                              >  1:4108138‑M2/45‑49 (MQ=255)
    gtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGc                                              >  1:2296528‑M2/45‑49 (MQ=255)
     tgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCa                                            >  1:936627‑M2/44‑50 (MQ=255)
     tgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCa                                            <  1:2002734‑M2/7‑1 (MQ=255)
       cggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCAtt                                          >  1:775023‑M2/42‑50 (MQ=255)
        ggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTg                                         >  1:3689684‑M2/41‑50 (MQ=255)
        ggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTg                                         <  1:1203556‑M2/10‑1 (MQ=255)
         gcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGt                                        <  1:500360‑M2/11‑1 (MQ=255)
          cacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTc                                       >  1:2957627‑M2/39‑50 (MQ=255)
             gggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTccc                                     >  1:5229039‑M2/36‑49 (MQ=255)
             gggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCa                                    <  1:4132259‑M2/15‑1 (MQ=255)
             gggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCa                                    <  1:2569130‑M2/15‑1 (MQ=255)
               gcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAAt                                  <  1:1305306‑M2/17‑1 (MQ=255)
                catttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAatt                                 >  1:3710041‑M2/33‑49 (MQ=255)
                catttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATc                                 <  1:224349‑M2/18‑1 (MQ=255)
                 atttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCg                                >  1:3539403‑M2/32‑50 (MQ=255)
                  tttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGt                               >  1:703711‑M2/31‑50 (MQ=255)
                  tttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGt                               >  1:3640274‑M2/31‑50 (MQ=255)
                     gcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCtt                            >  1:4973831‑M2/28‑50 (MQ=255)
                      cagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCtt                            <  1:2048600‑M2/23‑1 (MQ=255)
                      cagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTc                           <  1:4309627‑M2/24‑1 (MQ=255)
                           cttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGaa                      >  1:2567640‑M2/22‑50 (MQ=255)
                              gccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAAt                     <  1:5435731‑M2/30‑1 (MQ=255)
                              gccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATc                    <  1:1906342‑M2/31‑1 (MQ=255)
                              gccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATc                    <  1:2037861‑M2/31‑1 (MQ=255)
                              gccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCg                   >  1:4497060‑M2/19‑50 (MQ=255)
                              gccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCg                   >  1:3124853‑M2/19‑50 (MQ=255)
                               ccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgc                  <  1:280188‑M2/33‑1 (MQ=255)
                               ccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgc                  <  1:4473250‑M2/33‑1 (MQ=255)
                               ccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCg                   <  1:2650607‑M2/32‑1 (MQ=255)
                                cgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCg                   >  1:3439418‑M2/17‑48 (MQ=255)
                                cgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCg                   >  1:2289812‑M2/17‑48 (MQ=255)
                                cgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCg                   >  1:3130392‑M2/17‑48 (MQ=255)
                                 gaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcag                >  1:4214966‑M2/16‑50 (MQ=255)
                                 gaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcag                >  1:3350057‑M2/16‑50 (MQ=255)
                                  aaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcagc               <  1:366834‑M2/36‑1 (MQ=255)
                                  aaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcag                >  1:5026549‑M2/15‑49 (MQ=255)
                                   aggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcag                >  1:5057772‑M2/14‑48 (MQ=255)
                                      aagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAgcg           >  1:5255947‑M2/11‑50 (MQ=255)
                                      aagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAgcg           >  1:5138656‑M2/11‑50 (MQ=255)
                                       agattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAgcg           >  1:1360906‑M2/10‑49 (MQ=255)
                                        gattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAgcgc          >  1:3859131‑M2/9‑49 (MQ=255)
                                         attgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCaa        >  1:4631034‑M2/8‑50 (MQ=255)
                                            gagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATg      <  1:5520342‑M2/45‑1 (MQ=255)
                                             agtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGt    <  1:4691‑M2/47‑1 (MQ=255)
                                             agtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGt    >  1:4084121‑M2/4‑50 (MQ=255)
                                             agtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGt    <  1:354602‑M2/47‑1 (MQ=255)
                                             agtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGt    <  1:2998861‑M2/47‑1 (MQ=255)
                                             agtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGt    <  1:2164944‑M2/47‑1 (MQ=255)
                                               tTACGCCATTGTCCCAATCTTCTTCCTGAATCGCAGCAGCGCAATGGTGa  <  1:4668930‑M2/49‑1 (MQ=255)
                                                tACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGa  >  1:1016243/1‑49 (MQ=255)
                                                |                                                
TGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGA  >  NC_000913/3625632‑3625728

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