Open TanguyMuller opened 4 months ago
Note: could not parse as PED: sample.txt
It is not an error, only a warning (albeit admittedly a bit useless and confusing). Check if the GT fields have the correct ploidy for the right samples and chromosomes, then it worked.
Thank you for this response. But precisely because bcftools does not read the ped.file, the GT fields are diploids for the female individuals and don't give the correct ploidy. Then i don't realy understand why bcftools can't read this file.
The file sample.txt
you showed above is not a PED file. For the format definition see for example here https://gatk.broadinstitute.org/hc/en-us/articles/360035531972-PED-Pedigree-format
Regardless, if the ploidy is not determined and output correctly based on your current sample.txt
, that would indicate a bug. Any chance you could provide a small test case to reproduce the problem? A few lines from the mpileup
command, including the full headers, would suffice.
Hello everyone,
I am conducting research on the Pine Processionary Moth and I want to call variant on the sexual chromosome. In my data I have males that are homochromatic ZZ and females that are heterochromatic ZW. So I would like to make a variant calling on the chrZ specifying that the male individuals are diploid and female haploid. To do this I run the following command:
with variant_call.sh :
and this is the error I get
I am having difficulty understanding this because my file is in a specific format.
My sample.txt file :
and my ploidy.txt file is :
Could someone please help me ??