SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
bwa perfect match only guil Bioinformatics 4 08-30-2013 10:18 AM
Finding Perfect matches and mismatches bhargavpbk88 Bioinformatics 2 01-19-2012 08:00 AM
Perfect 454 FLX for sale (with FREE upgrade to long reads FLX+ available) 454FLXforsale Service Providers 0 01-18-2012 05:38 AM
Perfect 454 FLX for sale (with FREE upgrade to long reads FLX+ available ) 454FLXforsale 454 Pyrosequencing 0 01-18-2012 05:36 AM
Perfect match disagreement between bowtie and BLAT on human genome danielsbrewer Bioinformatics 3 03-10-2009 10:10 PM

Reply
 
Thread Tools
Old 10-06-2010, 06:05 AM   #1
NicoBxl
not just another member
 
Location: Belgium

Join Date: Aug 2010
Posts: 264
Default perfect alignment match

Hi,

I've a little problem with an alignment. I want to know to view the perfect alignment of a query to a database . Here's an example :

query

>query1
ATGC

subject :

>seq1
GATGC
>seq2
CGATGC
>seq3
ATGC

When I use Blast, the three alignment have the same e-value and score, but the best alignment for me is with seq3 ( same size )

is it a blast parameter to change ?

or an another program to use

Thanks,

Nicolas
NicoBxl is offline   Reply With Quote
Old 10-07-2010, 02:33 AM   #2
NicoBxl
not just another member
 
Location: Belgium

Join Date: Aug 2010
Posts: 264
Default

it's strange isn't it ?

nobody knows ?
NicoBxl is offline   Reply With Quote
Old 10-07-2010, 04:50 AM   #3
epigen
Senior Member
 
Location: Germany

Join Date: May 2010
Posts: 101
Default

Blast is a local alignment tool, so it will only consider the aligned part, in you example the ATGC, which is the same in all 3 sequences. There are 4 perfect matches, giving the same score for all 3. The e value is calcuated from the length of the match and the score - also the same. At no point the length of the whole target sequence is taken into account. That means to get your "best alignment" for the same size, you have to filter the results by target sequence length, probably by writing a small script. I don't know of any Blast parameter that could be useful.
epigen is offline   Reply With Quote
Old 10-07-2010, 05:28 AM   #4
NicoBxl
not just another member
 
Location: Belgium

Join Date: Aug 2010
Posts: 264
Default

ok I'll write a perl script

thanks
NicoBxl 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 03:30 AM.


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