Svim 基于long reads鉴定SV

SVIM 可基于long reads(pacbio, ONT,HIFI)进行call SV,deletion, insertion, inversion, tandem duplications, interspersed duplication and translocations.

githup:https://github.com/eldariont/svim
文章:SVIM: structural variant identification using mapped long reads

1、安装

conda create -n svim_env --channel bioconda svim

2、简单操练

所需数据:

  • long reads FASTA/FASTQ (压缩或uncompressed均可)
  • or long reads比对得到的bam文件, 需要\color{red}{sort,index}

在进行long reads比对的时候,可以选择 NGMLR
或者minimap2均可以。

SVIM主要包括4个步骤:

  • collect: 基于Long reads 检测SV
  • cluster:相同的SV进行合并
  • combine:来自不同基因组区域的cluster进行合并
  • genotype:确定基因型
svim alignment my_sample my_alignments.bam my_genome.fa

输出文件:

  • The log file: SVIM_{date}_{time}.log
  • The SV calls in VCF format: variants.vcf
  • The SV calls in BED format: candidates/candidates_*.bed
  • Intermediate signature clusters in BED format: signatures/*.bed

当然也可以进行适当的过滤,比如

## score > 10
bcftools view -i 'QUAL >= 10' variants.vcf'.
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容

  • 本文学习费章军老师文章Genome of Solanum pimpinellifolium provides in...
    斩毛毛阅读 7,690评论 1 23
  • 组学工具合集 整理了一些目前知道的与组学相关的工具,可能不全面,欢迎补充 基因组学 一代组装 Celera Ass...
    Boer223阅读 13,163评论 3 68
  • 本文转载自知乎 nanopor技术专栏[https://zhuanlan.zhihu.com/p/10239286...
    wo_monic阅读 14,410评论 0 25
  • 渐变的面目拼图要我怎么拼? 我是疲乏了还是投降了? 不是不允许自己坠落, 我没有滴水不进的保护膜。 就是害怕变得面...
    闷热当乘凉阅读 9,777评论 0 13
  • 夜莺2517阅读 127,806评论 1 9