BLINK(C版本) 进行 GWAS分析的教程

  • 0
  • 0
  • 星莓
  • 发布于 2023-12-27 11:31
  • 阅读 ( 6770 )

bam文件排序

通常进行SNP calling的之前我们需要对比对的bam文件进行排序、去重等工作,在这里介绍几种sort bam的方法。 1. 按照比对位置进行排序 sambamba比samtools快很多,samtools和sambamba默认按...

  • 0
  • 0
  • Ti Amo
  • 发布于 2023-12-27 11:05
  • 阅读 ( 3171 )

SURVIVOR工具进行结构变异结果文件的格式转换(bedtovcf)

结构变异结果文件格式转换工具——SURVIVOR

R语言中转换成dummy dummy-variables 虚拟变量

转换成dummy dummy-variables

  • 0
  • 0
  • omicsgene
  • 发布于 2023-12-21 15:11
  • 阅读 ( 2506 )

R中的错误处理函数

很多时候简单的try函数就能够增加脚本的鲁棒性,随着脚本逐渐复杂,一些较高级的错误处理函数就必不可少了,我们今天来介绍几个R中的 1. tryCatch() 函数 功能: tryCatch() 函数是一种更复杂...

  • 0
  • 0
  • xun
  • 发布于 2023-12-19 14:39
  • 阅读 ( 1580 )

awk同时读取两个文件

awk使用NR==FNR同时读取两个文件

ssh登陆报错

ssh登陆报错

  • 0
  • 0
  • 安生水
  • 发布于 2023-12-18 16:19
  • 阅读 ( 1858 )

下载AI字体

最近在使用AI软件对iTol生成的pdf文件进行修改。当AI导入文件后弹出提示窗口如下: “缺少字体。将使用替换字体显示受影响的文本。” 一般情况下可以忽略提示。 但这次的文件在忽略提示直接导...

  • 0
  • 0
  • rzx
  • 发布于 2023-12-18 15:40
  • 阅读 ( 1725 )

NCBI 数据下载

介绍如何在大型生物数据库NCBI上下载数据

  • 1
  • 0
  • Ti Amo
  • 发布于 2023-12-18 11:25
  • 阅读 ( 3321 )

GATK|Timeout of 30000 milliseconds was reached while trying to acquire a lock on file raw.vcf.tmp.idx.

解决办法:https://www.omicsclass.com/article/2194 这里还有一个需要做的步骤,就是先将raw.vcf.tmp.idx文件找到并删除,一般情况下这个文件是空的。 然后再在命令行中添加-disable_auto_in...

  • 0
  • 0
  • rzx
  • 发布于 2023-12-15 10:23
  • 阅读 ( 1686 )

37门付费精品课程免费送!(此链接微信中打开体验更好)

在这个科研内卷的年代,你没有好设备,没有好材料,没有充足经费,怎么立足?毕业、就业、发文、评职称,靠什么来实现?留给你破局的选择不多!!为此,我们组学大讲堂决定干一票大的!我们将永...

SURVIVOR合并单样品SV出错

在做SV、CNV变异检测时,会使用SURVIVOR将不用样品的vcf文件合并到一起,命令如下: ls *.max.vcf >sample_files && /share/work/biosoft/SURVIVOR/SURVIVOR/Debug/SURVIVOR merge...

  • 0
  • 0
  • 安生水
  • 发布于 2023-12-13 13:51
  • 阅读 ( 2798 )

InterProScan-蛋白质预测

InterProScan是InterPro数据库提供给用户的一个软件包,通过集成数据库鉴定蛋白质功能、所属的家族及其包含的结构域和位点

MarkDuplicates 过程中报错:Value was put into PairInfoMap more than once

RNA-seq分析进行SNP calling之前会对bam文件进行一系列处理,包括排序、去重、将落在外显子上的reads分离出来,取出N错误碱基,在去重(MarkDup)过程中出现报错: 是由于fq文件中出现了两次...

  • 0
  • 0
  • Ti Amo
  • 发布于 2023-12-08 14:33
  • 阅读 ( 2511 )

好用的perl模块推荐

Path::Tiny 是一个专注于简洁性和便利性的 Perl 模块,用于处理文件和目录操作。以下是 Path::Tiny 的一些主要功能 1,目录相关 在perl 里很多人会使用 mkdir path unless(-d path); 来创建...

  • 0
  • 0
  • xun
  • 发布于 2023-12-07 16:16
  • 阅读 ( 1393 )

docker配置国内源加速下载,linux版本

docker配置国内源加速下载,linux版本

  • 1
  • 0
  • omicsgene
  • 发布于 2023-12-07 11:14
  • 阅读 ( 2368 )

break报错

break报错

  • 0
  • 0
  • 安生水
  • 发布于 2023-12-06 14:47
  • 阅读 ( 1378 )

微生物多样性扩增子(16S/18S/ITS)课程更新v2.0

微生物多样性扩增子(16S/18S/ITS)课程更新v2.0

  • 0
  • 0
  • omicsgene
  • 发布于 2023-12-05 13:25
  • 阅读 ( 2496 )

Python模块导入的 路径索引

在使用python的时候经常需要调用库及module,有的时候分明已经import了,在运行的时候却收到 ModuleNotFoundError的报错。那么Python到底是从哪里、以什么样的顺序调用的呢? 在Python中有内建...

  • 0
  • 0
  • Ti Amo
  • 发布于 2023-12-05 10:17
  • 阅读 ( 1346 )

vcf过滤时GATK和bcftools的index报错及解决方式

gatk报错 An index is required but was not found for file drivingVariantFile。 bcftools报错 index: the file is not BGZF compressed, cannot index

  • 1
  • 0
  • 星莓
  • 发布于 2023-12-05 10:01
  • 阅读 ( 5403 )