Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Sequence size error #186

Open
onbio opened this issue Sep 24, 2024 · 2 comments
Open

Sequence size error #186

onbio opened this issue Sep 24, 2024 · 2 comments

Comments

@onbio
Copy link

onbio commented Sep 24, 2024

Hi,
I am trying to run AliTV for three plant genomes (same species but different sub-species). the genomes have been assembled using HifiAsm. As they are from the same species the sequences have long stretches of similar sequence. the genome size is about 1Gb.

Can you suggest a remedy for the below error? It seems that the long matching sequence is causing the LastZ program to break. One possible fix I can think of now is masking repeats, but it would be great if you could suggest a remedy.

INFO - Created temporary folder at '/tmp/dB3SQ7kwnl'
INFO - Starting alignment generation... (820 alignments required)
FAILURE: in add_segment()
table size (4,869,542,152 for 101,448,794 segments) exceeds allocation limit of 4,294,967,279;
consider raising scoring threshold (--hspthresh or --exact) or breaking your target sequence into smaller pieces
"/bin/lastz" unexpectedly returned exit value 1 at /AliTV-perl-interface/bin/../lib/AliTV/Alignment/lastz.pm line 107.

Thank you for your assistance

@onbio
Copy link
Author

onbio commented Sep 24, 2024

I did try running for just the first sequence/chromosome, but I still got the same error message.

@iimog
Copy link
Member

iimog commented Sep 24, 2024

This is a problem with lastz on large and very similar genomes. In addition to the suggestion in the message

consider raising scoring threshold (--hspthresh or --exact) or breaking your target sequence into smaller pieces

the discussion in lastz/lastz#53 might be helpful.

You can add lastz options to AliTV in the yml file, see the documentation

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants