在ggcoverage中绘制基因结构的注释图的三个问题

[wsweishang]

1）读取的gff文件需要包含哪些基因信息，才能正常绘制？ 原本用烟草参考基因组注释文件画图时会报错：

> basic.coverage + geom_gene(gtf.gr = gtf.gr)
Error in `dplyr::select()`:
! Can't subset columns that don't exist.
✖ Column `gene_type` doesn't exist.
Run `rlang::last_error()` to see where the error occurred.
> 

对比示例gtf文件，发现缺少一些基因注释信息，手动补充了一下，可以画出来了，但是读取的gff文件需要包含哪些基因信息，才能正常绘制？

2）目标基因有7个外显子，为何只显示六个？ 目标基因包含有7个外显子，为何只显示六个？图中的各形状分别代表哪些基因信息？此外是否有图例参考？

3）目标基因所在的基因组区间也不对 目标基因位于421068-432552，为何图上从420000之前就开始？

4）可否提供您的邮箱以便发送代码和输入文件测试？

[showteeth]

Thanks for your feedback!

E-mail