04. 病毒基因组组装

您所在的位置:网站首页 基因组测序从头组装到尾组装 04. 病毒基因组组装

04. 病毒基因组组装

2024-07-12 16:00:22| 来源: 网络整理| 查看: 265

病毒基因组组装¶

Mark Renton 更新于: 2020-09-14

与细菌一般可以培养从而获得高质量的基因组DNA,然后进行基因组测序不同。有些病毒难以培养(如诺如病毒),有些病毒培养生物安全等级要求较高,且时间长。 病毒基因组测序有很多方法,这里只针对二代高通量测序数据进行讨论。用二代测病毒基因组有不同的方法,比如对临床样本中已知病毒进行特异性PCR扩增后测序;有些则采用随机引物扩增后再测序,这种方法可以针对未知病毒序列,但由于不同样本(比如临床样本)中含有的干扰核酸较多(如人源DNA等),往往需要实验中进行特殊处理或者生物信息学分析时筛选。

因此我们只讨论特异性PCR扩增后(类似一代测序)用二代测序分析已知病毒的方法。和细菌单克隆培养物不同,临床样本即使PCR扩增后,

由于病毒基因组位点变异较明显,但是大片段的异位重组比较少,因此可以采用基于 mapping 的方法获得基因组。

首先要选择参考序列,由于病毒编译大,我们一般尽可能选择近源的毒株。我么那可以采用 denovo 的方法先获得一个大片段高覆盖的非污染序列,然后 blast 找到近源毒株基因组序列,将其作为参考基因组再使用 mapping 的方法获得测序毒株的基因组。

采用 freebayes¶ $ spades.py $ blast $ bwa index ref.fasta $ bwa mem $ samtools index $ samtools $ freebayes 采用 GATK4¶ $ trimmomatic # ID: # PL: Platform 测序平台,这里数据来自 illumina 平台 # PU: # LB: # SM: Sample 样本名称 $ bwa mem -t 40 -M -Y -R "@RG\tID:S1\tPL:illumina\tPU:\tLB:\tSM:" ref.fasta \ > $ samtools sort -@ 4 -m 16G -O bam -o .sorted.bam $ gatk MarkDuplicates -I -M -O $ samtools index sorted.markdup.bam $ gatk BaseRecalibrator -R ref.fasta -I sorted.markdup.bam \ > -O sorted.markdup.recal.table $ gatk ApplyBQSR --bqsr-recal-file sorted.markdup.recal.table -R ref.fasta > -I sorted.markdup.bam -O sorted.markdup.BQSR.bam $ samtools index $ gatk HaplotypeCaller -R ref.fasta -I sorted.markdup.BQSR.bam \ > -O HC.vcf.gz # SNP call $ gatk VariantRecalibrator # Indel call 使用 snippy¶ 2020-09-14


【本文地址】

公司简介

联系我们

今日新闻


点击排行

实验室常用的仪器、试剂和
说到实验室常用到的东西,主要就分为仪器、试剂和耗
不用再找了,全球10大实验
01、赛默飞世尔科技(热电)Thermo Fisher Scientif
三代水柜的量产巅峰T-72坦
作者:寞寒最近,西边闹腾挺大,本来小寞以为忙完这
通风柜跟实验室通风系统有
说到通风柜跟实验室通风,不少人都纠结二者到底是不
集消毒杀菌、烘干收纳为一
厨房是家里细菌较多的地方,潮湿的环境、没有完全密
实验室设备之全钢实验台如
全钢实验台是实验室家具中较为重要的家具之一,很多

推荐新闻


图片新闻

实验室药品柜的特性有哪些
实验室药品柜是实验室家具的重要组成部分之一,主要
小学科学实验中有哪些教学
计算机 计算器 一般 打孔器 打气筒 仪器车 显微镜
实验室各种仪器原理动图讲
1.紫外分光光谱UV分析原理:吸收紫外光能量,引起分
高中化学常见仪器及实验装
1、可加热仪器:2、计量仪器:(1)仪器A的名称:量
微生物操作主要设备和器具
今天盘点一下微生物操作主要设备和器具,别嫌我啰嗦
浅谈通风柜使用基本常识
 众所周知,通风柜功能中最主要的就是排气功能。在

专题文章

    CopyRight 2018-2019 实验室设备网 版权所有 win10的实时保护怎么永久关闭