获取单个基因的fst信息
摘要：#!/bin/bash echo "par1 is gene,par2 is fst_file" #awk '{(if ($3==gene) print $0}' $2>gene_info grep -i $1 gene_info|awk '{print $2,$3,$4}'>${1}_info w
阅读全文
从significant snp 的LD 范围内的SNP，进行finemap定位分析
摘要：#!/bin/bash #set -Eeuo pipefail if [ $# -lt 2 ] then echo "par1 is vcf,par2 is pvalue,sample_size=297" exit fi #vcftools --vcf $1 --freq --out ${1}maf
阅读全文
利用peak_snp，找其LD块，并找到其中的基因，然后将LD中的SNP信息变成finemap软件的输入格式
摘要：#!/bin/bash #LD PLOT #author lee echo"par1 is vcf,par2 is top_chr" Genome=/public/home/caisl/lee/genome/rice/msu.gff3 while read chr pos do pos_up=$((
阅读全文
