仅将此标记用于与生物信息学相关的编程相关问题。其他问题不属于此处,但可能是https://bioinformatics.stackexchange.com/上的主题。有关更多信息,请参阅标记维基。
我有一个包含三个字母氨基酸编码序列和RNA序列的字符串。我想按照它在字符串中出现的确切顺序提取氨基酸代码。 raw_seq ='''...
我需要一些关于使用heatmap.2生成热图的建议。我有一个基因表达值矩阵,它有15,616行和27列用于生成热图。问题是 ...
我想在我的基因组中绘制平均覆盖深度,染色体按顺序排列。我使用samtools计算了我的基因组每个位置的覆盖深度。我想要 ...
如何使用geom_label_repel避免绘图中标签的重叠?
我正试图用巨大的数据制作一个火山地块。在这里显示一些数据。 tab7
我正在尝试开发我的第一个小型shinyR应用程序,但它的UI不能与其功能正确集成。在运行app后浏览文件时,不执行该功能。 ...
我不擅长使用apply系列函数。我想使用apply系列函数而不是跟随类型的嵌套循环迭代9076 x 9076次以提高计算速度。最小的......
在biopython食谱中,我找不到如何实际运行clustalw。我已经完成了食谱上的内容,但它没有运行clustalw只是打印clustalw2 -infile = opuntia.fasta任何人都可以......
我正在尝试在openSUSE上安装samtools,我这样做了:cd htslib-1.2.1 ./configure make install工作正常。 bcftools-1.2 ./configure make install工作正常。对于samtools:cd samtools-1.2 ......
我有一个.fam,.bed和.bim文件,其中包含少数人的标记。我需要将其转换为VCF文件。有人可以帮助创建一个VCF文件。有没有可以做的开源工具......
Linkdatagen中的错误:在串联(。)或字符串中使用未启动的值$ chr
嗨,我试图使用linkdatagen,这是一个基于perl的工具。它需要一个vcf文件(使用SAMtools的mpileup)和hapmap注释文件(提供)。我按照指示但是......
在bash中评估表示单个管道命令的两个变量的最佳方法是什么?
我有一个函数产生,它确定文件是否存在,如果没有,它运行以下命令。当命令输出只是写入stdout时,这可以正常工作。但是在命令中......
处理12 Gb文件时遇到问题。谈到linux,我是新手。我希望有人能帮助我,任何建议都表示赞赏。我有一个名为phase_3.vcf的文件...
我正在尝试在shell / bash脚本中编写一个if-else条件,该脚本将用于许多不同的文件,因此它不适合某个结构。我有三个不同的文件,最多三个不同...
我正在使用bcftools共识从vcf文件中提取单倍型。给定输入文件:A.sorted.bam B.sorted.bam创建以下输出文件:A.hap1.fna A.hap2.fna B.hap1.fna B.hap2 ....
我写了一个准备一些文件的makefile。我创建ORIGINAL目录,然后我使用文件夹中的文件启动其他规则RDIR =。 RFILES:= $(通配符$(RDIR)/ * .vcf)OUTDIR = ...
下面是使用tabix从特定区域获取VCF文件的代码,然后使用vcftools中的“keep”选项对特定(欧洲)人口进行过滤。 ####如果[...选择特定人口]
这是确切的命令bcftools norm -f /path/hg19/ucsc.hg19.fasta -c s ./user1.vcf -o ../fixed/user2.vcf当我直接在shell中运行它时它可以正常工作。当我把它放入bash脚本时...
我有一个bash脚本接受参数,然后运行程序bcftools:$ infile = $ 1 $ outfile = $ 2 $ pop1 = $ 3 $ pop2 = $ 4我想使用pop1和pop2参数作为参数,所以在......