3
« on: February 01, 2024, 10:58:43 am »
I have run the find_pair program combined with analyze on a PDB file of the single stranded DNA of this sequence GCGCGAATAAACGCGC. The exact command that I run was this: find_pair AATAAA_dna.pdb | analyze and I got some output files from which I believe that AATAAA_dna.out is the main output with some base pair information. It contains a section like this that is somewhat similar to .ct file:
****************************************************************************
RMSD of the bases (----- for WC bp, + for isolated bp, x for helix change)
Strand I Strand II Helix
1 (0.022) ....>-:...1_:[.DG]G-----C[.DC]:..16_:-<.... (0.021) |
2 (0.024) ....>-:...2_:[.DC]C-----G[.DG]:..15_:-<.... (0.044) |
3 (0.042) ....>-:...3_:[.DG]G-----C[.DC]:..14_:-<.... (0.025) |
4 (0.025) ....>-:...4_:[.DC]C-----G[.DG]:..13_:-<.... (0.024) |
5 (0.024) ....>-:...5_:[.DG]G-----C[.DC]:..12_:-<.... (0.030) |
6 (0.029) ....>-:...8_:[.DT]T-**+-A[.DA]:..11_:-<.... (0.020) |
7 (0.030) ....>-:...7_:[.DA]A-*---A[.DA]:..10_:-<.... (0.023) |
Note: This structure contains 2[1] non-Watson-Crick base-pairs.
****************************************************************************
I would like to get the output in the .ct format like this:
1 G 0 2 16 1
2 C 1 3 15 2
3 G 2 4 14 3
4 C 3 5 13 4
5 G 4 6 12 5
6 A 5 7 0 6
7 A 6 8 0 7
8 T 7 9 0 8
9 A 8 10 0 9
10 A 9 11 0 10
11 A 10 12 0 11
12 C 11 13 5 12
13 G 12 14 4 13
14 C 13 15 3 14
15 G 14 16 2 15
16 C 15 0 1 16
or in dot-bracket format like this:
(((((......)))))