188 个回答
重测序额质控
需要提供CP162209.1.recode.vcf.gz截图看看是不是QD没有，正常GATK进行call SNP得vcf都有QD，vcftools分割的时候需要选择保留所有info才可以。可以在VariantFiltration之前用SplitIntervals进行染色体分割
回答于 2026-06-29 10:59
单倍型出现这种问题，怎么解决
hapSummary或 hapPrefix参数有问题，导致无法正确提取或匹配单倍型结果。我们脚本里会自动给hapPrefix赋值为“Hap”，所以应该是hapSummary出现了问题。hapSummary是对hapResult进行汇总整理得到的，通过你上面的这个log文件： Scanning file to determine attributes. File attributes: meta lines: 3374 header_lin...
回答于 2026-06-24 17:22
老师 我咱做比较基因组分析的时候，遇到下面的问题，我用其它的...
内存或系统资源被触发 OOM Killer 杀掉进程。如果是10.pep_aln.renamed这个结果有问题，那你可以ll aln_out/10.pep_aln.renamed看看这个比对结果的大小有没有问题，也可以grep ">" -c aln_out/10.pep_aln.renamed看看这个序列数目对不对。但是如果是一个比对结果有问题的话应该是报错，你显示的是killed更类似于资源不足...
回答于 2026-06-23 13:53
泛基因组 结构变异 注释
vcf基因组坐标区间采用了“倒位 断点上游 + 倒位区间 + 下游”的记录，ANNOVAR 使用 SV 的基因组坐标区间，你的 INV 被记录为覆盖从 GeneA 附近 到 GeneB 附近，注释按这个区间来。一个断点在 GeneA，另一个断点在 GeneB（或两者分别靠近两个基因），注释就会显示该 SV 影响 GeneA 和 GeneB。你可以检查一下断点是否在两个基...
回答于 2026-06-16 15:55
老师，我用perl $scriptdir/vcf2hmp.pl clean.vcf.gz hapmap.txt...
需要hapmap.txt.gz文件内容的截图，报错是格式问题
回答于 2026-05-25 09:59
运行 filterBAM_forHiC.pl hicxbly_aligned.REduced.paired_onl...
有报错信息吗？hicxbly_aligned.REduced.paired_only.bam没问题的话可能xblysample.clean.sam的头文件不完整，并且要求BAM 和 SAM 的染色体/contig 名称必须完全一致，可以从这两个角度检查一下
回答于 2026-05-15 17:59
老师，我做的是人的连续型变量，图片为什么会这样呢，可能是什么...
极端假阳性了。可能群体分层严重、存在亲缘关系没有较正、表型数据被当作协变量输入等原因都有可能。换一个模型也是这样的结果吗？
回答于 2026-05-15 15:45