Closed:Variant calling by scaffolds with bcftools
0
0
Entering edit mode
4.4 years ago
Dani • 0

Hi, I am trying to do variant calling using bcftools both for mpileup and calling. I have around 60 bam files that were mapped using bwa-mem to my reference genome that contains 92 scaffolds.

I tried doing it all at once but after 30 days it only went through 20 out 92 scaffolds. What I want to do now is to do my variant calling by each scaffold and then merge all of my vcf files...

I tried with -r but I get an error:

[mpileup] fail to load index for ind1.bam Failed to open -: unknown file type.

Has anyone done this using bcftools?

Thank you for your help :).

next-gen • 193 views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 2378 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6