nf-core / nf-core/rnavar

FILTER column empty and no annotation output

Open
#236 0 comments 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

bug
Dominant language
Nextflow
Stars
61
Forks
47
Avg merge
1d 1h
Merged PRs (30d)
1

Description

Description of the bug

Hi,

thanks for the code!

I am running Nextflow 25.04.6 and rnavar 1.2.0. The job completes, and in the results dir I see
pipeline_info
preprocessing
reports
samtools
seq2hla
variant_calling

The variant_calling dir reports a .haplotypecaller.combined.g.vcf.gz and .haplotypecaller.combined.g.vcf.gz.tbi for each sample

The FILTER column in the resulting VCFs is empty and the INFO column does not show any VEP annotation, as far as I can see.

Can you please let me know what I am missing ?

Thanks

Command used and terminal output
The command line is as follow
nextflow run nf-core/rnavar  -r 1.2.0 -c $inDIR/loc.config -profile apptainer \
    -work-dir $workDIR \
    --input $inDIR/samplesheet.csv \
    --outdir $outPATH/$outDIR \
    --fasta $GENOME \
    --fasta_fai $GENOMEIDX \
    --gtf $GTF_PATH \
    --read_length 101 \
    --aligner star \
    --star_twopass true \
    --star_index $STAR_INDEX \
    --remove_duplicates false \
    --no_intervals true \
    --generate_gvcf true \
    --known_indels $INDELS \
    --known_indels_tbi $INDELSTBI \
    --dbsnp $DBSNP \
    --dbsnp_tbi $DBSNPTBI \
    --snpeff_db "GRCh38.p13" \
    --tools seq2hla,bcfann,vep \
    --vep_cache $VEPCACHEPATH/.vep \
    --vep_genome "GRCh38" \
    --vep_species "homo_sapiens" \
    --vep_cache_version "110"
Relevant files

No response

System information

No response

Contributor guide

Open the contributing guide

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. Open a pull request that references the issue number.

Research direction

Reproduce the supplied Nextflow command and inspect the generated files under variant_calling, especially the combined VCF and its index. Trace the annotation and filtering steps for the listed tools, then confirm that the resulting VCF has populated FILTER and VEP INFO fields.

Written by the indexing model from the issue text.

Assessment

Domain
bioinformatics
Issue type
Bug
Difficulty
4/5
Estimated time
3-5 days
Activity status
Stale
Clarity
Needs clarification
Newbie friendliness
25/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.