找到约 15 条结果

问题 基因家族分析时,提取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...

文章 合并不同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...

问题 老师您好,我在安装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...

文章 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_...

问题 差异基因用p-vaule筛选可以吗?

...and true positives. Moreover, when the number of DEGs obtained is low, it makes sense to look at other methods of modeling gene count data or reducing the FDR cutoff, etc. RNA-Seq being an exploratory technique for hypothesis generation, must be accompanied by multiple testing correction to avoid ...

文章 染色体挂载:LACHESIS 安装

...sourceforge.net/projects/boost/files/boost/1.57.0/boost_1_57_0.tar.gz --no-check-certificate ./bootstrap.sh --prefix=/share/work/biosoft/boost/boost.1.57 --with-libraries=all ./b2 install ./bjam install

文章 biom格式的文件如何转化

...header-key taxonomys less feature_table_tax.tsv # Constructed from biom file #OTU ID FT-1 CKR-4 CK0-1 CKR-5 FTR-4 CK0-5 CKR-1 FTR-5 FT-2 FT-3 CKR-3 CK0-4 CK0-3 FTR-3 FTR-2 FTR-1 FT-5 CK0-2 FT-4 CKR-2 taxonomy 6f8a1344c343b000bbe6c0946ef7af3f 3.0...

问题 老师您好,我在运行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...

文章 shuf 随机抽取txt文件的行到新文件 linux命令

...ne  -n, --head-count=COUNT    output at most COUNT lines  -o, --output=FILE         write result to FILE instead of standard output      --random-source=FILE  get random bytes from FILE  -r, --repeat              output lines can be repeated  -z, --zero-terminated     line d...

文章 BUSCO评估基因组组装和完整性详解

...prok# or ignoring eukaryotes to save runtime, if compatible with your experimental goal.busco -m MODE -i INPUT -o OUTPUT --auto-lineage-euk# or ignoring non-eukaryotes to save runtime, if compatible with your experimental goal.  -i:基因组组装文件; -o:输...

文章 空间转录组数据读取

...置的函数直接读取了,代码如下 obj = Load10X_Spatial(data.dir = your_data_dir,filename = "filtered_feature_bc_matrix.h5",assay = assay) 这个datadir就是spaceranger跑出来的目录里有个子目录叫out,指定一下就行,后面的操作就和单细胞转录组差不多了,这个...

问题 我在做T2T,docker下载T2T镜像错误

... koq0g34v.mirror.aliyuncs.com: StatusCode: 403, This request is forbidden. Please proceed to https...

文章 bam文件转换fastq

...bamToFastq)Version: v2.25.0Summary: Convert BAM alignments to FASTQ files.Usage:   bamToFastq [OPTIONS] -i <BAM> -fq <FQ>Options:        -fq2    FASTQ for second end.  Used if BAM contains paired-end data.                BAM should b...