What the difference between last-split and single2_cov
0
1
Entering edit mode
3.1 years ago

I want to obtained a multiple alignment between some genomes by using TBA. First I used lastal of LAST to make pair-wise alignment of genomes. Then, I used single2_cov of TBA to find the returns only those alignments that are single coverage.

lastal -C2 -E0.05 ~/WORK/TBA/DB/A/A..db ~/DATABASE/D/D.fna \
            |tee|cat <(echo "##maf version=1") -|single_cov2 /dev/stdin > A.D.sing.maf

But, I found the single2_cov readed all data into memory, whose size beyond 100+ Gb. My server does not have such a large memory.

And I found the last-split finds a unique best alignment for each part of each read. It seems that the single2_cov and last-split have similar function. And latter uses less memory.

So can I use last-split instead of single2cov ?

Best wishes!

MSA TBA LAST • 493 views
ADD COMMENT

Login before adding your answer.

Traffic: 1025 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6