Extracting sequence less than N length from FASTA file
How do I extract sequences less than 800 bp from a fasta file, I used the following command to extract but I was not successful in doing so. I still have some sequences more than 800 bp in length.
command: $ awk 'BEGIN {RS = ">" ; ORS = ""} length ($2) <= 800 {print ">"$0}' input.fasta > output.fasta
• 1,894 views
•
link
written
by
atitparajuli2018 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Quast contig values not matching
written by jamie.pike
-
Catching error "QueryKey value not found in link output" when fetching taxonomy from NCBI
written by Frieda •
-
Extract sequences consecutively from a large fasta files and save in new fasta files
written by Kumar
-
Tutorial: Calculation of TSS Scores for ATAC-seq Data
written by Sael •
-
Forum: Linearize fasta files
written by finswimmer
-
grep/awk issue: can't extract sequences from a fasta file
written by benoit.kunath •
-
Extracting specific IDs + sequence from multifasta
written by maximilian.mayerhofer •
-
Calculate length of sequences from fasta file
written by KVC_bioinfo
-
Blastn: Blast short reads against genomes and retain only good hits
written by Tamandua •
-
interpretation Abyss stats
written by milanlove •
Hello atitparajuli2018!
Questions similar to yours can already be found at:
We have closed your question to allow us to keep similar content in the same thread.
If you disagree with this please tell us why in a reply below. We'll be happy to talk about it.
Cheers!