SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
finding hangovers in contig/scaffold alignments fahmida Bioinformatics 1 03-07-2013 05:30 AM
extracting predicted gene from scaffold: end position precedes start position amango Bioinformatics 1 09-02-2012 08:25 AM
Finding per gene GC content. chadn737 Bioinformatics 3 03-02-2012 02:43 AM
Filter SOLiD reads for gene of interest then assemble john1923 SOLiD 2 04-12-2011 03:47 AM
Gene Finding archananraja Bioinformatics 2 02-08-2011 07:25 AM

Reply
 
Thread Tools
Old 12-27-2013, 01:02 PM   #1
mvaliadi
Junior Member
 
Location: Germany

Join Date: Dec 2013
Posts: 2
Default Finding a gene of interest in a scaffold

Dear all,

I'm having my first go with working with a newly sequenced genome. I have a bacterial genome sequence in scaffolds assembled from 454 data. I want to quickly find out whether this organism contains some genes of interest (5 in total) without going through the whole annotation process for the entire genome.
I have reference sequences from my 5 interesting genes from closely related organisms. My first thought was to do an alignment but the scaffolds are too long and cutting them in fragments might cut through the sequence of interest.

Any alternatives? Is there an easy way to do this?
mvaliadi is offline   Reply With Quote
Old 12-27-2013, 03:25 PM   #2
boetsie
Senior Member
 
Location: NL, Leiden

Join Date: Feb 2010
Posts: 245
Default

You can align your five genes to the scaffolds with Nucmer (part of MUMmer);

http://mummer.sourceforge.net/
boetsie is offline   Reply With Quote
Old 12-28-2013, 02:28 AM   #3
mvaliadi
Junior Member
 
Location: Germany

Join Date: Dec 2013
Posts: 2
Default

Many thanks! I'll give that a go!
mvaliadi is offline   Reply With Quote
Old 12-28-2013, 06:47 PM   #4
maubp
Peter (Biopython etc)
 
Location: Dundee, Scotland, UK

Join Date: Jul 2009
Posts: 1,543
Default

Or make a BLAST database of your assembly, and search the five genes against that?
maubp is offline   Reply With Quote
Old 12-29-2013, 05:20 PM   #5
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,080
Default

BLAT should be the easiest and quickest solution in a simple case like this. http://genome.ucsc.edu/FAQ/FAQblat.html#blat3
GenoMax is offline   Reply With Quote
Reply

Tags
gene annotation, gene finding, genomics, scaffold

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 04:53 AM.


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