当 Salmon 构建 index 秒失败时应该反思哪里

操作环境:Win10 Ubuntu子系统

1. 下载参考转录组

下载传送门:http://ftp.ensembl.org/pub/release-76/fasta/homo_sapiens/cdna/Homo_sapiens.GRCh38.cdna.all.fa.gz

由于 wget 下载速度太慢,用了神奇的Motrix软件分分钟成功下载到本地。

就是来自技能树推文的安利之一,推文传送门:https://mp.weixin.qq.com/s/whx-n1ktT58WHhWAUa06dw


就是它👇

然后直接将 Homo_sapiens.GRCh38.cdna.all.fa.gz 拖到子系统的工作路径下。

2. Salmon 构建 index

假设已经装好了 Miniconda.

2.1 在小环境中安装 Salmon

conda create --name rnaseq
conda activate rnaseq
conda install -y salmon

2.2 构建 index

失败第一回

salmon index -t Homo_sapiens.GRCh38.cdna.all.fa.gz -i homo38_index

[Error reading from the FASTA/Q stream. Make sure the file is valid.]

难道是文件不完整吗,抱着这种想法甚至重新下载了1n遍。

失败第二回

看了别人构建好的 index, 里面有 .json 文件,难道是和网速有关?

开着战五渣的手机热点,依然 [Error reading from the FASTA/Q stream. Make sure the file is valid.] .

失败最终回/成功第一回

查看文件详细信息:

ls -lh

没有权限,毫无权限。

那么就来修改权限🤦‍♂️

chmod 777 Homo_sapiens.GRCh38.cdna.all.fa.gz

这是坠吼的

salmon index -t Homo_sapiens.GRCh38.cdna.all.fa.gz -i homo38_index

所以构建 index 秒失败应该反思哪里呢

  • 文件完整性,有必要
  • 网速,阔以,么得必要
  • 文件权限,非常必要,尤其用了除 wget 以外的操作下载

最后,向大家隆重推荐生信技能树的一系列干货!

  1. 生信技能树全球公益巡讲:https://mp.weixin.qq.com/s/E9ykuIbc-2Ja9HOY0bn_6g
  2. B站公益74小时生信工程师教学视频合辑:https://mp.weixin.qq.com/s/IyFK7l_WBAiUgqQi8O7Hxw
  3. 招学徒:https://mp.weixin.qq.com/s/KgbilzXnFjbKKunuw7NVfw
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

友情链接更多精彩内容