Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Is there a way to get transcription factors from gene list? Jluis Bioinformatics 6 12-14-2012 03:46 AM
How to get transcription factors from gene list ETHANol Bioinformatics 1 09-01-2012 10:01 AM
List of all COGs/KOs/gene functions leknif Bioinformatics 0 06-15-2011 04:27 AM
Gene list for calculating coverage nseh Bioinformatics 1 05-22-2011 07:31 AM
Is a Gene mapping tool available? ritzriya Bioinformatics 0 02-15-2011 01:45 AM

Thread Tools
Old 01-16-2013, 04:04 PM   #1
Junior Member
Location: Australia

Join Date: Feb 2012
Posts: 3
Default Tool to take mapping to gene list

Hello all,

I am relatively new to NGS and have just mapped a series of reads to a genome using Bowtie2. As I understand, the SAM/BAM output contains the information telling me where the reads map to the genome. I now want to:

a) transform this to a list of regions where reads map i.e., combine all the overlapping reads that map to one part of the genome into a starting position (bp), ending position (bp), and the number of reads that map to that region

b) cross-reference that list of regions against the genome annotation to identify whether that collection of reads maps to an exon, intron, etc. and, if an exon, extract the relevant gene and transcript id

I am sure there must be a tool for this - probably several! I am interested to know what people use? Any relative advantages/disadvantages, key tips or tricks, etc., would also be helpful.

Thanks for your help.


jarthur is offline   Reply With Quote
Old 01-16-2013, 04:31 PM   #2
Senior Member
Location: US

Join Date: Jan 2009
Posts: 392

I'm assuming this is RNA-seq data? Are you wanting simply a count of reads mapping to a feature?

If so, then there are a variety of tools that can accomplish this. The two most commonly used are probably BEDtools and HTSeq-count.
chadn737 is offline   Reply With Quote

gene_id, tool

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 03:31 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO