This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Lastz align sequence to self in blocks of 1kb

Hi all,

Sorry if this has been asked before, but I seem not to have found the answer.

I want to do a self-alignment of a 60Mb chromosome using a window size of 1kb.

First, the --self option does not run, so I opt for specifying the same sequence as the query and as the target, so I ended up with the command as follows:

lastz $inFile1[unmask] $inFile1[unmask] \ --notransition --step=20 --nogapped \ --format=general > $outFilePath"/"$sample"_identity_general.txt" The command ran correctly and I ended up with about 81GB size file as output.

Second, is the "--step" option specifying the window size (documentation mentioned seed length, is it the same as window size), and if not, how do I specify a 1kb window size for the alignment (I seem not to be able to locate this on the documentation page).

Thank you.

lastz chromosome self-align

0 answers

No answers yet.

Log in to answer this question.