美文网首页
水稻基因组变异注释

水稻基因组变异注释

作者: 花生学生信 | 来源:发表于2023-03-26 18:37 被阅读0次

太玄学了,仅供参考和日后借鉴。

运行命令:

###生成refgene 
~/tools/gtfToGenePred -genePredExt Oryza_sativa.IRGSP-1.0.50.gtf gcf_refGene.txt # *_refGene.txt,*

###fa
~/tools/annovar/retrieve_seq_from_fasta.pl --format refGene --seqfile IRGSP-1.0_genome.fasta gcf_refGene.txt --outfile gcf_refGeneMrna.fa
### 转annovar
 ~/tools/annovar/convert2annovar.pl -format vcf4 sample.vcf > sample.annovar
###注释命令
$> /public/home/fengting/tools/annovar/annotate_variation.pl -buildver gcf -geneanno -outfile sample.anno sv.annovar ricedb/

我看了下,只注释了后三条染色体,可能是染色体编号的原因,把前9条染色体名改一下:

sed -i "s/chr0/chr/g" sv.annovar 
okk,前9条染色体也出来了

接下来是提取相应区间的变异信息:

cat pav-gene-pos.txt | while read chr start end id
do


##awk '{if($3 == "'"$chr"'" && $3 == "gene" && $4 >= "'"${start}"'" && $5 <= "'"${end}"'")print } ' sample.anno.variant_function > tmp1.txt
awk '{if($3 == "'"$chr"'" && $4 >= "'"${start}"'" && $5 <= "'"${end}"'")print } ' sample.anno.variant_function > tmp1.txt

cat tmp1.txt |while read line1
do
    echo $line1 >> res/"${id}".csv
done

##awk '$1 == "1" && $3 == "gene" && $4 >= 10000 && $5 <= 200000 {print $0} ' Arabidopsis_thaliana.TAIR10.31.gff3 > out.gff

done
改成csv文件

借鉴,防遗忘笔记:
annovar基因(水稻)注释 - 知乎 (zhihu.com)
用annovar对vcf(SNP&INDEL)文件进行注释 - 简书 (jianshu.com)

相关文章

  • ANNOVAR软件使用

    ANNOVAR是一个高效的注释工具。它从不同的基因组功能注释基因变异检测(包括人类基因组hg18、hg19 、hg...

  • ANNOVAR 注释|自建数据库

    ANNOVAR|绵羊基因组变异注释 Gene-based Annotation ANNOVAR是由王凯老师编写的一...

  • RiceENCODE:水稻多元表观基因组数据库

    水稻是我国重要的粮食作物,也是生物学研究的重要模式植物。水稻基因组DNA顺式调控元件的注释和鉴定,对理解水稻基因表...

  • 使用bioawk对fasta文件进行序列分析

    最近在做水稻基因组注释,对于注释结果进行整理过滤,还在探索中。 evm注释结果 将结果转化成蛋白文件: 提取长度大...

  • annnovar——壹

    annovar介绍 ANNOVAR 是由王凯编写一款变异位点注释软件,利用最新信息功能注释从不同基因组检测到的遗传...

  • ANNOVAR gene-based annotation

    欢迎关注"生信修炼手册"! 通过基因相关注释,可以知道变异位点在基因组上的位置和对蛋白质编码的影响。在进行注释之前...

  • 变异检测与其相应的工具的选择

    什么是基因组变异 基因组变异是一个定义比较模糊的概念。 所谓的变异是相对于一个完美的“参考基因组”而言。但是其实完...

  • biostar handbook(十)|如何进行变异检测

    什么是基因组变异 基因组变异是一个定义比较模糊的概念. 所谓的变异是相对于一个完美的“参考基因组”而言。但是其实完...

  • 基础知识 — 基因组变异类型

    基因组变异是一个定义比较模糊的概念。 所谓的变异是相对于一个完美的“参考基因组”而言。但是其实完美的“参考基因组”...

  • 【基因组学工具】变异注释SnpEff

    写在前面在进行基因组学研究中,总避免不了对一些SNP和Indel进行分析,如何去注释它们的,就需要用到SnpEff...

网友评论

      本文标题:水稻基因组变异注释

      本文链接:https://www.haomeiwen.com/subject/xqxardtx.html