Skip to content
A VCF filter for WGS data to filter FPs without losing too many TPs
Shell Python Perl
Branch: master
Clone or download

Latest commit

Fetching latest commit…
Cannot retrieve the latest commit at this time.

Files

Permalink
Type Name Latest commit message Commit time
Failed to load latest commit information.
FPfilter-v1
conda_packages
LICENSE.txt
README.md

README.md

#FPfilter

  • FPfilter is a false-positive specific filter for vcf files from whole genome sequencing. It is proved to be more false-positive specific than GATK hard filter
  • ##Publications

    ##Test Datasets

  • The test vcf file "PE150-example.vcf" is provided in the script folder.
  • ##Setup

  • To run FPfilter, you need to have anaconda installed in your Linux system.
  • After that, you need to get FPfilter from binstar and create the running environment by:
  • conda create -c yuxiang fpfilter -n $env_name #(Recommendation: use FPF with version number as the env_name.)
  • Once it is setup correctly, it will tell you how to activate the environment. For example: source activate $env_name
  • ##Input Data Format

  • Regular vcf, such as from GATK.
  • ##How to Run

  • Before running FPfilter, you need to activate the FPfilter environment by:
  • source activate $env_name #(The one used in the setup step)
  • Running FPfilter is as simple as running a single command line by just typing "FPfilter -v vcf-file -p FPfilter_path".
  • [Details of all parameters: please see the help information in each FPfilter version for latest updates.]
  • *Note1: the output directory is the folder in which the vcf file is located. In order to avoid running error, the current working folder should be the folder in which the vcf file is located.*
  • *Note2: generally, the FPfilter_path is under path_of_miniconda2/envs/env_name/lib/FPfilter/. You can get the path_of_miniconda2 by "which FPfilter" and change "bin" to "lib" in the path.*
  • ##Parameters ###All the parameters are required.

  • -h help
  • -v The vcf file as input to be filtered *[No default value]
  • -p The path of FPfilter scripts *[No default value]
  • ##Output ###Final Output Stucture

  • ${input_vcf}_TP_after_filtered.vcf: the "TP" locations which passed the filter.
  • ${input_vcf}_filtered_FP.vcf :the "false positive" locations which were filtered.
  • You can’t perform that action at this time.