Question: How to obtain unique reads on Galaxy
0
gravatar for miriam115
4.1 years ago by
miriam1150
European Union
miriam1150 wrote:

Hi,

I am trying to use Galaxy for some NGS data. I aligned the data by using bowtie2. I can't figure out how to extract only the unique-mapped reads.

On the report of the alignment I can see this message:

"20552525 (13.02%) aligned 0 times

83720562 (53.02%) aligned exactly 1 time

53619781 (33.96%) aligned >1 times".

And I want only the 53.02%.

I am looking for any sort of help - on Galaxy.

Thanks

alignment • 1.4k views
ADD COMMENTlink modified 4.0 years ago by Biostar ♦♦ 20 • written 4.1 years ago by miriam1150

On Galaxy, both BAMtools (NGS:BAMtools and SAMtools (NGS:SAMtools) can filter on the bit flags. It looks like the BAMtools filter set is a little more intuitive. I think if you play with the combinations of filters, you should be able to get your 53.02%.

ADD REPLYlink written 4.1 years ago by Joseph Pearson450
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 2.3.0
Traffic: 1130 users visited in the last hour