Filtering 50 biggest contigs
Hi everyone,
After performing a novo assembly, with the followed command, i want to filter the biggest 50 contigs.
spades.py --untrusted-contigs lclav_genome.fa -1 randomnietnfectedFP.fastq.gz -2 randomnietinfectedRP.fastq.gz -t 2 -m 28 NINnovo --phred-offset 33
This creates a directory with the contigs in it. This file is sorted from biggest to smallest, and we want just the top 50 of these contigs. I have tried something with awk, but it is not working how i want it to. CAn someone help me out?
Thank you!
• 1,555 views
•
link
written
by
stan.aanhane •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
An alternative to haplotype phasing - assembly?
written by Jess •
-
spades.py Error 255
written by Scottie •
-
Rsubread align maximum nthreads
written by Rajesh •
-
Varying assembly statistics: Abyss & SPades
written by max_19
-
E coli Illumina paired assembly using a reference genome?
written by lg248 •
-
De-novo assembly using SPADES
written by ram.mainali •
-
Mauve on Ubuntu 16
written by marongiu.luigi
-
RNA sequencing depth for downstream proteomics
written by benaneely
-
Mauvealigner Error
written by HG
-
Abyss : Explain The K-Mer To Me!
written by madkitty
See past threads for inspiration:
How To Filter Multi Fasta By Length??
how to rearrange fasta file according to its length (add a filter for 50)
You could convert multiline fasta to single line using Multiline Fasta To Single Line Fasta and then extract the first 100 lines using
headwhich should extract the top 50 contigs for you.Don't forget to change them back to fasta format.