https://github.com/jhbadger/scripts
Raw File
Tip revision: 6ad694835e70be1c38b4cc13806c5b4361f669fc authored by Jonathan Badger on 12 January 2024, 20:33:42 UTC
report misses in virtualPCR
Tip revision: 6ad6948
sortFasta
#!/usr/bin/env ruby

require 'rubygems'
require'bio'
include Bio    

if (ARGV.size != 1)
  STDERR.printf("usage: %s seq.fasta\n", $0)
  exit(1)
end

file = ARGV.shift

seqs = Hash.new

FlatFile.new(FastaFormat, File.new(file)).each {|seq|
  seqs[seq.entry_id] = seq
}

seqs.keys.sort.each {|key|
  print seqs[key]
}
back to top