bug in prepare_reference.sh?

[onecarbon]

Hi

Thanks for the software!

I was unable to get prepare_reference.sh to finish without removing the extra "" in the sed calls. Is this meant to be there?

[KforKuma]

you are right, the original code returned me this:
sed: 无法读取 s/>HLA:HLA00037/>HLA:HLA00037 A*03:01:01:01 1098 bp/g: 没有那个文件或目录 sed: 无法读取 s/>HLA:HLA00132/>HLA:HLA00132 B*07:02:01:01 1089 bp/g: 没有那个文件或目录
but after delete the extra "" it works.

[agneantanaviciute]

Deleting the extra "" worked for me as well - thank you for the fix suggestion!