SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
any in silico enzyme digestion program for genome sunsnow86 Bioinformatics 5 10-08-2012 07:09 PM
.abi to fasta/fastq conversion script/program? AppleInformatics General 12 08-26-2012 10:17 PM
automated pipeline casava or FASTQ script/program sdj Illumina/Solexa 3 01-12-2012 06:02 AM
Program to relate SNPs to model genome annotation? pmiguel Bioinformatics 4 04-23-2010 06:33 AM
whole genome Bisulfite Sequence MAPping program wei Bioinformatics 0 08-07-2009 02:46 PM

Reply
 
Thread Tools
Old 08-13-2012, 12:34 PM   #1
mlafave
Junior Member
 
Location: Bethesda, MD

Join Date: Feb 2012
Posts: 5
Default Genome annotation program/script?

Hi, all - I'm looking for a way to annotate a list of positions in a genome. Basically, I have a list of positions that looks like this:

chr1 3026294
chr1 30175158
chr2 17717521

...and so on. I want to know:
1) the names of the gene(s) (if any) that overlap that position
2) the orientation of those genes, and
3) which part of the gene the listed position is (5' UTR, second exon, fourth intron, etc.).

I feel like someone must have made something to do this before, but I don't know where to look. Any ideas? Thanks!
mlafave is offline   Reply With Quote
Old 08-15-2012, 03:41 AM   #2
colindaven
Senior Member
 
Location: Germany

Join Date: Oct 2008
Posts: 415
Default

Try Bedtools.

You'll need an annotation of course.

Galaxy is also good - see the interval functions there.
colindaven is offline   Reply With Quote
Old 08-19-2012, 05:37 PM   #3
adaptivegenome
Super Moderator
 
Location: US

Join Date: Nov 2009
Posts: 437
Default

snpEff is good to, very simple to use
adaptivegenome is offline   Reply With Quote
Old 08-20-2012, 07:17 AM   #4
ishmael
Member
 
Location: NY, US

Join Date: Jul 2008
Posts: 17
Default

http://code.google.com/p/diffreps/
I would suggest diffReps. There is a script named region_analysis.pl.
It could read bed format files and generate annotated files fitting most your demands.
The output looks like:
Chrom Start End GName TName Strand TSS TES Feature D2TSS
chr1 2986341 2986690 Ust NM_001108458 - 2868195 3158231 Genebody 171715.5

Refseq and ensembl annotation surpported, but only human, mouse, rat genomes available now.
ishmael is offline   Reply With Quote
Old 08-20-2012, 06:29 PM   #5
adamyao
Member
 
Location: Taiwan

Join Date: Feb 2011
Posts: 19
Default

You can try VarioWatch (http://genepipe.ncgm.sinica.edu.tw/variowatch/main.do). It should be able to provide what you want with graphical output in real time. By using "Query by Bacth" can give you results for multiple positions. If you have thousands of positions to query then you can use MegaQuery.

VarioWatch only supports human genome.
adamyao is offline   Reply With Quote
Old 08-24-2012, 04:59 AM   #6
mlafave
Junior Member
 
Location: Bethesda, MD

Join Date: Feb 2012
Posts: 5
Default

Thanks for your help, everyone!
mlafave is offline   Reply With Quote
Reply

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 02:27 AM.


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