Merging mapping results
Hello
I think this is a silly question, nonetheless I need help with this issue :p
I have several bam files that I have obtained from some mapping procedures I have done with bwa, and I have retrieved the fastq sequences from these files.
So, I need to put all the fastq sequences in one file, but without repetitions. (If one read was a hit in all mappings I need to have it only once in the final file; if another appeared in just two mappings; again, I need it only once... you get the picture).
Thanks in advance.
1 answer
Just put them all in one file and filter out duplicates, see this thread:
• 0 views
•
link
written
by
Biomonika (Noolean)
Log in to answer this question.
More posts like this
-
Bowtie2 alignment on metagenomics - sum same sequence
written by rfour92 •
-
BWA MEM output directory
written by sarastrafella.ss •
-
Read Mapping with BWA mem with a bunch of files
written by zhaoenid •
-
How to loop over multiple files for BWA MEM in bash for loop?
written by DNAngel
-
Automating to convert multiple fastq files into one fastq file
written by zhou_1228 •
-
What kind of sequences I have, haploid or diploid?
written by Enrique López •
-
BWA mapping - only mapped the part of the reads, why?
written by gerberd1990 •
-
Contig alignment software
written by bpz •
-
Sequence extraction from a fastq file
written by bpz •
-
Sra-Toolkit 2.3.1 Fastq-Dump Multiplefiles At Once Windows Command Line
written by deepbiofever •
grep read names from fastq_file_1 into file_list, then grep read names from fastq_file_2 and add them to the same file_list... do it for all your samples, remove duplicates, and extract the reads with seqtk.