epo file format about glactools HOT 27 OPEN

from glactools.

Hello,
If I want to calculate β2 using an outgroup to call substitutions and ancestral alleles, can I just get epo files from outgroup vcf file, without common ancestor?
Thank you very much

from glactools.

Yes! You need to convert VCF to ACF and run replaceancestor I believe. Let me know if that works!

from glactools.

Thank you. I will try it !
If I use vcfm2acf to convert VCF to ACF , then I use usepopsrootanc to tell the root(outgroup), finally I use the replaceanc, the resulting file is the epo file that I need in betascan2 software?
Does this epo file include all the individuals I calculate later?

from glactools.

Hello, I have another question.
Is the file2.acf.gz in https://github.com/grenaud/glactools the result using the glactools with vcfm2acf --onlyGT --fai parameter?
I am not sure if using these parameters is correct?

Thanks

from glactools.

Hello,
the file1.acf.gz only have a outgroup individual:

when I use the "usepopsrootanc" parameter, I only tell the outgroup, then the error is as follows

I think it lack ancestral information, it is right? so the file1.acf.gz needs the root and ancestor information?

from glactools.

Hello,
What does "ChimpHumanAncestor " refer to specifically?

Thank you very much

from glactools.

Thank you. I will try it ! If I use vcfm2acf to convert VCF to ACF , then I use usepopsrootanc to tell the root(outgroup), finally I use the replaceanc, the resulting file is the epo file that I need in betascan2 software? Does this epo file include all the individuals I calculate later?

Would you describe seems like it will work. The EPO is only human+chimp+other great apes and the inferred ancestral states, there is a single human "individual" and that is the reference

from glactools.

Hello, I have another question. Is the file2.acf.gz in https://github.com/grenaud/glactools the result using the glactools with vcfm2acf --onlyGT --fai parameter? I am not sure if using these parameters is correct?

Thanks

I am not sure how this was generated but acf files can be generated in very different ways.

from glactools.

There is root and ancestor. In the case of humans, we can think of the chimps as the "root" as all humans are equitistant to the chimp. However there is a chimp-human ancestor, we do not have access to the sequence of its genome as it died out 6 million years ago. However we can make a probabilistic inference of what this genome had in terms of alleles. This is precisely what EPO does.

from glactools.

Hello,
Thank you for your reply!

1. May I ask if I have only root and ancestor, but no groups like other apes, is it OK to generate the follow file.acf.gz as my epo file for Betascan2？
   

2. After I get the epo file, I get an error adding the index file? How to create this index file?
   

3. The first line of the epo file seems to be different from the example you showed. Is it correct for my epo file?
   
   Thanks

from glactools.

How did you generate this file? Could you paste the command lines?

from glactools.

I think I misunderstood your meaning.
Is file4.acf.gz in your example (glactools replaceanc file2.acf.gz file3.acf.gz > file4.acf.gz) equal to the 1kg_chr6.AFRsubset.acf.gz in following code ( glactools vcfm2acf --onlyGT --epo all.epo.gz --fai human_g1k_v37.fasta.fai chr6.AFRsubset.recode.vcf > 1kg_chr6.AFRsubset.acf.gz)?
These two files are actually equivalent, right?

from glactools.

file4.acf.gz is just a placeholder name, it could be anything.

from glactools.

sorry，

1. glactools replaceanc file2.acf.gz <(glactools usepopsrootanc -u file1.acf.gz chimp ChimpHumanAncestor ) > file4.acf.gz
2. glactools vcfm2acf --onlyGT --epo all.epo.gz --fai human_g1k_v37.fasta.fai chr6.AFRsubset.recode.vcf > 1kg_chr6.AFRsubset.acf.gz

What I understand is that if there is no epo file, file4 generated according to step 1 and 1kg_chr6.AFRsubset.acf.gz generated using epo in 2 are both acf files containing outgroup information. steps 1 and 2 can generate the same file, and the steps are equivalent, is it right?

from glactools.

While I encourage people to use file descriptors, they are not great for debugging because you don't know which one fails.

Could you try to break down these commands into three separate ones, produce the output for each and inspect that using glactools view?

from glactools.

Yes, I used 3 commands to generate the file successfully.

Since my species is not human, I want to know if the output file of 1 is the equivalent of the output file of 2 ?

Thank you very much

from glactools.

Could you please paste the commands? I am not sure what is file 1 and file 2 here.

from glactools.

I mean these two ways above. Since my object is not human, there is no epo file, so I can only use the method of 1. I am not sure whether the file I generated is the input file of betascan2? the following is my output

from glactools.

Again, please no file descriptors, just single commands.

from glactools.

The following is code of my test file.
IN1 is outgroup, ZA1 is ancestor.
test.chr24.vcf is my pop1 file.

glactools vcfm2acf --onlyGT --fai $ref_genome outgroup_root_test.chr24.vcf > outgroup_root_test.chr24.acf.gz
glactools vcfm2acf --onlyGT --fai $ref_genome test.chr24.vcf > test.chr24.acf.gz
glactools usepopsrootanc outgroup_root_test.chr24.acf.gz IN1 ZA1 > outgroup_root_test_identify.chr24.acf.gz
glactools replaceanc test.chr24.acf.gz outgroup_root_test_identify.chr24.acf.gz > outgroup_root_test_merge.chr24.acf.gz

from glactools.

ok thank you! which command fails?

from glactools.

No error reported.
I'm just not sure if the outgroup_root_test_merge.chr24.acf.gz is the betascan2 input file?

from glactools.

ah no it is not! you need to use glactools acf2betascan

from glactools.

glactools acf2betascan --useroot outgroup_root_test_merge.chr24.acf.gz | gzip > outgroup_root_test.chr24.beta.txt.gz

Is it right to do ?

from glactools.

should be!

from glactools.

Ok! Thank you very much !

from glactools.