How can I create LASTdb for nr.gz data?
1
0
Entering edit mode
4.0 years ago
MAPK ★ 2.0k

I am trying to do long read alignments for my assembled reads with ncbi protein database. I was trying to use LAST tool available here (http://last.cbrc.jp/doc/last-tutorial.html) using nr.gz database using the command lastdb -cR01 nrdb nr.gz. However, it doesn't seem to work with .gz format. Is there a way I can extract all sequences from nr.gz into fasta format? Or if anyone has better suggestion on what I should be doing to perform this alignment? Thanks for your help in advance.

last alignement • 1.3k views
ADD COMMENT
1
Entering edit mode

You will need to gunzip nr.gz to get the fasta format file. Based on the tutorial that is what lastdb needs.

ADD REPLY
0
Entering edit mode

Thanks for the helpful comment.

ADD REPLY
2
Entering edit mode
4.0 years ago
h.mon 34k

You can stream with zcat:

zcat nr.gz | lastdb -cR01 nrdb -

You can also stream a compressed fasta file to lastal:

 zcat query.fa.gz | lastal nrdb > query.maf
ADD COMMENT

Login before adding your answer.

Traffic: 1465 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