pbs scripts

[aoife]

Hi,very new.hoping to incorporate python into my postgrad.

Basically I have 2,000 files.I want to write a script that says:

open each file in turn
for each file:
open this pbs script and run MUSCLE (a sequence alignment tool)
on each file
close this file
move on to next file.

any help would be great.
Aoife

 

[Simon Brunning]

Hi,very new.hoping to incorporate python into my postgrad.

Basically I have 2,000 files.I want to write a script that says:

open each file in turn

If they are in one directory, look at the glob module. If they are in
a bunch of sub-directories, see os.walk(), or <http://bit.ly/5Q5Qiv>.

for each file:
      open this pbs script and run MUSCLE (a sequence alignment tool)
on each file

Is MUSCLE a command-line tool? If so, see the subprocess module.

      close this file

Do you actually need to open the file, or just run a command on it?
Sounds like the latter to me.

      move on to next file.

Give it a go. Any problems, I'm sure we'd be happy to help.

 

[r0g]

Hi,very new.hoping to incorporate python into my postgrad.

Basically I have 2,000 files.I want to write a script that says:

open each file in turn
for each file:
open this pbs script and run MUSCLE (a sequence alignment tool)
on each file
close this file
move on to next file.

any help would be great.
Aoife

Hi Aoife,

import os

for each_pbs in os.listdir("/home/user/pbs_files/"): # loop thru dir
if each_pbs[-4:].upper() == ".PBS": # check extension
os.system("MUSCLE " + each_pbs) # call CLI prog


I don't think you need to open or close the files if their names (As
opposed to their contents) are just going to be parameters to this
MUSCLE program.

If you want it to look in subfolders replace 'os.listdir' with 'os.walk'.

Hope this helps, if not please clarify

Cheers,

Roger.