Entering edit mode
16 months ago
estpallag
▴
10
Hello, I need some help with SnpSift and fitler option. I am trying to run different filters but in every case the output is empty. No warnings, of fail run. all the works are completed and the .log file does not show errors. However, the .vcf file does not have any result.
This is an example of the script that I am trying:
#Load modules
module load apps/anaconda-4.7.12.tcl
#Activate conda enviroment for this analysis
eval "$(conda shell.bash hook)"
conda activate SnpSift
unset PYTHONPATH
#Filter by Chr17 region
java -jar /storage/users/epalominolago/Software/snpEff/SnpSift.jar filter "( CHROM = 'chr17') | ( POS > 63010000 ) & ( POS < 68200000 )" > SnpEff_Chr18.vcf
And this is the results:
> SnpSiftVersion="SnpSift 5.1d (build 2022-04-19 15:50), by Pablo Cingolani"
> SnpSiftCmd="SnpSift Filter '( CHROM = 'chr17') | ( POS > 63010000 ) & ( POS < 68200000 )'"
> FILTER=<ID=SnpSift,Description="SnpSift 5.1d (build 2022-04-19 15:50), by Pablo Cingolani, Expression used: ( CHROM = 'chr17') | ( POS > 63010000 ) & ( POS < 68200000 )"
as you can see it is empty wit any variant. I know I should have results because running it using VEP gave me the results.
any idea of what I am missing? thanks