找到约 15 条结果

文章 kraken2生成的biom文件处理

... metadata = table.metadata(axis='observation') new_metadata = [] for otu_id, meta in zip(otu_ids, metadata): if meta is not None and 'taxonomy' in meta: # Filter out unnamed taxonomic ranks taxonomy = [level for level in meta['taxonomy'] if not level.endswit...

问题 老师您好,重新启动docker时daemon is not running该怎么解决?

文章 numpy从文件直接读入得到数组对象

... np.loadtxt("iris.data", dtype = "S", delimiter=',', usecols = [4])print(b)for x in range(len(b)):    print(a[x], "is ", b[0])          #指定数据类型,导入数据t = np.dtype([('sepal_length', 'f'),('sepal_width','f'),('petal_length','f'),('petal_width','f'),('class','S20')])a = np...

文章 Biopython之序列输入

...ord迭代器。这使您可以执行以下操作: from Bio import SeqIO for record in SeqIO.parse("example.fasta", "fasta"): print(record.id) or using a handle: from Bio import SeqIO with open("example.fasta", "r") as handle: for record in SeqIO.parse(handle, "fasta"): print(rec...

问题 circos 分析

...rcos -outputfile circos  CONFIGURATION FILE ERROR  Could not parse value for configuration parameter [max] in block [plot]. The  location in the configuration is shown above.  If you are having trouble debugging this error, first read the best practices  tutorial for helpful tips that address m...

文章 Linux:echo命令详解

...现更复杂的输出格式控制   1.显示普通字符串: echo "It is a test" 这里的双引号完全可以省略,以下命令与上面实例效果一致: echo It is a test   2.显示转义字符 echo "\"It is a test\"" 结果将是: "It is a test" 同样,双引号也可以省...

问题 GATK CALL 变异后导入GenomicsDBImport 出现问题

...报错如下: GenomicsDBImport cannot use multiple VCF reader threads for initialization when the number of intervals is greater than 1. Falling back to serial VCF reader initialization。尝试了很多次,都不行。

问题 老师,我提取cds ,pep用基因的ID,作为序列的ID命令。运行出来"No such file or directory; can't open file ./.longest_isoform.gff3"但是我打开输出文件看是有数据的。是不是不用管No such file or direc...提示

下面是运行图和输出文件

问题 win10企业版本,专业版本安装docker失败 WSL 2 installation is incomplete.

问题 http://genome-e.ibrc.or.jp/home/bioinformatics-team/mutmap

您好,可以发一份这个软件吗?

问题 http://genome-e.ibrc.or.jp/home/bioinformatics-team/mutmap

您好,可以发一份这个软件给我吗?

问题 gwas分析群体画进化树

...命令行  WARNING: ignoring environment value of R_HOME --> Q&A for bioinformatics, please visit the website: https://www.omicsclass.com/ --> R beginners ? I suggest your  learning  R language: https://study.omicsclass.com/index Registered S3 methods overwritten by 'treeio':   me...

问题 mysql镜像启动问题

...示: 2023-04-24 13:44:46+00:00 [Note] [Entrypoint]: Entrypoint script for MySQL Server 5.7.42-1.el7 started. 2023-04-24 13:44:46+00:00 [ERROR] [Entrypoint]: mysqld failed while attempting to check config command was: mysqld --verbose --help --log-bin-index=/tmp/tmp.TNVtIp74qF mysqld: Can't rea...

问题 da

...组创建annovar库 gtfToGenePred -genePredExt $gtf unknown_refGene.txt retrieve_seq_from_fasta.pl --format refGene --seqfile $fa --outfile unknown_refGeneMrna.fa unknown_refGene.txt table_annovar.pl $data/goat.recode.vcf ./ -buildver unknown -out snp -remove -protocol refGene -operation g -nastr...

文章 LTR_HARVEST的并行软件LTR_HARVEST_parallel

...整性的指标之一。目前经常使用的计算 LAI 的软件是 LTR_retriever,他并非是一个单一的工具,而是通过整合多个LTR预测软件的结果,过滤其中的假阳性LTR-RT,得到高质量的LTR-RT库,因此需要配合LTRharvest, LTR_FINDER, MGEScan 3.0.0等软件...