()
| 53 | return output |
| 54 | |
| 55 | def main(): |
| 56 | if len(sys.argv) < 3: |
| 57 | print("Usage: python3 af2rechain.py pdbfilename.pdb 50 100 320") |
| 58 | return |
| 59 | |
| 60 | pdb_filename = sys.argv[1] |
| 61 | residue_breakpoints = list(map(int, sys.argv[2:])) |
| 62 | pdb_data = parse_pdb(pdb_filename) |
| 63 | rechained_data = rechain_pdb(pdb_data, residue_breakpoints) |
| 64 | |
| 65 | output_pdb_filename = pdb_filename.replace(".pdb","_rechain.pdb") |
| 66 | OUT=open(output_pdb_filename,'w') |
| 67 | |
| 68 | for line in rechained_data: |
| 69 | OUT.write(line) |
| 70 | |
| 71 | if __name__ == "__main__": |
| 72 | main() |
no test coverage detected