This is a test version of Biostars. For the public version, visit https://www.biostars.org.
LASTZ aligner for HaploMerger

Dear all,

I need to check repeat sequences in my genome assembly. For this, first I run RepeatMasker. Then, I want to run HaploMerger tool. But for HaploMerger, I have to run LASTZ tool. I checked the manual of it, but confused. How to use command tool for this?

Thank you

alignment assembly genome sequence
• 4,040 views • link
updated by Ram Ram • written by Mehmet Mehmet

You should also read the manual ;) what is it that confuses you (in detail)?

• 0 views • link
written by Michael Dondrup Michael Dondrup

Hi, I have already read many times as I mentioned in my message. For instance, my command was:

lastz target.fa (my genome file) no query file. I also used the action "multiple". I received the error Can't understand "multiple". If I don't put the action "multiple", I receive the error "contains more than one sequence". I want to align my genome file with itself.

• 0 views • link
updated by Ram Ram • written by Mehmet Mehmet

Yes I did, but I received this error

FAILURE: multiple action cannot be used with --lav

(lav has requirements on the order of alignments that would require additional

• 0 views • link
updated by Ram Ram • written by Mehmet Mehmet

I haven't tested lastz for a long time, but I didn't say --lav but --self

• 0 views • link
written by Michael Dondrup Michael Dondrup

1 answer

HaploMerger has batch scripts which run lastz for you. The HaploMerger manual (which comes with the distribution) details this.

• 0 views • link
updated by Ram Ram • written by Chris Fields Chris Fields

Log in to answer this question.

More posts like this