Hi everyone, I performed a local blastp with a query protein vs 5 genome (ORF data ). Now i want to extract all the top 5 sequence automatically and make a final top hit file through awk/perl/shell. Can any one please help me out how to do it.
# BLASTP 2.2.25 [Feb-01-2011]
# Query: gnl|ECOLI|EG10158-MONOMER ClpP 455901..456524 Escherichia coli K-12 substr. MG1655
# Database: mybd
# Fields: Query id, Subject id, % identity, alignment length, mismatches, gap openings, q. start, q. end, s. start, s. end, e-value, bit score
gnl|ECOLI|EG10158-MONOMER V161_10_length_174349_cov_20.0403_ID_19_63 100.00 207 0 0 1 207 1 207 9e-157 431
gnl|ECOLI|EG10158-MONOMER RS374_6_length_260428_cov_19.5118_ID_11_179 100.00 207 0 0 1 207 1 207 9e-157 431
gnl|ECOLI|EG10158-MONOMER H127_1_length_711606_cov_31.5057_ID_1_431 100.00 207 0 0 1 207 1 207 9e-157 431
gnl|ECOLI|EG10158-MONOMER sff_11_length_151555_cov_28.2432_ID_21_27 100.00 207 0 0 1 207 1 207 9e-157 431
gnl|ECOLI|EG10158-MONOMER H75_1_length_711552_cov_31.2637_ID_1_431 100.00 207 0 0 1 207 1 207 9e-157 431
gnl|ECOLI|EG10158-MONOMER H75_26_length_34692_cov_38.2655_ID_51_32 28.24 131 90 2 73 202 63 190 3e-11 60.8
gnl|ECOLI|EG10158-MONOMER RS374_42_length_11164_cov_52.9994_ID_83_6 28.24 131 90 2 73 202 6 ...
