script to filter contigs >1kb
1
0
Entering edit mode
9.9 years ago
krp0001 ▴ 40

Hi

I have been looking for a script that could filter out the contigs which are >1kb. my assembled file has the following stats. could some one help me with the script to get contigs >1kb. I have been through all the post for the similar script but couldn't find a right one. thanks in advance.

Contig statistics:
Total length:    49295072 (n:27000)
Max:    788230
N20:    149300 (n:40)
N50:    57668 (n:200)
N80:    10314 (n:808)
n > 1kb:    2485
Tot bp in contigs > 1kb:    45136108 (91.563124%)
genome Assembly • 5.8k views
ADD COMMENT
0
Entering edit mode

Another option, using the BBMap package:

reformat.sh in=contigs.fasta out=filtered.fasta minlength=1000

It's very fast. Also supports maxlength.

ADD REPLY
2
Entering edit mode
9.9 years ago
iraun 6.2k

Maybe this post can solve your issue.

ADD COMMENT
0
Entering edit mode

Thanks !

That really helped. perfect!

ADD REPLY
0
Entering edit mode

You're welcome :)

ADD REPLY

Login before adding your answer.

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