Filtration of paired reads with both length higher than 50 from bam file

Filtration of paired reads with both length higher than 50 from bam file

1

I understand that for filtering reads based on length i can simply use:

 samtools view -h file.bam | awk 'length($10) > 50 || $1 ~ /^@/' | samtools view -bS - > output.bam

But how can i select paired reads that both are longer than 50?


bam


awk


sed

• 127 views

updated 2 hours ago by

106k

written 5 hours ago by

▴

10

Read more here: Source link