Question: remove in mapping results the reads that aligned > 1 times
0
gravatar for cpozacarrion
19 months ago by
cpozacarrion • 0 wrote:

I would like to remove from the BAM file after mapping using bowtie2 tool, the reads that aligned > 1 times. Any know how remove these reads?. Thanks. César.

alignment bowtie galaxy • 1.2k views
ADD COMMENT • link • modified 19 months ago by Devon Ryan • 1.9k • written 19 months ago by cpozacarrion • 0
0
gravatar for Devon Ryan
19 months ago by
Devon Ryan • 1.9k
Germany
Devon Ryan • 1.9k wrote:

You can use the "Filter SAM or BAM, output SAM or BAM" tool under "NGS: Samtools". Just set the "Minimum MAPQ quality score" (the "quality" is redundant in that label...) to 2. This will remove slightly more than just those alignments that align >1 times, but the difference is minor (if you wanted to truly only remove those alignments counted as "aligned > 1 times" then you'd have to compare the AS and XS auxiliary tags...which is much more of a pain than it's worth)

ADD COMMENT • link written 19 months ago by Devon Ryan • 1.9k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 169 users visited in the last hour