Question: Getting Counts From Samtools
gravatar for Raygozak
8.7 years ago by
State College, PA, Penn State
Raygozak1.3k wrote:

I have the result from aligning illumina reads to a reference genome and i'm able to get the counts per mapped read, however it outputs chromosome name starting position , the end result i want is to replace this information with the actual gene/mrna name that's contained in a gff file for the reference sequence, i haven't found a tool that does this. does such thing exist?

samtools tophat rna • 2.9k views
ADD COMMENTlink modified 8.7 years ago by Sean Davis26k • written 8.7 years ago by Raygozak1.3k

You really need to clarify your question. It sounds like you want to get read count per annotation in the GFF file? What do you mean you already have the 'counts from a BAM files per read'?

ADD REPLYlink written 8.7 years ago by Damian Kao15k
gravatar for Sean Davis
8.7 years ago by
Sean Davis26k
National Institutes of Health, Bethesda, MD
Sean Davis26k wrote:

Just guessing what you want, but you'll want to take a look at, for example:

If you are an R/Bioconductor user, take a look at the GenomicFeatures and GenomicRanges packages.

ADD COMMENTlink written 8.7 years ago by Sean Davis26k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 733 users visited in the last hour