RefSeq IDの順番にソート

sort_refseq_id.pl

#!/usr/bin/perl

# RefSeqのIDを順番に並べ替える
#
# NM_123456.7
# NM_123456789.12
#
# (1:アルファベット部)_(2:整数部).(3:バージョン番号)
# 1 -> 2 -> 3 の優先順位でソートする
#
# 2023-08-21 Yuki Naito (@meso_cacase)

use warnings ;
use strict ;

print sort compare_refid <> ;

exit ;

sub compare_refid {
	my ($a_moltype, $a_id, $a_ver) = ($a =~ /([A-Z]+)_(\d+)(\.(\d+))?/) ? ($1, $2, $4) : () ;
	my ($b_moltype, $b_id, $b_ver) = ($b =~ /([A-Z]+)_(\d+)(\.(\d+))?/) ? ($1, $2, $4) : () ;
	return $a_moltype cmp $b_moltype || $a_id <=> $b_id || $a_ver <=> $b_ver ;
} ;