61
edits
Changes
From Genome Analysis Wiki
Jump to navigationJump to searchSeqShop: Variant Calling and Filtering for SNPs Practical, May 2015 (view source)
Revision as of 18:35, 19 August 2015
, 18:35, 19 August 2015→Filtered VCF: since the uncompressed vcf isn't created, use zless on the compressed
</div>
</div>
=== Prerequisite Tutorials ===
This tutorial builds on the alignment tutorial, if you have not already, please first run that tutorial: [[SeqShop:_Sequence_Mapping_and_Assembly_Practical|Alignment Tutorial]]
This tutorial builds on the alignment tutorial, if you have not already, please first run that tutorial: [[SeqShop:_Sequence_Mapping_and_Assembly_Practical|Alignment Tutorial]]
Let's look at the filtered sites file.
Let's look at the filtered sites file.
zless -S ${OUT}/vcfs/chr22/chr22.filtered.sites.vcf.gz
* Scroll down until you find some variants.
* Scroll down until you find some variants.
=== Running GotCloud Genotype Refinement ===
=== Running GotCloud Genotype Refinement ===
Since everything is setup, just run the following command (very similar to snpcall).
Since everything is setup, just run the following command (very similar to snpcall).
${GC}/gotcloud ldrefine --conf ${SS}/gotcloud.conf --numjobs 6 --region 22:36000000-37000000 --base_prefix ${SS} --outdir ${OUT}
${GC}/gotcloud ldrefine --conf ${SS}/gotcloud.conf --numjobs 8 --region 22:36000000-37000000 --base_prefix ${SS} --outdir ${OUT}
* Beagle will take about 1-3 minutes to complete
* Beagle will take about 1-3 minutes to complete
<li><code>1|1</code>; Homozygous Alt (C)</li>
<li><code>1|1</code>; Homozygous Alt (C)</li>
</ul>
</ul>
Note the '|' instead of '/'. The '|' indicates it is now phased.
</div>
</div>
</div>
</div>
