no gap-filling in "scaffold" mode?
- 主要语言
- 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