Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,761,993 Δ8,250 bp rhsAyibV rhsA, yibA, yibJ, yibG, yibS, yibV

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3761993 3770242 8250 43 [0] [0] 44 rhsA–yibV rhsA,yibA,yibJ,yibG,yibS,yibV

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 37619920 (0.000)41 (1.530) 33/94 0.0 100% intergenic (‑37/‑191) yibF/rhsA glutathione transferase‑like protein YibF/rhs element protein RhsA
?NC_000913 3770243 = 0 (0.000)coding (1137/1137 nt) yibH inner membrane protein YibH

TAATTAGTCCTTACAAGGAGCAAAGTTCCATAATCTCCCCCCTCCCCCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTAGCAG  >  NC_000913/3770195‑3770290
                                                |                                               
acgagtttcatgccgagtcctttgtgcgaggaaaaatatcagtatggct                                                 <  1:1457613‑M2/1‑1 (MQ=255)
   agtttcatgccgagtcctttgtgcgaggaaaaatatcagtatggcTCaa                                              >  1:1821039‑M2/46‑49 (MQ=255)
   agtttcatgccgagtcctttgtgcgaggaaaaatatcagtatggcTCaa                                              >  1:1650562‑M2/46‑49 (MQ=255)
    gtttcatgccgagtcctttgtgcgaggaaaaatatcagtatggcTCAAt                                             <  1:109991‑M2/5‑1 (MQ=255)
    gtttcatgccgagtcctttgtgcgaggaaaaatatcagtatggcTCAAt                                             >  1:1960459‑M2/45‑49 (MQ=255)
     tttcatgccgagtcctttgtgcgaggaaaaatatcagtatggcTCAAt                                             >  1:338507‑M2/44‑48 (MQ=255)
       tcatgccgagtcctttgtgcgaggaaaaatatcagtatggcTCAATGa                                           <  1:417642‑M2/7‑1 (MQ=255)
        catgccgagtcctttgtgcgaggaaaaatatcagtatggcTCAATGAt                                          <  1:269862‑M2/8‑1 (MQ=255)
         atgccgagtcctttgtgcgaggaaaaatatcagtatggcTCAATGATc                                         >  1:854359‑M2/40‑48 (MQ=255)
          tgccgagtcctttgtgcgaggaaaaatatcagtatggcTCAATGATcc                                        >  1:1675505‑M2/39‑48 (MQ=255)
             cgagtcctttgtgcgaggaaaaatatcagtatggcTCAATGATCCAAAt                                    <  1:1134098‑M2/14‑1 (MQ=255)
              gagtcctttgtgcgaggaaaaatatcagtatggcTCAATGATCCAAAt                                    >  1:1401579‑M2/35‑48 (MQ=255)
              gagtcctttgtgcgaggaaaaatatcagtatggcTCAATGATCCAAAt                                    >  1:497017‑M2/35‑48 (MQ=255)
                gtcctttgtgcgaggaaaaatatcagtatggcTCAATGATCCAAATaaa                                 <  1:1578349‑M2/17‑1 (MQ=255)
                   ctttgtgcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGa                               >  1:70265‑M2/30‑48 (MQ=255)
                   ctttgtgcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGa                               >  1:225828‑M2/30‑48 (MQ=255)
                   ctttgtgcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGa                               <  1:2039765‑M2/19‑1 (MQ=255)
                      tgtgcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATa                             <  1:1286253‑M2/21‑1 (MQ=255)
                        tgcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATg                          <  1:143412‑M2/24‑1 (MQ=255)
                         gcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCa                        >  1:1161200‑M2/24‑49 (MQ=255)
                         gcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCa                        >  1:1705723‑M2/24‑49 (MQ=255)
                          cgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCAt                       <  1:880030‑M2/27‑1 (MQ=255)
                          cgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCAt                       >  1:1750503‑M2/23‑49 (MQ=255)
                          cgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCAt                       <  1:1239158‑M2/27‑1 (MQ=255)
                            aggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATc                      >  1:1436324‑M2/21‑48 (MQ=255)
                            aggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATc                      >  1:665326‑M2/21‑48 (MQ=255)
                              gaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAg                   >  1:1854056‑M2/19‑49 (MQ=255)
                              gaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAg                   <  1:1332230‑M2/31‑1 (MQ=255)
                              gaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAg                   >  1:461646‑M2/19‑49 (MQ=255)
                                  aatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTgg               >  1:194271‑M2/15‑49 (MQ=255)
                                  aatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTgg               >  1:758218‑M2/15‑49 (MQ=255)
                                   atatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGt              >  1:1807228‑M2/14‑49 (MQ=255)
                                     atcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCa            >  1:1722464‑M2/12‑49 (MQ=255)
                                      tcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCAt           <  1:715862‑M2/39‑1 (MQ=255)
                                        agtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTc         >  1:1195342‑M2/9‑49 (MQ=255)
                                        agtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTc         <  1:119202‑M2/41‑1 (MQ=255)
                                          tatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCtt       <  1:1838995‑M2/43‑1 (MQ=255)
                                            tggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTag     <  1:1230367‑M2/45‑1 (MQ=255)
                                            tggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTa      <  1:212668‑M2/44‑1 (MQ=255)
                                              gcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTagca   <  1:1285627‑M2/47‑1 (MQ=255)
                                              gcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTag     >  1:313593‑M2/3‑47 (MQ=255)
                                               cTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTagcag  >  1:1454354‑M2/2‑49 (MQ=255)
                                                tcAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTagcag  >  1:273417/1‑48 (MQ=255)
                                                tcAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTagcag  >  1:2003888/1‑48 (MQ=255)
                                                |                                               
TAATTAGTCCTTACAAGGAGCAAAGTTCCATAATCTCCCCCCTCCCCCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTAGCAG  >  NC_000913/3770195‑3770290

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