How to filter both reads of paired end reads by AS when only 1 of the 2 reads has an AS below 50
I am trying to filter reads with an AS below 50. I found that I could use sambamba for this with the command:
sambamba view -F '[AS]' >= 50 /path/to/input/file > /path/to/output/file
The problem with this command is that if a read has a mate with an AS above that threshold, that read won't be deleted and the flag of the mate with no pair still says that it is read is paired. Is there a command with sambamba to filter both pairs if 1 has a lower AS than the threshold or can I filter unpaired reads with a wrong flag with another tool?
Kind Regards, Jelle
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Clarification regarding SAM flags "mate reverse strand" (flag 16/0x10) and "read reverse strand" (f…
written by kalavattam •
-
Importing read pairs from BAM with specific R1 and R2 start/end positions for pair
written by Emily Wasson •
-
Bamtools stats and samtools flagstat produce different results and confusing my sambamba output.
written by DNAngel
-
Error during pre-processing of RNA-seq reads for variant calling using Opossum/0.2
written by UDAY.AGRI123 •
-
Filtering BAM file to include PE reads that did not map, even though their mate did
written by nute11a •
-
Fixing inconsistent flags (strand info set for only one read) between RNA-seq read pairs with samto…
written by ophilia30092 •
-
How would you filter PCR duplicates for merged paired-end reads
written by hohoku •
-
Samtools stats "inward oriented pairs" - what flags does this correspond to?
written by roryC •
-
Samtools mpileup vcf file
written by AW
-
Doing genome assembly using SPAdes with mate-pair library only on public Galaxy server
written by sentausa
