#Algorithms for computing parsimonious #evolutionary scenarios for genome evolution: https://bmcevolbiol.biomedcentral.com/articles/10.1186/1471-2148-3-2
#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