Results for "fasta"

Wire posts

  • Sort FASTA sequence by length: cat f1.fasta| awk '{if($0 ~ /^>/){printf "%st", $0}else{print length($0)"t"$0}}'| sort -k2n,2| awk '{print $1"n"$3}' > sequence.fasta #Sort #Fasta #Trick #Onliner

    3185 days ago

  • 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

    3158 days ago

  • Remove the fasta description and just keep the identifier. perl -p -i -e 's/>(.+?) .+/>$1/g' sample1.fa #Perl #Oneliner #Trick #Tips #Fasta #Identifier #Remove

    3158 days ago

  • Delete space from FASTA header: perl -pe 's/>\S+\K.+//' < file.fasta > new_file.fasta #PerlOneLiner #Trick #Perl #Space #Fasta

    2976 days ago

  • #BAM file to fasta #conversion with #Samtool: samtools view filename.bam | \ awk '{OFS="\t"; print ">"$1"\n"$10}' - > filename.fasta

    2962 days ago

  • Rename the fasta name perl -ane 'if(/\>/){$a++; @foo = split /\s+/;$_ = join '_', @foo; print "$_\n"}else{print;}' your.fasta > new.fa #PerlTrick #PerlOneliner #Perl

    2786 days ago

  • Extract fasta sequence by Ids: perl -ne 'if(/^>(\S+)/){$c=grep{/^$1$/}qw(id1 id2)}print if $c' fasta.file #Perl #PerlOneliner #PerlTrick #Extract #Ids #Fasta

    2785 days ago

  • Extract sequence by using Ids file (assuming one ids in a line): perl -ne 'if(/^>(\S+)/){$c=$i{$1}}$c?print:chomp;$i{$_}=1 if @ARGV' ids.file fasta.file #PerlTricks #PerlOneliner #Perl #Extract #Ids #

    2785 days ago

  • A perl script for mapping fasta formated sequences to multiple reference sequences using one of several alignment programs. http://seq.crg.es/download/software/Miro/doc/run_Mapping.html #Perl #Script

    2767 days ago

  • Extract a fasta sequence from multifasta file: samtools faidx multitest.fa and then samtools faidx multitest.fa fasta_id > out.fa #Samtools #Fasta #Extract #Multifasta

    2745 days ago