求助mRNA ID 批量+1问题

你好，我在分析转录组数据是出现如下问题：

查资料发现你也遇到类似问题，看老师给出的解决办法：mRNA ID和基因ID相同了，你把所有的mRNAID后面批量加一个.1，请问你这个问题解决了吗？

0 条评论
分类：转录组

默认排序时间排序

2 个回答

omicsgene - 生物信息 2020-11-12 17:56

擅长：重测序,遗传进化,转录组,GWAS

可以自己写脚本修改，

或者看看这个修改方法：https://www.omicsclass.com/article/816

Sirius 2020-11-13 15:00

老师，利用这个脚本修改，awk -F"[\t=]" 'BEGIN{OFS="\t"}{if($3=="mRNA"){print $1,$2,"gene",$4,$5,$6,$7,$8,$9"="$10 ; print $1,$2,$3,$4,$5,$6,$7,$8,$9"="$10"Parent="$10;}else{print $0}}' PO.gene.gff >PO.gene.gff1

可否提示一下这个脚本修改那一部分可以实现mRNAID后面批量加一个.1

求助mRNA ID 批量+1问题

2 个回答

相似问题