Question: remove in mapping results the reads that aligned > 1 times
0
gravatar for cpozacarrion
19 months ago by
cpozacarrion0 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 COMMENTlink modified 19 months ago by Devon Ryan1.9k • written 19 months ago by cpozacarrion0
0
gravatar for Devon Ryan
19 months ago by
Devon Ryan1.9k
Germany
Devon Ryan1.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 COMMENTlink written 19 months ago by Devon Ryan1.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