Skip to content
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

ValueError: Unable to infer reference genome #192

Open
Shuhui-d opened this issue Dec 29, 2020 · 1 comment
Open

ValueError: Unable to infer reference genome #192

Shuhui-d opened this issue Dec 29, 2020 · 1 comment

Comments

@Shuhui-d
Copy link

Shuhui-d commented Dec 29, 2020

When I run the code as below

vaxrank
--vcf /mydata.vcf
--bam /mydata.bam
--vaccine-peptide-length 25
--mhc-predictor netmhc
--mhc-alleles H2-Kb,H2-Db
--padding-around-mutation 5
--output-ascii-report vaccine-peptides.txt
--output-pdf-report vaccine-peptides.pdf
I got the Error saying:
ValueError: Unable to infer reference genome for /mydata.vcf

mydata was the somatic variants called from hg19, do I run pyensembl install --release 75 --species human before this.
I used the physical computer with Win10 system and installed VMwareWorkstation (Centos8).
Should the vaxrank work like the codes above? Did I do any thing wrong?
Looking forward to your reply. Thank you for helping!

@leldershaw
Copy link

I had this issue too, but I looked at vaxrank --help which suggests adding --genome hg19 to your command.

--genome GENOME
What reference assembly your variant coordinates are using. Examples: 'hg19', 'GRCh38', or 'mm9'. This argument is ignored for MAF files, since each row includes the reference. For VCF files, this is used if specified, and otherwise is guessed from the header. For variants specfied on the commandline with --variant, this option is required.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants