./pseq ex1.vcf v-view --geno --vmeta --gmeta

./pseq ex1.vcf v-view --geno --vmeta --gmeta --mask indiv.ex=P002,P004 include=" VM != 2"

./pseq ~/mini/proj1 new-project --vcf ex1.vcf ex2.vcf --metameta ~pshaun/pseq/example/mini/proj1.meta

./pseq proj1 load-vcf

./pseq proj1 vardb-summary

./pseq proj1 v-view --geno

## In R

source("~pshaun/pseq/R/pseq.R")

pseq.project("proj1")

k <- var.fetch()

length(k)

str(k[[1]])

x.consensus.core( k , "REF" )

x.consensus.genotype( k )

x.consensus.genotype( k , "GM" )