https://github.com/jhbadger/scripts
Tip revision: 6ad694835e70be1c38b4cc13806c5b4361f669fc authored by Jonathan Badger on 12 January 2024, 20:33:42 UTC
report misses in virtualPCR
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]
}