I need to extract fastq seq. from a BAM file for other experiments. But the fastq file should only contain unique fastq seq. [Here unique fastq means if there are multiple copies of a seq. in BAM file then in the fastq file I only want a single copy of that seq. ]
•
written
by
ved_vyas •
2 answers
The operation that you ask about is called deduplication.
But first, you'd need to extract the sequences which can be slightly more complicated for paired-end reads, you would need to collate then extract, see this post
https://lh3.github.io/2021/07/06/remapping-an-aligned-bam
Deduplication:
If you want to make sure that you get each read only once, even if it aligned to multiple places, and don't care if two separate reads have identical sequences, you can filter your bam file to remove secondary alignments. Each read will have exactly one primary alignment.
Log in to answer this question.
-
how to extract unique snps in a vcf file by comparing with multiple vcf files
written by nikhil •
-
Is there any tutorial page that tells you how to extract read count from illumina single-cell RNAse…
written by Simon Ahn •
-
extract gene sequence from multiple bam files in order to do Ka/Ks test
written by annaA •
-
Separating each cell in fastq file
written by sidrah.maryam •
-
Extract Sequences from Multiple sequence alignment containing special character like gap?
written by MSRS
-
scRNAseq: Split FASTQ into cell specific FASTQ Files
written by niklas.lang •
-
Getting consensus fasta sequence for a chromosome from a bam file
written by nagarsaggi •
-
tool for a pattern of repeats
written by 2nelly
-
Ensembl Hg38 dna, dna_rm, and dna_sm
written by bharata1803
-
Extracting from tophat outputs reads pairs and splice-junctions with a single best match
written by trakhtenberg
