问题描述
我查询了很多地方,发现我研究的物种没有组装完整的染色体,因此它的染色体名称不是standard/typical,以下是我的.ded12文件。解决方法,忽略染色体名称
gene_anno =readTranscriptFeatures("genes.bed12",remove.unusual=FALSE, up.flank = 2000,
down.flank = 200)
加上remove.unusual=FALSE即可。