Filtration of paired reads with both length higher than 50 from bam file
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?
• 127 views
Read more here: Source link