Question

How to extract from a fasta file, all the protein sequences that have more than 800 aa?

0

Entering edit mode

6.1 years ago

ayixon • 0

How to extract from a file all the protein sequences that have more than 800 aa using linux terminal?

sequence Assembly • 1.5k views

ADD COMMENT • link updated 6.1 years ago by cpad0112 21k • written 6.1 years ago by ayixon • 0

0

Entering edit mode

Solved. Thank you very much. The perl script worked fine.

ADD REPLY • link 6.1 years ago by ayixon • 0

2

Entering edit mode

6.1 years ago

cpad0112 21k

For future reference, to filter a multifasta file for sequences with more than 800 aa,

$ seqkit seq -m 800 test.fa

-m = minimum length of the sequence to be printed

ADD COMMENT • link 6.1 years ago by cpad0112 21k

score 2 · Accepted Answer · 2018-04-20

2

Entering edit mode

6.1 years ago

GenoMax 142k

ADD COMMENT • link 6.1 years ago by GenoMax 142k