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 30 [0] [0] 28 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)27 (1.010) 24/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

CGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAA  >  NC_000913/3761943‑3762041
                                                 |                                                 
cgacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC                                                   <  1:1086938‑M1/50‑1 (MQ=255)
cgacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC                                                   <  1:896642‑M1/50‑1 (MQ=255)
   cgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtca                                                <  1:549252‑M1/50‑4 (MQ=255)
     aGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaat                                              >  1:1294400‑M1/1‑45 (MQ=255)
      gTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatg                                             <  1:2059098‑M1/50‑7 (MQ=255)
       tttCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatga                                            <  1:557987‑M1/50‑8 (MQ=255)
          cATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatc                                          <  1:2206528‑M1/49‑10 (MQ=255)
              ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaa                                      >  1:2063321‑M1/1‑36 (MQ=255)
                   tCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                                <  1:1794103‑M1/50‑20 (MQ=255)
                    ccTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                                >  1:810753‑M1/1‑30 (MQ=255)
                       ttGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                            >  1:1285583‑M1/1‑27 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatg                           >  1:1581449‑M1/1‑26 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                            <  1:1182548‑M1/49‑24 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataa                             <  1:1899321‑M1/48‑23 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataa                             <  1:506609‑M1/48‑23 (MQ=255)
                         gtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatg                           <  1:543276‑M1/49‑25 (MQ=255)
                         gtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatg                           >  1:882813‑M1/1‑25 (MQ=255)
                          tgCGAGGAAAAATATCATTATGGCtcaatgatccaaataaagataatgc                          <  1:2006149‑M1/49‑26 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                       >  1:921926‑M1/1‑22 (MQ=255)
                             gAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcc                      <  1:2240814‑M1/50‑30 (MQ=255)
                             gAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                       >  1:2149361‑M1/1‑21 (MQ=255)
                                gAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccagc                   <  1:483139‑M1/50‑33 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                 >  1:1046361‑M1/1‑16 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                 >  1:73874‑M1/1‑16 (MQ=255)
                                   aaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctgg                >  1:939073‑M1/1‑15 (MQ=255)
                                            tATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         >  1:1013834‑M1/1‑6 (MQ=255)
                                              tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag      >  1:2045536‑M1/1‑4 (MQ=255)
                                              tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag      >  1:1903057‑M1/1‑4 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagca    <  1:90885‑M1/49‑48 (MQ=255)
                                                 ctcaatgatccaaataaagataatgcatccagctggtcattcttagcagc  <  1:323495‑M1/50‑50 (MQ=255)
                                                 |                                                 
CGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAA  >  NC_000913/3761943‑3762041

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: