Мне нужно перебрать пары входных файлов, R1 и R2 , которые будут обрабатываться вместе, и к каждому из них будет добавлено расширение ".trim".
У меня ограниченные знания Python и glob, которые хорошо работали для аналогичных задач, используя парные входы, но один выход. Я полностью озадачен тем, почему это не работает, хотя я уверен, что это легко исправить.
#!/usr/bin/python
import glob
import os
files = glob.glob("ATAC*R1*.fastq.gz")
for ifile in files:
os.system("cutadapt -a CTGTCTCTTATACACATCT -A CTGTCTCTTATACACATCT -a AGATGTGTATAAGAGA -o %s.trim -p %s.trim %s %s" % (ifile, ifile.replace("R1","R2"), ifile, ifile.replace("R1","R2"))
Error:
line 7
^
SyntaxError: invalid syntax
Это должно работать следующим образом:
ATAC2-1_R1_.fastq.gz and ATAC2-1_R1_.fastq.gz -> ATAC2-1_R1_.fastq.gz.trim and ATAC2-1_R1_.fastq.gz.trim.
Есть несколько пар R1 и R2, которые я хотел бы перебрать. Спасибо за помощь!