merge paired-end sequencing
Heys,
I'm creating a pipeline to do the SNP calling of whole-genome sequencing for a non-model organism and I have a doubt regarding the merging of the paired-end data. When should it occur? I already did all the pipeline (quality filtering, mapping, sorting, remove duplicates, and SNP calling with GATK), but in GATK best practises I never saw anything about merging the paired-end files. Could you light me on this?
Thanks, Gabriel
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Overlapping reads, how do we counter them while variant calling, haplotype calling, does it affect …
written by ayeraselvan •
-
Seeking Opinions on Pipeline for Detecting Somatic Variants and Rare Germline Mutations
written by George •
-
Merge single end and paired end illumina raw sequences for SNP calling
written by prashantwaiker •
-
problem with a for loop
written by gubrins
-
quality check of bam files
written by gubrins
-
ctDNA WES duplication rates from NovaSeq 6000 and HiSeq 4000
written by graeme.thorn
-
How do I chose variant quality filters for a VCF file?
written by James Reeve
-
Query related to GATK
written by Bioinfonext
-
workflow for variant calling using GATK in multiple paired end reads
written by reza
-
Contig/Scaffold Length- Variant Calling
written by williamjohn360
The alignment takes care of the paired-end data, beyond that there should be no specific step. Is there a specific concern you have, do the reads overlap each other?
Thanks for your reply, I think my reads don't overlap, I was just afraid I missed a step, as rad-seq's softwares (as ipyrad) does a merging in the final outputs. Then, if I aligned with bwa my pair-end data that's all I should have done? Thanks again!