Changes between Version 13 and Version 14 of SOPs/vcf_manipulation


Ignore:
Timestamp:
10/26/18 12:14:46 (6 years ago)
Author:
gbell
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • SOPs/vcf_manipulation

    v13 v14  
    101101Extract fields from a VCF file to a TXT, tab separated format with [http://snpeff.sourceforge.net/SnpSift.html#Extract snpSift extractFields] function
    102102{{{
    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.vcf
     103# 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):
     104java -jar SnpSift.jar extractFields foo.vcf CHROM POS REF ALT "GEN[*].ADF" "GEN[*].ADR" > Foo.ADF.ADR.vcf
    105105
    106 # Similar command, but extract and split all PL fields
    107 java -jar /nfs/BaRC/apps/snpEff/Version_4/SnpSift.jar extractFields - CHROM POS REF ALT "GEN[*].PL[*]" foo.vcf > PL_values.vcf
     106# Similar command, but extract and split all ADF and ADR fields
     107java -jar SnpSift.jar extractFields foo.vcf CHROM POS REF ALT "GEN[*].ADF[*]" "GEN[*].ADR[*]" > Foo.ADF.ADR.split.vcf
    108108}}}