Basic statistics for a given hapmap file. Including genotype number, missing rate, het rate, maf etc.
Usage:
perl basicStatHapmap.pl test1.hapmap ./test1.output
test1.hapmap
test1.output.site.stat, test1.output.ind.stat
Extract hapmapSNPs genotype from the interested SNP set, when user provides a list of SNP positions in rice genome.
Usage:
perl extractTarHapmapPos.pl test2.snplist test2.hapmap > test2.output
test2.snplist, test2.hapmap
test2.output
Extract flanking sequences from the interested SNP set, when user provides a list of SNP positions in rice genome. (for primer design)
Usage:
perl extractSNPflankSeq.pl ../SR4Rdatabase/IRGSP-1.0_genome.fasta test3.snplist 50 >
test3.output
IRGSP-1.0_genome.fasta, test3.snplist
test3.output
flank-seq-length=[50/100/150/200/250]
Convert a provided hapmap file to ped format.
Usage:
perl hapmap2plink.pl test4.hapmap ./test4.output
test4.hapmap
test4.output.ped, test4.output.map
Convert a provided hapmap file to vcf format.
Usage:
perl hapmap2vcf.pl test5.hapmap > test5.vcf
test5.hapmap
test5.vcf
Convert a provided hapmap file to bgl format.
Usage:
perl hapmap2beagle.pl test6.hapmap > test6.bgl
test6.hapmap
test6.bgl
Imputation of missing genotype and phasing for a provided hapmap file (Please performed on each chromosome separately).
Usage:
perl ../1.6.hapmap2beagle/hapmap2beagle.pl test7.hapmap > test7.bgl
java -Xmx215180m -XX:MaxPermSize=215180m -jar beagle.jar unphased=test7.bgl missing=N
out=test7.bgl
perl beaglephasedgz2hapmap.pl test7.hapmap test7.bgl.test7.bgl.phased.gz
>test7.output.hapmap
test7.hapmap
test7.output.hapmap
missing=N