ultimatesource / denovogear

A program to detect denovo-variants using next-generation sequencing data.
http://www.nature.com/nmeth/journal/v10/n10/full/nmeth.2611.html
GNU General Public License v3.0
49 stars 25 forks source link

[E::bcf_hdr_add_sample] Duplicated sample name '.' Aborted (core dumped) #302

Closed Saranya-Balachandran closed 5 years ago

Saranya-Balachandran commented 5 years ago

Hi, On running dng dnm auto --vcf vcf_f --ped ped_f -w out.vcf, I get [E::bcf_hdr_add_sample] Duplicated sample name '.' Aborted (core dumped) On querying the sample names using bcftools I get three distinct sample names. The bcf file has DP, PL record for each sample individually.

Kindly help in understanding this issue. I'm runing in the denovogear developer mode

Regards, Saranya

Saranya-Balachandran commented 5 years ago

On running on the test data generated by make test in human_trio receiving the same error dng dnm auto --ped human_trio/ped/trio.ped --vcf human_trio/bcftools/trio.call.vcf

reedacartwright commented 5 years ago

The sample data you used is for dng call so the error is expected.

It appears you are using the wrong pedigree format for dng dnm. Can you update to the format sorted by dng dnm?

On Thu, Jul 18, 2019, 07:54 Saranya-Balachandran notifications@github.com wrote:

On running on the test data generated by make test in human_trio receiving the same error dng dnm auto --ped human_trio/ped/trio.ped --vcf human_trio/bcftools/trio.call.vcf

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/ultimatesource/denovogear/issues/302?email_source=notifications&email_token=AAAY4HWEIU36RA7XNGROL4DQABKXRA5CNFSM4ID76SO2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOD2IHRDY#issuecomment-512784527, or mute the thread https://github.com/notifications/unsubscribe-auth/AAAY4HU5GMSQGKMX6S5XFHTQABKXRANCNFSM4ID76SOQ .

Saranya-Balachandran commented 5 years ago

Please find below my ped format

PEDNG v1.0

ID Dad Mom Sex Samples

SID93046 . . 2 = SID93044 . . 1 = SID93041 SID93044 SID93046 1 =

reedacartwright commented 5 years ago

That's the format recognized by dng-call. To use dng-dnm follow the instructions here:

https://github.com/ultimatesource/denovogear/blob/develop/README.md#about-sampleped

On Thu, Jul 18, 2019, 09:43 Saranya-Balachandran notifications@github.com wrote:

Please find below my ped format

PEDNG v1.0

ID Dad Mom Sex Samples

SID93046 . . 2 = SID93044 . . 1 = SID93041 SID93044 SID93046 1 =

— You are receiving this because you commented. Reply to this email directly, view it on GitHub https://github.com/ultimatesource/denovogear/issues/302?email_source=notifications&email_token=AAAY4HQFRCVRG4IFYNDTVNDQABXPHA5CNFSM4ID76SO2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOD2IQRDY#issuecomment-512821391, or mute the thread https://github.com/notifications/unsubscribe-auth/AAAY4HSHAMOJF2KHW6UPF3LQABXPHANCNFSM4ID76SOQ .

Saranya-Balachandran commented 5 years ago

Thank you!