使用以下命令时出现的报错:

cellranger count --fastqs=${fastq_dir} \
        --transcriptome=${ref_dir}/Human_index \
        --localcores=30 \
        --sample=file1,file2 \
        --id=disease

报错原因是我的fa文件命名方式为:

file1_1.fastq.gz, file1_2.fastq.gz
file2_1.fastq.gz, file2_2.fastq.gz

但是 cellranger 要求的 fastq 标准命名方式为:

[Sample Name]_S1L00[Lane Number]_[Read Type]_001.fastq.gz
Where Read Type is one of:
I1: Sample index read (optional)
I2: Sample index read (optional)
R1: Read 1
R2: Read 2

因此只需要把file1_1.fastq.gz, file1_2.fastq.gz,file2_1.fastq.gz, file2_2.fastq.gz 改为 cellranger 要求的 fastq 标准命名方式即可。

修改以后变为:
file1_S1_L001_R1_001.fastq.gz, file1_S1_L001_R2_001.fastq.gz
file2_S1_L001_R1_001.fastq.gz, file2_S1_L001_R2_001.fastq.gz

重命名后再跑以下脚本就没有问题了

cellranger count --fastqs=${fastq_dir} \
        --transcriptome=${ref_dir}/Human_index \
        --localcores=30 \
        --sample=file1,file2 \
        --id=disease

原文地址:http://www.cnblogs.com/chenwenyan/p/16924743.html

1. 本站所有资源来源于用户上传和网络,如有侵权请邮件联系站长! 2. 分享目的仅供大家学习和交流,请务用于商业用途! 3. 如果你也有好源码或者教程,可以到用户中心发布,分享有积分奖励和额外收入! 4. 本站提供的源码、模板、插件等等其他资源,都不包含技术服务请大家谅解! 5. 如有链接无法下载、失效或广告,请联系管理员处理! 6. 本站资源售价只是赞助,收取费用仅维持本站的日常运营所需! 7. 如遇到加密压缩包,默认解压密码为"gltf",如遇到无法解压的请联系管理员! 8. 因为资源和程序源码均为可复制品,所以不支持任何理由的退款兑现,请斟酌后支付下载 声明:如果标题没有注明"已测试"或者"测试可用"等字样的资源源码均未经过站长测试.特别注意没有标注的源码不保证任何可用性