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 72 [0] [0] 73 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)68 (0.970) 40/66 0.0 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

GTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAG  >  NC_000913/3625646‑3625714
                                  |                                  
ggcatttgcagtgcttgccgaaggaagattgagtt                                    <  1:6795626‑M2/1‑1 (MQ=255)
 gcatttgcagtgcttgccgaaggaagattgagtTa                                   <  1:4856848‑M2/2‑1 (MQ=255)
 gcatttgcagtgcttgccgaaggaagattgagtTa                                   <  1:6438076‑M2/2‑1 (MQ=255)
 gcatttgcagtgcttgccgaaggaagattgagtTa                                   <  1:660975‑M2/2‑1 (MQ=255)
  catttgcagtgcttgccgaaggaagattgagtTAc                                  <  1:7452006‑M2/3‑1 (MQ=255)
    tttgcagtgcttgccgaaggaagattgagtTACGc                                >  1:1067600‑M2/31‑35 (MQ=255)
    tttgcagtgcttgccgaaggaagattgagtTACGc                                <  1:6057709‑M2/5‑1 (MQ=255)
    tttgcagtgcttgccgaaggaagattgagtTACGc                                >  1:7764165‑M2/31‑35 (MQ=255)
    tttgcagtgcttgccgaaggaagattgagtTACGc                                >  1:6349422‑M2/31‑35 (MQ=255)
    tttgcagtgcttgccgaaggaagattgagtTACGc                                >  1:413778‑M2/31‑35 (MQ=255)
    tttgcagtgcttgccgaaggaagattgagtTACGc                                <  1:6919763‑M2/5‑1 (MQ=255)
    tttgcagtgcttgccgaaggaagattgagtTACGc                                >  1:6925606‑M2/31‑35 (MQ=255)
     ttgcagtgcttgccgaaggaagattgagtTACGcc                               <  1:7655779‑M2/6‑1 (MQ=255)
      tgcagtgcttgccgaaggaagattgagtTACGCCa                              <  1:3242472‑M2/7‑1 (MQ=255)
      tgcagtgcttgccgaaggaagattgagtTACGCCa                              <  1:2669013‑M2/7‑1 (MQ=255)
      tgcagtgcttgccgaaggaagattgagtTACGCCa                              >  1:2279105‑M2/29‑35 (MQ=255)
       gcagtgcttgccgaaggaggattgagtTACGCCAt                             >  1:2760748‑M2/28‑35 (MQ=255)
       gcagtgcttgccgaaggaagattgagtTACGCCAt                             <  1:5605649‑M2/8‑1 (MQ=255)
        cagtgcttgccgaaggaagattgagtTACGCCAtt                            <  1:6058953‑M2/9‑1 (MQ=255)
        cagtgcttgccgaaggaagattgagtTACGCCAtt                            >  1:6168352‑M2/27‑35 (MQ=255)
         agtgcttgccgaaggaagattgagtTACGCCATTg                           >  1:5721371‑M2/26‑35 (MQ=255)
         agtgcttgccgaaggaagattgagtTACGCCATTg                           <  1:4889412‑M2/10‑1 (MQ=255)
         agtgcttgccgaaggaagattgagtTACGCCATTg                           <  1:1234718‑M2/10‑1 (MQ=255)
         agtgcttgccgaaggaagattgagtTACGCCATTg                           >  1:3455079‑M2/26‑35 (MQ=255)
         agtgcttgccgaaggaagattgagtTACGCCATTg                           >  1:7115391‑M2/26‑35 (MQ=255)
          gtgcttgccgaaggaagattgagtTACGCCATTGt                          <  1:2983863‑M2/11‑1 (MQ=255)
           tgcttgccgaaggaagattgagtTACGCCATTGTc                         <  1:3826138‑M2/12‑1 (MQ=255)
           tgcttgccgaaggaagattgagtTACGCCATTGTc                         >  1:1865502‑M2/24‑35 (MQ=255)
            gcttgccgaaggaagattgagtTACGCCATTGTcc                        <  1:2517355‑M2/13‑1 (MQ=255)
             cttgccgaaggaagattgagtTACGCCATTGTccc                       >  1:1340872‑M2/22‑35 (MQ=255)
             cttgccgaaggaagattgagtTACGCCATTGTccc                       <  1:7238495‑M2/14‑1 (MQ=255)
             cttgccgaaggaagattgagtTACGCCATTGTccc                       >  1:1217594‑M2/22‑35 (MQ=255)
             cttgccgaaggaagattgagtTACGCCATTGTccc                       >  1:524902‑M2/22‑35 (MQ=255)
             cttgccgaaggaagattgagtTACGCCATTGTccc                       >  1:5467316‑M2/22‑35 (MQ=255)
             cttgccgaaggaagattgagtTACGCCATTGTccc                       >  1:699248‑M2/22‑35 (MQ=255)
                gccgaaggaagattgagtTACGCCATTGTCCCAAt                    <  1:3958309‑M2/17‑1 (MQ=255)
                gccgaaggaagattgagtTACGCCATTGTCCCAAt                    <  1:6204024‑M2/17‑1 (MQ=255)
                gccgaaggaagattgagtTACGCCATTGTCCCAAt                    >  1:6509812‑M2/19‑35 (MQ=255)
                gccgaaggaagattgagtTACGCCATTGTCCCAAt                    <  1:113445‑M2/17‑1 (MQ=255)
                   gaaggaagattgagtTACGCCATTGTCCCAATCGt                 >  1:4308474‑M2/16‑35 (MQ=255)
                    aaggaagattgagtTACGCCATTGTCCCAATCGTc                >  1:2183532‑M2/15‑35 (MQ=255)
                     aggaagattgagtTACGCCATTGTCCCAATCGTCt               >  1:3371206‑M2/14‑35 (MQ=255)
                      ggaagattgagtTACGCCATTGTCCCAATCGTCtt              >  1:848167‑M2/13‑35 (MQ=255)
                      ggaagattgagtTACGCCATTGTCCCAATCGTCtt              >  1:5836060‑M2/13‑35 (MQ=255)
                      ggaagattgagtTACGCCATTGTCCCAATCGTCtt              <  1:4451485‑M2/23‑1 (MQ=255)
                      ggaagattgagtTACGCCATTGTCCCAATCGTCtt              <  1:2197697‑M2/23‑1 (MQ=255)
                      ggaagattgagtTACGCCATTGTCCCAATCGCCtt              >  1:7108645‑M2/13‑35 (MQ=255)
                       gaagattgagtTACGCCATTGTCCCAATCGTCTTc             <  1:1136273‑M2/24‑1 (MQ=255)
                       gaagattgagtTACGCCATTGTCCCAATCGTCTTc             >  1:2872256‑M2/12‑35 (MQ=255)
                        aagattgagtTACGCCATTGTCCCAATCGTCTTcc            >  1:1070585‑M2/11‑35 (MQ=255)
                        aagattgagtTACGCCATTGTCCCAATCGTCTTcc            >  1:438197‑M2/11‑35 (MQ=255)
                         agattgagtTACGCCATTGTCCCAATCGTCTTCCt           <  1:5798042‑M2/26‑1 (MQ=255)
                         agattgagtTACGCCATTGTCCCAATCGTCTTCCt           <  1:1832281‑M2/26‑1 (MQ=255)
                         agattgagtTACGCCATTGTCCCAATCGTCTTCCt           >  1:6916947‑M2/10‑35 (MQ=255)
                         agattgagtTACGCCATTGTCCCAATCGTCTTCCt           >  1:3680409‑M2/10‑35 (MQ=255)
                         agattgagtTACGCCATTGTCCCAATCGTCTTCCt           >  1:6423418‑M2/10‑35 (MQ=255)
                         agattgagtTACGCCATTGTCCCAATCGTCTTCCt           <  1:4381039‑M2/26‑1 (MQ=255)
                           attgagtTACGCCATTGTCCCAATCGTCTTCCTGa         <  1:2958901‑M2/28‑1 (MQ=255)
                           attgagtTACGCCATTGTCCCAATCGTCTTCCTGa         >  1:4992365‑M2/8‑35 (MQ=255)
                            ttgagtTACGCCATTGTCCCAATCGTCTTCCTGaa        >  1:6179824‑M2/7‑35 (MQ=255)
                            ttgagtTACGCCATTGTCCCAATCGTCTTCCTGaa        <  1:5620790‑M2/29‑1 (MQ=255)
                            ttgagtTACGCCATTGTCCCAATCGTCTTCCTGaa        <  1:2122264‑M2/29‑1 (MQ=255)
                            ttgagtTACGCCATTGTCCCAATCGTCTTCCTGaa        >  1:4057699‑M2/7‑35 (MQ=255)
                             tgagtTACGCCATTGTCCCAATCGTCTTCCTGAAt       <  1:7388220‑M2/30‑1 (MQ=255)
                             tgagtTACGCCATTGTCCCAATCGTCTTCCTGAAt       <  1:2797295‑M2/30‑1 (MQ=255)
                              gagtTACGCCATTGTCCCAATCGTCTTCCTGAATc      >  1:2233457‑M2/5‑35 (MQ=255)
                               agtTACGCCATTGTCCCAATCGTCTTCCTGAATCg     >  1:7089617‑M2/4‑35 (MQ=255)
                                gtTACGCCATTGTCCCAATCGTCTTCCTGAATCGc    >  1:4938886‑M2/3‑35 (MQ=255)
                                 tTACGCCATTGTCCCAATCGTCTTCCTGAATCGCa   >  1:6262238‑M2/2‑35 (MQ=255)
                                 tTACGCCATTGTCCCAATCGTCTTCCTGAATCGCa   >  1:3885148‑M2/2‑35 (MQ=255)
                                  tACGCCATTGTCCCAATCGTCTTCCTGAATCgcag  <  1:2793774/35‑1 (MQ=255)
                                  tACGCCATTGTCCCAATCGTCTTCCTGAATCgcag  <  1:4450604/35‑1 (MQ=255)
                                  tACGCCATTGTCCCAATCGTCTTCCTGAATCgcag  <  1:7546999/35‑1 (MQ=255)
                                  |                                  
GTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAG  >  NC_000913/3625646‑3625714

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