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

GACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCA  >  NC_000913/3761944‑3762040
                                                |                                                
gacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC                                                  >  1:1048323‑M1/1‑49 (MQ=255)
 acgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                                 <  1:1457613‑M1/49‑2 (MQ=255)
    aGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaa                                              >  1:1821039‑M1/1‑45 (MQ=255)
    aGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaa                                              >  1:1650562‑M1/1‑45 (MQ=255)
     gTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaat                                             <  1:109991‑M1/49‑6 (MQ=255)
     gTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaat                                             >  1:1960459‑M1/1‑44 (MQ=255)
      tttCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaat                                             >  1:338507‑M1/1‑43 (MQ=255)
        tCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatga                                           <  1:417642‑M1/48‑8 (MQ=255)
         cATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgat                                          <  1:269862‑M1/48‑9 (MQ=255)
          aTGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatc                                         >  1:854359‑M1/1‑39 (MQ=255)
           tGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatcc                                        >  1:1675505‑M1/1‑38 (MQ=255)
              cGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaat                                    <  1:1134098‑M1/49‑15 (MQ=255)
               gAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaat                                    >  1:497017‑M1/1‑34 (MQ=255)
               gAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaat                                    >  1:1401579‑M1/1‑34 (MQ=255)
                 gTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaa                                 <  1:1578349‑M1/49‑18 (MQ=255)
                    cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                               >  1:225828‑M1/1‑29 (MQ=255)
                    cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                               <  1:2039765‑M1/48‑20 (MQ=255)
                    cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                               >  1:70265‑M1/1‑29 (MQ=255)
                       tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagata                             <  1:1286253‑M1/47‑22 (MQ=255)
                         tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatg                          <  1:143412‑M1/48‑25 (MQ=255)
                          gCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgca                        >  1:1161200‑M1/1‑23 (MQ=255)
                          gCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgca                        >  1:1705723‑M1/1‑23 (MQ=255)
                           cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcat                       <  1:880030‑M1/49‑28 (MQ=255)
                           cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcat                       >  1:1750503‑M1/1‑22 (MQ=255)
                           cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcat                       <  1:1239158‑M1/49‑28 (MQ=255)
                             aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                      >  1:1436324‑M1/1‑20 (MQ=255)
                             aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                      >  1:665326‑M1/1‑20 (MQ=255)
                               gAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccag                   >  1:1854056‑M1/1‑18 (MQ=255)
                               gAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccag                   <  1:1332230‑M1/49‑32 (MQ=255)
                               gAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccag                   >  1:461646‑M1/1‑18 (MQ=255)
                                   aaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctgg               >  1:194271‑M1/1‑14 (MQ=255)
                                   aaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctgg               >  1:758218‑M1/1‑14 (MQ=255)
                                    atatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt              >  1:1807228‑M1/1‑13 (MQ=255)
                                      atCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtca            >  1:1722464‑M1/1‑11 (MQ=255)
                                       tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcat           <  1:715862‑M1/49‑40 (MQ=255)
                                         aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattc         <  1:119202‑M1/49‑42 (MQ=255)
                                         aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattc         >  1:1195342‑M1/1‑8 (MQ=255)
                                           tATGGCtcaatgatccaaataaagataatgcatccagctggtcattctt       <  1:1838995‑M1/49‑44 (MQ=255)
                                             tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag     <  1:1230367‑M1/49‑46 (MQ=255)
                                             tGGCtcaatgatccaaataaagataatgcatccagctggtcattctta      <  1:212668‑M1/48‑45 (MQ=255)
                                               gCtcaatgatccaaataaagataatgcatccagctggtcattcttagca   <  1:1285627‑M1/49‑48 (MQ=255)
                                               gCtcaatgatccaaataaagataatgcatccagctggtcattcttag     >  1:313593‑M1/1‑2 (MQ=255)
                                                ctcaatgatccaaataaagataatgcatccagctggtcattcttagcag  >  1:1454354‑M1/1‑1 (MQ=255)
                                                |                                                
GACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCA  >  NC_000913/3761944‑3762040

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.