使用二代和三代WGS数据进行基因组组装

本教程中,我们的目标是使用第二代和第三代全基因组测序 (WGS) 数据组装细菌基因组。我们将以此为例来探讨 WGS 数据分析,并探讨测序技术之间的差异。

生物信息

使用NART通过读取分类进行长扩增子分析

NART设计用于基于图谱的纳米孔扩增子(实时)分析,例如 16S rRNA 基因。NART由NART(Nanopore Amplicon Real-Time entry)和 NAWF(Nanopore Amplicon snakemake WorkFlow entry)组成。通过基于映射的策略提供从基础调用读取到最终计数矩阵的(实时)端到端解决方案。

生物信息

使用LACA从长扩增子中从头挑选OTU

LACA是用于长扩增子一致性分析(例如 16S rRNA 基因扩增子分析)的可重复且可扩展的工作流程。它使用用snakemake管理工作流程以及conda来管理环境。

生物信息

16S rRNA基因扩增子分析

在此工作流程中,介绍了 Qiime2 和 R 中 16S rRNA 基因扩增子数据分析的主要步骤。本教程是为哥本哈根大学食品科学系的 MAC 2023 课程准备的。尽管这些步骤是为 Oxford Nanopore Tech (ONT) 测序设计的,但也在 Ilumina 短读长上进行了测试。

生物信息

用VueJS和Fastapi通过websocket实现进度条追踪

简介在我们的应用小程序中,我们是前后端分离的。前端页面只负责渲染,而后端需要处理数据。但是如果遇到数据量很大的情况下,我们处理起来就很缓慢,如果我们想通过 AJAX 的方法追踪后台数据变化的进度,需要用到轮询的方案,这个是非常消耗资源的。这里我们用 VueJS 和 Fastapi 的小例子演示前端传递数据,后台用 10 秒处理数据并实时反应进度给前台的实现。

网站搭建

生信小白如何用vuejs简单创建一个交互式网页

什么是 Vue.jsVueJS 是一个渐进式的前端框架,所谓渐进式的意思就是你可以用它快速完成原型创作,然后在此基础上逐步完善。他可以足够简单,也可以足够完善,那么对于新手小白来说,这简直就是福利!

网站搭建

Julia语言编写Needleman Wunsch全局比对算法

输入是两个字符串,输出是对齐后的两个字符串。

编程语言

Julia语言模仿BAM文件的pileup类似操作

简介使用过pysam和samtools的小伙伴肯定了解 pileup的操作,如果把 BAM 文件看作表格的话,那么通常我们是按行去解析它的 record,进而获得一些信息,例如比对到哪条染色体,比对的开始位置和结束位置等. 另一种情况下,我们想要按照列去循环解析,得到这个列上的具体信息,典型的就是这个列上比对序列的碱基是什么?比对序列的位置是什么?以及是 Match or Mismatch or indel 等。那么,该操作就需要引入pileup操作了。

编程语言

Python编写拆分Barcode的脚本,并用Codon编译为Native code

摘要成对的 reads 中,read_2 的开头包含两份 barcode 序列,分别长 10bp,中间有一段固定长度为 15bp 的序列分割,例如 1ATCTATGACATGTTACGTTAACTCCNATCTATCACTTAGCGCTGNCCCTGTCCTCTACACTCCACCCCCTCCCCACCAGACTAAACAACGCCCTTTCCCC 该序列中ATTTATGACA及AATCTATCAA为 barcode 序列。要注意,barcode 因为测序的原因存在一定的错配,需要对其有一定的容纳。

编程语言

单细胞数据如何绘制stacked violin?

Python 的Scanpy包和Seurat包一样,是单细胞数据处理的利器,其中,Scanpy中有一种堆积的小提琴图,可以很好的展示 marker 的表达情况,但是在Seurat中并没有内置命令。因此,我自己尝试提取数据并用ggplot2包来画该图。 首先来展示以下画图的成果,如图

生物信息

站内搜索

没有找到内容!