Error with parsing alleles depth

[glaudelrio]

Error with parsing alleles depth

[btmartin721]

Hi.

The error is occurring because the ipyrad and stacks VCF files have a field that has the allele depths in it, but your GATK 4 file seems to either not have the allele depths (AD) field or it might be in a different location or have a different format than ipyrad and stacks. vcf2bgc.py hasn't been tested with GATK 4 VCF files, and VCF files have a lot of variant formats.

If you could send me your VCF file, I can see if it would be possible to add support for GATK 4 VCF files. If it has the Allele Depth (AD) field in the FORMAT tag, then it should be possible to add it. See the example record below to see where you can find the AD field. You can email me your VCF file if you want and I'll see if I can add support for that format.

#CHROM  POS  ID  REF  ALT  QUAL  FILTER  INFO  FORMAT  NA12878 [other samples...]
20  10001019   .  T  G  364.77  .  [CLIPPED]  GT:AD:DP:GQ:PL  0/1:18,15:33:99:393,0,480

[glaudelrio]

Hi Bradley! Thank you so much for your help! Here is one of the VCF files that I am trying to use. I was able to work with the genind though, but then there was another error when running combine_bgc_output. I posted the last error a few minutes ago. I really appreciate all the help, and sorry for all the inconvenience. Best wishes! chr22_ldna_test.recode.vcf <https://drive.google.com/file/d/1EmTDVITmULYL5sCiohgloCCLus4zCFlq/view?usp=drive_web>

…

On Mon, May 2, 2022 at 4:35 PM Bradley Martin ***@***.***> wrote: Hi. The error is occurring because the ipyrad and stacks VCF files have a field that has the allele depths in it, but your GATK 4 file seems to either not have the allele depths (AD) field or it might be in a different location or have a different format than ipyrad and stacks. vcf2bgc.py hasn't been tested with GATK 4 VCF files, and VCF files have a lot of variant formats. If you could send me your VCF file, I can see if it would be possible to add support for GATK 4 VCF files. If it has the Allele Depth (AD) field in the FORMAT tag, then it should be possible to add it. See the example record below to see where you can find the AD field. You can email me your VCF file if you want and I'll see if I can add support for that format. #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT NA12878 [other samples...] 20 10001019 . T G 364.77 . [CLIPPED] GT:AD:DP:GQ:PL 0/1:18,15:33:99:393,0,480 — Reply to this email directly, view it on GitHub <#19 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/ACQPWYLTL4T3BXACBWWLY4DVIBDDXANCNFSM5U2PAGPQ> . You are receiving this because you authored the thread.Message ID: ***@***.***>

-- *Glaucia Del-Rio* *www.gdel-rio.com <http://www.gdel-rio.com>* Museum of Natural Science Louisiana State University

[btmartin721]

Hi,

I was not able to reproduce the error you received with the VCF file you sent me. Note that I subset the number of sites to the same number that it printed above, and I am not 100% sure that I assigned the same individuals to P1, P2, and admixed populations. But it seemed to run all the way through. Is the VCF file you ran above any different than the one you sent me?

-Bradley

[glaudelrio]

Hi Bradley, I tried running absolutely the same vcf (I prefer working with the vcf to bgc instead of geeing) and I have the same problem again. Do you think there is anything wrong with my installation? I am not familiar with docker.

Thank you so much!

[btmartin721]

Again, I apologize for my delayed reply. I think I know what the issue is. I added support for stacks recently, and I think I need to rebuild the docker image. So in other words, the vcf2bgc.py file in the ClineHelpR/scripts directory is correct, but I haven't added the updated file to the docker image.

[btmartin721]

the reason I mentioned stacks is because I think that format also works with the VCF file generated by GATK4.

[btmartin721]

Hi. I fixed the issue with the Docker version of the vcf2bgc.py script. It should work now. I'm closing the issue, but if it still doesn't work let me know.

You should run:

sudo docker pull btmartin721/clinehelpr:latest

to pull the new docker image.