OthroFinder-物种同源基因分析工作原理

Orthofinder是做直系同源基因分析时常用软件,除此之外,还可以使用OrthoMCL,今天主要介绍Orthofinder软件的工作流程。 1. Orthofinder发表的工作流程如上图:     1.1 (a)将所有物种的...

  • 0
  • 0
  • rzx
  • 发布于 2023-06-15 16:15
  • 阅读 ( 2400 )

提取 VCF 文件中的基因型信息

VCF 是一种常用的生物信息学文件格式,用于存储基因组中的遗传变异数据,特别是单核苷酸多态性(SNP)和小插入/缺失(Indel)等变异。 VCF 文件通常由文本格式组成,可以用文本编辑器进行查看...

  • 0
  • 0
  • 星莓
  • 发布于 2023-06-15 15:59
  • 阅读 ( 2821 )

csvtk,A powerful software for handling CSV files

CSVTK,即CSV工具包,是一个用于处理逗号分隔值(CSV)文件的实用命令行程序。由于其简单性、灵活性和效率,对于经常处理CSV文件的人来说,特别是在数据科学、生物信息学以及任何需要数据分析的...

  • 0
  • 0
  • xun
  • 发布于 2023-06-08 09:54
  • 阅读 ( 534 )

运行命令makeblastdb建库时报错,生成pdb-lock 空文件。

运行命令makeblastdb -in all.pep.fa -dbtype prot -title all.pep.fa 时卡住报错,生成all.pep.fa.pdb-lock 空文件。 解决:makeblastdb 工具版本问题,从2.13.0版本改成2.6.0版本,就可以...

  • 0
  • 0
  • rzx
  • 发布于 2023-06-05 15:20
  • 阅读 ( 675 )

bedtools 对bed 文件进行排序

bedtools 对bed 文件进行排序

  • 0
  • 0
  • 安生水
  • 发布于 2023-06-05 11:35
  • 阅读 ( 1233 )

blastn使用教程及结果详解

BLAST (Basic Local Alignment Search Tool) 是我们常用的短序列比对工具,直接输入fasta格式的序列文件就可进行比对

  • 0
  • 0
  • 星莓
  • 发布于 2023-05-29 14:49
  • 阅读 ( 2814 )

计算RNASeq 的Power值

计算RNASeq 的Power值

  • 0
  • 0
  • 安生水
  • 发布于 2023-05-29 13:25
  • 阅读 ( 851 )

构建一致性序列

构建一致性序列

  • 0
  • 0
  • 安生水
  • 发布于 2023-05-19 11:33
  • 阅读 ( 3093 )

python2运行lefse软件的时候报错

使用python2运行lefse软件的时候报错这是因为当前LEfSe依赖的R包rpy2的版本与你的R版本不兼容,你可以尝试更新一个rpy2的小版本,比如我更新了 pip2 install rpy2==2.8.6 随后就可以直接使用...

  • 0
  • 0
  • xun
  • 发布于 2023-05-12 18:06
  • 阅读 ( 390 )

Minimap2安装与使用

minimap2是生信大牛Heng Li在2018年发表的三代序列比对工具。与传统的bwa相比,Minimap2可以用于三代测序, 也支持 splicing awared 比对;与一些传统三代比对工具相比,Minimap2 速度非常快,...

  • 0
  • 1
  • 星莓
  • 发布于 2023-05-12 11:07
  • 阅读 ( 4739 )

MNP标记设计中DP的计算方法

MNP标记设计中DP的算法

  • 0
  • 0
  • 安生水
  • 发布于 2023-05-08 16:34
  • 阅读 ( 587 )

MNP标记设计

MNP标记设计

  • 0
  • 0
  • 安生水
  • 发布于 2023-05-08 16:22
  • 阅读 ( 3827 )

很大的压缩文件计算行数

这个其实很多办法都可以做到,基本就是靠zcat后接awk,sed,或者wc命令,不过速度差强人意,比如这个 awk应该是这几个里最快的,但一个3个G左右的文件仍然运行了接近2分钟,考虑到还有几个30个...

  • 0
  • 0
  • xun
  • 发布于 2023-05-05 15:22
  • 阅读 ( 423 )

根据fastq ID列表,输出指定ID对应的双端或单端fastq文件

如果某fq_clean文件的其中一端出现了错误,我们手里还持有他的原始数据,那我们就可以用以下方法处理1,首先提取clean文件另一端的id,我用了python脚本 import gzipimport argparsedef extrac...

  • 0
  • 0
  • xun
  • 发布于 2023-05-05 13:31
  • 阅读 ( 430 )

三代数据SV检测软件安装与使用——cuteSV

三代测序在检测基因组结构变异方面有着很大的优势,但是由于数据分析算法、软件还处于不断开发更新中,当前还没有公认推荐的软件。这里为大家推荐一款国人开发的三代数据SV检测软件—cuteSV。该...

  • 0
  • 0
  • 星莓
  • 发布于 2023-04-28 14:14
  • 阅读 ( 965 )

使用vcftools对vcf文件做滑窗处理

使用vcftools对vcf文件做滑窗处理

  • 0
  • 0
  • 安生水
  • 发布于 2023-04-26 10:13
  • 阅读 ( 895 )

统计每个窗口内SNP变异位点数量

统计每个窗口内SNP变异位点数量

  • 0
  • 0
  • 安生水
  • 发布于 2023-04-07 16:20
  • 阅读 ( 1197 )

linux使用ossutil下载阿里云数据

linux使用ossutil下载阿里云数据

  • 0
  • 2
  • 星莓
  • 发布于 2023-03-31 16:30
  • 阅读 ( 1746 )

muscle新版,旧版使用的一些差别

muscle新版,旧版使用的一些差别

从fasta基因组中提取反向互补序列

从fasta基因组中提取反向互补序列,借助bedtools工具。 bedtools getfasta -fi Dlong_asm_chr.fasta -bed DlNIP.bed -s -fo DlNIP.bed.fa -fi 基因组文件     -bed 基因位置 共6列:【染色...

  • 0
  • 0
  • rzx
  • 发布于 2023-02-27 14:22
  • 阅读 ( 1018 )