6 GI:50845406 NM_031444.2 snoRNA3 -9 Box D except for snoRNA4 Query 3 CTGGAGTCAAGGCT 16 |||||||||||||| Sbjct 1297 CTGGAGTCAAGGCT 1284 Homo sapiens chromosome 22 open reading frame 13 (C22orf13), mRNA. http://www.ncbi.nlm.nih.gov/sites/entrez?cmd=Retrieve&db=nucleotide&dopt=GenBank&RID=UEDK8KR6016&log%24=nuclalign&blast_rank=13&list_uids=50845406 5 GI:38327560 NM_006282.2 snoRNA3 -9 Box D except for snoRNA4 Query 5 GGAGTCAAGGCTAC 18 |||||||||||||| Sbjct 5129 GGAGTCAAGGCTAC 5116 Homo sapiens serine/threonine kinase 4 (STK4), mRNA http://www.ncbi.nlm.nih.gov/sites/entrez?cmd=Retrieve&db=nucleotide&dopt=GenBank&RID=UEDK8KR6016&log%24=nuclalign&blast_rank=14&list_uids=38327560 3 GI:91982771 NM_001040105.1 snoRNA 10 Query 4 TGGAGTCAAT 13 |||||||||| Sbjct 4854 TGGAGTCAAT 4845 Homo sapiens mucin 17, cell surface associated (MUC17), mRNA. http://www.ncbi.nlm.nih.gov/sites/entrez?cmd=Retrieve&db=nucleotide&dopt=GenBank&RID=UDU305DZ01N&log%24=nuclalign&blast_rank=97&list_uids=91982771 3 GI:154448895 NM_001100162.1 snoRNA 25, 26 and 27 Query 2 CCTGGAGTCGAGTG 15 |||||||||||||| Sbjct 146 CCTGGAGTCGAGTG 133 Homo sapiens exportin 7 (XPO7), transcript variant 3, mRNA. http://www.ncbi.nlm.nih.gov/sites/entrez?cmd=Retrieve&db=nucleotide&dopt=GenBank&RID=UDW41RSS01S&log%24=nuclalign&blast_rank=2&list_uids=154448895 31 GI:153945877 NM_002458.1 snoRNA 25, 26 and 27 Query 3 CTGGAGTCGAGTG 15 ||||||||||||| Sbjct 6818 CTGGAGTCGAGTG 6806 Query 3 CTGGAGTCGAGTG 15 ||||||||||||| Sbjct 8489 CTGGAGTCGAGTG 8477 Query 3 CTGGAGTCGAGTG 15 ||||||||||||| Sbjct 10589 CTGGAGTCGAGTG 10577 Query 3 CTGGAGTCGAGTG 15 ||||||||||||| Sbjct 12260 CTGGAGTCGAGTG 12248 Homo sapiens mucin 5B, oligomeric mucus/gel-forming (MUC5B), mRNA. http://www.ncbi.nlm.nih.gov/sites/entrez?cmd=Retrieve&db=nucleotide&dopt=GenBank&RID=UDW41RSS01S&log%24=nuclalign&blast_rank=9&list_uids=153945877 #### 'snoRNA 25,26 and 27'=>[ { GI=>'GI:15444889', NM=>'NM_001100162.1', exon=>'3', seq=>[ Query 2 CCTGGAGTCGAGTG 15 |||||||||||||| Sbjct 146 CCTGGAGTCGAGTG 133 ], geneName =>'Homo sapiens mucin 17, cell surface associated (MUC17), mRNA.', weblink=>'http://......' }, { GI=>'GI:153945877', NM=>'NM_002458.1' exon=>'31', seq=>[ #more than one seq ] }, ], 'snoRNA3 -9 Box D except for snoRNA4'=>[ #more than record once again ....... ] #### use strict; use Data::Dump qw[ pp ]; my %records; until(eof(DATA)){ chomp(my $snoRNA = ); push @{$records{$snoRNA}},{} ; my $seqs = 1; my $line = ; if( $line =~ m[(\d+) different hits] ) { $seqs = $1; chomp( $records{ $snoRNA }[ -1 ]{ exon } = ); } else { chomp( $records{ $snoRNA }[ -1 ]{ exon } = $line ); } chomp( $records{ $snoRNA }[ -1 ]{ GeneID } = ); chomp( $records{ $snoRNA }[ -1 ]{ NM_ID } = ); for( 1 .. $seqs ) { chomp( my $query = ); scalar (); chomp( my $sbjct = ); push @{ $records{ $snoRNA }[ -1 ]{ seqs } }, { $query => $sbjct }; } chomp( $records{ $snoRNA }[ -1 ]{ gene_name } = ); chomp( $records{ $snoRNA }[ -1 ]{ web_link } = ); } pp \%records; __DATA__ snoRNA 25, 26 and 27 2 GI:142387131 NM_006299.3 Query 2 CCTGGAGTCGAGT 14 ||||||||||||| Sbjct 371 CCTGGAGTCGAGT 359 Homo sapiens zinc finger protein 193 (ZNF193), mRNA. http://www.ncbi.nlm.nih.gov/sites/entrez?cmd=Retrieve&db=nucleotide&dopt=GenBank&RID=UDW41RSS01S&log%24=nuclalign&blast_rank=11&list_uids=142387131 snoRNA 25, 26 and 27 1 NM_001005236.3 GI:256773198 Query 3 CTGGAGTCGAGTGTCT 18 |||||| ||||||||| Sbjct 168 CTGGAGACGAGTGTCT 153 Homo sapiens olfactory receptor, family 1, subfamily L, member 1 (OR1L1), mRNA. http://www.ncbi.nlm.nih.gov/ snoRNA 25, 26 and 27 4 different hits 31 GI:153945877 NM_002458.1 Query 3 CTGGAGTCGAGTG 15 ||||||||||||| Sbjct 6818 CTGGAGTCGAGTG 6806 Query 3 CTGGAGTCGAGTG 15 ||||||||||||| Sbjct 8489 CTGGAGTCGAGTG 8477 Query 3 CTGGAGTCGAGTG 15 ||||||||||||| Sbjct 10589 CTGGAGTCGAGTG 10577 Query 3 CTGGAGTCGAGTG 15 ||||||||||||| Sbjct 12260 CTGGAGTCGAGTG 12248 Homo sapiens mucin 5B, oligomeric mucus/gel-forming (MUC5B), mRNA. http://www.ncbi.nlm.nih.gov/sites/entrez?cmd=Retrieve&db=nucleotide&dopt=GenBank&RID=UDW41RSS01S&log%24=nuclalign&blast_rank=9&list_uids=153945877 ####