Skip to content

Instantly share code, notes, and snippets.

@naphthalene

naphthalene/format_protein.sh

Created January 31, 2016 16:48
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save naphthalene/9459d3824dc3eb6d0d4a to your computer and use it in GitHub Desktop.
Save naphthalene/9459d3824dc3eb6d0d4a to your computer and use it in GitHub Desktop.
Download ZIP
Raw
format_protein.sh
efetch -db protein "ybfE" -format gpc -mode xml | xtract -insd complete CDS coded_by | awk '{match($0,"([A-Z0-9\\.]+)\\t([A-Z0-9\\.]+):([0-9]+)\\.\\.([0-9]+)",a)}END{print a[1], a[2], a[3], a[4]}' | tee output.txt
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment