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

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

查看全部 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具体是什么意思呢
全基因组关联分析结果与模型选择
根据已知序列查找其第三代基因组编号
请问hemmer是否可以使用转录组搜索HMM库?
seurat中如何通过坐标/细胞名删掉细胞
知道一批基因的具体位置,如何批量的从基因组中取出基因序列?
请教一个问题,一个转录因子TF只会调控位于同一条链的基因,还是有可能调控反义链的基因?
非已知转录因子预测
无对照的转录组数据如何寻找差异表达基因
转录组和代谢组联合分析的方法有那些?

推荐内容

L1插入位点的百分比密度分布
sc-ATAC数据质控
Ancestry_hmm
去除批次效应选定靶标细胞群后如何进行亚群分析?
3D基因组里compartment里一般是包含好多TAD的,但这图为什么compartment数量比TAD多这么多呢
infercnv运行报错
SPSS报错
染色体号是罗马数字怎么写sh循环
qpAdm结果解读
All Rights Reserved Powered BY WeCenter V4.1.0 © 2025