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 42 [0] [0] 40 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)40 (1.100) 32/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

GACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTC  >  NC_000913/3761944‑3762038
                                                |                                              
gacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC                                                >  1:1147102‑M1/1‑49 (MQ=255)
gacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC                                                >  1:1159326‑M1/1‑49 (MQ=255)
 acgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                               >  1:829809‑M1/1‑48 (MQ=255)
 acgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                               >  1:913699‑M1/1‑48 (MQ=255)
 acgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                               <  1:2671036‑M1/49‑2 (MQ=255)
  cgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtc                                              <  1:1063726‑M1/49‑3 (MQ=255)
   gaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtca                                             <  1:1184533‑M1/49‑4 (MQ=255)
    aGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaa                                            >  1:2741958‑M1/1‑45 (MQ=255)
    aGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaa                                            >  1:2577250‑M1/1‑45 (MQ=255)
     gTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaat                                           <  1:2726386‑M1/49‑6 (MQ=255)
      tttCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaat                                           >  1:757469‑M1/1‑43 (MQ=255)
      tttCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaat                                           >  1:1397666‑M1/1‑43 (MQ=255)
          aTGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatc                                       >  1:1752667‑M1/1‑39 (MQ=255)
             ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaa                                   <  1:2599623‑M1/49‑14 (MQ=255)
             ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaa                                    <  1:2432653‑M1/48‑13 (MQ=255)
             ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaa                                    <  1:885153‑M1/48‑13 (MQ=255)
              cGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaat                                  <  1:2255750‑M1/49‑15 (MQ=255)
              cGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaat                                  <  1:1756416‑M1/49‑15 (MQ=255)
               gAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaat                                  >  1:236823‑M1/1‑34 (MQ=255)
                aGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                 >  1:635542‑M1/1‑33 (MQ=255)
                  tCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                 >  1:2324214‑M1/1‑31 (MQ=255)
                    cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                             <  1:2011414‑M1/48‑20 (MQ=255)
                    cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                             <  1:2749849‑M1/48‑20 (MQ=255)
                      ttGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataa                          >  1:1389682‑M1/1‑27 (MQ=255)
                        gtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                         <  1:712484‑M1/48‑24 (MQ=255)
                         tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatg                        >  1:1468911‑M1/1‑24 (MQ=255)
                             aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcc                   <  1:2348789‑M1/49‑30 (MQ=255)
                              ggAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                    >  1:1198036‑M1/1‑19 (MQ=255)
                               gAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccag                 >  1:2243797‑M1/1‑18 (MQ=255)
                                 aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct               >  1:828841‑M1/1‑16 (MQ=255)
                                  aaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg              <  1:1361229‑M1/49‑35 (MQ=255)
                                  aaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg              <  1:1195889‑M1/49‑35 (MQ=255)
                                   aaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctgg             >  1:1361677‑M1/1‑14 (MQ=255)
                                   aaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctgg             >  1:267660‑M1/1‑14 (MQ=255)
                                      atCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtca          >  1:1194004‑M1/1‑11 (MQ=255)
                                       tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcat         >  1:633385‑M1/1‑10 (MQ=255)
                                         aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattc       >  1:738502‑M1/1‑8 (MQ=255)
                                         aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcatt        <  1:2078164‑M1/48‑41 (MQ=255)
                                          gTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct      >  1:1019312‑M1/1‑7 (MQ=255)
                                            aTGGCtcaatgatccaaataaagataatgcatccagctggtcattct      >  1:1051132‑M1/1‑5 (MQ=255)
                                               gCtcaatgatccaaataaagataatgcatccagctggtcattcttagc  <  1:1124194‑M1/48‑47 (MQ=255)
                                               gCtcaatgatccaaataaagataatgcatccagctggtcattcttag   >  1:22964‑M1/1‑2 (MQ=255)
                                                |                                              
GACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTC  >  NC_000913/3761944‑3762038

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