Question: How do I find unmapped reads in my sam file?
0
gravatar for blur
3.4 years ago by
blur110
European Union
blur110 wrote:

Hi,

I want to see why my unmapped reads did not map to the reference - I wanted to pull out a few of the unmapped reads to blast them against the reference. Any idea how I do so? Thanks

sam • 2.8k views
ADD COMMENTlink modified 3.4 years ago by vivekbhr540 • written 3.4 years ago by blur110

hi, i am working on differential gene expression analysis on zebrafish datasets, i am working with GPL 14664 agilent platform but ther is annotation package available. i used limma package and i got results. but how to convert into gene symbols.

ADD REPLYlink written 3.4 years ago by rajasekargutha30
1
gravatar for vivekbhr
3.4 years ago by
vivekbhr540
Germany
vivekbhr540 wrote:

Use samtools. Just do samtools view -f 4 <bamfile> ..

You can extract all kinds of reads by specifying the flag No. after samtools -f . Get the relevant flags here

ADD COMMENTlink written 3.4 years ago by vivekbhr540

so, first turn it into bam, right?

ADD REPLYlink written 3.4 years ago by blur110

It does not matter if the file is SAM or BAM formatted. Samtools view command is able to read both files. ;)

ADD REPLYlink modified 3.4 years ago • written 3.4 years ago by glihm610

thanks, that was super helpful! :)

ADD REPLYlink written 3.4 years ago by blur110
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: 1076 users visited in the last hour