BOL

Results for "awk"

Wire posts

• • Jit@jit.aber

  

  #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

  2327 days ago

• • Jit@jit.aber

  

  samtools view -f 256 Input.bam | awk '$6 ~/S/ && $7 == "=" {print $0}' > Secondary_clipped.sam

  2287 days ago

• • Jit@jit.aber

  

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

  2274 days ago

• • Jit@jit.aber

  

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

  2267 days ago

• • Rahul Nayak@rahul

  

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

  2228 days ago

• • Rahul Nayak@rahul

  

  #Split #multifasta $ cat hg18.fa | awk '{ if (substr($0, 1, 1)==">") {filename=(substr($0,2) ".fa")} print $0 > filename }'

  2210 days ago

• • Manisha Mishra@manisha

  

  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

  2202 days ago

• • Neelam Jha@neelam

  

  #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

  2146 days ago

• • Jit@jit.aber

  

  #Zero coverage regions ~/Tools/bedtools2/bin/bedtools genomecov -ibam PvsH_aln_sorted.bam -bga | awk '$4==0' > ZERO_coverage_regions.txt

  2137 days ago

• • Jit@jit.aber

  

  #Multiline Fasta To #SingleLine Fasta $ awk '/^>/ {printf("\n%s\n",$0);next; } { printf("%s",$0);} END {printf("\n");}' < file.fa

  2070 days ago

• « Previous
• 1
• 2
• 3
• 4
• 5
• Next »