#Remove a #sequence by id from multifasta: cat vaga.fa | awk '{if (substr($0,1) == ">scaffold_1 1087316 bp") censor=1; else if (substr($0,1,1) == ">") censor=0; if (censor==0) print $0}' > fixed.fasta
get the fasta genes using gff3 and genome.fa sequence: ~/Tools/bedtools2/bin/bedtools getfasta -fi vaga.fa -bed Adineta_vaga.v2.gff3 -fo Adineta_vaga_genes.fa
#Convert #fasta file to uppercase/lowercase, without altering the annotation line $ awk '{ if ($0 !~ />/) {print toupper($0)} else {print $0} }' name.fasta
#Convert a #multi-line #fasta to a #singleline #fasta $ awk '!/^>/ { printf "%s", $0; n = "\n" } /^>/ { print n $0; n = "" } END { printf "%s", n }' sample1.fa > sample1_singleline.fa