Getting all unmapped reads

Question: Getting all unmapped reads

17 months ago by

mm4184 • 20

mm4184 • 20 wrote:

I have a fastqsanger dataset from an RNAseq experiment, and I am trying to remove contaminating human sequences. I can run BWA for Illumina against the human genome and get BAM files of the mapping, but I am trying to get a fastqsanger file with only the unmapped reads. I have been able get this as a standard output option in other programs, but I can't seem to find out how to do this here. I also tried RNA Star, but it is listed as "new" and it looks to me like no reference genomes are available to it yet.

Thanks for your help.

--Michael

rna-seq bwa • 1.4k views

ADD COMMENT • link •

modified 17 months ago by Guy Reeves • 1.0k • written 17 months ago by mm4184 • 20

17 months ago by

Jennifer Hillman Jackson ♦ 25k

United States

Jennifer Hillman Jackson ♦ 25k wrote:

Hello,

All BWA tools at http://usegalaxy.org do not offer an option to output fastq datasets of unmapped reads directly. Instead, use tools to filter the mapped BAM result dataset to only include unmapped reads then convert SAM/BAM to fastq. There are a few tools to do the filtering and multiple filtering criteria can be set - including filters to capture low-quality mapping results.

One SAM/BAM filtering tool is in the Picard group along with a tool to convert SAM to FASTQ (it also accepts BAM input). SamTools and BAMtools also have BAM/SAM filtering tools. When using these tools be sure to sort the BAM/SAM dataset first, as a distinct step: https://galaxyproject.org/support/sort-your-inputs/

Some of the newer wrappers for mapping tools do include this option, including RNASTAR. RNASTAR is pending indexes, so use a custom genome for now: https://galaxyproject.org/learn/custom-genomes/

The same fastq data generally not appropriate for both BWA and RNASTAR but this may be why you want to capture the unmapped reads in your search for contamination.

Hope this helps! Jen, Galaxy team

ADD COMMENT • link modified 17 months ago • written 17 months ago by Jennifer Hillman Jackson ♦ 25k

17 months ago by

Guy Reeves • 1.0k

Germany

Guy Reeves • 1.0k wrote:

You may want to look at this post which asks a related question https://biostar.usegalaxy.org/p/21898/

ADD COMMENT • link written 17 months ago by Guy Reeves • 1.0k

Please log in to add an answer.

Similar posts • Search »