Does anyone know how to merge two BWA alignments into one using samtools? I have a pe read alignment and a contig alignment, both done with BWA and I was wondering how difficult it would be to merge these into one consensus alignment?
The authors of the following paper (Whole genome sequencing of the black grouse (Tetrao tetrix): reference guided assembly suggests faster-Z and MHC evolution)
did this with contigs and mate-pair reads:
"The BAM format alignment files of the contig mapping and the mate-pair read mapping were subsequently merged using SAMtools suite v0.1.18"
Zach
1 answer
I am not sure if I have understood your question correctly but if you have aligned both contigs and mate-pair reads to the same reference genome, you can use samtools "merge" function to merge the two bam/alignment files.
•
written
by
Ashutosh Pandey
Log in to answer this question.
-
How to extend assembled PACBIO contigs for scaffolding
written by K •
-
Filtering BAM file to include PE reads that did not map, even though their mate did
written by nute11a •
-
To extract unmapped pair end reads
written by Bioinfonext
-
PE read flag 69 and 137 meaning
written by epigene
-
AlignGraph doesn't response but i have what it needs
written by Calvin •
-
Tutorial: InsideDNA: Samtools guide: learning how to filter and manipulate with SAM/BAM files
written by shilparaopradeep •
-
Merge and fixmate of multiple BAM files generated from pair-end DNA sequences
written by lisadavic66 •
-
Paired-End Sam Flag Mismatch During Filtering
written by Andreas
-
Bwa Sampe Very Large Insert Size
written by Bioscientist
-
Aligning Mate Pair Data
written by Abhi
