该问题已被锁定!
4
关注
1890
浏览

转录本坐标转换成基因组坐标

查看全部 2 个回答

chatGPT机器人 机器人 用户来自于: IANA
2023-06-17 20:20
首先,将转录本坐标转换成基因组坐标需要考虑基因组中的剪切变异和多个转录本的存在。您的代码中使用了GenomicFeatures包中的mapFromTranscripts函数进行坐标转换,这是正确的方法。 然而,您提到转换后得到的坐标与基因组上的碱基不对应的问题可能是因为您的代码中没有考虑转录本的方向性。您可以使用strand参数指定转录本的方向,例如: ```map2genome <- mapFromTranscripts(gr, transcripts, strand = strand(gr))``` 另外,您提到转换后得到的基因组坐标与RMbase数据库中的修饰信息做overlap得到的数量很少,可能是因为RMbase数据库的修饰信息是基于特定转录本的,而您的代码中使用的是所有转录本的坐标进行转换。您可以尝试只使用特定转录本的坐标进行转换,例如: ```gr <- GRanges(seqnames = "ENST00000281150", ranges = IRanges(start = 123456, end = 123457), strand = "+")``` 其中seqnames参数指定需要转换的转录本,ranges参数指定需要转换的m6A修饰位置。 希望这些建议能帮助您解决问题。

关于作者

问题动态

发布时间
2023-06-17 20:20
更新时间
2023-06-19 15:43
关注人数
4 人关注

相关问题

已经有不同材料的全基因组二代测序数据,查找是否存在某个已知基因?
参考基因组添加外源基因序列进行比对
基因组文件中chrUn具体是什么意思呢
耐药细菌的比较基因组分析流程和机制探索
绵羊转录组测序后差异表达基因太多了
seurat中如何通过坐标/细胞名删掉细胞
3D基因组里,如何分清compartment、TAD和chromatin loop
annotatePeak peak基因组注释
stringtie 得到的gtf 通过DESeq2分析后stringtieID 如何转换成esmbleID
无参基因组分析GO富集到了这一步不会做

推荐内容

Ancestry_hmm
链特异性文库(mRNA/lncRNA/circRNA)如何将RNA类型分开?
关于基因间的相关性分析
SPSS报错
不同测序平台的scRNAseq数据集整合(10✖️、smart-seq)
关于ceRNA网络构建的后续分析有哪些?
使用Tracking Tumor ImmunoPhenotype(TIP)网站分析TCGA的BLCA_tpm数据
3D基因组里compartment里一般是包含好多TAD的,但这图为什么compartment数量比TAD多这么多呢
HHblits使用相同的数据库,得出的结果明显差于HHpred网页版的结果。该如何使用HHsuite命令行得到与网页端类似的效果
All Rights Reserved Powered BY WeCenter V4.1.0 © 2025