BSA,重测序完整篇

重要声明:许多工具、文档是需要爬墙才能看到的,链接无法打开时,可能是被墙了。
BSA综述2016年发表
0.批量转移文件的方法
find ./bwa -name '*2447-20*'|xargs -i mv {} ./bwa/first-2447
把bwa目录的所有文件名包含2447-20的文件全部转移到./bwa/first-2447文件夹中。

流程图


BSA分析流程图

一般流程

  1. 第一部分内容是从测序数据下机数据到变异的SNP检出
    GATK4的使用

2.BSA分析
SHORE-map分析法
SHORE-MAP安装
QTL-seq法
Annovar 软件注释流程介绍
玉米基因组下载

  1. shore-map使用
SHOREmap convert --marker 119-8.vcf --folder path/to/folder -runid int
SHOREmap extract --chrsizes chromsize.txt --folder ../SHOREmap_analysis --marker 11_converted_variant.txt --consen 11_converted_consen.txt -verbose
SHOREmap backcross --chrsizes chromsize.txt --marker ../convert/11_converted_variant.txt --consen extracted_consensus_0.txt --folder ../BC_analysis -plot-bc --marker-score 40 --marker-freq 0.0 --min-converage 10 --max-coverage 80 -bg ../convert/12_converted_variant.txt  --bg-cov 1 --bg-freq 0.0 --bg-score 1 -non-EMS --cluster 1 --marker-hit 1 -verbose
SHOREmap annotate --chrsizes chromsize.txt --folder ../BC_analysis/ann --snp ../convert/11_converted_variant.txt --chrom 2 --start 1 --end 4000000 --genome ../../TAIR10_chr_all.fas --gff ../../TAIR10_GFF3_genes.gff
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容

  • 注意:版本不同,命令会不一致。一定要用对应的版本。 需要登录google账户 GATK3的下载地址[https:/...
    wo_monic阅读 16,921评论 12 57
  • 全基因组重测序是对已知基因组序列的物种进行不同个体的基因组测序,并在此基础上对个体或群体进行差异性分析。所以首先我...
    lizg阅读 12,610评论 1 35
  • 刘小泽写于18.7.18表型与基因型相关联(基因定位)是遗传学上的重要问题。 什么是重测序Resequencing...
    刘小泽阅读 23,228评论 5 47
  • wes定义: 全外显子组测序,是利用目标序列捕获技术, 将全基因组编码基因外显子区域的DNA捕获并富集后,进行高通...
    凤凰_0949阅读 4,511评论 0 7
  • 从知道这个词开始就对这个词没有太多的好感,一直不清楚大家口中的三观到底是哪三观,世界观?人生观?价值观?那宇宙观呢...
    长马阅读 338评论 0 7