Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

奇怪的BUG:报错信息solving row 120798: 'allow.nonnarrowing' is FALSE and the supplied end (16581) is > refwidth #32

Closed
b-niu opened this issue Feb 12, 2022 · 2 comments

Comments

@b-niu
Copy link

b-niu commented Feb 12, 2022

Dear @ShixiangWang ,

在运行sigflow时,遇到了一个报错信息:
image

查询得知,一般该报错来自于使用了错误的参考基因组,导致坐标出界。例如:
PoisonAlien/maftools#85

但是我仔细检查了使用的基因组,确认bwa比对和sigflow使用的都是hg19,且报错的这一行row 120798的坐标信息确实位于基因组内。

根据 https://www.genome.ucsc.edu/goldenpath/help/hg19.chrom.sizes:
chr15 size = 102531392
row 120798是
image

image

请问您有什么线索可以让我找出错误所在吗?
谢谢

@b-niu
Copy link
Author

b-niu commented Feb 12, 2022

更新一些信息:

目前已经定位,报错位于sig_tally()函数。将chrM染色体的变异移除后,程序可以成功运行。
仔细检查了下,未发现chrM上有变异的位点大于chrM的染色体size。希望这些信息能有帮助。

@ShixiangWang
Copy link
Owner

@b-niu 谢谢report。这个可能是基因组版本不完全兼容造成的。

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants