找到约 15 条结果

问题 您好,我扩增子分类注释报错了,麻烦您帮我看看

sh: line 1:    78 Segmentation fault      uclust --input "otus.fa" --id 0.8 --rev --maxaccepts 3 --allhits --libonly --lib "/work/database/SILVA_132_QIIME_release/rep_set/rep_set_16S_only/97/silva_132_97_16S.fna" --uc "/work/tmp/UclustConsensusTaxonAssigner_al463L.uc" > "/work/tmp/tmpFpr35JV...

问题 MEME就运行脚本的时候,其中的maxsize 指最大输入文件的字符数,maximum dataset size in characters,这个是指多少b吗,6000000是指文件最大6M,单位是什么?时间的最大值单位是秒吧?,motif的个数最小长度一般是多少?有什么参考的吗,最大长度是多少,这个长度的单位是什么,多少b,还是多少kb?

/biosoft/meme/meme-v4.12.0/bin/meme ../WRKY_pep_confirmed.fa -protein -oc ./ -nostatus -time 18000 -maxsize 6000000 -mod anr -nmotifs 10 -minw 6 -maxw 100

文章 trimAL | 多序列比对中剪切掉不保守部分的参数

... |剪切结果序列长度最低不能低于原序列长度的80% trimal -in aligned.fa -out aligned_trimed.fa -gt 0.8 -st 0.001 -cons 80 有了这三个参数,就可以对比对结果中一些不保守的列进行剪切/删除。只保留保守的序列部分。 除此之外还有一些自动...

问题 基因家族分析时获取基因与mRNA的对应关系

错误提示:Use of uninitialized value $pid in exists at script/mRNAid_to_geneid.pl line 40, <IN> line 188156 gff文件内容有什么问题

文章 IF=4.2 | 不同西瓜品种果实品质的转录组代谢组联合分析

...Scientia Horticulturae杂志发表了一篇题为“Combined transcriptome and metabolome analysis reveals the mechanism of fruit quality formation in different watermelon (Citrullus lanatus) cultivars”的转录组代谢组联合分析文章。” 01 — 研究背景 西瓜作为夏季最受...

问题 bioconductor常用包 安装

...quired to build R packages but is not currently installed. Please download and install the appropriate version of Rtools before proceeding:https://cran.rstudio.com/bin/windows/Rtools/Installing package into ‘C:/Users/Lenovo/Documents/R/win-library/3.6’(as ‘lib’ is unspecified)试开URL’htt...

文章 快速了解转录调控家族的新贵-circRNAs

...用这个好用的软件CIRI。 5 参考文献 Yuan Gao, Jinfeng Wang, and Fangqing Zhao.CIRI: an efficient and unbiased algorithm for de novo circular RNA identification.Genome Biol. 2015; 16(1): 4. Chen G, Cui J, Jin B al.Genome-Wide Identification of Circular RNAs in Arabidopsis thaliana.Front...

问题 误: 无法分配大小为7.4 Mb的矢量 Killed after 9 m

...’ The following object is masked from ‘package:base’:     expand.grid 载入需要的程辑包:IRanges 载入程辑包:‘IRanges’ The following object is masked from ‘package:grDevices’:     windows 载入需要的程辑包:GenomeInfoDb 载入需要的程...

文章 echo输出不换行

...可以用echo命令的-n选项,作用是输出是不追加换行。 #!/bin/bash/for i in `cat sample.txt`do        echo -n "cat ">> run.sh        for j in `cat primerid.txt`        do                echo -n " $i/${i}_$j.fa">> run.sh        done        echo...

问题 顺势作用原件功能评估

...www.omicsclass.com/article/1154  Fig4中描述的比如顺势原件 MYB binding site involved in flavonoid biosynthetic genes regulation; LTR, cis-acting element involved in low-temperature responsiveness。请问老师如何评估这些元件的功能,文中也没有引用参考。我觉得评...

问题 老师您好!

...  directory = DataDirectory,files.per.chunk=6, method='client') Downloading data for project TCGA-CHOL 试开URL’https://gdc.cancer.gov/files/public/file/gdc-client_v1.6.0_Windows_x64-py3.7_0.zip' Content type 'application/zip' length 16207515 bytes (15.5 MB) downloaded 440 KB Error in unz...

文章 Scanpy数据结构:AnnData

... AnnData object with n_obs × n_vars = 6050 × 32285 obs: 'in_tissue', 'array_row', 'array_col', 'gender', 'age', 'tissue', 'Strain', 'n_genes_by_counts', 'log1p_n_genes_by_counts', 'total_counts', 'log1p_total_counts', 'pct_counts_in_top_50_genes', 'pct_counts_in_top_100_genes', 'pct...

文章 扩增子序列直接物种注释代码

...代码如下 import re with open("results_clean_batch1.txt", "r", encoding="utf-8") as file: text = file.read() matches = re.findall(r'(Query=.*?)(?=Query=|$)', text, re.DOTALL) results = [] for match in matches: query_id = re.search(r'Query= (\S+)', match).group(1) species_matches =...

文章 Uniprot蛋白数据库介绍及使用详解!

...titute)、SIB(the Swiss Institute of Bioinformatics)、PIR(Protein Information Resource)三大数据库的资源。 EBI( European Bioinformatics Institute): 欧洲生物信息学研究所(EMBL-EBI)是欧洲生命科学旗舰实验室EMBL的一部分。位于英国剑桥欣克...

文章 你了解SAM和BAM文件吗

...s2的reads数量properly paired:正确配对的reads数量with itself and mate mapped:一对reads均比对上的reads数singletons:只有单条reads比对上的reads数以上计数均以reads条数计,一对reads计为两条。 还可以通过以下命令快速查看flag值所对应的...