nixCraft

Hi all!

I am working on Linux Platform.
I am in very much need to separate out my files according to the information given in another file.I am giving you the format of each file here and than I will desribe what output actually I need.

This is my File no. 1:


Column starting from > mark represent the following:
>Serial number of nucleotides
----->Residue number as in PDB
---------->Residue name
------------->Chain ID
-------------->Serial number of the paired base
------------------->Residue number of the paired base as in PDB
------------------------>Residue name of the paired base
--------------------------->Chain ID of the paired base
----------------------------->Base pair type
---------------------------------->Base pair indicator
------------------------------------->E-value indicating base pair deformation
----------------------------------------->Equivalent information for other pairs
1 601 U R 71 672 A R W:WC BP 0.01
2 602 C R 70 671 G R W:WC BP 0.08
3 603 C R 69 670 G R W:WC BP 0.04
4 604 G R 68 669 C R W:WC BP 0.04
5 605 U R 67 668 G R W:WC BP 0.15
6 606 G R 66 667 C R W:WC BP 0.05
7 607 A R 65 666 U R W:WC BP 0.04
8 608 U R 14 614 A R W:HT BP 0.24
9 609 A R 23 623 A R H:HT BP 0.02
10 610 G R 25 625 U R W:WC BP 0.15
11 611 U R 24 624 A R W:WC BP 0.08
12 612 U R 23 623 A R W:WC BP 0.06
13 613 U R 22 622 G R W:WT BP 0.09
14 614 A R 8 608 U R H:WT BP 0.24
15 615 A R 47 648 U R W:WT BP 0.03
16 616 U R 58 659 U R 2:2T BP 0.59
17 617 G R
18 618 G R 55 656 C R W:WC BP 0.14
19 619 U R
20 620 C R
21 621 A R
22 622 G R 13 613 U R W:WT BP 0.09 46 646 A R H:+T TP 0.15
23 623 A R 9 609 A R H:HT BP 0.02 12 612 U R W:WC TP 0.06
24 624 A R 11 611 U R W:WC BP 0.08
25 625 U R 10 610 G R W:WC BP 0.15
26 626 G R 44 644 A R W:WC BP 0.18
27 627 G R 43 643 C R W:WC BP 0.03
28 628 G R 42 642 C R W:WC BP 0.03
29 629 C R 41 641 G R W:WC BP 0.00
30 630 G R 40 640 U R W:WC BP 0.04
31 631 C R 39 639 G R W:WC BP 0.04

and... here is file 2:

ATOM 3948 O3P U R 601 87.443 56.290 -33.341 1.00 68.37 1ASY4244
ATOM 3949 P U R 601 86.329 55.822 -34.357 1.00 69.36 1ASY4245
ATOM 3950 O1P U R 601 85.172 55.252 -33.617 1.00 70.90 1ASY4246
ATOM 3951 O2P U R 601 86.854 54.828 -35.356 1.00 71.50 1ASY4247
ATOM 3952 O5* U R 601 85.845 57.094 -35.175 1.00 63.25 1ASY4248
ATOM 3953 C5* U R 601 85.151 56.948 -36.433 1.00 55.85 1ASY4249
ATOM 3954 C4* U R 601 84.563 58.263 -36.846 1.00 51.42 1ASY4250
ATOM 3955 O4* U R 601 83.911 58.892 -35.718 1.00 49.27 1ASY4251
ATOM 3956 C3* U R 601 85.698 59.228 -37.191 1.00 49.91 1ASY4252
ATOM 3957 O3* U R 601 85.992 59.082 -38.595 1.00 51.35 1ASY4253
ATOM 3958 C2* U R 601 85.046 60.587 -36.933 1.00 47.11 1ASY4254
ATOM 3959 O2* U R 601 84.247 60.953 -38.052 1.00 46.28 1ASY4255
ATOM 3960 C1* U R 601 84.163 60.295 -35.727 1.00 46.00 1ASY4256
ATOM 3961 N1 U R 601 84.731 60.719 -34.426 1.00 43.00 1ASY4257
ATOM 3962 C2 U R 601 84.878 62.078 -34.176 1.00 42.61 1ASY4258
ATOM 3963 O2 U R 601 84.605 62.926 -35.029 1.00 47.96 1ASY4259
ATOM 3964 N3 U R 601 85.298 62.429 -32.927 1.00 39.71 1ASY4260
ATOM 3965 C4 U R 601 85.614 61.581 -31.911 1.00 40.20 1ASY4261
ATOM 3966 O4 U R 601 85.992 62.005 -30.809 1.00 39.70 1ASY4262
ATOM 3967 C5 U R 601 85.446 60.193 -32.230 1.00 41.95 1ASY4263
ATOM 3968 C6 U R 601 85.025 59.813 -33.444 1.00 42.69 1ASY4264
ATOM 3969 P C R 602 87.506 58.907 -39.120 1.00 58.39 1ASY4265
ATOM 3970 O1P C R 602 88.347 58.249 -38.078 1.00 57.73 1ASY4266


Now, this is just the part of my file.
So, I am interseted to fetch out the whole part from file 2 which is in series as it is there in file1...example 601,602,603.In another words, I m looking for the whole part from file 2: as example if there is "601 U R " in file 1, I need everything from file 2 which comprises of 601 U R and other things like
ATOM 3948 O3P U R 601 87.443 56.290 -33.341 1.00 68.37 1ASY4244.
PLEASE SOLVE THIS PROBLEM FOR ME BY WRITING A CODE AS THERE IS WHOLE LOT OF FILES LIKE THIS.

Thnks in advance.