-
Notifications
You must be signed in to change notification settings - Fork 1
Description
I was unable to input the reference genome from local in despite of using fasta or fna.gz or xip file. Once complains about the format, then name, then, matching names etc ...
Could you please indicate which is the format of the input file with the reference genome and how to input it from config.yaml or from samples.csv. My reference genome input file is zip and contains two files, fna and gbff files.
Thanks in adavance,
Joan Pons
Error in rule index_reference:
jobid: 6
input: results/GCA_037176745.1/data/genome/GCA_037176745.1.fna
output: results/GCA_037176745.1/data/genome/GCA_037176745.1.fna.sa, results/GCA_037176745.1/data/genome/GCA_037176745.1.fna.pac, results/GCA_037176745.1/data/genome/GCA_037176745.1.fna.bwt, results/GCA_037176745.1/data/genome/GCA_037176745.1.fna.ann, results/GCA_037176745.1/data/genome/GCA_037176745.1.fna.amb, results/GCA_037176745.1/data/genome/GCA_037176745.1.fna.fai, results/GCA_037176745.1/data/genome/GCA_037176745.1.dict
log: logs/GCA_037176745.1/index_ref/log.txt (check log file(s) for error details)
conda-env: /grass/jpp/posigenome_project01/.snakemake/conda/f2c4e92d92a9365985051436885c333f_
shell:
bwa index results/GCA_037176745.1/data/genome/GCA_037176745.1.fna 2> logs/GCA_037176745.1/index_ref/log.txt
samtools faidx results/GCA_037176745.1/data/genome/GCA_037176745.1.fna --output results/GCA_037176745.1/data/genome/GCA_037176745.1.fna.fai >> logs/GCA_037176745.1/index_ref/log.txt
samtools dict results/GCA_037176745.1/data/genome/GCA_037176745.1.fna -o results/GCA_037176745.1/data/genome/GCA_037176745.1.dict >> logs/GCA_037176745.1/index_ref/log.txt 2>&1
(one of the commands exited with non-zero exit code; note that snakemake uses bash strict mode!)