Extract sequences by their ID from a fasta file. perl -ne 'if(/^>(\S+)/){$c=$i{$1}}$c?print:chomp;$i{$_}=1 if @ARGV' ids.txt sample1.fa #Extract #Sequences #IDs #Fasta #Perl #Trick #Oneliner
#Mappy provides a convenient #interface to #minimap2, a fast and accurate #C program to #align #genomic and transcribe nucleotide #sequences. https://pypi.org/project/mappy/