freebayes not reading the chunks of regions file
0
0
Entering edit mode
3.6 years ago
evafinegan • 0

Hello,

I am using freebayes for SNP calling using:

RUN=${SLURM_ARRAY_TASK_ID}
INPUT_DIR=/path/
OUTPUT_DIR=/path/
INPUT=$(ls $INPUT_DIR/*.txt | sed -n ${RUN}p)
SAMPLE=$(basename "$INPUT" .txt)

freebayes  -f genome.fasta \
          -L list_of_bam.txt \
          --targets  ${INPUT_DIR}/myfile${RUN}.txt > ${OUTPUT_DIR}/${SAMPLE}_.vcf

However, it gives error:

RUN #1 with sample myfile00
bed targets file is not open

I think it is not reading the subset files (total 50) that I made using split. Can you please help me figure out what is my mistake in the array. Thank you!

snp • 1000 views
ADD COMMENT
0
Entering edit mode

Where are you specifying the region?

ADD REPLY

Login before adding your answer.

Traffic: 2438 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6