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 32 [0] [0] 34 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)32 (1.070) 29/96 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

GACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3761944‑3761992
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTAGCAG  >  NC_000913/3770243‑3770290
                                                                                                 
GACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCT                                                 >  1:2638844/1‑50
 ACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCTC                                                <  1:1654134/50‑1
 ACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCT                                                 >  1:1322750/1‑49
  CGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCTCA                                               <  1:1526618/50‑1
   GAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCTCAA                                              >  1:1470146/1‑50
       TTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCTCAATGAT                                          >  1:1370085/1‑50
          ATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCTCAATGATCCA                                       <  1:107884/50‑1
             CCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAA                                     <  1:1901577/49‑1
              CGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATA                                   <  1:2520894/50‑1
               GAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAA                                  <  1:1882503/50‑1
                    CTTTGTGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAAAGAT                              <  1:687101/49‑1
                     TTTGTGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAAAGATAA                            >  1:1712327/1‑50
                       TGTGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAAAGATAAT                           <  1:2329414/49‑1
                       TGTGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAAAGATAAT                           <  1:2190445/49‑1
                        GTGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAAAGATAATG                          >  1:2685631/1‑49
                        GTGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAAAGATAATG                          <  1:458808/49‑1
                         TGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAAAGATAATGCA                        >  1:2209453/1‑50
                         TGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAAAGATAATGC                         <  1:634250/49‑1
                           CGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAAAGATAATGCATC                      >  1:2140051/1‑50
                            GAGGAAAAATATCAGTATGGCTCAATGATCCAAATAAAGATAATGCATC                      >  1:1800927/1‑49
                             AGGAAAAATATCAGTATGGCTCAATGATCCAAATAAAGATAATGCATC                      >  1:2316379/1‑48
                               GAAAAATATCAGTATGGCTCAATGATCCAAATAAAGATAATGCATCCAGC                  <  1:274636/50‑1
                                     TATCAGTATGGCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCA            >  1:571484/1‑50
                                      ATCAGTATGGCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCAT           >  1:2700735/1‑50
                                         AGTATGGCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTC         <  1:1007628/49‑1
                                          GTATGGCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTT       <  1:2335728/50‑1
                                          GTATGGCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTT       <  1:932256/50‑1
                                           TATGGCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTA      <  1:1355831/50‑1
                                             TGGCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTAGC    <  1:1055637/50‑1
                                               GCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTAGCAG  >  1:9485/1‑50
                                               GCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTAGCAG  >  1:1045005/1‑50
                                               GCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTAGCA   <  1:1646770/49‑1
                                                                                                 
GACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3761944‑3761992
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTAGCAG  >  NC_000913/3770243‑3770290

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.