修复Clojure脚本转换FASTA至PanSN-spec时的符号解析错误

错误原因 #

这个错误是因为你的脚本将FASTA序列头的起始符号>和后续的contig名称（如Chr01）当成了一个未定义的Clojure符号，而非字符串处理。Clojure会尝试解析该符号，但找不到对应的定义，因此抛出错误。

修复步骤 #

1. 正确处理FASTA头字符串 #

读取FASTA文件时，需将以>开头的整行作为字符串读取，先剥离开头的>，再提取原contig名称，最后拼接成符合PanSN-spec规范的新头：

(defn process-header [header-line sample-name haplotype-id]
  (let [original-contig (subs header-line 1) ; 移除开头的">"
        new-header-content (str sample-name "#" haplotype-id "#" original-contig)]
    (str ">" new-header-content))) ; 新头需重新添加">"

2. 完整的FASTA转换逻辑 #

确保整个文件处理流程中，所有序列头相关内容都以字符串形式处理，避免符号解析错误：

(require '[clojure.java.io :as io])

(defn convert-to-pansn [input-path output-path sample haplotype]
  (with-open [reader (io/reader input-path)
              writer (io/writer output-path)]
    (doseq [line (line-seq reader)]
      (if (.startsWith line ">")
        ;; 处理序列头
        (.write writer (str (process-header line sample haplotype) "\n"))
        ;; 序列行直接写入
        (.write writer (str line "\n"))))))

;; 调用示例：替换为你的参数
(convert-to-pansn "input.fasta" "output_pansn.fasta" "SampleX" "HapA")

3. 验证结果 #

运行脚本后，检查输出FASTA的序列头是否符合规范：

原头示例：>Chr01
转换后示例：>SampleX#HapA#Chr01

内容的提问来源于stack exchange，提问作者user3523406