我對python相當陌生,而且我的python腳本(split_fasta.py)有問題。這裏是我的問題的一個例子:在python中使用變量作爲新文件的名稱的一部分
list = ["1.fasta", "2.fasta", "3.fasta"]
for file in list:
contents = open(file, "r")
for line in contents:
if line[0] == ">":
new_file = open(file + "_chromosome.fasta", "w")
new_file.write(line)
我已經離開了程序的底部部分,因爲它不需要。我的問題是,當我在同一個direcoty我的fasta123文件運行這個程序,它的偉大工程:
蟒蛇split_fasta.py * .fasta
但是,如果我在不同的目錄我希望程序輸出新文件(例如1.fasta_chromsome.fasta)到我當前的目錄中......它不會:
python /home/bin/split_fasta.py/home/data /*.fasta
這仍會在與fasta文件相同的目錄中創建新文件。在這裏,我敢肯定,問題是這一行:
new_file = open(file + "_chromosome.fasta", "w")
因爲如果我把它改成這樣:
new_file = open("seq" + "_chromosome.fasta", "w")
它創造了我的當前目錄下的輸出文件。
我希望這對你們中的一些人有意義,我可以得到一些建議。
您應該結賬BioPython –