Я бегу в очень странном случае ... Вчера я написал небольшой сценарий. Его цель - проверить одно условие в файле на основе другого файла. Это сработало как задумано. Но с этого утра это не так. Я ничего не изменил, насколько мне известно. код не выдает никакой ошибки. Я думаю, что виновником является glob.glob
for file in glob.glob('*private.vcf.gz'):
seen = False
vcf = VCF(file)
print("test")
if not (file == "controlH.g.vcf.gz" or file == "output.g.vcf.gz"):
sample_name = file.split('.')[0]
out = "{}.FalsePositiveRefCallPurged.vcf".format(sample_name)
w = Writer(out, vcf)
for v in vcf:
seen = False
ref = VCF('output.g.vcf.gz')
for r in ref:
if seen :
break
if not seen:
if v.CHROM == r.CHROM:
if v.start == r.start or v.start > r.start and v.start < r.end:
if r.FILTER == "RefCall":
continue#print(str(v))
else:
w.write_record(v)
seen = True
w.close()
Действительно, при простом запуске
glob.glob('.*private.vcf.gz')
я получаю пустой список. Вот вывод bash
ls *.private.vcf.gz
D2A1.private.vcf.gz D3A1.private.vcf.gz D5B3.private.vcf.gz H2C3.private.vcf.gz H4C2.private.vcf.gz
D2B3.private.vcf.gz D4A3.private.vcf.gz H2A3.private.vcf.gz H4A4.private.vcf.gz H5A3.private.vcf.gz
Так что я уверен, что файлы есть ... Я действительно не понимаю, почему вдруг glob.glob не может найти их.
Любая помощь будет принята с благодарностью, спасибо