Question: Extract Raw Counts From Sam File
0
gravatar for mad.cichlids
6.1 years ago by
mad.cichlids100
Texas
mad.cichlids100 wrote:

Hi, All

I was performing a de novo assembly of RNA seq data. I aligned my reads to the assembly for each of my samples and generated sam file for each, my question is how am i to extract the raw reads from the sam files? I know couple of tools such as htseq can extract reads with the reference genome, but how about the tools without any reference genome to perform the task of extracting raw reads?

Thanks a lot for any suggestions.

Best,

rnaseq sam counts • 3.4k views
ADD COMMENTlink modified 6.1 years ago by Rm7.9k • written 6.1 years ago by mad.cichlids100
1
gravatar for Gabriel R.
6.1 years ago by
Gabriel R.2.6k
Center for Geogenetik Københavns Universitet
Gabriel R.2.6k wrote:

Sam or bam ?

Do you just want to see your reads + quality ?

If sam, just use

cut -f 10,11 file.sam

if bam

samtools view file.bam | cut -f 10,11

Now if you want your data fasta or fastq, use this hacked version of samtools which produces fastq/a : https://github.com/udo-stenzel/samtools-patched

ADD COMMENTlink written 6.1 years ago by Gabriel R.2.6k

Thank you so much! Sorry, i could have stated it more clearly. I meant to ask how many raw reads aligned to each of my assembled contigs, so that i can compare the gene expression differences later. the 10, 11 are respective sequence and quality score right?

ADD REPLYlink modified 6.1 years ago • written 6.1 years ago by mad.cichlids100
0
gravatar for Rm
6.1 years ago by
Rm7.9k
Danville, PA
Rm7.9k wrote:

use eXpress pipeline

ADD COMMENTlink written 6.1 years ago by Rm7.9k

awesome, really appreciate it!

ADD REPLYlink written 6.1 years ago by mad.cichlids100
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: 1316 users visited in the last hour