acf2betascan empty output about glactools HOT 6 CLOSED

Hi josieparis,
Is the ancestral state defined?

Gabriel

from glactools.

Hey Gabriel,

Thanks for getting back to me so quick. I have the ancestral state defined in the vcf file like so:
chr1 4620 chr1_4620 A G 5527.83 PASS AA=A
chr1 4827 chr1_4827 G A 38420.2 PASS AA=G
chr1 4947 chr1_4947 T G 6433.52 PASS AA=T
chr1 5182 chr1_5182 A C 5021.24 PASS AA=A

It also does not produce an output when I use --useroot

from glactools.

Ok the conversion to ACF cannot process custom tags for now. You could add it as a feature request.

For now you could either:

1. if this is hg18, you could use our per-parsed EPO files during acf conversion
2. otherwise, get a VCF of the ancestral state and use replaceanc+usepopsrootanc to use it as ancestral
3. or specify --fold but then you lose which allele is ancestral and derived, you just get minor/major.

from glactools.

ah gotya!

yeah unfortunately these are custom vcfs. I'll add as a feature request

in the meantime, I'll probably just fold them.
thanks for your help! (and glactools - very useful)

from glactools.

from glactools.

Ok for now, may I close this issue? tor. 16. maj 2019 16.22 skrev josieparis [email protected]:
…
ah gotya! yeah unfortunately these are custom vcfs. I'll add as a feature request in the meantime, I'll probably just fold them. thanks for your help! (and glactools - very useful) — You are receiving this because you commented. Reply to this email directly, view it on GitHub <#6?email_source=notifications&email_token=AAQRNI6WPGDDCESLHIGZI4LPVVNXXA5CNFSM4HNFH3KKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGODVRZHAY#issuecomment-493065091>, or mute the thread https://github.com/notifications/unsubscribe-auth/AAQRNI3BKQZ4FB2OM6JL23TPVVNXXANCNFSM4HNFH3KA .

hello grenaud! your reply was helpful. now, I just have outgroup vcf file (other mammals not human), could you tell me how can I transform vcf file to epo file?

from glactools.