Changes

From Genome Analysis Wiki
Jump to navigationJump to search
6 bytes removed ,  21:45, 26 January 2015
Line 169: Line 169:  
SHAPEIT can be used to phase a small number of samples (a reference panel required) as well as a large number of samples (NO reference panel required). The reference panels and genetic map files are available for download on their [https://mathgen.stats.ox.ac.uk/genetics_software/shapeit/shapeit.html home-page].
 
SHAPEIT can be used to phase a small number of samples (a reference panel required) as well as a large number of samples (NO reference panel required). The reference panels and genetic map files are available for download on their [https://mathgen.stats.ox.ac.uk/genetics_software/shapeit/shapeit.html home-page].
   −
* The following example shows a typical SHAPEIT command line to phase a LARGE number (>200) of GWAS samples (<code>Gwas.chr20.Unphased.vcf</code> is the quality control GWAS data set in VCF format).
+
* The following example shows a typical SHAPEIT command line to phase a LARGE number (>200) of GWAS samples (<code>Gwas.chr20.Unphased.vcf</code> is the quality controlled GWAS data set in VCF format).
    
  shapeit -V Gwas.chr20.Unphased.vcf \
 
  shapeit -V Gwas.chr20.Unphased.vcf \
Line 175: Line 175:  
         -O Gwas.Chr20.Phased.Output
 
         -O Gwas.Chr20.Phased.Output
   −
* The following example shows a typical SHAPEIT command line to phase a SMALL number (<200) of GWAS samples (<code>Gwas.chr20.Unphased.vcf</code> is the quality control GWAS data set in VCF format).
+
* The following example shows a typical SHAPEIT command line to phase a SMALL number (<200) of GWAS samples (<code>Gwas.chr20.Unphased.vcf</code> is the quality controlled GWAS data set in VCF format).
   −
  ## The following step splits out variants mis-aligned between the reference and target/gwas panel
+
  ## The following step splits out variants mis-aligned between the reference and gwas panel
 
  shapeit -check \
 
  shapeit -check \
 
         -V Gwas.chr20.Unphased.vcf\
 
         -V Gwas.chr20.Unphased.vcf\
Line 183: Line 183:  
         --input-ref reference.haplotypes.gz reference.legend.gz reference.sample \
 
         --input-ref reference.haplotypes.gz reference.legend.gz reference.sample \
 
         --output-log gwas.alignments
 
         --output-log gwas.alignments
  ## The following step phases target/gwas panel using the reference panel while excluding the markers found in the step above.
+
 +
  ## The following step phases gwas panel using the reference panel while excluding the markers found in the step above.
 
  shapeit -B gwas \
 
  shapeit -B gwas \
 
         -V Gwas.chr20.Unphased.vcf \
 
         -V Gwas.chr20.Unphased.vcf \
487

edits

Navigation menu