Can't find the vcf_stats file（MergeVcfWideStats.merged_bed_file）, how to generate it?

[jingydz]

Bug Report

$ Rscript plot_sv_vcf_distribs.R -N $( cat 4388.nyuwa.after.change.lst | sort | uniq | wc -l ) -S SV_colors.txt nvwaCHBCHS.4388.genotype.vcffilter.vcf.gz.gz.stats plotQC_vcfwide_output/ Error in read.table(INFILE, comment.char = "", sep = "\t", header = T, : more columns than column names Execution halted

Affected module(s) or script(s)

plot_sv_vcf_distribs.R

Description

I want to display the length distribution of SV, and I found the specific code as follows:

Plot VCF-wide distributions

/opt/sv-pipeline/scripts/vcf_qc/plot_sv_vcf_distribs.R
-N $( cat ~{samples_list} | sort | uniq | wc -l )
-S /opt/sv-pipeline/scripts/vcf_qc/SV_colors.txt
~{vcf_stats}
plotQC_vcfwide_output/

Prep outputs

tar -czvf ~{prefix}.plotQC_vcfwide_output.tar.gz
plotQC_vcfwide_output

Question

But I'm curious what file the 'vcf_stats' should be? The source code shows 'vcf_stats=MergeVcfWideStats.merged_bed_file' How did you obtain this file?