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

Support non-human variant analysis? #124

Open
yige-luo opened this issue Jul 17, 2021 · 0 comments
Open

Support non-human variant analysis? #124

yige-luo opened this issue Jul 17, 2021 · 0 comments

Comments

@yige-luo
Copy link

Dear Xiaowei,

I want to use rvtests to analyze variants data in a non-human data, whose CHROM field in the .vcf file is labeled in a different way.

 #CHROM  POS     ID      REF     ALT     QUAL    FILTER  INFO    FORMAT  LC000001_14_I_LC000001_14_I     LC000001_17_F_LC000001_17_F     LC000004_01_I_LC00000
Scaffold_2      78417   Scaffold_2:78417        ATATTCATTG      A       .       .       PR      GT      0/0     0/0     0/0     0/0     0/1     0/1     0/0
Scaffold_2      80513   Scaffold_2:80513        TAAAAGAC        T       .       .       PR      GT      0/0     0/0     0/0     0/1     0/0     0/0     0/0
Scaffold_4      65529   Scaffold_4:65529        C       A       .       .       PR      GT      0/0     0/0     0/0     0/0     0/0     0/0     0/1     0/0
Scaffold_4      71567   Scaffold_4:71567        C       T       .       .       PR      GT      0/0     0/0     0/0     0/0     0/0     0/0     0/0     0/0

I tried to use vcf2kinship function to get the kinship matrix, but kept encountering this error

vcf2kinship --inVcf my.vcf.gz --bn --out output

Total [ 2330000 ] VCF records ( 0 a[INFO]al/X-PATotal [ 2413517 ] VCF records have been processed.
[ERROR] There are not enough variants to create kinship matrix.

I guess it may have something to do with the naming of the CHROM field, since none of my vcf records are deemed as "autosomal" by the program. I'm wondering if rvtests currently supports non-human variant analysis?

Sincerely,
Yige

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

1 participant