Samir
  • Member for 3 years, 5 months
  • Last seen more than a month ago
2 answers
18 votes
9k views
Obtaining uniquely mapped reads from BWA mem alignment
4 votes

For much faster and stable$ implementation, sambamba can handle this using following one-liner: sambamba view -t 12 -h -f bam -F "mapping_quality >= 1 and not (unmapped or secondary_alignment) and ...

View answer