文章 微生物多样性分类柱状图按照分组进行排序脚本升级；tax_bar_plot.r

...一行：skip=1df<-read.table(opt$input,sep="\t",skip=1,header = TRUE,check.names=FALSE,comment.char="@")#处理第一列物种分类名字，没有物种名字的全部归为Otherspe.names<-str_split(df$`#OTU ID`,";",simplify=T)spe.name<-str_replace(spe.names[,ncol(spe.names)],"^\\w__","")...

文章 测序数据文件中MD5是个啥 ？不要无视我的存在！

...载该软件后，对下载后的文件用专门的软件（如Windows MD5 Check等）做一次MD5校验，以确保我们获得的文件与该站点提供的文件为同一文件，而不是被篡改植入病毒的。 讲到这里，大家应该对MD5有了一个大概的了解。我们回到文...

文章 蛋白质比对结果快速分类taxid

... int taxid; UT_hash_handle hh; }; struct thread_arg { char input_filename[256]; char output_filename[256]; }; void add_accession(struct my_struct **s, char *accession, int taxid) { struct my_struct *new_item; HASH_FIND_STR(*s, accession, new_item); // Search to see if accessi...

问题 老师，这个警告问题大不大

Wrote 92128018 bytes to primary GFM file: Amphitetranychus_viennensis.genome.1.ht2Wrote 35712720 bytes to secondary GFM file: Amphitetranychus_viennensis.genome.2.ht2Re-opening _in1 and _in2 as input streamsReturning from GFM constructorWarning: a local graph exploded (offset: 1239040, length: 57344...

文章 人类癌症数据库突变注释

...btype of cancer. For more updated cosmic, see instructions below on how to make them.icgc21: International Cancer Genome Consortium version 21 mutations.nci60: NCI-60 human tumor cell line panel exome sequencing allele frequency data

文章 BiocManager::install("seqinr")时报错

...ge ‘Rcpp’ ... ** package ‘Rcpp’ successfully unpacked and MD5 sums checked ** using staged installation ** libs g++ -std=gnu++14 -I"/share/work/biosoft/R/R-v4.2.2/lib64/R/include" -DNDEBUG -I../inst/include/ -I/share/work/biosoft/zlib/include -I/share/work/biosoft/bzip2/bzip2-v1.0.6/include...

文章 R 包安装报错：configure: error: gdal-config not found or not executable.

...ackage ‘sf’ ... ** package ‘sf’ successfully unpacked and MD5 sums checked ** using staged installation configure: CC: /share/work/biosoft/gcc/gcc-v6.4.0/bin/gcc configure: CXX: g++ -std=gnu++11 checking for gdal-config... no no configure: error: gdal-config not found or not executable. ERRO...

问题 KaKs_Calculator2 安装和使用求助, make 后产生错误，请问该怎么解决呢？

问题 基因家族分析时，提取cds与pep序列报错

...因组注释文件。 /share/work/biosoft/TransDecoder/latest/util/gff3_file_to_proteins.pl --gff3 $species.longest_isoform.gff3 --fasta $genome --seqType CDS >$species.cds.fa#提取pep序列/share/work/biosoft/TransDecoder/latest/util/gff3_file_to_proteins.pl --gff3 $species.longest_isofor...

话题 check.names

问题 老师您好，我在安装MCScanX时，按照手册修改好了三个文件，make之后出现了错误

文章 perl自动获取网页上的信息

...d store response at var $resp     my $resp = $ua->get($url);     # check response     if ($resp->is_success) {         # get response content(html source code)         my $content = $resp->decoded_content;         # use Regex get page title from $content         if ( ...

问题 老师您好，我在运行megahit对200多个环境微生物宏基因组数据进行组装时，megahit报错Exit code -9

...536 reads to search for mercy k-mers' 2025-05-26 08:27:54 - Error occurs, please refer to /data2/temp/megahit/log for detail 2025-05-26 08:27:54 - Command: /home/anaconda3/envs/megahit/bin/megahit_core seq2sdbg --host_mem 973841920819 --mem_flag 1 --output_prefix /data2/temp/megahit/tmp/k21/21 --n...

文章 合并不同sample的vcf文件，通过bcftools

...ge.vcf.gz -c, --stdout write on standard output, keep original files unchanged -f, --force overwrite files without asking -@, --threads INT number of compression threads to use [1]   2. 对生成的vcf.gz进行index： bcftools index [options] <in.bcf>|<i...

文章 blast报错

...会有如下错误提示： USAGE  makeblastdb [-h] [-help] [-in input_file] [-input_type type]    -dbtype molecule_type [-title database_title] [-parse_seqids]    [-hash_index] [-mask_data mask_data_files] [-mask_id mask_algo_ids]    [-mask_desc mask_algo_descriptions] [-gi_mask]    [-gi_...