Results for "fasta"

Wire posts

  • #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

    2331 days ago

  • 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

    2298 days ago

  • #GFA to #FASTA file conversion. awk '/^S/{print ">"$2"\n"$3}' in.gfa | fold > out.fa

    2278 days ago

  • #GFA to #Fasta file conversion. $ grep '^S' example.gfa |awk '{print ">"$2; print $3}' > example.fa

    2271 days ago

  • #Convert #fasta file to uppercase/lowercase, without altering the annotation line $ awk '{ if ($0 !~ />/) {print toupper($0)} else {print $0} }' name.fasta

    2232 days ago

  • Convert lowercase letter to N:$ sed -e '/^>/! s/[[:lower:]]/N/g' in.fasta > out.fasta #Convert #Lowercase #Fasta #N

    2215 days ago

  • Extract fasta with Ids. $ perl -ne 'if(/^>(\S+)/){$c=grep{/^$1$/}qw(id1 id2)}print if $c' fasta.file #extract #fasta #sequence #ids

    2206 days ago

  • Remove the sequence by fasta ids. awk 'BEGIN{while((getline<"ids.txt")>0)l[">"$1]=1}/^>/{f=!l[$1]}f' seq.fa #Remove #Delete #ids #fasta

    2206 days ago

  • #Convert a #multi-line #fasta to a #singleline #fasta $ awk '!/^>/ { printf "%s", $0; n = "\n" } /^>/ { print n $0; n = "" } END { printf "%s", n }' sample1.fa > sample1_singleline.fa

    2151 days ago

  • Convert fasta to 2bit $ faToTwoBit myAssembly.fa myAssembly.2bit #Fasta #2bit #Convert

    2179 days ago