Combine haploid fastq's into diploid "individuals"
0
0
Entering edit mode
7 weeks ago

I have fastq files, where each file is a this sequence of a distinct haploid individual. I need to run these through GATK as though they were diploids, in order to use a software which takes only a VCF with diploid samples as an input. This post: How to merge two haploid samples (vcf, or g.vcf) into a pseudo-diploid? suggests merging bam files, and some post-merge processing, but couldn't I just merge pairs of my fastq's to get two haplotypes into one file? eg:

> cat read1_indv1.fq.gz read1_indv2.fq.gz > read1_combined.fq.gz 
> cat read2_indv1.fq.gz read2_indv2.fq.gz > read2_combined.fq.gz

Then do all the sam bam GATK stuff after the fact?

Thanks

GATK pop-gen haplodiploid • 225 views
ADD COMMENT
0
Entering edit mode

Are you absolutely sure this is necessary? What happens if you just put one fastq through the software? Will the software really refuse to run if it doesn't detect any mixed loci?

ADD REPLY
0
Entering edit mode

Based on my conversation with the author, merging to diploid individuals would be the easiest route.

ADD REPLY

Login before adding your answer.

Traffic: 2675 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