-
Notifications
You must be signed in to change notification settings - Fork 216
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
VCF2MAF and allele frequency #344
Comments
anyone? |
In my VCF, the AF field is in the FORMAT column, so I could get VAF with |
@tanghaibao AF in the FORMAT column represents the Back to the original issue, I am also missing GnomAD AF in my MAF file. Would be interested if anyone has figured out the issue and a solution. |
Hi @berguner I have been using the Nextflow pipeline here (dev branch) https://github.com/FriederikeHanssen/vcftomaf/tree/dev To do the conversion it automatically includes all VEP annotations present in the info column |
I am converting a VCF to MAF using the VCF2MAF program, and importing that into cbioportal. I notice that the "af" from the INFO field is not making to the MAF output file.
Snippet from VCF: (the "af" value is no where to be found in maf file)
##INFO=<ID=af,Number=1,Type=Float,Description="Alternate allele frequency">
#CHROM POS ID REF ALT QUAL FILTER INFO
chr17 27419372 . C A . PASS af=0.31410256;cds_syntax=5176G>T;cosmic_status=UNKNOWN;depth=468;effect=MISSENSE;gene_name=MYO18A;protein_syntax=A1726S;transcript_name=NM_078471
Perl VCF2MAF command:
perl vcf2maf.pl --input-vcf /home/jmf432/test_vcf_files/ORD-1543927-02.vcf --output-maf /home/jmf432/test_vcf_files/ORD-1543927-02-JF2.maf --ref-fasta /home/jmf432/.vep/homo_sapiens/109_GRCh38/Homo_sapiens.GRCh38.dna.toplevel.fa.gz --vep-path /home/jmf432/vep/ensembl-vep-release-109.3 --ncbi-build GRCh37 --verbose --vep-overwrite --any-allele --retain-info af
Snippet from the output MAF file (all the af fields are blank):
DOMAINS AF AFR_AF AMR_AF ASN_AF EAS_AF EUR_AF SAS_AF AA_AF EA_AF
I was expecting to see the "af" value from the VCF file in the MAF output file, but it is not.
Can someone please explain why and correct me if I have a misunderstanding ?****
The text was updated successfully, but these errors were encountered: