CITE-Seq-Count修剪Fastq文件/运行逗号分隔列表

我正在使用CITE-seq-Count程序来分析单元格哈希/ cite-seq数据。我们在多个泳道(1-4)上对AdT / HTO进行了测序。我尝试将fastq.gz文件与cat函数合并,并按如下所示运行它们:

CITE-seq-Count 
-R1 ./adt_read1/merged_adt_S2_L005_R1_001.fastq.gz 
-R2 ./adt_read2/merged_adt_S2_L005_R2_001.fastq.gz ... (tags,cells,etc).

,我收到错误:./ adt_read2 / merged_adt_S2_L005_R2_001.fastq.gz中的序列长度不一致。请以相同的长度修剪所有序列。退出应用程序。

我看到在1.4.3版本的CITE-seq-Count中可以使用多个通道(而不是合并.fastq.gz文件)

Inputs:
  Required input files.

  -R1 READ1_PATH,--read1 READ1_PATH
                        The path of Read1 in gz format,or a comma-separated list of paths to all Read1 files in gz format (E.g. A1.fq.gz,B1.fq,gz,...
  -R2 READ2_PATH,--read2 READ2_PATH
                        The path of Read2 in gz format,or a comma-separated list of paths to all Read2 files in gz format (E.g. A2.fq.gz,B2.fq,...

但是当我尝试使用下面的代码运行时,它将失败。

CITE-seq-Count
-R1 ./adt_read1/ADT703_S2_L001_R1_001.fastq.gz,./adt_read1/ADT703_S2_L002_R1_001.fastq.gz,./adt_read1/ADT703_S2_L003_R1_001.fastq.gz,./adt_read1/ADT703_S2_L004_R1_001.fastq.gz 

-R2 ./adt_read2/ADT703_S2_L001_R2_001.fastq.gz,./adt_read2/ADT703_S2_L002_R2_001.fastq.gz,./adt_read2/ADT703_S2_L003_R2_001.fastq.gz,./adt_read2/ADT703_S2_L004_R2_001.fastq.gz 

我收到一条错误消息,说这些参数未使用。有没有一种方法可以将所有L001-4车道合并为一个车道?我的逗号分隔列表不正确吗?

谢谢!

felixwoo80 回答:CITE-Seq-Count修剪Fastq文件/运行逗号分隔列表

尝试此命令。要解决此命令,请执行以下操作:删除逗号(foo,bar-> bar,bar)等空白,并在以下位置添加\(反斜杠无尾随空格)跨多行的命令的每一行的末尾。

CITE-seq-Count \
-R1 ./adt_read1/ADT703_S2_L001_R1_001.fastq.gz,./adt_read1/ADT703_S2_L002_R1_001.fastq.gz,./adt_read1/ADT703_S2_L003_R1_001.fastq.gz,./adt_read1/ADT703_S2_L004_R1_001.fastq.gz \
-R2 ./adt_read2/ADT703_S2_L001_R2_001.fastq.gz,./adt_read2/ADT703_S2_L002_R2_001.fastq.gz,./adt_read2/ADT703_S2_L003_R2_001.fastq.gz,./adt_read2/ADT703_S2_L004_R2_001.fastq.gz
本文链接:https://www.f2er.com/3108549.html

大家都在问