a-ludi / a-ludi/dentist

no gap-filling in "scaffold" mode?

未关闭
#33 11 条评论 0 个 reaction 已指派 0 人 在 GitHub 查看
主要语言
D
星标
52
派生
6
PR 合并指标
30 天内没有已合并 PR

描述

Hi again!

I have now run dentist successfully a few times and tested the different join policies, using as input the raw contig assembly or a scaffold assembly, scaffolded with `LRscaf`.

I get the best N50 running dentist with `join-policy: scaffolds` on the already scaffolded assembly (13.5 Mb). However, the final `gap-closed.fasta` contains almost as many `N`s (186183) as the input assembly `reference.fasta` (186185).

When running with `join-policy: contigs`, neither the reference.fasta nor the gap-closed.fasta contain any `N`s.

Is it intended that gaps are not closed in the `join-policy:scaffold`-mode? Do I have to run dentist a second time with `join-policy: scaffoldGaps` to actually glose the gaps in the 13.5 Mb assembly?

As far as I can tell, all dentist runs finished without any errors.

Here are some stats of the different assemblies:

### `join-policy: scaffolds`
```
file format type num_seqs sum_len min_len avg_len max_len Q1 Q2 Q3 sum_gap N50 Q20(%) Q30(%)
gap-closed.fasta FASTA DNA 117 340,348,369 1,843 2,908,960.4 28,227,288 9,395 77,974 2,243,295 0 13,556,940 0 0
reference.fasta FASTA DNA 122 340,337,486 940 2,789,651.5 28,227,288 8,609 89,797 2,243,295 0 12,859,076 0 0
```

### `join-policy: contigs`
```
file format type num_seqs sum_len min_len avg_len max_len Q1 Q2 Q3 sum_gap N50 Q20(%) Q30(%)
gap-closed.fasta FASTA DNA 161 303,401,106 1,069 1,884,478.9 28,227,288 7,516 31,549 988,561 0 10,239,100 0 0
reference.fasta FASTA DNA 171 303,403,113 520 1,774,287.2 25,232,339 7,314 31,549 971,039.5 0 9,931,962 0 0
```

贡献指南

这个仓库没有索引到贡献指南

调研方向

该 issue 涉及 'scaffold' 模式下的间隙填充行为。检查代码库中的 join-policy 逻辑,它可能位于主处理模块或配置处理器中。查找 'scaffold' 和 'scaffoldGaps' 模式,以了解关闭间隙的流程。使用 scaffolded input 运行测试,以复现 N-count 问题,并验证间隙是否正在被处理。

由索引模型根据 Issue 内容生成。

评估

领域
bioinformatics
Issue 类型
缺陷
难度
3/5
预计耗时
1-2 天
活跃度
停滞
描述清晰度
描述清楚
新手友好度
45/100

把新 issue 发到你的邮箱

精选适合新手参与的 GitHub issue 摘要。