Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

VCF file reading problem #9

Open
Caseyb4ob opened this issue Feb 9, 2022 · 1 comment
Open

VCF file reading problem #9

Caseyb4ob opened this issue Feb 9, 2022 · 1 comment

Comments

@Caseyb4ob
Copy link

Hi,

I am trying to use DNAscan, I do not have access to run sudo commands but I was able to accomplish do everything in these commands without access for sudo (except for tff-dejavu).

I am getting the folowing error:
"WARNING: Fastq format input data is required if you want to perform the alignment stage
WARNING: Using input vcf as variant file. Do not provide vcf file if you wish to perform variant calling
Hard filtering of SNV/indel variants is being performed...
bcftools filter -i " QUAL > 1 & QUAL / INFO/AO > 10 & SAF > 0 & SAR > 0 & RPR > 1 & RPL > 1 " All_data/ALS57_CKDN200006030-1A_HLTMNDSXY_L3_filtered_snps.vcf | bgzip -c > All_data/DNAscan/SETUP/DNAscanOutput_snps/ALS57_CKDN200006030-1A_HLTMNDSXY_L3/sample_sorted_filtered.vcf.gz ; tabix -fp vcf All_data/DNAscan/SETUP/DNAscanOutput_snps/ALS57_CKDN200006030-1A_HLTMNDSXY_L3/sample_sorted_filtered.vcf.gz
WARNING: All_data/DNAscan/SETUP/DNAscanOutput_snps/ALS57_CKDN200006030-1A_HLTMNDSXY_L3/sample_sorted_filtered.vcf.gz only contains the header and no data, therefore DNAscan will now terminate."

The commands I used to call are:
"DNAscan.py -in $VCF/${names[${SLURM_ARRAY_TASK_ID}]}_filtered_InDel.vcf -format vcf -reference hg38 -alsgenescanner -out $OUT
DNAscan.py -in $VCF/${names[${SLURM_ARRAY_TASK_ID}]}_filtered_snps.vcf -format vcf -reference hg38 -alsgenescanner -out $OUT2"

I set up Annovar and all other dependecies to my knowledge. Do you have any troubleshooting suggestions?
Also, I checked the path to the vcfs and it is correct, the vcf is not empty, but the file that gets copied into the output folder is.

@AlfredoKCL
Copy link
Member

Hello,
we are in the process of updating DNAscan. I think it is best if you try the new version directly https://github.com/hevmarriott/DNAscanv2 . Let me know if you prefer to receive support for the older version instead
BW
Alfredo

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants