Question: How do I filter a BAM file for proper pairs WITHOUT samtools
0
2
Devon Ryan ♦ 94k wrote:
In deepTools:
alignmentSieve --samFlagInclude 2 -b input.bam -o output.bam
You can use multiple threads if you like. Of course this uses pysam, which uses htslib, which is part of samtools, but most tools will use something related to samtools somewhere...it's a standard tool after all.
2
Pierre Lindenbaum ♦ 126k wrote:
using samjdk http://lindenb.github.io/jvarkit/SamJdk.html
java -jar dist/samjdk.jar -e 'return record.getProperPairFlag();' input.bam
Please log in to add an answer.
Use of this site constitutes acceptance of our User
Agreement
and Privacy
Policy.
Powered by Biostar
version 2.3.0
Traffic: 1526 users visited in the last hour
May I ask you why samtools is not allowed?
We have some conflicts with samtools
That's incredibly unusual, it's likely you should just fix your environment then.
I would normally agree, but we have reasons for not using samtools in this instance