Skip to content

KevinZ的小窝

Menu
  • Home
  • Categories
Menu

使用Delly进行结构变异分析

Posted on 2026年 2月 9日 by KevinZhou

github官网:https://github.com/dellytools/delly
hg38.excl 下载:https://github.com/dellytools/delly/tree/1ac07d0420e8a111076949ac8aab6361eb625e60/excludeTemplates

Germline short-read SV calling

# SV calling is done by sample for high-coverage genomes or in small batches for low-coverage genomes
delly call -g hg38.fa -o s1.bcf -x hg38.excl sample1.bam

# Merge SV sites into a unified site list
delly merge -o sites.bcf s1.bcf s2.bcf ... sN.bcf

# Genotype this merged SV site list across all samples. This can be run in parallel for each sample.
delly call -g hg38.fa -v sites.bcf -o s1.geno.bcf -x hg38.excl s1.bam
delly call -g hg38.fa -v sites.bcf -o sN.geno.bcf -x hg38.excl sN.bam

# Merge all genotyped samples to get a single VCF/BCF using bcftools merge
bcftools merge -m id -O b -o merged.bcf s1.geno.bcf s2.geno.bcf ... sN.geno.bcf

# Apply the germline SV filter which requires at least 20 unrelated samples
delly filter -f germline -o germline.bcf merged.bcf

# 构建批量运行代码

cnv calling

delly cnv -g example/ref.fa -m example/map.fa.gz -c out.cov.gz -o cnv.bcf example/sr.bam
2026 年 2 月
一 二 三 四 五 六 日
 1
2345678
9101112131415
16171819202122
232425262728  
« 12 月    

俺家的猫~

胖达~

© 2026 KevinZ的小窝 |

粤ICP备2023017690号

|

粤公网安备 44010402003004号