Getting Counts From Samtools
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?
• 4,007 views
•
link
updated
by
Sean Davis
•
written
by
Raygozak
•
written
by
Raygozak
1 answer
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.
• 2 views
•
link
written
by
Sean Davis
Log in to answer this question.
More posts like this
-
How does kallisto handle multi mapped reads?
written by bioinfo •
-
Graphing NGS data
written by n.joshi17 •
-
How to quantify gene expression when stringtie considers two adjacent genes as single gene?
written by kousi31 •
-
Mapping SNP position to gff file
written by shachigahoimbi •
-
TCGA SNP arrays
written by joan.frigola •
-
IGV genome browser problems
written by Marvin
-
Counting sense reads in bacterial paired-end RNA-seq data
written by biotech
-
Find The Number Of Simple Repeats On The Y Chromosome
written by zhaolin20042008 •
-
Should Rna-Seq Reads Map To One Strand Of Cdna Reference?
written by chjiao3456 •
-
Calculating Ds/Dn Ratio Through Contig Snp Annotation
written by jackuser1979
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'?