script to filter contigs >1kb

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

• 6,446 views • link

updated 4.3 years ago by Ram

• written 11.5 years ago by krp0001 •

Another option, using the BBMap package:

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

It's very fast. Also supports maxlength.

• 0 views • link

updated 4.3 years ago by Ram

• written 11.5 years ago by Brian Bushnell

1 answer

Maybe this post can solve your issue.

• 0 views • link

updated 4.3 years ago by Ram

• written 11.5 years ago by iraun

Thanks !

That really helped. perfect!

• 0 views • link

updated 6.8 years ago by Ram

• written 11.5 years ago by krp0001 •

You're welcome :)

• 0 views • link

updated 4.3 years ago by Ram

• written 11.5 years ago by iraun

More posts like this

RNA de novo assembly - blasts - KEGG - GO
written 6.4 years ago by dimitrischat

Hello, I am a phd candidate to bioninformatics and with (almost) 0 guidance. Seeking help here.. I was asked to do a de novo RNA …
Assembly statistics: ABySS & BBMap
written 6.8 years ago by el97004 •

Hi! I have noticed some differences in resulting assembly statistics from Abyss and BBMap stats.sh and was wondering if anyone knew why. For example, this …
Problem with "split_and_run_sparc.sh" from DBG2OLC pipeline
written 9.9 years ago by Josué Barrera •

Hi everybody! I'm having a problem in the consensus stage of the DBG2OLC pipeline. I'm using the script "split_and_run_sparc.sh" to obtain the "final_assembly.fasta" file from …
Using the ABySS pipeline to scaffold a DISCOVAR de novo assembly
written 10.7 years ago by dbrowne.up •

Hi all, I have a DISCOVAR de novo assembly that was made with 2x250 bp paired end reads with a fragment size of ~800. Since …
Cleaning up trinity assemblies
written 10.9 years ago by dpearton •

Hi, I have performed de novo assemblies on mRNA using trinity with PE illumina data (125bp PE). Trinity gives a large number of contigs (around …
Scaffolding with SSPACE returned scaffold the same as initial contigs set
written 10.9 years ago by pbigbig

Hi all, I have used Minia to assemble a contigs set from my paired end reads, as Minia instructed that it doesn't use pairing information …
i got two different trinity.fasta by using two version of trinity
written 11.2 years ago by Kurban

Hi guys, we have RNA-seq data sequenced of an insect in 2012, and assembled them by using one of the Trinity 2011 versions at the …
Fragmented ABySS assembly
written 11.2 years ago by zgayk •

Hello, We are trying to assemble the genome of the common loon, and I have used abyss (v. 1.5.2) to produce de novo assemblies with …
How To Interpret Abyss Stats
written 14.0 years ago by lin.barnum •

I ran abyss on a dataset and obtained the following stats. What do the various columns mean? <pre><code>n n:200 n:N50 min N80 N50 N20 max …
Abyss : Explain The K-Mer To Me!
written 14.2 years ago by madkitty

Hey guys, I ran a de novo assembly with ABYSS and I'm struggling to find the right value of k-mer. The average read length …