最近工作很忙，忙的喘不上气的感觉，对人生又多了很多感慨，不管未来怎样，还是要坚定的走下去。 从一个做DM的朋友哪里了解到sas程序员这个工作，开始也没有太在意,直到两年后，和...
发简信
IP属地：安徽
最近工作很忙，忙的喘不上气的感觉，对人生又多了很多感慨，不管未来怎样，还是要坚定的走下去。 从一个做DM的朋友哪里了解到sas程序员这个工作，开始也没有太在意,直到两年后，和...
@LRivers 用法是 python get_ASprofile_ref_hdrs.py path/species.genome.fa species，如果是人类可以输入 python get_ASprofile_ref_hdrs.py Homo_sapiens.GRCh38.99.dna.primary_assembly.fa Homo_Sapiens
哪里可以获得.hdrs文件？extract-as 简介 该脚本需要两个参数： 转录本对应的 gtf 文件 基因组长度统计文件，后缀为 .hdrs，内容如下：>chr1 /len=249250621 /n...
报错 解决办法 如果需要导入的文件较少，可以手动转换成 .csv 或 .txt 格式，再用相应的函数导入； 可以试试 read. table() 后记 虽然问题解决了，但是始...
gffcompare 这个报错真让我 yue 了！ 报错 不对，怎么报错了？我前面也用了这个命令，能成功运行啊？ 经过 那有错也还得解决啊，还能咋办？ 把“老三套”独立解决报...
已有的工具如circos, 只能绘制弯的染色体，又如R包chromoMap，IdeoViz，karyploteR，ggbio和在线工具Idiographica，虽然能绘制直的...
对于strand-specific的RNA-seq而言，我们必须得知道它是哪一种建库方式，才能进行后续的定量分析。 stringtie: kallisto: 现在比较常用的方...
前言 这次给大家带来的是16年发表在NATURE PROTOCOLS上面的一篇处理RNA-seq数据的文章：Transcript-level expression analy...
报错内容 报错原因 内存不够，清点内存就行了。 亲测成功！ 写在后面 这个报错信息显示了四句话，检索一句两句，会检索出很多内容，有以下几个方面： 单端数据用成了双端的代码，或...
Learn Git Branching 总结 链接：https://learngitbranching.js.org 主要 基础篇 1.Git Commit 2.Git Br...
注意 python 中，一个需要特别注意的点。 gzip.open() 默认的 mode 参数是 rb，即使设置 r，无论是否有额外的 b ，都会以 rb 的模式读入。如果需...
报错 解决办法 对于内存问题（如果你有一些很长的 reads，nanopore 肯定会触发这些问题），那么你可以尝试添加 -t 2 （甚至大于 2）来强制 fastqc 拥有...
extract-as 简介 该脚本需要两个参数： 转录本对应的 gtf 文件 基因组长度统计文件，后缀为 .hdrs，内容如下：>chr1 /len=249250621 /n...
今天在使用 conda activate 环境的时候，出现了如下报错： 在 github 上有个同样的问题，链接为：https://github.com/conda/cond...
今天在运行 bwa mem 命令时，出现如下错误提示： 报错是因为 ref.fa 文件里存在着空格、空行，一般是出现在> 行。 这其实是一个很小的问题，但是很容易忽略，特别是...
