Movatterモバイル変換


[0]ホーム

URL:


Skip to content

Navigation Menu

Search code, repositories, users, issues, pull requests...

Provide feedback

We read every piece of feedback, and take your input very seriously.

Saved searches

Use saved searches to filter your results more quickly

Sign up

traying to make a maf from an already annotated vcf #354

Open
@javierAPC

Description

@javierAPC

Hi im traying to get some visual summaries from a vcf i annotated with VEP, and for that i want to use maftools but firts i have to get a maf file, so i use this comand:perl vcf2maf.pl --input-vcf /home/victor/ensembl-vep/test_ann.vcf --output-maf tests/test.maf --tumor-id TUMOUR --normal-id NORMAL --inhibit-vep.
which gives me this errors

[.......][W::fai_get_val] Reference chr10:101362827-101362833 not found in FASTA file, returning empty sequence[faidx] Failed to fetch sequence in chr10:101362827-101362833ERROR: You're either using an outdated samtools, or --ref-fasta is not the same genome build as your --input-vcf. at vcf2maf.pl line 410.

the thing is that for VEP i use the cache option (it didn't let me connect to the database) and the GRCh38 version of the human genome i have is in.vep/homo_sapiens/111_GRCh38, so i changed the vcf2maf.pl to point at that directory where i store Homo_sapiens.GRCh38.dna.toplevel.fa.gz.
The vcf where also align to GRCh38.

Also before the test commands give this error
ERROR: Provided --ref-fasta is missing or empty: /home/victor/.vep/homo_sapiens/102_GRCh37
but with the new reference give me this error:

[.......]WARNING: Reference allele G at 7:116412043 doesn't match A (flanking bps: CAC) from reference FASTA: /home/victor/.vep/homo_sapiens/111_GRCh38/Homo_sapiens.GRCh38.dna.toplevel.fa.gzSTATUS: Running VEP and writing to: tests/test.vep.vcfERROR: Cannot find VEP script under: /home/victor/miniconda3/bin

Uptade:
cheking the reference in the header of the vcf, teh path isreference=/nfs-1-c1/work/d1/0034c2fad733d8de3f89c93ed42117/reference_files/genome.fa, again it should also be GRCh38

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions


      [8]ページ先頭

      ©2009-2025 Movatter.jp