标签：gff3 基因组 overlap gmap 注释 intersect Bedtools gene bedtools
目录问题
原问题来自：How to calculate overlapping genes between two genome annotation versions？
其实可分为两个问题：
- 一是我组装了一个新的基因组，做了多个注释版本，如何比较它们的feature？比如gene
- 二是我组装了一个新的参考基因组，并做了注释，想和其他已有的同物种参考基因组比较，如何寻找共有和特有的基因（或其他feature）？
思路
第一个问题是比较好解决的，使用bedtools即可。
bedtools比较gff、bed、bam的方法类似，具体可参考这篇教程：
bedtools求overlap
要比较gene，可先从gff中提取gene后再进行比较。或者比较所有feature后再筛选也行。
# 将所有overlap 区域成对输出
bedtools intersect -a A.gene.gff3 -b B.gene.gff3 -wa -wb >gene_wa_wb.out
#只要A中的这段区域与B中区域有交集，就输出，而且overlap几次，就输出几次
bedtools intersect -a A.gene.gff3 -b B.gene.gff3 -wa >gene_wa.out
#除了输出A中的overlap区域外，还会输出B中的整个区间
bedtools intersect -a A.gene.gff3 -b B.gene.gff3 -wb >gene_wb.out
#统计A中每个区域与B overlap的次数
bedtools intersect -a A.gene.gff3 -b B.gene.gff3 -c >gene_overlap.count
#只输出A中没有与B overlap的区域
bedtools intersect -a A.gene.gff3 -b B.gene.gff3 -v >gene_nonoverlap.count
bedtools intersect -a B.gene.gff3 -b A.gene.gff3 -v >gene_msu_uniq.count
第二个问题需要用比对软件，如gmap进行比对，建立两个基因组的联系，得到gff文件。再利用bedtools比较。
/gmap/bin/gmap_build -D ./ -d A A.fa
/gmap/bin/gmap -D ./ -t 30 -d A -f gff3_gene ../B.cdna > B.gff3
最后的结果要注意，feature不是一一对应的，有一对多，多对一，unique等情况。
标签：gff3,基因组,overlap,gmap,注释,intersect,Bedtools,gene,bedtools 来源： https://www.cnblogs.com/jessepeng/p/14823527.html
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。