Skip to content

Instantly share code, notes, and snippets.

@yuifu

yuifu/R one-liner for calculating seq lengths in FASTA

Created February 24, 2016 11:14
Show Gist options
  • Star 1 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save yuifu/c769564d37e938f00633 to your computer and use it in GitHub Desktop.
Save yuifu/c769564d37e938f00633 to your computer and use it in GitHub Desktop.
Download ZIP
Raw
R one-liner for calculating seq lengths in FASTA
R --slave -e 'library(Biostrings); seqs <- readDNAStringSet("test.fa"); write.table(data.frame(name=names(seqs), length=width(seqs)),stdout(), quote=F, row.names=F)' > hoge.txt
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment