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 31 [0] [0] 32 [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)31 (1.060) 24/92 0.1 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

TAATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCAT  >  NC_000913/325364‑325458
                                               |                                               
tctaccagaaaattcacggcatcggctagggtgataccggcatcgata                                                 <  1:570442‑M2/1‑1 (MQ=255)
  taccagaaaattcacggcatcggctagggtgataccggcatcgatAGTc                                              >  1:2023856‑M2/46‑49 (MQ=255)
    ccagaaaattcacggcatcggctagggtgataccggcatcgatAGt                                               <  1:376646‑M2/3‑1 (MQ=255)
    ccagaaaattcacggcatcggctagggtgataccggcatcgatAGt                                               <  1:281946‑M2/3‑1 (MQ=255)
    ccagaaaattcacggcatcggctagggtgataccggcatcgatAGt                                               <  1:857484‑M2/3‑1 (MQ=255)
     cagaaaattcacggcatcggctagggtgataccggcatcgatAGTCTg                                            >  1:1334359‑M2/43‑48 (MQ=255)
      agaaaattcacggcatcggctagggtgataccggcatcgatAGTCTGTa                                          >  1:361080‑M2/42‑49 (MQ=255)
      agaaaattcacggcatcggctagggtgataccggcatcgatAGTCTGTa                                          >  1:184628‑M2/42‑49 (MQ=255)
         aaattcacggcatcggctagggtgataccggcatcgatAGTCTGTAt                                         >  1:1955439‑M2/39‑47 (MQ=255)
         aaattcacggcatcggctagggtgataccggcatcgatAGTCTGTAt                                         >  1:782753‑M2/39‑47 (MQ=255)
          aattcacggcatcggctagggtgataccggcatcgatAGTCTGTATCAg                                      <  1:1893618‑M2/12‑1 (MQ=255)
              cacggcatcggctagggtgataccggcatcgatAGTCTGTATCAGGAAt                                  >  1:1800483‑M2/34‑49 (MQ=255)
                cggcatcggctagggtgataccggcatcgatAGTCTGTATCAGGAATGt                                <  1:397478‑M2/18‑1 (MQ=255)
                cggcatcggctagggtgataccggcatcgatAGTCTGTATCAGGAATGt                                <  1:705937‑M2/18‑1 (MQ=255)
                  gcatcggctagggtgataccggcatcgatAGTCTGTATCAGGAATGTTc                              <  1:741986‑M2/20‑1 (MQ=255)
                    atcggctagggtgataccggcatcgatAGTCTGTATCAGGAATGTTCgg                            <  1:1257103‑M2/22‑1 (MQ=255)
                    atcggctagggtgataccggcatcgatAGTCTGTATCAGGAATGTTCgg                            >  1:1385022‑M2/28‑49 (MQ=255)
                        gctagggtgataccggcatcgatAGTCTGTATCAGGAATGTTCGGGtt                         >  1:2292222‑M2/24‑48 (MQ=255)
                        gctagggtgataccggcatcgatAGTCTGTATCAGGAATGTTCGGGTTa                        <  1:1340560‑M2/26‑1 (MQ=255)
                         ctagggtgataccggcatcgatAGTCTGTATCAGGAATGTTCGGGTTaa                       <  1:1796226‑M2/27‑1 (MQ=255)
                           agggtgataccggcatcgatAGTCTGTATCAGGAATGTTCGGGTTaaa                      >  1:834330‑M2/21‑48 (MQ=255)
                            gggtgataccggcatcgatAGTCTGTATCAGGAATGTTCGGGTTAAata                    <  1:1331095‑M2/30‑1 (MQ=255)
                             ggtgataccggcatcgatAGTCTGTATCAGGAATGTTCGGGTTAAatat                   <  1:766916‑M2/31‑1 (MQ=255)
                              gtgataccggcatcgatAGTCTGTATCAGGAATGTTCGGGTTAAatat                   >  1:954100‑M2/18‑48 (MQ=255)
                              gtgataccggcatcgatAGTCTGTATCAGGAATGTTCGGGTTAAatat                   >  1:1861743‑M2/18‑48 (MQ=255)
                                gataccggcatcgatAGTCTGTATCAGGAATGTTCGGGTTAAATATCAg                >  1:62060‑M2/16‑49 (MQ=255)
                                    ccggcatcgatAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCaaa            <  1:44542‑M2/38‑1 (MQ=255)
                                    ccggcatcgatAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCaaa            <  1:1623370‑M2/38‑1 (MQ=255)
                                    ccggcatcgatAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCaa             <  1:1988409‑M2/37‑1 (MQ=255)
                                           cgatAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCg     >  1:231009‑M2/5‑49 (MQ=255)
                                            gatAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCg     <  1:1453281‑M2/45‑1 (MQ=255)
                                               aGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGcat  >  1:33398/1‑48 (MQ=255)
                                               |                                               
TAATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCAT  >  NC_000913/325364‑325458

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