Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Forum
Planet
All »
Home
Log In
Question: Get multimapped reads from STAR in Galaxy.
0
gravatar for luisafermercadom
8 months ago by
luisafermercadom • 0
luisafermercadom • 0 wrote:

Hi,

Using RNA STAR in Galaxy to map my yeast read back to the genome: (single end reads), I'm getting a huge number of multi-mapped reads (40-60%). I want to get those multimapped reads, so I can look at them ...

How do I write out a multi-mapped reads file from Galaxy? (We do not have an instance here.)

Thank you!
fastq rna-seq sam filter bam • 401 views
ADD COMMENT • link •
modified 21 days ago by sf.rocha • 0 • written 8 months ago by luisafermercadom • 0
1
gravatar for Jennifer Hillman Jackson
8 months ago by
Jennifer Hillman Jackson ♦ 25k
United States
Jennifer Hillman Jackson ♦ 25k wrote:

Hello,

BAM datasets can be filtered with either of the tools below. If you want the fastq sequences, those can be extracted from the filtered BAM dataset.

Filter BAM:

  • NGS: SAMtools > Filter SAM or BAM, output SAM or BAM
  • NGS: BamTools > Filter BAM datasets on a variety of attributes

BAM/SAM to Fastq:

  • NGS: Picard > SamToFastq extract reads and qualities from SAM/BAM dataset and convert to fastq

Thanks! Jen, Galaxy team
ADD COMMENT • link written 8 months ago by Jennifer Hillman Jackson ♦ 25k
0
gravatar for luisafermercadom
8 months ago by
luisafermercadom • 0
luisafermercadom • 0 wrote:

Thank you,

I could easily separate the multi-mapped reads using NGS: BamTools > Filter BAM datasets on a variety of attributes
ADD COMMENT • link written 8 months ago by luisafermercadom • 0
0
gravatar for sf.rocha
21 days ago by
sf.rocha • 0
sf.rocha • 0 wrote:

I'm having issues with using the NGS: BamTools > Filter BAM datasets on a variety of attributes. I'm using as input a BAM file obtained with HISAT2 and I want to filter for only the reads that have only 1 mapping hit. For this I filtered on: Select BAM property to filter on > tag Filter on a particular tag > NH:=1

In the 'filtered bam' output file I get this as a summary: { "filters": [ { "tag": "NH:=1", "id": "1" } ] }

and if I visualize the data by clicking on the 'eye' icon, I only have the header.

Am I doing anything wrong?
ADD COMMENT • link written 21 days ago by sf.rocha • 0

Hi - Use this term instead: NH:1

Thanks! Jen, Galaxy team

ADD REPLY • link written 20 days ago by Jennifer Hillman Jackson ♦ 25k
Please log in to add an answer.

Similar posts • Search »

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