checkv的基本使用

最新推荐文章于 2024-05-17 16:41:20 发布

筱贺学生信

最新推荐文章于 2024-05-17 16:41:20 发布

阅读量2.3k

点赞数

分类专栏：生信文章标签： linux

本文链接：https://blog.csdn.net/m0_70166478/article/details/126723834

版权

生信专栏收录该内容

17 篇文章 0 订阅

订阅专栏

1.checkv运行的的步骤

A: Remove host contamination

B:Estime genome completeness

C:Predict closed genomes

D:Summarieze quality

2.下载与安装checkv

checkv基于conda的环境

conda install -c conda-forge -c bioconda checkv
#数据库下载（自动）
checkv download_database ./
#数据库下载（手动）
wget https://portal.nersc.gov/CheckV/checkv-db-v1.0.tar.gz
tar -zxvf checkv-db-v1.0.tar.gz
export CHECKVDB=/path/to/checkv-db-v1.0

3.checkv的基本使用

一步法

checkv end_to_end input_file.fna output_fiel -t -16

分步法

checkv contamination input_file.fna output_directory -t 16
checkv completeness input_file.fna output_directory -t 16
checkv complete_genomes input_file.fna output_directory
checkv quality_summary input_file.fna output_directory

4.checkv的运行的过程

CheckV v0.8.1: contamination
[1/8] Reading database info...
[2/8] Reading genome info...
[3/8] Calling genes with Prodigal...
[4/8] Reading gene info...
[5/8] Running hmmsearch...
[6/8] Annotating genes...
[7/8] Identifying host regions...
[8/8] Writing results...
Run time: 128.65 seconds
Peak mem: 0.12 GB

CheckV v0.8.1: completeness
[1/8] Skipping gene calling...
[2/8] Initializing queries and database...
[3/8] Running DIAMOND blastp search...
[4/8] Computing AAI...
[5/8] Running AAI based completeness estimation...
[6/8] Running HMM based completeness estimation...
[7/8] Determining genome copy number...
[8/8] Writing results...
Run time: 28.43 seconds
Peak mem: 0.25 GB

CheckV v0.8.1: complete_genomes
[1/7] Reading input sequences...
[2/7] Finding complete proviruses...
[3/7] Finding direct/inverted terminal repeats...
[4/7] Filtering terminal repeats...
[5/7] Checking genome for completeness...
[6/7] Checking genome for large duplications...
[7/7] Writing results...
Run time: 0.13 seconds
Peak mem: 0.25 GB

CheckV v0.8.1: quality_summary
[1/6] Reading input sequences...
[2/6] Reading results from contamination module...
[3/6] Reading results from completeness module...
[4/6] Reading results from complete genomes module...
[5/6] Classifying contigs into quality tiers...
[6/6] Writing results...
Run time: 0.04 seconds
Peak mem: 0.25 GB