vcfisec.o: vcfisec.c $(htslib_vcf_h) $(htslib_synced_bcf_reader_h) $(htslib_vcfutils_h) $(bcftools_h) $(filter_h) vcfisec.o: vcfisec.c $(htslib_vcf_h) $(htslib_synced_bcf_reader_h) $(htslib_vcfutils_h) $(htslib_hts_os_h) $(bcftools_h) $(filter_h) vcfmerge.o: vcfmer...
bcftools isec 示例:# 创建两个文件的交集和补集将其保存在 dir文件夹<br>bcftools isec -p dir A.vcf.gz B.vcf.gz<br> # 对A文件(INFO/MAF>=0.01)和B文件(INFO/dbSNP)进行过滤,但不对C文件进行过滤;在 dir 文件夹生成 3 个文件的交集位点,这些位点至少在 2 个输入文件中出现 bcftools isec -e'...
$ bcftools isec a.vcf.gz b.vcf.gz-p dir 默认参数就是取交集,更多高级用法请参考帮助文档。 5. stats stats命令用于统计VCF文件的基本信息。比如,突变个数、突变类型的个数、转换颠换个数、测序深度、Indel长度。还可以利用plot-vcfstats进行可视化处理。用法如下 ...
convert convert VCF/BCF files to different formats and back isec intersections of VCF/BCF files merge merge VCF/BCF files files from non-overlapping sample sets norm left-align and normalize indels plugin user-defined plugins query transform VCF/BCF into user-defined formats reheader modify VCF/BC...
# First,compress the VCF using bgzip,then index the gzipped VCF >bgzip first.vcf >tabix -p vcf first.vcf.gz >bcftools isec first.vcf.gz second.vcf.gz -p folder 0000.vcf # records private to first.vcf.gz 0001.vcf # records private to second.vcf.gz ...
-p Combine @PG headers with colliding IDs [alter IDs to be distinct] -s VALUE Override random seed -b FILE List of input BAM filenames,one per line [null] --input-fmt-option OPT[=VAL] Specify a single input file format option in the form ...