Entering edit mode
5.9 years ago
Janey
▴
30
Hi
I use the following script to extract sequence from fasta file by samtools. How can I write similar command to remove the sequences by sequence IDs list from fasta file by samtools.
cut -c 2- ID.text | xargs -n 1 samtools faidx in.fasta > out.fasta
Thanks for your help