Question: remove in mapping results the reads that aligned > 1 times
0
gravatar for cpozacarrion
9 weeks 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 • 114 views
ADD COMMENTlink modified 9 weeks ago by Devon Ryan1.5k • written 9 weeks ago by cpozacarrion0
0
gravatar for Devon Ryan
9 weeks ago by
Devon Ryan1.5k
Germany
Devon Ryan1.5k 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 9 weeks ago by Devon Ryan1.5k
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: 88 users visited in the last hour