SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Lastz help Esther Bioinformatics 23 08-24-2015 01:16 PM
lastz and yasra install problems LVAndrews Bioinformatics 4 07-13-2014 09:35 AM
LASTZ installation as1438 Bioinformatics 14 05-15-2014 07:28 AM
MIRA failure ajthomas Bioinformatics 1 01-29-2014 09:56 AM
lastz parameters Vasilis General 0 12-17-2012 01:00 AM

Reply
 
Thread Tools
Old 03-25-2014, 08:19 PM   #1
Hyunmin
Member
 
Location: Seoul, Korea

Join Date: Feb 2012
Posts: 13
Default Using lastz (FAILURE: in load_fasta_sequence)

This is my lastz command
$ lastz hg19/chrAll.fa[multiple] merged.fsa

Error
FAILURE: in load_fasta_sequence for hg19/chrAll.fa, sequence length 2,084,766,314+115,169,878 exceeds maximum (2,147,483,637)

I got the error in above.

What's wrong thing?

Thanks,
Hyunmin is offline   Reply With Quote
Old 03-26-2014, 03:10 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,077
Default

Have you tried to do the comparison choosing a smaller set of chromosomes from hg19 at a time and checking if that works?
GenoMax is offline   Reply With Quote
Old 03-28-2014, 04:56 AM   #3
Bob-Harris
Member
 
Location: USA

Join Date: Mar 2014
Posts: 13
Default

Quote:
Originally Posted by Hyunmin View Post
This is my lastz command
$ lastz hg19/chrAll.fa[multiple] merged.fsa
FAILURE: in load_fasta_sequence for hg19/chrAll.fa, sequence length 2,084,766,314+115,169,878 exceeds maximum (2,147,483,637)
The short answer is you need to use lastz_32 instead of lastz.

The long answer… the default lastz build can only handle a target up to 2Gbp combined length. This was a conscious design choice (circa 2004), as there is some efficiency gain from limiting the variables that track sequence positions to 31 bits. The lastz_32 build relaxes that constraint so that it can handle a target up to 4Gbp combined length.

The efficiency gain was more important on machines 10 years ago than it is now, but for backward compatibility I am keeping that as the default behavior. I'll (probably) change this error report on future releases so that it tells the user to consider using lastz_32.

The issue is discussed in the readme file here:
http://www.bx.psu.edu/~rsharris/last...v_whole_genome

Bob H
(lastz author)
Bob-Harris is offline   Reply With Quote
Old 05-14-2014, 02:10 AM   #4
Hyunmin
Member
 
Location: Seoul, Korea

Join Date: Feb 2012
Posts: 13
Default

Quote:
Originally Posted by Bob-Harris View Post
The short answer is you need to use lastz_32 instead of lastz.

The long answer… the default lastz build can only handle a target up to 2Gbp combined length. This was a conscious design choice (circa 2004), as there is some efficiency gain from limiting the variables that track sequence positions to 31 bits. The lastz_32 build relaxes that constraint so that it can handle a target up to 4Gbp combined length.

The efficiency gain was more important on machines 10 years ago than it is now, but for backward compatibility I am keeping that as the default behavior. I'll (probably) change this error report on future releases so that it tells the user to consider using lastz_32.

The issue is discussed in the readme file here:
http://www.bx.psu.edu/~rsharris/last...v_whole_genome

Bob H
(lastz author)
sorry, my late thread..

I solve the problem with upper.
Thanks!
Hyunmin is offline   Reply With Quote
Old 07-21-2016, 06:18 AM   #5
dan
wiki wiki
 
Location: Cambridge, England

Join Date: Jul 2008
Posts: 266
Default

I'm trying lastz-1.03.73 from here: http://www.bx.psu.edu/~rsharris/lastz/newer/

But I don't see lastz_32... Is this now the default behaviour?

I just tried 'make lastz_32' and it appeared :-D


Thanks,
__________________
Homepage: Dan Bolser
MetaBase the database of biological databases.

Last edited by dan; 07-21-2016 at 06:20 AM.
dan is offline   Reply With Quote
Reply

Tags
blastz, lastz, multiple alignment

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 08:15 PM.


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