Changes between Version 19 and Version 20 of SOPs/variant_calling
- Timestamp:
- 07/18/17 08:14:47 (8 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
SOPs/variant_calling
v19 v20 58 58 == Call raw variants with mpileup+bcftools == 59 59 60 Call variants (one sample vs. reference) with samtools' mpileup+bcftools (see the [http:// samtools.sourceforge.net/mpileup.shtml samtools' variant calling page] for more details).60 Call variants (one sample vs. reference) with samtools' mpileup+bcftools (see the [http://www.htslib.org/workflow/#mapping_to_variant samtools' variant calling workflow] for more details). 61 61 In our experience, "-B" (disable BAQ) or "-E" (recalculate BAQ) works better than the default method, which can remove some obvious variants. 62 62 {{{ … … 71 71 == Add filters and annotations to raw variants == 72 72 73 This step uses vcf-annotate from the [http ://vcftools.sourceforge.net/docs.html VCFtools suite]73 This step uses vcf-annotate from the [https://vcftools.github.io/perl_module.html VCFtools suite] 74 74 75 75 Annotate variants by adding tags ("filters" but all variants are kept) to each variant, using all default filters.