Perl script to extract a sequence from multifasta with range !
...mp; next unless /\w/; s/>$//gs; my @chunk = split /\n/; my $header = shift @chunk; my $seqlen = length join "", @chunk; print ">...2081 days ago
Perl script to run awk inside perl
...r while () { chomp; # strip record separator @Fld = split(' ', $_, -1); if ($awk) { $Fld[(1)-1] = ''; $Fld[(2)-1] = ''; print join($,,@Fld); } print join($,,@Fld) if $f1 . . $txt; }1178 days ago
bash script to extract sequence by ids !
Use a Perl one-liner, grep and seqtk subseq to extract the desired fasta sequences: # Create test input: cat > in.fasta BGI_novel_T016313 Solyc03g025570.2.1 TTC...838 days ago