Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 263,514 Δ61,897 bp [ykfI]ykgH 71 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 263514 325410 61897 47 [0] [0] 47 [ykfI]–ykgH [ykfI],yafW,ykfH,ykfG,yafX,ykfF,ykfB,yafY,yafZ,ykfA,perR,insN,insI‑1,insN,eyeA,ykfC,insH‑1,mmuP,mmuM,afuC,afuB,insB‑2,insA‑2,yagB,yagA,yagE,yagF,yagG,yagH,xynR,argF,ykgS,insB‑3,insA‑3,yagJ,yagK,yagL,yagM,ykgV,yagN,intF,yagP,paoD,paoC,paoB,paoA,yagU,ykgJ,ecpE,ecpD,ecpC,ecpB,ecpA,ecpR,ykgL,ykgO,ykgM,ykgR,ykgP,eaeH,insE‑1,insF‑1,ykgA,rclC,rclB,rclA,rclR,ykgE,ykgF,ykgG,ykgH

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2635130 (0.000)43 (1.210) 33/92 0.0 100% coding (157/342 nt) ykfI CP4‑6 prophage; toxin of the YkfI‑YafW toxin‑antitoxin system
?NC_000913 325411 = 0 (0.000)intergenic (‑47/+166) ykgH/betA uncharacterized protein YkgH/choline dehydrogenase

ATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCAT  >  NC_000913/325368‑325458
                                           |                                               
ccagaaaattcacggcatcggctagggtgataccggcatcgatAGTc                                              <  1:702741‑M2/4‑1 (MQ=255)
  agaaaattcacggcatcggctagggtgataccggcatcgatAGTCTGTa                                          <  1:1160678‑M2/8‑1 (MQ=255)
   gaaaattcacggcatcggctagggtgataccggcatcgatAGTCTGTAt                                         <  1:236850‑M2/9‑1 (MQ=255)
    aaaattcacggcatcggctagggtgataccggcatcgatAGTCTGTATc                                        <  1:2325396‑M2/10‑1 (MQ=255)
     aaattcacggcatcggctagggtgataccggcatcgatAGTCTGTAt                                         >  1:1315464‑M2/39‑47 (MQ=255)
     aaattcacggcatcggctagggtgataccggcatcgatAGTCTGTATCa                                       <  1:620249‑M2/11‑1 (MQ=255)
         tcacggcatcggctagggtgataccggcatcgatAGTCTGTATCAGGaa                                   >  1:69566‑M2/35‑49 (MQ=255)
           acggcatcggctagggtgataccggcatcgatAGTCTGTATCAGGAATg                                 >  1:2239997‑M2/33‑49 (MQ=255)
           acggcatcggctagggtgataccggcatcgatAGTCTGTATCAGGAATg                                 >  1:2174053‑M2/33‑49 (MQ=255)
            cggcatcggctagggtgataccggcatcgatAGTCTGTATCAGGAATGt                                <  1:561775‑M2/18‑1 (MQ=255)
             ggcatcggctagggtgataccggcatcgatAGTCTGTATCAGGAATGtt                               >  1:1264451‑M2/31‑49 (MQ=255)
             ggcatcggctagggtgataccggcatcgatAGTCTGTATCAGGAATGtt                               >  1:854034‑M2/31‑49 (MQ=255)
             ggcatcggctagggtgataccggcatcgatAGTCTGTATCAGGAATGtt                               <  1:15895‑M2/19‑1 (MQ=255)
              gcatcggctagggtgataccggcatcgatAGTCTGTATCAGGAATGTTc                              <  1:2712796‑M2/20‑1 (MQ=255)
              gcatcggctagggtgataccggcatcgatAGTCTGTATCAGGAATGTTc                              <  1:2119635‑M2/20‑1 (MQ=255)
                 tcggctagggtgataccggcatcgatAGTCTGTATCAGGAATGTtcgag                           <  1:1784254‑M2/23‑3 (MQ=255)
                 tcggctagggtgataccggcatcgatAGTCTGTATCAGGAATGTTCggg                           <  1:1784269‑M2/23‑1 (MQ=255)
                  cggctagggtgataccggcatcgatAGTCTGTATCAGGAATGTTCGGGt                          >  1:1536952‑M2/26‑49 (MQ=255)
                   ggctagggtgataccggcatcgatAGTCTGTATCAGGAATGTTCGGGt                          >  1:2566991‑M2/25‑48 (MQ=255)
                    gctagggtgataccggcatcgatAGTCTGTATCAGGAATGTTCGGGtt                         >  1:367569‑M2/24‑48 (MQ=255)
                    gctagggtgataccggcatcgatAGTCTGTATCAGGAATGTTCGGGTTa                        <  1:75624‑M2/26‑1 (MQ=255)
                     ctagggtgataccggcatcgatAGTCTGTATCAGGAATGTTCGGGTTaa                       <  1:1498687‑M2/27‑1 (MQ=255)
                     ctagggtgataccggcatcgatAGTCTGTATCAGGAATGTTCGGGTTa                        >  1:2499598‑M2/23‑48 (MQ=255)
                      tagggtgataccggcatcgatAGTCTGTATCAGGAATGTTCGGGTTaaa                      >  1:279777‑M2/22‑49 (MQ=255)
                        gggtgataccggcatcgatAGTCTGTATCAGGAATGTTCGGGTTAAata                    <  1:1942538‑M2/30‑1 (MQ=255)
                        gggtgataccggcatcgatAGTCTGTATCAGGAATGTTCGGGTTAAata                    <  1:1762866‑M2/30‑1 (MQ=255)
                        gggtgataccggcatcgatAGTCTGTATCAGGAATGTTCGGGTTAAat                     <  1:2470758‑M2/29‑1 (MQ=255)
                          gtgataccggcatcgatAGTCTGTATCAGGAATGTTAGGGGGAAATATc                  <  1:789074‑M2/32‑1 (MQ=255)
                           tgataccggcatcgatAGTCTGTATCAGGAATGTTCGGGTTAAATATCa                 <  1:1485768‑M2/33‑1 (MQ=255)
                             ataccggcatcgatAGTCTGTATCAGGAATGTTCGGGTTAAATATCAg                >  1:1316193‑M2/15‑48 (MQ=255)
                                 cggcatcgatAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCaaaa           <  1:881266‑M2/39‑1 (MQ=255)
                                  ggcatcgatAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCaaaaa          <  1:2844512‑M2/40‑1 (MQ=255)
                                  ggcatcgatAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCaaaaa          <  1:347269‑M2/40‑1 (MQ=255)
                                  ggcatcgatAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCaaa            >  1:561542‑M2/10‑47 (MQ=255)
                                   gcatcgatAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAg         <  1:943895‑M2/41‑1 (MQ=255)
                                     atcgatAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGcc       <  1:1052679‑M2/43‑1 (MQ=255)
                                     atcgatAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGcc       <  1:552141‑M2/43‑1 (MQ=255)
                                     atcgatAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGcc       <  1:56061‑M2/43‑1 (MQ=255)
                                      tcgatAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGccc      <  1:2232345‑M2/44‑1 (MQ=255)
                                      tcgatAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGccc      <  1:403519‑M2/44‑1 (MQ=255)
                                      tcgatAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGccc      <  1:1392729‑M2/44‑1 (MQ=255)
                                        gatAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCg     >  1:2768778‑M2/4‑48 (MQ=255)
                                        gatAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCg     <  1:2728335‑M2/45‑1 (MQ=255)
                                          tAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCAt  <  1:1762774‑M2/48‑1 (MQ=255)
                                          tAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCAt  >  1:586614‑M2/2‑49 (MQ=255)
                                          tAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCAt  >  1:2138489‑M2/2‑49 (MQ=255)
                                           aGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGcat  >  1:2039296/1‑48 (MQ=255)
                                           |                                               
ATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCAT  >  NC_000913/325368‑325458

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.