Skip to content

Instantly share code, notes, and snippets.

@yfu
Last active August 29, 2015 14:09
Show Gist options
  • Select an option

  • Save yfu/230ad83dd8f125b2928b to your computer and use it in GitHub Desktop.

Select an option

Save yfu/230ad83dd8f125b2928b to your computer and use it in GitHub Desktop.
Run Multiz for SOAPdenovo and Velvet assemblies
cat /data/fuy2/cl/results/2014-11-06/Hi5_soapdenovo_69mer.fa | tr ' ' '_' | gawk '{ if($1 ~ "^>") { match($1, "^>(.+)", array); printf "\n" array[1] "\t"} else { printf $1 } }' | grep -E -v '^$' | awk '{ print ">sd69:" $1 ":1:+:" length($2); print $2 }' > sd69
cat /data/fuy2/cl/results/2014-11-06/Hi5_velvet_69mer.fa | gawk '{ if($1 ~ "^>") { match($1, ">NODE_([0-9]+).+", array); printf "\n" array[1] "\t" } else{printf $1} }' | grep -v -E '^$' | awk '{ print ">vv:" $1 ":1:+:"length($2); print $2 }' > vv69
all_bz - '(sd69 vv69)' > all_bz.log
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment