#! perl -slw use strict; my( $ID, $query, $off ); while( ) { if( m[^>] ) { ( $ID ) = (split '\|', $_)[ 1 ]; next; } if( m[^Query] ) { ( $query ) = m[^Query\s+(\d+)]; my $top = substr $_, 15; my $pipes = substr , 15; my $bot = substr , 15; my $p = 0; while( $p = 1+index $pipes, ' ', $p ) { printf "%20s :(%d) %1s/%1s\n", $ID, $query+$p, substr( $top, $p-1, 1 ), substr( $bot, $p-1, 1 ); } } } __DATA__ #### >lcl|14079 ref|NC_000009.11|:4900000-5300000 Homo sapiens chromosome 9, GRCh37 primary reference assembly Length=400001 Score = 270 bits (146), Expect = 2e-74 Identities = 148/149 (99%), Gaps = 0/149 (0%) Strand=Plus/Minus Query 1 TGGGCAAGGACTTCATGTCTAAAACACCAAAAGCAATGGCAACAAAAGCCAAAATTGACA 60 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Sbjct 48784 TGGGCAAGGACTTCATGTCTAAAACACCAAAAGCAATGGCAACAAAAGCCAAAATTGACA 48725 Query 61 AATGGGATCTAATTAAACTAAAGAGCTTCTGCACAGCAAAAGAAACTACCATCAGAGTGA 120 |||||||||||||| ||||||||||||||||||||||||||||||||||||||||||||| Sbjct 48724 AATGGGATCTAATTCAACTAAAGAGCTTCTGCACAGCAAAAGAAACTACCATCAGAGTGA 48665 Query 121 ACAGGCAACCTACAGAATGGGAGAACATT 149 ||||||||||||||||||||||||||||| Sbjct 48664 ACAGGCAACCTACAGAATGGGAGAACATT 48636 >lcl|14080 ref|NC_000009.11|:4900000-5300000 Homo sapiens chromosome 9, GRCh37 primary reference assembly Length=400001 Score = 270 bits (146), Expect = 2e-74 Identities = 148/149 (99%), Gaps = 0/149 (0%) Strand=Plus/Minus Query 1 TGGGCAAGGACTTCATGTCTAAAACACCAAAAGCAATGGCAACAAAAGCCAAAATTGACA 60 |||||||||||||| ||||||||||||||||||||||||||||||||||||||||||||| Sbjct 48784 TGGGCAAGGACTTCATGTCTAAAACACCAAAAGCAATGGCAACAAAAGCCAAAATTGACA 48725 Query 61 AATGGGATCTAATTAAACTAAAGAGCTTCTGCACAGCAAAAGAAACTACCATCAGAGTGA 120 ||||||||||||||||||||||||||||| |||||||||||||||||||||||||||||| Sbjct 48724 AATGGGATCTAATTCAACTAAAGAGCTTCTGCACAGCAAAAGAAACTACCATCAGAGTGA 48665 Query 121 ACAGGCAACCTACAGAATGGGAGAACATT 149 ||| |||| |||||| ||||||||||| Sbjct 48664 ACAGGCAACCTACAGAATGGGAGAACATT 48636 >lcl|14081 ref|NC_000009.11|:4900000-5300000 Homo sapiens chromosome 9, GRCh37 primary reference assembly Length=400001 Score = 270 bits (146), Expect = 2e-74 Identities = 148/149 (99%), Gaps = 0/149 (0%) Strand=Plus/Minus Query 1 TGGGCAAGGACTTCATGTCTAAAACACCAAAAGCAATGGCAACAAAAGCCAAAATTGACA 60 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Sbjct 48784 TGGGCAAGGACTTCATGTCTAAAACACCAAAAGCAATGGCAACAAAAGCCAAAATTGACA 48725 Query 61 AATGGGATCTAATTAAACTAAAGAGCTTCTGCACAGCAAAAGAAACTACCATCAGAGTGA 120 |||||||||||||| ||||||||||||||||||||||||||||||||||||||||||||| Sbjct 48724 AATGGGATCTAATTCAACTAAAGAGCTTCTGCACAGCAAAAGAAACTACCATCAGAGTGA 48665 Query 121 ACAGGCAACCTACAGAATGGGAGAACATT 149 ||||||||||||||||||||||||||||| Sbjct 48664 ACAGGCAACCTACAGAATGGGAGAACATT 48636 #### c:\test\blast>..\853819.pl 14079 ref :(75) A/C 14080 ref :(15) A/A 14080 ref :(90) T/T 14080 ref :(124) G/G 14080 ref :(129) C/C 14080 ref :(136) A/A 14081 ref :(75) A/C