How to obtain unique reads on Galaxy
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
• 2,705 views
•
link
updated
by
Biostar
•
written
by
miriam115 •
•
written
by
miriam115 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Why Is My Uniquely Mapped Read Count Higher in Samtools Compared to Bowtie2?
written by daffodil •
-
Request to validate my RNAseq workflow
written by umn_bist •
-
Extract Sequences From [Gtf File] + [Genome Fasta File], combine matching XLOCs
written by Tom Koch
-
bowtie2 error in mapping pair ended sample
written by kkj3083 •
-
Please help me with bowtie2 output and samtools report
written by seta
-
Paired bowtie2: 91% overall alignment but 0.01% concordant - should I be worried?
written by Bilal Akil •
-
R: remove rows by a list of rownames
written by silvi_free88 •
-
Unaligned Reads From Bowtie2
written by deepthithomaskannan
-
Does Bowtie2 Preform Unique Mapping
written by Raghav
-
Extraction Of Desired Alignments From Sam File
written by Raghav
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%.