笑死了,我真的笑死了分析的时候发现怎么宏基因组数据这么小,怎么双端300bp,怎么拼接后这么烂才发现用的16s数据,无语无语从新来正好检查一遍,不得不说AI真好用,写脚本很专...
笑死了,我真的笑死了分析的时候发现怎么宏基因组数据这么小,怎么双端300bp,怎么拼接后这么烂才发现用的16s数据,无语无语从新来正好检查一遍,不得不说AI真好用,写脚本很专...
困难重重的拼接,按照软件教程,一直报错 随后咨询了朋友,seqtk排除文件损坏,排除内存不足,最后发现是环境不匹配。用下面代码就成功了 随后使用脚本 因为上一步每一个样本在一...
1、软件配置KneadData 2、宿主基因组下载 3、脚本跑
step1:aspera下载 step2:质量控制FastQC + Trimmomatic:经典组合 拿一个数据为例 fastqc后会生成一个html网页文件,和一个fast...
以前写过的方法在;X9信息列因为有些基因没有基因名而错位更新
激活环境 conda activate cj 如何看有啥软件已经安装 conda list aspera下载数据 记得联网(悄悄地) 查找ascp钥匙openssh在哪 .....
@cf7de4c43d4d 结果不是输出为一个列表了吗
R获取KEGG/GO中某一通路所有基因1、通过KEGGREST包来探索KEGG数据库的12大代谢通路 通过观察KEGG官网 :https://www.genome.jp/kegg/pathway.html[htt...
1、通过KEGGREST包来探索KEGG数据库的12大代谢通路 通过观察KEGG官网 :https://www.genome.jp/kegg/pathway.html[htt...
!!!计算方法是以前初师兄根本文献扒的,后续把文章doi附上!!!KEGG基因list,利用R可以很快扒下来,我后面再更新,百度一下你就知道 准备好转录组表达矩阵,KEGG基...
以下列数据为例子,需要count矩阵并包含基因长度、(没有的话从gtf提取,前面我写过教程) count矩阵读入 count to TPM count to FPKM FPK...
以大设施anlei01环境为例 1、配置miniconda软件 2、设置环境 建立一个RNA-seq分析环境,这样是为了和其他环境区分开,以免本环境损坏造成其他环境也异常,也...
把sortbam文件全部上传到越哥服务器上给服务器上安装了bedgraphtobw,直接用就行了传文件使用的命令是 1、合并bam 2、bam转bw 3、callpeak
跟着官网[https://yezhengstat.github.io/CUTTag_tutorial/]学习CUT&Tag数据分析,参考:protocol[https://w...
在学校服务器上conda安装软件后 提示安装成功但是bowtie2报错缺失库文件 由于我是非root用户,权限受限在两位老师@exps和@wandering大佬的指导下解决了...
1、trim_galore去接头及低质量reads 使用nohup挂后台后可以用jobs -l查看 跑完以后在clean文件夹下,会有val_1和val_2生成,同时有一个r...
@大王_bf6a 大家的服务器不一样,缺什么东西的时候sudo没办法看,给我难受坏了
cellphoneDB呕心沥血下载安装教程我安装了一天之后快哭了,然后去100块买了个老师指导生信自学要么哭要么花钱,就是我这个垃圾1、安装https://github.com/ventolab/CellphoneD...
我的count的行名是转录本id就不行
CellPhoneDB的使用CellPhoneDB主要分成method、plot、query和database4个模块。 安装成功CellPhoneDB后,要先下载CellphoneDB databas...
从cellphonedb得到以下四个文件,用中间两个画图 画图