contig or scaffold GFF 文件通过AGP文件转换成 染色体 gff文件

contig GFF 文件通过AGP文件转换成 染色体 gff文件


注释基因组文件 contig.gff  或者 scaffold.gff  to  chromosome gff

python3  -m jcvi.formats.chain fromagp genome.final.agp  contig.fa  genome.fa


得到:genome.final.chain  再用 liftover转换:http://hgdownload.cse.ucsc.edu/admin/exe/linux.x86_64/liftOver


liftOver -gff contig.gff3 genome.final.chain genome.gff3 unmapped


The arguments are:

  • contig.gff3: scaffold annotations (INPUT)
  • contig_to_genome.chain: translations (INPUT)
  • genome.gff3: genome annotations (OUTPUT)
  • unmapped: things that fail (OUTPUT)


或者用

CrossMap.py gff genome.final.chain contig.gff genome.gff
  • 发表于 2022-11-17 21:14
  • 阅读 ( 1666 )
  • 分类:基因组学

你可能感兴趣的文章

相关问题

0 条评论

请先 登录 后评论
omicsgene
omicsgene

生物信息

698 篇文章

作家榜 »

  1. omicsgene 698 文章
  2. 安生水 347 文章
  3. Daitoue 167 文章
  4. 生物女学霸 120 文章
  5. xun 82 文章
  6. 红橙子 78 文章
  7. rzx 74 文章
  8. CORNERSTONE 72 文章