Changes between Version 13 and Version 14 of SOPs/vcf_manipulation
- Timestamp:
- 10/26/18 12:14:46 (6 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
SOPs/vcf_manipulation
v13 v14 101 101 Extract fields from a VCF file to a TXT, tab separated format with [http://snpeff.sourceforge.net/SnpSift.html#Extract snpSift extractFields] function 102 102 {{{ 103 # Extract chromosome, position, reference base(s), alternate bases(s), all the samples' DP4 (number of high-quality ref-fwd, ref-reverse, alt-fwd,alt-reverse bases):104 java -jar /nfs/BaRC/apps/snpEff/Version_4/SnpSift.jar extractFields - CHROM POS REF ALT GEN[*].DP4 foo.vcf > DP4.vcf103 # Extract chromosome, position, reference base(s), alternate bases(s), all the samples' ADF and ADR (number of high-quality ref-fwd, alt-fwd, ref-reverse, and alt-reverse bases): 104 java -jar SnpSift.jar extractFields foo.vcf CHROM POS REF ALT "GEN[*].ADF" "GEN[*].ADR" > Foo.ADF.ADR.vcf 105 105 106 # Similar command, but extract and split all PLfields107 java -jar /nfs/BaRC/apps/snpEff/Version_4/SnpSift.jar extractFields - CHROM POS REF ALT "GEN[*].PL[*]" foo.vcf > PL_values.vcf106 # Similar command, but extract and split all ADF and ADR fields 107 java -jar SnpSift.jar extractFields foo.vcf CHROM POS REF ALT "GEN[*].ADF[*]" "GEN[*].ADR[*]" > Foo.ADF.ADR.split.vcf 108 108 }}}