很久不见。错过了很多的询问和回复,很抱歉。这段时间在整理我的第一篇paper,大部分涉及的东西不太具有可分享性(其实也是应该记录的)—————————————————————...
很久不见。错过了很多的询问和回复,很抱歉。这段时间在整理我的第一篇paper,大部分涉及的东西不太具有可分享性(其实也是应该记录的)—————————————————————...
最近在学习当下很火的环状DNA(eccDNA)的建库技术以及生信分析,特此记录。 No.1 去低质量&去接头 这一步需要注意的是低质量碱基的quality cutoff, 2...
记录今天学会的一个小知识:统计测序文件的reads长度分布 由于某些原因,需要统计reads长度分布对建库方法进行质控,对测序数据进行可用与否的评估。 有两个文库:质检差别比...
这一笔记的主要目的是记录学习和实践过程,所写皆尝试过为可用操作,看到的各位小伙伴有更好的操作欢迎随时交流~ 1. 创建 R 4.0以上的conda 环境[1]: 2. 安装...
发现年纪越大越容易不记得事情,以前记在本子上,收藏夹里的回忆都因为找不到本子、换了电脑而丢失,再想回去看看时,几番寻找都再得不到。凡此过往,皆为序章,铭记于心。脑子记不住,就...
今天要记录一个很实用的实验技巧:如题,怎么利用低浓度的质粒包好慢病毒吗? 慢病毒已经成为目前生命科学领域一个很常用的进行基因表达调控的实验方法了。目前也很多公司商业化包装能达...
课题需要某一篇已发表文章的基于10x Genomics平台的scATAC-seq数据,于是乎,开始了踩坑指南。 No.1 数据下载 众所周知,公共数据的下载一般主要有GEO和...
写得很好,且非常有帮助的文章!!!但是博主,按照这个流程,我发下一个问题: fasterq-dump转SRA生成的fastq文件可能reads的名字会被更改,对于PE测序不太友好。在fastq-dump里也是需要用到--defline-seq 这个参数来指定reads命名,这样才可行,找了一圈没在fasterq-dump中发现相应参数呢
生信 | scATAC-Seq数据获取(超快sra转fq)部分代码参考自10x genomics官网[what-is-cell-ranger-atac],但绝大多数代码融入了我自己的思考 1、上游分析(软件安装 + 数据获取 + c...