Beginner’s #guide to #comparative #bacterial genome #analysis using #next-generation sequence data https://microbialinformaticsj.biomedcentral.com/articles/10.1186/2042-5783-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
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