bioinformatics 相关问题

仅将此标记用于与生物信息学相关的编程相关问题。其他问题不属于此处,但可能是https://bioinformatics.stackexchange.com/上的主题。有关更多信息,请参阅标记维基。

无法访问jar文件abra2.jar

$ java -Xmx8G -jar abra2.jar --in Sorted_Output.bam --out Sorted_Output.abra.bam --ref GCA_genomic.fna --threads 8 --tmpdir tmpdir > abra.log 错误:无法访问 jarfile abra2.jar 我不...

回答 1 投票 0

Julia 中用于解决 Rosalind 问题“开放阅读框架”的对象类型操作问题

代码上下文:我正在尝试解决 Rosalind 问题“开放阅读框架”(https://rosalind.info/problems/orf/)。我使用的方法是将每个 ORF 结果存储在一个变量中

回答 1 投票 0

从一个 fastq 文件中获取在另一个 fastq 文件中找不到的读取内容

我想要一个由 F1.fastq 中但不在 F2.fastq 中的 fastq 读取组成的 fastq 文件。 此处提供以下代码作为解决方案。 所需读取数 = [] 读取数组 = []

回答 2 投票 0

获取F1.fastq中不在F2.fastq中的fastq读取

我想要一个由 F1.fastq 中但不在 F2.fastq 中的 fastq 读取组成的 fastq 文件。 此处提供以下代码作为解决方案。 所需读取数 = [] 读取数组 = []

回答 1 投票 0

错误:!未能收集惰性表。由 `db_collect()` 中的错误引起 - 在 R 中使用 biomaRt 包

我目前正在使用 R 开发一个生物信息学项目,在尝试使用 biomaRt 包时遇到错误。安装包并将其加载到 R 中后,我尝试选择...

回答 1 投票 0

将字符串和数据框转换为另一个数据框

我有一个单行数据框,如下所示: 捐赠者治疗时间点 MK434_016 WT5 ST002_50uM 6小时 和一个看起来像这样的字符串: [1]“

回答 1 投票 0

如何使用 Perl 搜索字符串中生物信息学主题的不同变体?

我有一个程序输出,其中有不同变体的一个串联重复。是否可以(在字符串中)搜索主题并告诉程序查找最大“3”不匹配的所有变体/

回答 3 投票 0

尝试安装 BioGeoBears 包依赖项 rexpokit 时出错,不再在 CRAN 上且无法安装

我正在尝试安装 R 包 rexpokit,它是 BioGeoBears 的依赖项。该软件包的 CRAN 版本刚刚于 2023 年 10 月 10 日被删除。 我尝试从ar下载并安装...

回答 1 投票 0

如何使用Bio.PDB.Selection仅展开蛋白质原子?

从 Bio.PDB 导入 PDBParser 从 Bio.PDB 导入选择 Structure = PDBParser().get_struct('4GBX', '4GBX.pdb') # 加载你的分子 atom_list = Selection.unfold_entities(结构[0]['E'], '...

回答 1 投票 0

Biopython 代码给出断言错误

biophoton 中的以下代码给了我一个断言错误:- 从 Bio 导入 SeqFeature 开始 = SeqFeature.ExactPosition(5) 结束 = SeqFeature.BetweenPosition(25, 左=8, 右=28).to_location()

回答 1 投票 0

BioPython:KEGG REST 不断报告 HTTP 错误 403:禁止

我正在尝试使用Bio.KEGG中的BioPython的REST模块来查询KEGG数据库以获取某些化合物的名称和分子式,使用化合物化学识别号(CID),例如

回答 2 投票 0

通过匹配其他文件中的顺序来打印行

文件1.txt SRR1071717 SRR1079830 SRR1081765 SRR1085402 SRR1086236 SRR1092208 SRR1093930 SRR1097296 SRR1099957 SRR1120296 文件2.txt SRR1120296,女 SRR1099957,男 SRR1097296,女 SRR1093930,女性...

回答 1 投票 0

如何向第二列作为目录名称的列添加标题

我有多个目录中包含 ID 和计数的计数文件(对于每个加入 SRRXXXXXX)。我想使用 bash 循环在每个文件中添加标头“gene_id”和 SRRabcdXXX。 目录

回答 1 投票 0

如何安装 Bioconda 软件包

这是一个常见问题的通用表述。我正在尝试从 Bioconda 安装软件包 foo。我尝试: conda 安装-c bioconda foo 然而,这会导致无法满足的依赖...

回答 1 投票 0

无法在 Mac 上安装 Bioconda 软件包

第一次使用命令行,我尝试在带有 M1 芯片的 Mac 环境中安装 Bioconda 软件包(fastp 和 Bowtie1),但我不断收到相同的错误: 收集包元数据 (

回答 2 投票 0

查找具有非零元素的最大尺寸子矩阵

我有一个矩阵,其中行是基因名称,列是样本名称。该矩阵的元素为0或1(1表示该样本中表达的基因,0表示不表达)。我想找到最大子集...

回答 1 投票 0

使用 AWK 将两个不同的 gsub 结果输出到两列

我正在尝试优化 awk 代码以在我的 VCF 文件中生成新列。 这是输入(最小工作示例): #CHROM POS 基因参考替代定性过滤器信息格式 1 2 3 ...

回答 1 投票 0

SAM 对齐:提取查询序列中的特定区域及其 CIGAR 字符串中的封闭部分

我需要对已进行全局比对的DNA序列的给定区域执行局部比对,并更新全局CIGAR字符串的相应部分。 步骤如下...

回答 1 投票 0

Google Cloud 平台 - vertex AI 无法获取 a100 GPU 实例

所以我一直收到这个错误。 区域“projects/xxxxx/zones/Asia-northeast3-a”没有足够的资源来满足该请求。尝试不同的区域,或者稍后再试...

回答 1 投票 0

snakemake 无法在参数中获取 lambda 返回

我有一个snakemake参数,例如: def get_papreDE(config, super_group): gtf_config = [] 对于sample.keys()中的sample_name: group_data=样本[样本名称][0] 超级组数据=...

回答 1 投票 0

© www.soinside.com 2019 - 2024. All rights reserved.