pax_global_header 0000666 0000000 0000000 00000000064 13331057401 0014507 g ustar 00root root 0000000 0000000 52 comment=99c96e52d046db580495e88bde6083d63e641040
vcftools-0.1.16/ 0000775 0000000 0000000 00000000000 13331057401 0013433 5 ustar 00root root 0000000 0000000 vcftools-0.1.16/.gitignore 0000664 0000000 0000000 00000000664 13331057401 0015431 0 ustar 00root root 0000000 0000000 # created by autoscan
/autom4te.cache
/autoscan.log
/configure.scan
# created by autoreconf -fi (includes autoconf, automake, libtool)
Makefile.in
/aclocal.m4
/autom4te.cache
/compile
/config.guess
/config.h.in*
/config.sub
/configure
/depcomp
/install-sh
/ltmain.sh
/m4
/missing
# created by configure
Makefile
/config.h
/config.log
/config.status
/libtool
/stamp-h1
.deps
# created by make
*.la
*.lo
*.o
.libs
vcftools
*.tar.*
vcftools-0.1.16/.tarball-version 0000664 0000000 0000000 00000000007 13331057401 0016535 0 ustar 00root root 0000000 0000000 0.1.16
vcftools-0.1.16/LICENSE 0000664 0000000 0000000 00000016744 13331057401 0014454 0 ustar 00root root 0000000 0000000 GNU LESSER GENERAL PUBLIC LICENSE
Version 3, 29 June 2007
Copyright (C) 2007 Free Software Foundation, Inc.
Everyone is permitted to copy and distribute verbatim copies
of this license document, but changing it is not allowed.
This version of the GNU Lesser General Public License incorporates
the terms and conditions of version 3 of the GNU General Public
License, supplemented by the additional permissions listed below.
0. Additional Definitions.
As used herein, "this License" refers to version 3 of the GNU Lesser
General Public License, and the "GNU GPL" refers to version 3 of the GNU
General Public License.
"The Library" refers to a covered work governed by this License,
other than an Application or a Combined Work as defined below.
An "Application" is any work that makes use of an interface provided
by the Library, but which is not otherwise based on the Library.
Defining a subclass of a class defined by the Library is deemed a mode
of using an interface provided by the Library.
A "Combined Work" is a work produced by combining or linking an
Application with the Library. The particular version of the Library
with which the Combined Work was made is also called the "Linked
Version".
The "Minimal Corresponding Source" for a Combined Work means the
Corresponding Source for the Combined Work, excluding any source code
for portions of the Combined Work that, considered in isolation, are
based on the Application, and not on the Linked Version.
The "Corresponding Application Code" for a Combined Work means the
object code and/or source code for the Application, including any data
and utility programs needed for reproducing the Combined Work from the
Application, but excluding the System Libraries of the Combined Work.
1. Exception to Section 3 of the GNU GPL.
You may convey a covered work under sections 3 and 4 of this License
without being bound by section 3 of the GNU GPL.
2. Conveying Modified Versions.
If you modify a copy of the Library, and, in your modifications, a
facility refers to a function or data to be supplied by an Application
that uses the facility (other than as an argument passed when the
facility is invoked), then you may convey a copy of the modified
version:
a) under this License, provided that you make a good faith effort to
ensure that, in the event an Application does not supply the
function or data, the facility still operates, and performs
whatever part of its purpose remains meaningful, or
b) under the GNU GPL, with none of the additional permissions of
this License applicable to that copy.
3. Object Code Incorporating Material from Library Header Files.
The object code form of an Application may incorporate material from
a header file that is part of the Library. You may convey such object
code under terms of your choice, provided that, if the incorporated
material is not limited to numerical parameters, data structure
layouts and accessors, or small macros, inline functions and templates
(ten or fewer lines in length), you do both of the following:
a) Give prominent notice with each copy of the object code that the
Library is used in it and that the Library and its use are
covered by this License.
b) Accompany the object code with a copy of the GNU GPL and this license
document.
4. Combined Works.
You may convey a Combined Work under terms of your choice that,
taken together, effectively do not restrict modification of the
portions of the Library contained in the Combined Work and reverse
engineering for debugging such modifications, if you also do each of
the following:
a) Give prominent notice with each copy of the Combined Work that
the Library is used in it and that the Library and its use are
covered by this License.
b) Accompany the Combined Work with a copy of the GNU GPL and this license
document.
c) For a Combined Work that displays copyright notices during
execution, include the copyright notice for the Library among
these notices, as well as a reference directing the user to the
copies of the GNU GPL and this license document.
d) Do one of the following:
0) Convey the Minimal Corresponding Source under the terms of this
License, and the Corresponding Application Code in a form
suitable for, and under terms that permit, the user to
recombine or relink the Application with a modified version of
the Linked Version to produce a modified Combined Work, in the
manner specified by section 6 of the GNU GPL for conveying
Corresponding Source.
1) Use a suitable shared library mechanism for linking with the
Library. A suitable mechanism is one that (a) uses at run time
a copy of the Library already present on the user's computer
system, and (b) will operate properly with a modified version
of the Library that is interface-compatible with the Linked
Version.
e) Provide Installation Information, but only if you would otherwise
be required to provide such information under section 6 of the
GNU GPL, and only to the extent that such information is
necessary to install and execute a modified version of the
Combined Work produced by recombining or relinking the
Application with a modified version of the Linked Version. (If
you use option 4d0, the Installation Information must accompany
the Minimal Corresponding Source and Corresponding Application
Code. If you use option 4d1, you must provide the Installation
Information in the manner specified by section 6 of the GNU GPL
for conveying Corresponding Source.)
5. Combined Libraries.
You may place library facilities that are a work based on the
Library side by side in a single library together with other library
facilities that are not Applications and are not covered by this
License, and convey such a combined library under terms of your
choice, if you do both of the following:
a) Accompany the combined library with a copy of the same work based
on the Library, uncombined with any other library facilities,
conveyed under the terms of this License.
b) Give prominent notice with the combined library that part of it
is a work based on the Library, and explaining where to find the
accompanying uncombined form of the same work.
6. Revised Versions of the GNU Lesser General Public License.
The Free Software Foundation may publish revised and/or new versions
of the GNU Lesser General Public License from time to time. Such new
versions will be similar in spirit to the present version, but may
differ in detail to address new problems or concerns.
Each version is given a distinguishing version number. If the
Library as you received it specifies that a certain numbered version
of the GNU Lesser General Public License "or any later version"
applies to it, you have the option of following the terms and
conditions either of that published version or of any later version
published by the Free Software Foundation. If the Library as you
received it does not specify a version number of the GNU Lesser
General Public License, you may choose any version of the GNU Lesser
General Public License ever published by the Free Software Foundation.
If the Library as you received it specifies that a proxy can decide
whether future versions of the GNU Lesser General Public License shall
apply, that proxy's public statement of acceptance of any version is
permanent authorization for you to choose that version for the
Library.
vcftools-0.1.16/Makefile.am 0000664 0000000 0000000 00000000404 13331057401 0015465 0 ustar 00root root 0000000 0000000 SUBDIRS = src
EXTRA_DIST = LICENSE README.md examples
# Create a '.tarball-version' file containing the version string
# when creating a distribution tarball
# (not needed when using git repository)
dist-hook:
@echo $(VERSION) > $(distdir)/.tarball-version
vcftools-0.1.16/README.md 0000664 0000000 0000000 00000003610 13331057401 0014712 0 ustar 00root root 0000000 0000000 # VCFtools
A set of tools written in Perl and C++ for working with [VCF files](https://samtools.github.io/hts-specs/VCFv4.2.pdf), such as those generated by the
[1000 Genomes Project](http://www.1000genomes.org/).
Project website: https://vcftools.github.io/
License
-------
The program package is released under the GNU Lesser General Public License version 3.0
(LGPLv3). See the `LICENSE` file for the complete LGPL license text.
Credits
-------
- Adam Auton (Binary Executable)
- Petr Danecek (Perl Module)
- Anthony Marcketta (Binary Executable)
Building VCFtools
-----------------
General help about the building process's configuration step can be acquired via:
```
./configure --help
```
### Build from Release Tarball
```
./configure
make
make install
```
You may need `sudo` permissions to run `make install`.
### Build from GitHub
```
git clone https://github.com/vcftools/vcftools.git
cd vcftools
./autogen.sh
./configure
make
make install
```
You many need `sudo` permissions to run `make install`.
Documentation
-------------
Documentation and usage examples can be found here:
https://vcftools.github.io/examples.html
A manual page is also available. If prefix is set to `/usr` or if `MANPATH` points to
`$prefix/share/man`, you can access the manual page via:
```
man vcftools
```
Getting Help
------------
The best way to get help regarding VCFtools is to email the mailing list:
vcftools-help@lists.sourceforge.net
Citation
--------
If you make use of VCFtools in your research, we would appreciate a citation of the following paper:
> **The Variant Call Format and VCFtools**, Petr Danecek, Adam Auton, Goncalo Abecasis, Cornelis
> A. Albers, Eric Banks, Mark A. DePristo, Robert Handsaker, Gerton Lunter, Gabor Marth, Stephen
> T. Sherry, Gilean McVean, Richard Durbin and 1000 Genomes Project Analysis Group,
> **Bioinformatics**, 2011 http://dx.doi.org/10.1093/bioinformatics/btr330
vcftools-0.1.16/autogen.sh 0000775 0000000 0000000 00000000032 13331057401 0015427 0 ustar 00root root 0000000 0000000 #!/bin/sh
autoreconf -fi
vcftools-0.1.16/build-aux/ 0000775 0000000 0000000 00000000000 13331057401 0015325 5 ustar 00root root 0000000 0000000 vcftools-0.1.16/build-aux/git-version-gen 0000775 0000000 0000000 00000017571 13331057401 0020303 0 ustar 00root root 0000000 0000000 #!/bin/sh
# Print a version string.
scriptversion=2012-12-31.23; # UTC
# Copyright (C) 2007, 2008, 2009, 2010, 2011, 2012, 2013, 2014, 2015
# Free Software Foundation, Inc.
#
# This program is free software: you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation; either version 3 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program. If not, see .
# This script is derived from GIT-VERSION-GEN from GIT: http://git.or.cz/.
# It may be run two ways:
# - from a git repository in which the "git describe" command below
# produces useful output (thus requiring at least one signed tag)
# - from a non-git-repo directory containing a .tarball-version file, which
# presumes this script is invoked like "./git-version-gen .tarball-version".
# In order to use intra-version strings in your project, you will need two
# separate generated version string files:
#
# .tarball-version - present only in a distribution tarball, and not in
# a checked-out repository. Created with contents that were learned at
# the last time autoconf was run, and used by git-version-gen. Must not
# be present in either $(srcdir) or $(builddir) for git-version-gen to
# give accurate answers during normal development with a checked out tree,
# but must be present in a tarball when there is no version control system.
# Therefore, it cannot be used in any dependencies. GNUmakefile has
# hooks to force a reconfigure at distribution time to get the value
# correct, without penalizing normal development with extra reconfigures.
#
# .version - present in a checked-out repository and in a distribution
# tarball. Usable in dependencies, particularly for files that don't
# want to depend on config.h but do want to track version changes.
# Delete this file prior to any autoconf run where you want to rebuild
# files to pick up a version string change; and leave it stale to
# minimize rebuild time after unrelated changes to configure sources.
#
# As with any generated file in a VC'd directory, you should add
# /.version to .gitignore, so that you don't accidentally commit it.
# .tarball-version is never generated in a VC'd directory, so needn't
# be listed there.
#
# Use the following line in your configure.ac, so that $(VERSION) will
# automatically be up-to-date each time configure is run (and note that
# since configure.ac no longer includes a version string, Makefile rules
# should not depend on configure.ac for version updates).
#
# AC_INIT([GNU project],
# m4_esyscmd([build-aux/git-version-gen .tarball-version]),
# [bug-project@example])
#
# Then use the following lines in your Makefile.am, so that .version
# will be present for dependencies, and so that .version and
# .tarball-version will exist in distribution tarballs.
#
# EXTRA_DIST = $(top_srcdir)/.version
# BUILT_SOURCES = $(top_srcdir)/.version
# $(top_srcdir)/.version:
# echo $(VERSION) > $@-t && mv $@-t $@
# dist-hook:
# echo $(VERSION) > $(distdir)/.tarball-version
me=$0
version="git-version-gen $scriptversion
Copyright 2011 Free Software Foundation, Inc.
There is NO warranty. You may redistribute this software
under the terms of the GNU General Public License.
For more information about these matters, see the files named COPYING."
usage="\
Usage: $me [OPTION]... \$srcdir/.tarball-version [TAG-NORMALIZATION-SED-SCRIPT]
Print a version string.
Options:
--prefix prefix of git tags (default 'v')
--fallback fallback version to use if \"git --version\" fails
--help display this help and exit
--version output version information and exit
Running without arguments will suffice in most cases."
prefix=v
fallback=
while test $# -gt 0; do
case $1 in
--help) echo "$usage"; exit 0;;
--version) echo "$version"; exit 0;;
--prefix) shift; prefix="$1";;
--fallback) shift; fallback="$1";;
-*)
echo "$0: Unknown option '$1'." >&2
echo "$0: Try '--help' for more information." >&2
exit 1;;
*)
if test "x$tarball_version_file" = x; then
tarball_version_file="$1"
elif test "x$tag_sed_script" = x; then
tag_sed_script="$1"
else
echo "$0: extra non-option argument '$1'." >&2
exit 1
fi;;
esac
shift
done
if test "x$tarball_version_file" = x; then
echo "$usage"
exit 1
fi
tag_sed_script="${tag_sed_script:-s/x/x/}"
nl='
'
# Avoid meddling by environment variable of the same name.
v=
v_from_git=
# First see if there is a tarball-only version file.
# then try "git describe", then default.
if test -f $tarball_version_file
then
v=`cat $tarball_version_file` || v=
case $v in
*$nl*) v= ;; # reject multi-line output
[0-9]*) ;;
*) v= ;;
esac
test "x$v" = x \
&& echo "$0: WARNING: $tarball_version_file is missing or damaged" 1>&2
fi
if test "x$v" != x
then
: # use $v
# Otherwise, if there is at least one git commit involving the working
# directory, and "git describe" output looks sensible, use that to
# derive a version string.
elif test "`git log -1 --pretty=format:x . 2>&1`" = x \
&& v=`git describe --abbrev=4 --match="$prefix*" HEAD 2>/dev/null \
|| git describe --abbrev=4 HEAD 2>/dev/null` \
&& v=`printf '%s\n' "$v" | sed "$tag_sed_script"` \
&& case $v in
$prefix[0-9]*) ;;
*) (exit 1) ;;
esac
then
# Is this a new git that lists number of commits since the last
# tag or the previous older version that did not?
# Newer: v6.10-77-g0f8faeb
# Older: v6.10-g0f8faeb
case $v in
*-*-*) : git describe is okay three part flavor ;;
*-*)
: git describe is older two part flavor
# Recreate the number of commits and rewrite such that the
# result is the same as if we were using the newer version
# of git describe.
vtag=`echo "$v" | sed 's/-.*//'`
commit_list=`git rev-list "$vtag"..HEAD 2>/dev/null` \
|| { commit_list=failed;
echo "$0: WARNING: git rev-list failed" 1>&2; }
numcommits=`echo "$commit_list" | wc -l`
v=`echo "$v" | sed "s/\(.*\)-\(.*\)/\1-$numcommits-\2/"`;
test "$commit_list" = failed && v=UNKNOWN
;;
esac
# Change the first '-' to a '.', so version-comparing tools work properly.
# Remove the "g" in git describe's output string, to save a byte.
v=`echo "$v" | sed 's/-/./;s/\(.*\)-g/\1-/'`;
v_from_git=1
elif test "x$fallback" = x || git --version >/dev/null 2>&1; then
v=UNKNOWN
else
v=$fallback
fi
v=`echo "$v" |sed "s/^$prefix//"`
# Test whether to append the "-dirty" suffix only if the version
# string we're using came from git. I.e., skip the test if it's "UNKNOWN"
# or if it came from .tarball-version.
if test "x$v_from_git" != x; then
# Don't declare a version "dirty" merely because a time stamp has changed.
git update-index --refresh > /dev/null 2>&1
dirty=`exec 2>/dev/null;git diff-index --name-only HEAD` || dirty=
case "$dirty" in
'') ;;
*) # Append the suffix only if there isn't one already.
case $v in
*-dirty) ;;
*) v="$v-dirty" ;;
esac ;;
esac
fi
# Omit the trailing newline, so that m4_esyscmd can use the result directly.
echo "$v" | tr -d "$nl"
# Local variables:
# eval: (add-hook 'write-file-hooks 'time-stamp)
# time-stamp-start: "scriptversion="
# time-stamp-format: "%:y-%02m-%02d.%02H"
# time-stamp-time-zone: "UTC"
# time-stamp-end: "; # UTC"
# End:
vcftools-0.1.16/configure.ac 0000664 0000000 0000000 00000004051 13331057401 0015721 0 ustar 00root root 0000000 0000000 # -*- Autoconf -*-
# Process this file with autoconf to produce a configure script.
AC_PREREQ([2.63])
AC_INIT([vcftools],
[m4_esyscmd([build-aux/git-version-gen .tarball-version])],
[https://github.com/vcftools/vcftools/issues])
AC_CONFIG_SRCDIR([src/cpp/vcftools.cpp])
AC_CONFIG_HEADERS([config.h])
# Automake invocation.
AM_INIT_AUTOMAKE([foreign])
# Checks for programs.
AC_PROG_CXX
AC_PROG_CC
AC_PROG_CPP
# Checks for perl.
AC_PATH_PROGS([PERL], [perl] [perl5], [false])
AS_IF([test "x$PERL" = "xfalse"],[
AC_MSG_ERROR([Perl not found; check your \$PATH.])
])
pmdir_relative_path=`\
$PERL -MConfig \
-wle '($_ = $Config{installsitelib})
=~ s!^\Q$Config{siteprefix}/!!; \
print'`
AC_ARG_WITH(
[pmdir],
AS_HELP_STRING(
[--with-pmdir=DIR],
[install Perl modules in DIR]),
[PMDIR=${withval}],
[PMDIR="$pmdir_relative_path"])
AC_SUBST([PMDIR])
# Checks for libraries.
PKG_CHECK_MODULES([ZLIB], [zlib])
# Checks for header files.
AC_CHECK_HEADERS([arpa/inet.h fcntl.h limits.h netdb.h stdint.h stdlib.h string.h sys/socket.h unistd.h])
# Checks for typedefs, structures, and compiler characteristics.
AC_HEADER_STDBOOL
AC_C_INLINE
AC_TYPE_INT16_T
AC_TYPE_INT32_T
AC_TYPE_INT64_T
AC_TYPE_INT8_T
AC_TYPE_OFF_T
AC_TYPE_SIZE_T
AC_TYPE_SSIZE_T
AC_TYPE_UINT16_T
AC_TYPE_UINT32_T
AC_TYPE_UINT8_T
# Checks for operating system services or capabilities.
AC_SYS_LARGEFILE
# Checks for library functions.
AC_FUNC_ERROR_AT_LINE
AC_FUNC_MALLOC
AC_FUNC_REALLOC
AC_CHECK_FUNCS([gethostbyaddr gethostbyname memset pow select socket sqrt strchr strdup strerror strstr strtol])
# Optional features.
AC_ARG_ENABLE([pca],
AS_HELP_STRING([--enable-pca], [enable PCA feature]),
[pca=${enableval}],
[pca=no])
AS_IF([test "x$pca" = "xyes"],[
AC_CHECK_LIB([lapack], [dgeev_])
])
# Generate output.
AC_CONFIG_FILES([Makefile
src/Makefile
src/cpp/Makefile
src/perl/Makefile])
AC_OUTPUT
vcftools-0.1.16/examples/ 0000775 0000000 0000000 00000000000 13331057401 0015251 5 ustar 00root root 0000000 0000000 vcftools-0.1.16/examples/annotate-test.vcf 0000664 0000000 0000000 00000003225 13331057401 0020541 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.0
##INFO=
##FORMAT=
##FORMAT=
##FORMAT=
##FILTER=
#CHROM POS ID REF ALT QUAL FILTER INFO
1 100 . GTTT G 1806 q10 DP=35
1 104 . C . 1792 PASS DP=32
1 105 . C T 246 PASS DP=10
1 106 . C A 246 PASS DP=10
2 107 . C . 1806 q10 DP=35
2 108 . C . 1792 PASS DP=32
2 109 . C . 628 q10 DP=21
2 110 . C G 1016 PASS DP=22
2 111 . C G 727 PASS DP=30
2 112 . C G 246 PASS DP=10
2 113 . C . 246 PASS DP=10
2 114 . T . 246 PASS DP=10
2 115 . T . 246 PASS DP=10
2 116 . T . 246 PASS DP=10
2 117 . T A 246 PASS DP=10
2 118 . T C 246 PASS DP=10
2 119 . TAAA T 246 PASS DP=10
2 124 . TA T 246 PASS DP=10
2 128 . T TA 246 PASS DP=10
2 130 . C A 246 PASS DP=10
2 131 . T A 246 PASS DP=10
2 132 . T A 246 PASS DP=10
2 133 . T A 246 PASS DP=10
2 134 . T A 246 PASS DP=10
2 135 . T C 246 PASS DP=10
2 136 . TT T 246 PASS DP=10;AF=0.1
2 138 . TT T 246 PASS DP=10;AF=0.2
2 140 . TT T 246 PASS DP=10;AF=0.1
17 12412 . CAGAGAGAGA CAGAGAGAGAGA 74.8 . INDEL;DP=5388;AF1=0.006576;CI95=0.005525,0.01105;DP4=2077,2367,21,22;MQ=47;FQ=74.8;PV4=0.88,1,0.34,0.021
17 12427 . G A 999 . DP=5557;AF1=0.06028;CI95=0.04972,0.07182;DP4=2461,2689,106,74;MQ=47;FQ=999;PV4=0.0038,1,2.6e-12,1
17 69284 . G A 14.6 . DP=3946;AF1=0.003468;CI95=0.002762,0.008287;DP4=1529,2177,7,9;MQ=44;FQ=14.6;PV4=1,0.035,0.098,1
17 69293 . GTTTCATTTC GTTTCTTTTCATTTC 999 . INDEL;DP=3568;AF1=0.1295;CI95=0.1077,0.1547;DP4=1014,1238,118,121;MQ=44;FQ=999;PV4=0.22,1,9.4e-54,1
vcftools-0.1.16/examples/annotate.out 0000664 0000000 0000000 00000002540 13331057401 0017614 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.0
##FORMAT=
##FORMAT=
##FORMAT=
##FILTER=
##INFO=
##INFO=
##INFO=
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT A
1 100 . GTTT G 1806 q10 DP=5;GN=gene1;HM2 GT:GQ:DP 0/1:409:35
1 110 . C T,G 1792 PASS DP=6,6 GT:GQ:DP 0/1:245:32
1 110 . CAAA C 1792 PASS DP=6,6 GT:GQ:DP 0/1:245:32
1 120 . GA G 628 q10 DP=21 GT:GQ:DP 1/1:21:21
1 130 . G T 1016 PASS DP=7,7;HM2=, GT:GQ:DP 0/1:212:22
1 130 . GAA GG 1016 PASS DP=7,7;HM2=, GT:GQ:DP 0/1:212:22
1 140 . GT G 727 PASS DP=8 GT:GQ:DP 0/1:150:30
1 150 . TAAAA TA,T 246 PASS DP=9 GT:GQ:DP 1/2:12:10
1 160 . TAAAA TA,T 246 PASS DP=10 GT:GQ:DP 1/2:12:10
2 100 . GTTT G 1806 q10 DP=35 GT:GQ:DP 0/1:409:35
2 110 . CAAA C 1792 PASS GN=gene2;HM2 GT:GQ:DP 0/1:245:32
2 120 . GA G 628 q10 GN=gene2;HM2 GT:GQ:DP 1/1:21:21
2 130 . GAA G 1016 PASS GN=gene2;HM2 GT:GQ:DP 0/1:212:22
2 140 . GT G 727 PASS GN=gene2;HM2 GT:GQ:DP 0/1:150:30
2 150 . TAAAA TA,T 246 PASS GN=gene2;HM2 GT:GQ:DP 1/2:12:10
2 160 . TAAAA TA,TC,T 246 PASS DP=11;GN=gene3 GT:GQ:DP 0/2:12:10
vcftools-0.1.16/examples/annotate.txt 0000664 0000000 0000000 00000000453 13331057401 0017625 0 ustar 00root root 0000000 0000000 100 100 1 id1_100 . . HM2 gene1 5
110 110 1 id1_110 CAAA C,CA 0 . 6
110 110 1 id2_110 C T 0 . 6
130 130 1 id1_130 G T HM2 . 7
130 130 1 id2_130 GAA GG HM2 . 7
140 140 1 id1_140 GT G 0 . 8
150 150 1 id1_150 TAAAA T 0 . 9
110 150 2 id2_110_150 CAAA C HM2 gene2 .
160 160 2 id2_160 TAAAA TC 0 gene3 11
vcftools-0.1.16/examples/annotate2.out 0000664 0000000 0000000 00000006043 13331057401 0017700 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.0
##INFO=
##FORMAT=
##FORMAT=
##FORMAT=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
#CHROM POS ID REF ALT QUAL FILTER INFO
1 100 . GTTT G 1806 q10;MaxDP DP=35
1 104 . C . 1792 PASS DP=32
1 105 . C T 246 SnpGap;SnpCluster DP=10
1 106 . C A 246 SnpGap;SnpCluster DP=10
2 107 . C . 1806 q10;MaxDP DP=35
2 108 . C . 1792 PASS DP=32
2 109 . C . 628 q10 DP=21
2 110 . C G 1016 SnpGap;SnpCluster DP=22
2 111 . C G 727 SnpGap;SnpCluster DP=30
2 112 . C G 246 SnpGap;SnpCluster DP=10
2 113 . C . 246 PASS DP=10
2 114 . T . 246 PASS DP=10
2 115 . T . 246 PASS DP=10
2 116 . T . 246 PASS DP=10
2 117 . T A 246 SnpGap;SnpCluster DP=10
2 118 . T C 246 SnpGap;SnpCluster DP=10
2 119 . TAAA T 246 SnpCluster DP=10
2 124 . TA T 246 GapWin DP=10
2 128 . T TA 246 SnpCluster DP=10
2 130 . C A 246 SnpGap;SnpCluster DP=10
2 131 . T A 246 SnpGap;SnpCluster DP=10
2 132 . T A 246 SnpGap;SnpCluster DP=10
2 133 . T A 246 SnpGap;SnpCluster DP=10
2 134 . T A 246 SnpGap;SnpCluster DP=10
2 135 . T C 246 SnpGap;SnpCluster DP=10
2 136 . TT T 246 GapWin;SnpCluster DP=10;AF=0.1
2 138 . TT T 246 SnpCluster DP=10;AF=0.2
2 140 . TT T 246 GapWin;SnpCluster DP=10;AF=0.1
17 12412 . CAGAGAGAGA CAGAGAGAGAGA 74.8 MaxDP INDEL;DP=5388;AF1=0.006576;CI95=0.005525,0.01105;DP4=2077,2367,21,22;MQ=47;FQ=74.8;PV4=0.88,1,0.34,0.021
17 12427 . G A 999 MaxDP;SnpGap DP=5557;AF1=0.06028;CI95=0.04972,0.07182;DP4=2461,2689,106,74;MQ=47;FQ=999;PV4=0.0038,1,2.6e-12,1
17 69284 . G A 14.6 MaxDP;SnpGap DP=3946;AF1=0.003468;CI95=0.002762,0.008287;DP4=1529,2177,7,9;MQ=44;FQ=14.6;PV4=1,0.035,0.098,1
17 69293 . GTTTCATTTC GTTTCTTTTCATTTC 999 MaxDP INDEL;DP=3568;AF1=0.1295;CI95=0.1077,0.1547;DP4=1014,1238,118,121;MQ=44;FQ=999;PV4=0.22,1,9.4e-54,1
vcftools-0.1.16/examples/annotate3.out 0000664 0000000 0000000 00000002562 13331057401 0017703 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.0
##FORMAT=
##FORMAT=
##FORMAT=
##FILTER=
##INFO=
##INFO=
##INFO=
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT A
1 100 id1_100 GTTT G 1806 q10 DP=5;GN=gene1;HM2 GT:GQ:DP 0/1:409:35
1 110 id2_110 C T,G 1792 PASS DP=6 GT:GQ:DP 0/1:245:32
1 110 id1_110 CAAA C 1792 PASS DP=6 GT:GQ:DP 0/1:245:32
1 120 . GA G 628 q10 DP=21 GT:GQ:DP 1/1:21:21
1 130 id1_130 G T 1016 PASS DP=7;HM2 GT:GQ:DP 0/1:212:22
1 130 id2_130 GAA GG 1016 PASS DP=7;HM2 GT:GQ:DP 0/1:212:22
1 140 id1_140 GT G 727 PASS DP=8 GT:GQ:DP 0/1:150:30
1 150 id1_150 TAAAA TA,T 246 PASS DP=9 GT:GQ:DP 1/2:12:10
1 160 . TAAAA TA,T 246 PASS DP=10 GT:GQ:DP 1/2:12:10
2 100 . GTTT G 1806 q10 DP=35 GT:GQ:DP 0/1:409:35
2 110 id2_110_150 CAAA C 1792 PASS GN=gene2;HM2 GT:GQ:DP 0/1:245:32
2 120 . GA G 628 q10 DP=21 GT:GQ:DP 1/1:21:21
2 130 . GAA G 1016 PASS DP=22 GT:GQ:DP 0/1:212:22
2 140 . GT G 727 PASS DP=30 GT:GQ:DP 0/1:150:30
2 150 . TAAAA TA,T 246 PASS DP=10 GT:GQ:DP 1/2:12:10
2 160 id2_160 TAAAA TA,TC,T 246 PASS DP=11;GN=gene3 GT:GQ:DP 0/2:12:10
vcftools-0.1.16/examples/cmp-test-a-3.3.vcf 0000664 0000000 0000000 00000001012 13331057401 0020216 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv3.3
##INFO=DP,1,Integer,"Total Depth"
##FORMAT=GT,1,String,"Genotype"
##FORMAT=GQ,1,Integer,"Genotype Quality"
##FORMAT=DP,1,Integer,"Read Depth"
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT A B
1 100100 . G C 0 . DP=1 GT:GQ:DP 0|1:40:1 0/1:40:1
1 100200 . G C 0 . DP=1 GT:GQ:DP 0|1:40:1 0/0:40:1
1 100300 . G C 0 . DP=1 GT:GQ:DP 1/1:40:1 ./.:40:1
1 100400 . C G,T 35 . DP=1 GT:GQ:DP 1/1:41:1 0/2:40:1
1 100500 . A G 0 . DP=1 GT:GQ:DP 1/1:40:1 0/0:40:1
1 100600 . C G 0 . DP=1 GT:GQ:DP 1/1:40:1 0/0:40:1
vcftools-0.1.16/examples/cmp-test-a.vcf 0000664 0000000 0000000 00000001176 13331057401 0017730 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.0
##INFO=
##FORMAT=
##FORMAT=
##FORMAT=
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT A B
1 100100 . G C 0 . DP=1 GT:GQ:DP 0|1:40:1 0/1:40:1
1 100200 . G C 0 . DP=1 GT:GQ:DP 0|1:40:1 0/0:40:1
1 100300 . G C 0 . DP=1 GT:GQ:DP 1/1:40:1 ./.:40:1
1 100400 . C G,T 35 . DP=1 GT:GQ:DP 1/1:41:1 0/2:40:1
1 100500 . A G 0 . DP=1 GT:GQ:DP 1/1:40:1 0/0:40:1
1 100600 . C G 0 . DP=1 GT:GQ:DP 1/1:40:1 0/0:40:1
vcftools-0.1.16/examples/cmp-test-b-3.3.vcf 0000664 0000000 0000000 00000001010 13331057401 0020215 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv3.3
##INFO=DP,1,Integer,"Total Depth"
##FORMAT=GT,1,String,"Genotype"
##FORMAT=GQ,1,Integer,"Genotype Quality"
##FORMAT=DP,1,Integer,"Read Depth"
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT A B
1 100100 . G C 0 . DP=1 GT:GQ:DP 0|1:40:1 1/0:40:1
1 100200 . G C 0 . DP=1 GT:GQ:DP 1|0:40:1 0/0:40:1
1 100300 . G C 0 . DP=1 GT:GQ:DP 1/1:40:1 0/0:40:1
1 100400 . C G 35 . DP=1 GT:GQ:DP 1/1:41:1 0/0:40:1
1 100500 . A G 0 . DP=1 GT:GQ:DP 1/1:40:1 0/0:40:1
1 100600 . C G 0 . DP=1 GT:GQ:DP 1/1:40:1 0/0:40:1
vcftools-0.1.16/examples/cmp-test-b.vcf 0000664 0000000 0000000 00000001174 13331057401 0017727 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.0
##INFO=
##FORMAT=
##FORMAT=
##FORMAT=
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT A B
1 100100 . G C 0 . DP=1 GT:GQ:DP 0|1:40:1 1/0:40:1
1 100200 . G C 0 . DP=1 GT:GQ:DP 1|0:40:1 0/0:40:1
1 100300 . G C 0 . DP=1 GT:GQ:DP 1/1:40:1 0/0:40:1
1 100400 . C G 35 . DP=1 GT:GQ:DP 1/1:41:1 0/0:40:1
1 100500 . A G 0 . DP=1 GT:GQ:DP 1/1:40:1 0/0:40:1
1 100600 . C G 0 . DP=1 GT:GQ:DP 1/1:40:1 0/0:40:1
vcftools-0.1.16/examples/cmp-test.out 0000664 0000000 0000000 00000005111 13331057401 0017534 0 ustar 00root root 0000000 0000000 # This file was generated by vcf-compare.
#
#VN 'Venn-Diagram Numbers'. Use `grep ^VN | cut -f 2-` to extract this part.
#VN The columns are:
#VN 1 .. number of sites unique to this particular combination of files
#VN 2- .. combination of files and space-separated number, a fraction of sites in the file
VN 6 cmp-test-a.vcf.gz (100.0%) cmp-test-b.vcf.gz (100.0%)
#SN Summary Numbers. Use `grep ^SN | cut -f 2-` to extract this part.
SN Number of REF matches: 6
SN Number of ALT matches: 5
SN Number of REF mismatches: 0
SN Number of ALT mismatches: 1
SN Number of samples in GT comparison: 2
#GS Genotype Comparison Summary. Use `grep ^GS | cut -f 2-` to extract this part.
#GS The columns are:
#GS 1 .. variant type
#GS 2 .. number of mismatches
#GS 3 .. number of matches
#GS 4 .. discordance
GS hom_RR 0 3 0.00%
GS het_RA 1 3 25.00%
GS hom_AA 0 4 0.00%
GS het_AA 0 0 0.00%
SN Non-reference Discordance Rate (NDR): 12.50
SN Summary: NDR 12.50, RR 0.00, RA 25.00, AA 0.00
#GC Genotype Comparison. Use `grep ^GC | cut -f 2-` to extract this part.
#GC The columns are:
#GC 1 .. Sample
#GC 2-6 .. Gtype mismatches: total hom_RR hom_AA het_RA het_AA
#GC 7-9 .. Gtype lost: total het_RA het_AA
#GC 10-14 .. Gtype gained: total hom_RR hom_AA het_RA het_AA
#GC 15-17 .. Phase lost: total het_RA het_AA
#GC 18 .. Phase gained
#GC 19-23 .. Matching sites: total hom_RR hom_AA het_RA het_AA
#GC 24 .. Phased matches: het_RA
#GC 25 .. Misphased matches: het_RA
GC A 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 6 0 4 2 0 2 1
GC B 1 0 0 1 0 0 0 0 1 1 0 0 0 0 0 0 0 4 3 0 1 0 0 0
#AF Number of matching and mismatching genotypes vs non-ref allele frequency. Use `^AF | cut -f 2-` to extract this part.
#AF The columns are:
#AF 1 .. Non-ref allele count
#AF 2 .. Hom(RR) matches
#AF 3 .. Het(RA) matches
#AF 4 .. Hom(AA) matches
#AF 5 .. Het(AA) matches
#AF 6 .. Hom(RR) mismatches
#AF 7 .. Het(RA) mismatches
#AF 8 .. Hom(AA) mismatches
#AF 9 .. Het(AA) mismatches
AF 0.25 1 1 0 0 0 0 0 0
AF 0.50 2 2 2 0 0 0 0 0
AF 0.75 0 0 1 0 0 1 0 0
AF 1.00 0 0 1 0 0 0 0 0
#DP Counts by depth. Use `grep ^DP | cut -f 2-` to extract this part.
#DP The columns are:
#DP 1 .. depth
#DP 2 .. RR matches
#DP 3 .. RA matches
#DP 4 .. AA matches
#DP 5 .. RR -> RA mismatches
#DP 6 .. RR -> AA mismatches
#DP 7 .. RA -> RR mismatches
#DP 8 .. RA -> AA mismatches
#DP 9 .. AA -> RR mismatches
#DP 10 .. AA -> RA mismatches
DP 1 3 3 4 0 0 1 0 0 0
vcftools-0.1.16/examples/concat-a.vcf 0000664 0000000 0000000 00000002234 13331057401 0017437 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.0
##INFO=
##FORMAT=
##FORMAT=
##FORMAT=
##FILTER=
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT A
1 100 . GTTT G 1806 q10 DP=35 GT:GQ:DP 0/1:409:35
1 110 . C T,G 1792 PASS DP=32 GT:GQ:DP 0/1:245:32
1 110 . CAAA C 1792 PASS DP=32 GT:GQ:DP 0/1:245:32
1 120 . GA G 628 q10 DP=21 GT:GQ:DP 1/1:21:21
1 130 . G T 1016 PASS DP=22 GT:GQ:DP 0/1:212:22
1 130 . GAA GG 1016 PASS DP=22 GT:GQ:DP 0/1:212:22
1 140 . GT G 727 PASS DP=30 GT:GQ:DP 0/1:150:30
1 150 . TAAAA TA,T 246 PASS DP=10 GT:GQ:DP 1/2:12:10
1 160 . TAAAA TA,T 246 PASS DP=10 GT:GQ:DP 1/2:12:10
2 100 . GTTT G 1806 q10 DP=35 GT:GQ:DP 0/1:409:35
2 110 . CAAA C 1792 PASS DP=32 GT:GQ:DP 0/1:245:32
2 120 . GA G 628 q10 DP=21 GT:GQ:DP 1/1:21:21
2 130 . GAA G 1016 PASS DP=22 GT:GQ:DP 0/1:212:22
2 140 . GT G 727 PASS DP=30 GT:GQ:DP 0/1:150:30
2 150 . TAAAA TA,T 246 PASS DP=10 GT:GQ:DP 1/2:12:10
2 160 . TAAAA TA,TC,T 246 PASS DP=10 GT:GQ:DP 0/2:12:10
vcftools-0.1.16/examples/concat-b.vcf 0000664 0000000 0000000 00000001242 13331057401 0017436 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.0
##INFO=
##FORMAT=
##FORMAT=
##FORMAT=
##FILTER=
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT A
1 141 . GTTT G 1806 q10 DP=35 GT:GQ:DP 0/1:409:35
1 151 . CAAA C 1792 PASS DP=32 GT:GQ:DP 0/1:245:32
1 161 . GA G 628 q10 DP=21 GT:GQ:DP 1/1:21:21
1 171 . GAA G 1016 PASS DP=22 GT:GQ:DP 0/1:212:22
1 181 . GT G 727 PASS DP=30 GT:GQ:DP 0/1:150:30
1 191 . TAAAA TA,T 246 PASS DP=10 GT:GQ:DP 1/2:12:10
vcftools-0.1.16/examples/concat-c.vcf 0000664 0000000 0000000 00000001714 13331057401 0017443 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.0
##INFO=
##FORMAT=
##FORMAT=
##FORMAT=
##FILTER=
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT A
2 142 . GTTT G 1806 q10 DP=35 GT:GQ:DP 0/1:409:35
2 152 . CAAA C 1792 PASS DP=32 GT:GQ:DP 0/1:245:32
2 162 . GA G 628 q10 DP=21 GT:GQ:DP 1/1:21:21
2 172 . GAA G 1016 PASS DP=22 GT:GQ:DP 0/1:212:22
2 182 . GT G 727 PASS DP=30 GT:GQ:DP 0/1:150:30
2 192 . TAAAA TA,T 246 PASS DP=10 GT:GQ:DP 1/2:12:10
1 142 . GTTT G 1806 q10 DP=35 GT:GQ:DP 0/1:409:35
1 152 . CAAA C 1792 PASS DP=32 GT:GQ:DP 0/1:245:32
1 162 . GA G 628 q10 DP=21 GT:GQ:DP 1/1:21:21
1 172 . GAA G 1016 PASS DP=22 GT:GQ:DP 0/1:212:22
1 182 . GT G 727 PASS DP=30 GT:GQ:DP 0/1:150:30
1 192 . TAAAA TA,T 246 PASS DP=10 GT:GQ:DP 1/2:12:10
vcftools-0.1.16/examples/concat.out 0000664 0000000 0000000 00000004032 13331057401 0017250 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.0
##INFO=
##FORMAT=
##FORMAT=
##FORMAT=
##FILTER=
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT A
1 100 . GTTT G 1806 q10 DP=35 GT:GQ:DP 0/1:409:35
1 110 . C T,G 1792 PASS DP=32 GT:GQ:DP 0/1:245:32
1 110 . CAAA C 1792 PASS DP=32 GT:GQ:DP 0/1:245:32
1 120 . GA G 628 q10 DP=21 GT:GQ:DP 1/1:21:21
1 130 . G T 1016 PASS DP=22 GT:GQ:DP 0/1:212:22
1 130 . GAA GG 1016 PASS DP=22 GT:GQ:DP 0/1:212:22
1 140 . GT G 727 PASS DP=30 GT:GQ:DP 0/1:150:30
1 141 . GTTT G 1806 q10 DP=35 GT:GQ:DP 0/1:409:35
1 142 . GTTT G 1806 q10 DP=35 GT:GQ:DP 0/1:409:35
1 150 . TAAAA TA,T 246 PASS DP=10 GT:GQ:DP 1/2:12:10
1 151 . CAAA C 1792 PASS DP=32 GT:GQ:DP 0/1:245:32
1 152 . CAAA C 1792 PASS DP=32 GT:GQ:DP 0/1:245:32
1 160 . TAAAA TA,T 246 PASS DP=10 GT:GQ:DP 1/2:12:10
1 161 . GA G 628 q10 DP=21 GT:GQ:DP 1/1:21:21
1 162 . GA G 628 q10 DP=21 GT:GQ:DP 1/1:21:21
1 171 . GAA G 1016 PASS DP=22 GT:GQ:DP 0/1:212:22
1 172 . GAA G 1016 PASS DP=22 GT:GQ:DP 0/1:212:22
1 181 . GT G 727 PASS DP=30 GT:GQ:DP 0/1:150:30
1 182 . GT G 727 PASS DP=30 GT:GQ:DP 0/1:150:30
1 191 . TAAAA TA,T 246 PASS DP=10 GT:GQ:DP 1/2:12:10
1 192 . TAAAA TA,T 246 PASS DP=10 GT:GQ:DP 1/2:12:10
2 100 . GTTT G 1806 q10 DP=35 GT:GQ:DP 0/1:409:35
2 110 . CAAA C 1792 PASS DP=32 GT:GQ:DP 0/1:245:32
2 120 . GA G 628 q10 DP=21 GT:GQ:DP 1/1:21:21
2 130 . GAA G 1016 PASS DP=22 GT:GQ:DP 0/1:212:22
2 140 . GT G 727 PASS DP=30 GT:GQ:DP 0/1:150:30
2 142 . GTTT G 1806 q10 DP=35 GT:GQ:DP 0/1:409:35
2 150 . TAAAA TA,T 246 PASS DP=10 GT:GQ:DP 1/2:12:10
2 152 . CAAA C 1792 PASS DP=32 GT:GQ:DP 0/1:245:32
2 160 . TAAAA TA,TC,T 246 PASS DP=10 GT:GQ:DP 0/2:12:10
2 162 . GA G 628 q10 DP=21 GT:GQ:DP 1/1:21:21
2 172 . GAA G 1016 PASS DP=22 GT:GQ:DP 0/1:212:22
2 182 . GT G 727 PASS DP=30 GT:GQ:DP 0/1:150:30
2 192 . TAAAA TA,T 246 PASS DP=10 GT:GQ:DP 1/2:12:10
vcftools-0.1.16/examples/consensus.fa 0000664 0000000 0000000 00000002016 13331057401 0017600 0 ustar 00root root 0000000 0000000 >1:1-500
ATACCATATGTGACTTATAAAAAAGAACATAACCTACGTATCAACTAAAGTGGTTGTTTG
CAGAAAAGGAAGACTTAAAAAGAGTCAGTACTAACCTACATAATATATACAATGTTCATT
AAATAATAAAATGAGCTCATCATACTTAGGTCATCATAAATATATCTGAAATTCACAAAT
ATTGATCAAATGGTAAAATAGACAAGTAGATTTTAATAGGTTAAACAATTACTGATTCTC
TTGAAAGAATAAATTTAATATGAGACCTATTTCATTATAATGAACTCACAAATTAGAAAC
TTCACACTGGGGGCTGGAGAGATGGCTCAGTAGTTAAGAACACTGACTGCTCTTCTGAAG
GTCCTGAGTTCAAATCCCAGCAACCACATGGTGACTTACAACCATCTGTAATGACATCTG
ATGCCCTCTGGTGTGTCTGAAGACAGCTACAGTGTACTTACATAAAATAATAAATAAATC
TTTAAAAACAAAAAAAAAGAA
>2:1-500
GAAGATCTTTTCCTTATTAAGGATCTGAAGCTCTGTAGATTTGTATTCTATTAAACATGG
AGAGATTAGTGATTTTCCATATTCTTTAAGTCATTTTAGAGTAATGTGTTCTTAAGATAA
ATCAGAAAAACAAAAACTTGTGCTTTCCTGTTTGAAAAACAAACAGCTGTGGGGAATGGT
GTCGGGACAGCCTTTTTATAAAATTTTTCTAAATAATGTTGAGGCTTTGATACGTCAAAG
TTATATTTCAAATGGAATCACTTAGACCTCGTTTCTGAGTGTCAATGGCCATATTGGGGA
TTTGCTGCTGCCAATGACAGCACACCCTGGGAATGCCCCAACTACTTACTACAAAGCAGT
GTTACATGGAGAAGATCTTCAAGAGTCTTTTTGCTAGATCTTTCCTTGGCTTTTGATGTG
ACTCCTCTCAATAAAATCCACAGTAATATAGTGAGTGGTCTCCTGCTCCAAACCAGTATT
TCAGACACAGTTAATCCAGAC
vcftools-0.1.16/examples/consensus.out 0000664 0000000 0000000 00000002015 13331057401 0020020 0 ustar 00root root 0000000 0000000 >1:1-500
ATAC*ATAT*TG*T***ATAAAAAAGAACATAACCTACGTATCAACTAAAGTGGTTGTTT
G*AGAAAAGGAAGACTTAAAAAGAGTCAGTACTAACCTACATAATATATACAATGTTCAT
TAAATAATAAAATGAGCTCATCATACTTAGGTCATCATAAATATATCTGAAATTCACAAA
TATTGATCAAATGGTAAAATAGACAAGTAGATTTTAATAGGTTAAACAATTACTGATTCT
CTTGAAAGAATAAATTTAATATGAGACCTATTTCATTATAATGAACTCACAAATTAGAAA
CTTCACACTGGGGGCTGGAGAGATGGCTCAGTAGTTAAGAACACTGACTGCTCTTCTGAA
GGTCCTGAGTTCAAATCCCAGCAACCACATGGTGACTTACAACCATCTGTAATGACATCT
GATGCCCTCTGGTGTGTCTGAAGACAGCTACAGTGTACTTACATAAAATAATAAATAAAT
CTTTAAAAACAAAAAAAAAGAA
>2:1-500
GAAGATCTTTTCCTTATTAAGGATCTGAAGCTCTGTAGATTTGTATTCTATTAAACATGG
A*ATTAGTGATTTTCCATATTCTTTAAGTCATTTTAGAGTAATGTGTTCTTAAGATAAAT
CAGAAAAACAAAAACTTGTGCTTTCCTGTTTGAAAAACAAACAGCTGTGGGGAATGGTGT
CGGGACAGCCTTTTTATAAAATTTTTCTAAATAATGTTGAGGCTTTGATACGTCAAAGTT
ATATTTCAAATGGAATCACTTAGACCTCGTTTCTGAGTGTCAATGGCCATATTGGGGATT
TGCTGCTGCCAATGACAGCACACCCTGGGAATGCCCCAACTACTTACTACAAAGCAGTGT
TACATGGAGAAGATCTTCAAGAGTCTTTTTGCTAGATCTTTCCTTGGCTTTTGATGTGAC
TCCTCTCAATAAAATCCACAGTAATATAGTGAGTGGTCTCCTGCTCCAAACCAGTATT*C
AGACACAGTTAATCCAGAC
vcftools-0.1.16/examples/consensus.out2 0000664 0000000 0000000 00000002015 13331057401 0020102 0 ustar 00root root 0000000 0000000 >1:1-500
ATAC*ATATGTG*T***ATAAAAAAGAACATAACCTACGTATCAACTAAAGTGGTTGTTT
G*AGAAAAGGAAGACTTAAAAAGAGTCAGTACTAACCTACATAATATATACAATGTTCAT
TAAATAATAAAATGAGCTCATCATACTTAGGTCATCATAAATATATCTGAAATTCACAAA
TATTGATCAAATGGTAAAATAGACAAGTAGATTTTAATAGGTTAAACAATTACTGATTCT
CTTGAAAGAATAAATTTAATATGAGACCTATTTCATTATAATGAACTCACAAATTAGAAA
CTTCACACTGGGGGCTGGAGAGATGGCTCAGTAGTTAAGAACACTGACTGCTCTTCTGAA
GGTCCTGAGTTCAAATCCCAGCAACCACATGGTGACTTACAACCATCTGTAATGACATCT
GATGCCCTCTGGTGTGTCTGAAGACAGCTACAGTGTACTTACATAAAATAATAAATAAAT
CTTTAAAAACAAAAAAAAAGAA
>2:1-500
GAAGATCTTTTCCTTATTAAGGATCTGAAGCTCTGTAGATTTGTATTCTATTAAACATGG
A*ATTAGTGATTTTCCATATTCTTTAAGTCATTTTAGAGTAATGTGTTCTTAAGATAAAT
CAGAAAAACAAAAACTTGTGCTTTCCTGTTTGAAAAACAAACAGCTGTGGGGAATGGTGT
CGGGACAGCCTTTTTATAAAATTTTTCTAAATAATGTTGAGGCTTTGATACGTCAAAGTT
ATATTTCAAATGGAATCACTTAGACCTCGTTTCTGAGTGTCAATGGCCATATTGGGGATT
TGCTGCTGCCAATGACAGCACACCCTGGGAATGCCCCAACTACTTACTACAAAGCAGTGT
TACATGGAGAAGATCTTCAAGAGTCTTTTTGCTAGATCTTTCCTTGGCTTTTGATGTGAC
TCCTCTCAATAAAATCCACAGTAATATAGTGAGTGGTCTCCTGCTCCAAACCAGTATT@C
AGACACAGTTAATCCAGAC
vcftools-0.1.16/examples/consensus.vcf 0000664 0000000 0000000 00000000423 13331057401 0017770 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.1
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT NA001
1 5 . C * . PASS . GT 0/1
1 10 . G * . PASS . GT 0/0
1 12 . GACT G* . PASS . GT 0/1
1 16 . T T*** . PASS . GT 1/1
1 61 . C * . PASS . GT 1/1
2 61 . AGAG A* . PASS . GT 0/1
2 481 . T *,@ . PASS . GT 0/2
vcftools-0.1.16/examples/contrast.out 0000664 0000000 0000000 00000015447 13331057401 0017652 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.1
##samtoolsVersion=0.1.18-r572
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##FORMAT=
##FORMAT=
##FORMAT=
##FORMAT=
##FORMAT=
##FORMAT=
##FORMAT=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##INFO=
##INFO=
##FILTER=
##INFO=
##INFO=
##INFO=
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT A B C D
1 10250 . A C 61 MinMQ DP=271;VDB=0.0265;AF1=0.125;AC1=1;DP4=87,78,18,9;MQ=17;FQ=61;PV4=0.21,1,1,0.1;AN=8;AC=1;NOVELAL=D;NOVELTY=255 GT:DP:SP:GQ 0/0:60:0:99 0/0:32:5:53 0/0:50:2:83 0/1:50:3:62
1 10352 . TACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCC TAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCC 253 MinMQ INDEL;DP=413;VDB=0.0226;AF1=0.8598;AC1=7;DP4=7,17,13,44;MQ=15;FQ=4.35;PV4=0.58,1,1,0.0055;AN=8;AC=7;NOVELAL=D;NOVELTY=6 GT:PL:DP:SP:GQ 1/1:67,6,0:18:2:11 1/1:14,7,0:12:0:12 1/1:111,22,0:23:0:26 0/1:83,0,22:28:2:18
1 17538 . C A 64 MinMQ DP=393;VDB=0.0314;AF1=0.125;AC1=1;DP4=138,205,17,27;MQ=28;FQ=64;PV4=0.87,0.32,1,1;AN=8;AC=1;NOVELAL=D;NOVELTY=29 GT:PL:DP:SP:GQ 0/0:0,152,255:148:1:99 0/0:0,29,227:72:4:34 0/0:0,71,255:86:6:76 0/1:70,0,226:81:4:65
1 28563 . A G 999 MinMQ DP=124;VDB=0.0072;AF1=1;AC1=8;DP4=22,31,27,39;MQ=18;FQ=-3.67;PV4=1,1,1,1;AN=8;AC=8;NOVELAL=D;NOVELTY=1 GT:PL:DP:SP:GQ 1/1:191,6,0:41:1:14 1/1:90,2,0:24:0:11 1/1:213,20,0:31:4:28 1/1:104,0,1:23:0:8
1 28590 . TT TTGGT 116 MinMQ INDEL;DP=112;VDB=0.0233;AF1=0.3933;AC1=3;DP4=5,46,10,16;MQ=19;FQ=54.6;PV4=0.005,1,1,0.00097;AN=8;AC=3;NOVELTY=9;NOVELGT=D GT:PL:DP:SP:GQ 0/1:80,0,2:23:10:8 0/1:9,0,9:15:15:9 0/1:51,0,26:21:2:31 0/0:0,17,39:18:5:16
1 55085 . T A 149 MinMQ DP=190;VDB=0.0199;AF1=0.3891;AC1=3;DP4=73,61,13,39;MQ=25;FQ=149;PV4=0.0003,0.35,0.01,1;AN=8;AC=3;NOVELTY=7;NOVELGT=D GT:PL:DP:SP:GQ 0/1:79,0,161:48:4:80 0/1:9,0,146:22:13:10 0/1:68,0,250:49:12:69 0/0:0,7,228:67:12:7
1 58176 . G A 94.7 MinMQ DP=93;VDB=0.0330;AF1=0.3746;AC1=3;DP4=51,13,15,9;MQ=17;FQ=94.7;PV4=0.11,0.0027,1,1;AN=8;AC=3;NOVELTY=18;NOVELGT=D GT:PL:DP:SP:GQ 0/1:30,0,23:22:0:26 0/1:18,0,15:12:7:17 0/1:55,0,102:29:2:56 0/0:0,42,114:25:9:41
1 66507 . T A 999 PASS DP=202;VDB=0.0385;AF1=0.626;AC1=5;DP4=25,14,63,82;MQ=42;FQ=999;PV4=0.03,0.023,1,0.0014;AN=8;AC=5;NOVELTY=20;NOVELGT=D GT:PL:DP:SP:GQ 0/1:255,0,205:42:7:99 0/1:255,0,20:37:12:21 0/1:255,0,155:57:4:99 1/1:255,72,0:48:0:71
1 66521 . TATATAATATA TATATAATATAATATA 999 PASS INDEL;DP=200;VDB=0.0384;AF1=0.3747;AC1=3;DP4=61,75,25,12;MQ=43;FQ=999;PV4=0.016,1,3.8e-20,0.38;AN=8;AC=3;NOVELTY=25;NOVELGT=D GT:PL:DP:SP:GQ 0/1:233,0,255:40:7:99 0/1:25,0,255:32:16:26 0/1:178,0,255:56:3:99 0/0:0,75,255:45:3:74
1 73841 . C T 999 PASS DP=182;VDB=0.0366;AF1=0.3748;AC1=3;DP4=50,64,12,26;MQ=30;FQ=999;PV4=0.25,1.6e-10,0.084,1;AN=8;AC=3;NOVELTY=28;NOVELGT=D GT:PL:DP:SP:GQ 0/1:95,0,255:33:3:96 0/1:174,0,204:27:9:99 0/1:28,0,255:53:17:29 0/0:0,64,255:39:6:63
X 2 . A G 89 PASS DP=304;VDB=0.0327;AF1=0.125;AC1=1;DP4=99,171,10,15;MQ=30;FQ=89;PV4=0.83,0.3,0.035,1;AN=8;AC=1;NOVELTY=255;NOVELGT=D GT:PL:DP 0/1:95,0,255:11 0/0:0,72,255:11 0/0:0,101,255:11 1:255,0:11
X 3 . A G 89 PASS DP=304;VDB=0.0327;AF1=0.125;AC1=1;DP4=99,171,10,15;MQ=30;FQ=89;PV4=0.83,0.3,0.035,1;AN=8;AC=1;NOVELAL=D;NOVELTY=255 GT:PL:DP 0/0:0,95,255:11 0/0:0,72,255:11 0/0:0,101,255:11 1:255,0:11
vcftools-0.1.16/examples/contrast.vcf 0000664 0000000 0000000 00000066641 13331057401 0017623 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.1
##samtoolsVersion=0.1.18-r572
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##FORMAT=
##FORMAT=
##FORMAT=
##FORMAT=
##FORMAT=
##FORMAT=
##FORMAT=
##source_20120424.1=vcf-annotate(r735) --fill-AC-AN -f +
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##INFO=
##INFO=
##source_20120424.2=vcf-annotate(r735) --fill-AC-AN -f +
##FILTER=
##source_20120710.1=vcf-annotate(r761) -f q=30 mpileup-v1/merged.filt.vcf.gz
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT A B C D
1 10177 . A C 37 MinMQ DP=495;VDB=0.0168;AF1=0.1596;AC1=1;DP4=167,82,52,21;MQ=13;FQ=37;PV4=0.57,1,1,1;AN=8;AC=1 GT:PL:DP:SP:GQ 0/0:0,0,15:101:2:5 0/0:0,36,89:51:5:40 0/0:0,79,103:85:1:83 0/1:41,0,31:85:16:36
1 10250 . A C 61 MinMQ DP=271;VDB=0.0265;AF1=0.125;AC1=1;DP4=87,78,18,9;MQ=17;FQ=61;PV4=0.21,1,1,0.1;AN=8;AC=1 GT:DP:SP:GQ 0/0:60:0:99 0/0:32:5:53 0/0:50:2:83 0/1:50:3:62
1 10257 . A C 31.9 MinMQ DP=400;VDB=0.0245;AF1=0.2404;AC1=2;DP4=93,100,26,10;MQ=16;FQ=31.9;PV4=0.01,1,1,0.013;AN=8;AC=2 GT:PL:DP:SP:GQ 0/0:0,93,197:65:3:95 0/1:13,0,92:41:9:11 0/0:0,91,128:59:0:93 0/1:27,0,70:64:14:25
1 10329 . ACCCC ACCC 26.4 MinMQ INDEL;DP=315;VDB=0.0160;AF1=0.2047;AC1=2;DP4=2,42,9,16;MQ=17;FQ=29.3;PV4=0.0011,1,0.061,1;AN=8;AC=1 GT:PL:DP:SP:GQ 0/0:0,4,68:22:10:7 0/0:2,0,16:7:0:3 0/0:0,15,61:18:0:17 0/1:46,0,34:22:7:40
1 10352 . TACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCC TAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCC 253 MinMQ INDEL;DP=413;VDB=0.0226;AF1=0.8598;AC1=7;DP4=7,17,13,44;MQ=15;FQ=4.35;PV4=0.58,1,1,0.0055;AN=8;AC=7 GT:PL:DP:SP:GQ 1/1:67,6,0:18:2:11 1/1:14,7,0:12:0:12 1/1:111,22,0:23:0:26 0/1:83,0,22:28:2:18
1 10492 . C T 999 PASS DP=213;VDB=0.0102;AF1=0.375;AC1=3;DP4=84,74,34,19;MQ=32;FQ=999;PV4=0.2,0.11,0.057,0.13;AN=8;AC=3 GT:PL:DP:SP:GQ 0/1:85,0,255:57:3:86 0/0:0,123,255:41:0:99 0/1:255,0,255:47:0:99 0/1:114,0,255:66:4:99
1 10583 . G A 20 PASS DP=134;VDB=0.0071;AF1=0.1242;AC1=1;DP4=78,41,6,6;MQ=32;FQ=20;PV4=0.35,0.29,0.052,1;AN=8;AC=1 GT:PL:DP:SP:GQ 0/1:26,0,227:40:8:21 0/0:0,35,255:21:0:40 0/0:0,108,255:36:0:99 0/0:0,21,255:34:0:26
1 10797 . CAGA CAGAGA 90.4 MinMQ INDEL;DP=37;VDB=0.0243;AF1=0.2819;AC1=2;DP4=7,3,0,6;MQ=29;FQ=93.3;PV4=0.011,1,1,1;AN=8;AC=2 GT:PL:DP:SP:GQ 0/0:0,9,104:3:0:10 0/0:0,6,39:2:0:8 0/1:59,0,36:5:0:43 0/1:56,0,32:6:4:39
1 10821 . T A 49.7 MinMQ DP=9;VDB=0.0091;AF1=1;AC1=5;DP4=1,3,0,4;MQ=12;FQ=7.75;PV4=1,1,1,1;AN=8;AC=6 GT:PL:DP:SP:GQ 1/1:42,9,0:3:0:10 0/1:0,3,4:1:0:2 1/1:12,1,0:2:0:4 0/1:0,6,8:2:0:2
1 14907 . A G 999 MinMQ DP=461;VDB=0.0384;AF1=0.5;G3=8.874e-45,1,8.011e-40;HWE=0.0185;AC1=4;DP4=101,122,129,102;MQ=25;FQ=999;PV4=0.031,0.011,1,1;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:225,0,255:133:0:99 0/1:213,0,225:91:20:99 0/1:255,0,188:104:0:99 0/1:255,0,208:126:4:99
1 14930 . A G 999 MinMQ DP=502;VDB=0.0393;AF1=0.5;G3=1.282e-48,1,7.866e-46;HWE=0.0185;AC1=4;DP4=117,121,135,111;MQ=28;FQ=999;PV4=0.24,0.02,0.42,1;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:255,0,255:150:0:99 0/1:232,0,255:84:9:99 0/1:255,0,250:114:0:99 0/1:255,0,218:136:4:99
1 15118 . A G 196 MinMQ DP=408;VDB=0.0389;AF1=0.4995;G3=4.894e-09,1,2.035e-09;HWE=0.0193;AC1=4;DP4=79,107,98,101;MQ=13;FQ=196;PV4=0.19,1,0.16,1;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:42,0,97:110:1:45 0/1:14,0,34:82:7:17 0/1:54,0,93:90:2:57 0/1:92,0,15:103:1:18
1 15211 . T G 999 MinMQ DP=381;VDB=0.0374;AF1=0.6993;AC1=6;DP4=52,44,122,137;MQ=17;FQ=156;PV4=0.28,0.31,1,1;AN=8;AC=6 GT:PL:DP:SP:GQ 0/1:146,0,101:114:5:99 1/1:77,1,0:67:2:4 0/1:121,0,61:78:7:60 1/1:192,89,0:96:11:90
1 15274 . A T 999 MinMQ DP=229;VDB=0.0313;AF1=1;AC1=8;DP4=0,0,99,120;MQ=11;FQ=-92.5;AN=8;AC=8 GT:PL:DP:SP:GQ 1/1:83,114,0:54:0:99 1/1:82,108,0:48:0:99 1/1:84,105,0:47:0:99 1/1:112,175,0:70:0:99
1 15820 . G T 90.6 MinMQ DP=149;VDB=0.0252;AF1=0.374;AC1=3;DP4=24,68,15,40;MQ=17;FQ=90.6;PV4=1,1,2e-07,1;AN=8;AC=3 GT:PL:DP:SP:GQ 0/1:40,0,124:48:1:41 0/0:0,27,153:33:2:26 1/1:65,25,0:15:3:20 0/0:0,65,195:51:5:64
1 15903 . GCC GCCC 158 MinMQ INDEL;DP=14;VDB=0.0182;AF1=1;AC1=8;DP4=0,0,0,7;MQ=29;FQ=-18.2;AN=8;AC=8 GT:PL:DP:SP:GQ 1/1:29,3,0:1:0:13 1/1:44,6,0:2:0:16 1/1:29,3,0:1:0:13 1/1:72,9,0:3:0:19
1 16103 . T G 24.3 PASS DP=110;VDB=0.0122;AF1=0.2119;AC1=2;DP4=49,26,23,2;MQ=31;FQ=24.3;PV4=0.01,1,6.9e-13,0.33;AN=8;AC=1 GT:PL:DP:SP:GQ 0/0:0,2,228:39:11:5 0/0:0,7,189:15:6:10 0/0:0,0,218:25:0:4 0/1:26,0,192:21:5:24
1 16378 . T C 999 MinMQ DP=587;VDB=0.0267;AF1=0.5;G3=9.954e-26,1,3.125e-18;HWE=0.0185;AC1=4;DP4=128,75,245,120;MQ=18;FQ=999;PV4=0.36,0.23,0.024,1;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:127,0,150:194:9:99 0/1:118,0,139:108:4:99 0/1:156,0,80:130:3:83 0/1:166,0,148:136:1:99
1 16495 . G C 97.7 MinMQ DP=644;VDB=0.0239;AF1=0.2493;AC1=2;DP4=226,252,67,87;MQ=19;FQ=97.7;PV4=0.46,0.14,1,1;AN=8;AC=2 GT:PL:DP:SP:GQ 0/0:0,41,126:190:5:43 0/0:0,17,168:115:1:19 0/1:19,0,132:166:1:17 0/1:87,0,176:161:4:85
1 16534 . C T 264 MinMQ DP=516;VDB=0.0397;AF1=0.5;G3=3.737e-14,1,2.067e-30;HWE=0.0185;AC1=4;DP4=129,149,109,113;MQ=14;FQ=264;PV4=0.59,0.34,1,0.0011;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:97,0,110:151:4:99 0/1:38,0,92:115:2:41 0/1:50,0,120:118:5:53 0/1:85,0,158:116:5:88
1 16571 . G A 120 MinMQ DP=435;VDB=0.0388;AF1=0.4998;G3=1.594e-10,1,7.561e-11;HWE=0.0189;AC1=4;DP4=94,134,84,109;MQ=10;FQ=120;PV4=0.69,0.018,1,1;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:33,0,20:123:0:23 0/1:42,0,24:95:2:27 0/1:18,0,43:107:0:21 0/1:33,0,70:96:4:36
1 17538 . C A 64 MinMQ DP=393;VDB=0.0314;AF1=0.125;AC1=1;DP4=138,205,17,27;MQ=28;FQ=64;PV4=0.87,0.32,1,1;AN=8;AC=1 GT:PL:DP:SP:GQ 0/0:0,152,255:148:1:99 0/0:0,29,227:72:4:34 0/0:0,71,255:86:6:76 0/1:70,0,226:81:4:65
1 20144 . G A 98.2 MinMQ DP=304;VDB=0.0356;AF1=0.4851;G3=4.916e-07,1,8.127e-35;HWE=0.0213;AC1=4;DP4=91,122,40,43;MQ=15;FQ=98.2;PV4=0.44,0.0094,1,1;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:6,0,72:94:6:9 0/1:32,0,80:44:13:35 0/1:28,0,112:81:9:31 0/1:38,0,62:77:0:41
1 28558 . C T 164 MinMQ DP=142;VDB=0.0026;AF1=0.4529;G3=1.465e-06,1,4.392e-30;HWE=0.0307;AC1=4;DP4=38,62,18,20;MQ=17;FQ=164;PV4=0.34,1,1,1;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:0,0,104:44:9:4 0/1:27,0,28:32:3:27 0/1:77,0,113:35:5:79 0/1:64,0,31:27:0:35
1 28563 . A G 999 MinMQ DP=124;VDB=0.0072;AF1=1;AC1=8;DP4=22,31,27,39;MQ=18;FQ=-3.67;PV4=1,1,1,1;AN=8;AC=8 GT:PL:DP:SP:GQ 1/1:191,6,0:41:1:14 1/1:90,2,0:24:0:11 1/1:213,20,0:31:4:28 1/1:104,0,1:23:0:8
1 28590 . TT TTGGT 116 MinMQ INDEL;DP=112;VDB=0.0233;AF1=0.3933;AC1=3;DP4=5,46,10,16;MQ=19;FQ=54.6;PV4=0.005,1,1,0.00097;AN=8;AC=3 GT:PL:DP:SP:GQ 0/1:80,0,2:23:10:8 0/1:9,0,9:15:15:9 0/1:51,0,26:21:2:31 0/0:0,17,39:18:5:16
1 30867 . CCTCTCTCTCTCTCTCTCTCTCTCTC CCTCTCTCTCTCTCTCTCTCTC 999 PASS INDEL;DP=229;VDB=0.0320;AF1=0.5;G3=4.953e-17,1,5e-52;HWE=0.0185;AC1=4;DP4=56,66,27,32;MQ=37;FQ=999;PV4=1,1,1,1;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:211,0,255:47:1:99 0/1:74,0,255:20:2:77 0/1:255,0,255:70:3:99 0/1:176,0,255:44:0:99
1 30923 . G T 999 PASS DP=107;VDB=0.0022;AF1=1;AC1=8;DP4=0,0,47,50;MQ=37;FQ=-36;AN=8;AC=8 GT:PL:DP:SP:GQ 1/1:255,75,0:25:0:99 1/1:221,30,0:10:0:72 1/1:255,117,0:39:0:99 1/1:255,69,0:23:0:99
1 40639 . CTTTTTTTTTTTTTTTTTTT CTTTTTTTTTTTTTTTT 118 PASS INDEL;DP=72;VDB=0.0379;AF1=1;AC1=8;DP4=0,0,14,2;MQ=33;FQ=-18.8;AN=8;AC=8 GT:PL:DP:SP:GQ 1/1:77,30,0:11:0:40 1/1:16,3,0:1:0:14 1/1:16,3,0:1:0:14 1/1:25,6,0:3:0:16
1 46633 . T A 46.4 MinMQ DP=169;VDB=0.0275;AF1=0.1322;AC1=1;DP4=67,81,9,9;MQ=15;FQ=46.4;PV4=0.8,0.5,1,1;AN=8;AC=1 GT:PL:DP:SP:GQ 0/1:52,0,60:47:0:47 0/0:0,7,71:30:0:12 0/0:0,114,146:38:0:99 0/0:0,154,179:51:0:99
1 49298 . T C 999 MinMQ DP=124;VDB=0.0376;AF1=1;AC1=8;DP4=17,14,49,36;MQ=24;FQ=-3.76;PV4=0.83,1,1,0.49;AN=8;AC=8 GT:PL:DP:SP:GQ 1/1:130,0,4:20:4:6 1/1:127,8,0:12:3:17 1/1:247,27,0:45:26:36 1/1:252,60,0:39:12:69
1 51803 . T C 999 MinMQ DP=88;VDB=0.0284;AF1=1;AC1=8;DP4=9,30,20,25;MQ=15;FQ=-3.64;PV4=0.065,0.21,1,1;AN=8;AC=8 GT:PL:DP:SP:GQ 1/1:45,2,0:22:11:10 1/1:60,5,0:6:3:13 1/1:105,1,0:30:2:9 1/1:153,6,0:26:2:14
1 51898 . C A 22.2 PASS DP=128;VDB=0.0230;AF1=0.1272;AC1=1;DP4=56,50,14,2;MQ=41;FQ=22.2;PV4=0.013,0.069,9.8e-17,1;AN=8;AC=1 GT:PL:DP:SP:GQ 0/1:28,0,255:23:13:23 0/0:0,35,255:19:3:40 0/0:0,59,255:44:0:64 0/0:0,11,255:36:9:16
1 51928 . G A 54.1 PASS DP=149;VDB=0.0311;AF1=0.1269;AC1=1;DP4=67,52,22,5;MQ=41;FQ=54.1;PV4=0.017,0.0073,7.3e-34,0.37;AN=8;AC=1 GT:PL:DP:SP:GQ 0/1:60,0,255:29:13:55 0/0:0,27,255:19:0:32 0/0:0,30,255:51:2:35 0/0:0,13,255:47:11:18
1 52058 . G C 17.5 PASS DP=132;VDB=0.0277;AF1=0.2308;AC1=2;DP4=55,57,15,2;MQ=35;FQ=17.5;PV4=0.0031,0.036,0.00091,0.03;AN=8;AC=2 GT:PL:DP:SP:GQ 0/1:13,0,178:20:7:11 0/0:0,60,255:20:0:62 0/1:12,0,255:51:8:10 0/0:0,15,255:38:9:17
1 52238 . T G 999 PASS DP=138;VDB=0.0125;AF1=1;AC1=8;DP4=0,0,65,60;MQ=37;FQ=-42;AN=8;AC=8 GT:PL:DP:SP:GQ 1/1:255,63,0:21:0:99 1/1:255,36,0:12:0:84 1/1:255,166,0:55:0:99 1/1:255,111,0:37:0:99
1 54586 . T C 51.1 PASS DP=116;VDB=0.0136;AF1=0.236;AC1=2;DP4=47,45,14,7;MQ=36;FQ=51.1;PV4=0.23,1,6.1e-11,1;AN=8;AC=2 GT:PL:DP:SP:GQ 0/1:48,0,186:24:0:46 0/0:0,26,252:14:0:28 0/1:11,0,255:32:7:9 0/0:0,40,255:43:9:42
1 54676 . C T 999 PASS DP=143;VDB=0.0244;AF1=0.4969;G3=1.224e-08,1,4.851e-96;HWE=0.0191;AC1=4;DP4=54,48,19,21;MQ=40;FQ=999;PV4=0.58,0.47,8.6e-13,0.2;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:82,0,233:23:7:85 0/1:121,0,237:20:2:99 0/1:173,0,255:49:6:99 0/1:13,0,255:50:2:16
1 54712 . TTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTT TTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTT 423 PASS INDEL;DP=161;VDB=0.0367;AF1=0.8839;AC1=7;DP4=1,0,4,12;MQ=45;FQ=-9.11;PV4=0.29,1,1,0.06;AN=8;AC=7 GT:PL:DP:SP:GQ 0/1:0,3,13:1:0:5 1/1:125,9,0:3:0:14 1/1:67,6,0:2:0:11 1/1:255,33,0:11:0:38
1 54753 . T G 61.5 PASS DP=177;VDB=0.0130;AF1=0.2019;AC1=2;DP4=48,82,1,5;MQ=40;FQ=61.5;PV4=0.42,1,0.27,1;AN=8;AC=2 GT:PL:DP:SP:GQ 0/1:5,0,237:20:3:4 0/1:63,0,193:16:0:60 0/0:0,160,255:53:0:99 0/0:0,141,255:47:0:99
1 54844 . G A 999 MinMQ DP=172;VDB=0.0254;AF1=0.4999;G3=4.104e-12,1,1.068e-33;HWE=0.0185;AC1=4;DP4=70,44,38,18;MQ=20;FQ=999;PV4=0.5,0.27,1,0.29;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:88,0,103:45:5:91 0/1:31,0,120:25:0:34 0/1:97,0,124:58:2:99 0/1:49,0,181:42:5:52
1 55085 . T A 149 MinMQ DP=190;VDB=0.0199;AF1=0.3891;AC1=3;DP4=73,61,13,39;MQ=25;FQ=149;PV4=0.0003,0.35,0.01,1;AN=8;AC=3 GT:PL:DP:SP:GQ 0/1:79,0,161:48:4:80 0/1:9,0,146:22:13:10 0/1:68,0,250:49:12:69 0/0:0,7,228:67:12:7
1 55164 . C A 999 MinMQ DP=96;VDB=0.0334;AF1=1;AC1=8;DP4=0,0,54,35;MQ=23;FQ=-36;AN=8;AC=8 GT:PL:DP:SP:GQ 1/1:163,60,0:20:0:99 1/1:124,30,0:10:0:72 1/1:198,69,0:23:0:99 1/1:203,108,0:36:0:99
1 55926 . T C 999 MinMQ DP=56;VDB=0.0269;AF1=1;AC1=8;DP4=0,0,23,32;MQ=14;FQ=-11.9;AN=8;AC=8 GT:PL:DP:SP:GQ 1/1:114,30,0:10:0:48 1/1:8,6,0:2:0:24 1/1:130,63,0:21:0:81 1/1:78,66,0:22:0:84
1 57376 . C T 16 MinMQ DP=143;VDB=0.0237;AF1=0.1883;AC1=2;DP4=70,55,2,13;MQ=28;FQ=16;PV4=0.002,0.034,0.001,0.027;AN=8;AC=2 GT:PL:DP:SP:GQ 0/1:4,0,201:26:8:3 0/1:18,0,188:19:10:15 0/0:0,26,246:38:15:29 0/0:0,154,255:57:0:99
1 57856 . T A 999 MinMQ DP=191;VDB=0.0263;AF1=0.5;G3=5.154e-23,1,1.244e-12;HWE=0.0185;AC1=4;DP4=58,51,29,51;MQ=21;FQ=999;PV4=0.027,0.028,1,1;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:204,0,154:54:13:99 0/1:121,0,52:29:9:55 0/1:142,0,175:58:2:99 0/1:104,0,198:48:0:99
1 57952 . A C 118 MinMQ DP=30;VDB=0.0356;AF1=0.8939;AC1=7;DP4=1,1,8,19;MQ=10;FQ=7.51;PV4=0.53,0.23,0.0064,0.46;AN=8;AC=7 GT:PL:DP:SP:GQ 1/1:36,21,0:7:0:27 1/1:17,6,0:2:0:12 1/1:40,33,0:11:0:39 0/1:29,0,12:9:4:7
1 58176 . G A 94.7 MinMQ DP=93;VDB=0.0330;AF1=0.3746;AC1=3;DP4=51,13,15,9;MQ=17;FQ=94.7;PV4=0.11,0.0027,1,1;AN=8;AC=3 GT:PL:DP:SP:GQ 0/1:30,0,23:22:0:26 0/1:18,0,15:12:7:17 0/1:55,0,102:29:2:56 0/0:0,42,114:25:9:41
1 58211 . A G 999 MinMQ DP=46;VDB=0.0332;AF1=1;AC1=8;DP4=0,0,30,15;MQ=22;FQ=-11.8;AN=8;AC=8 GT:PL:DP:SP:GQ 1/1:113,30,0:10:0:47 1/1:7,6,0:2:0:23 1/1:148,60,0:20:0:77 1/1:143,39,0:13:0:56
1 58771 . T C 999 MinMQ DP=263;VDB=0.0270;AF1=0.375;AC1=3;DP4=99,85,29,44;MQ=20;FQ=999;PV4=0.053,0.46,1,1;AN=8;AC=3 GT:PL:DP:SP:GQ 0/1:95,0,124:60:2:96 0/0:0,123,209:41:0:99 0/1:179,0,208:101:6:99 0/1:94,0,169:55:6:95
1 58866 . C G 119 MinMQ DP=233;VDB=0.0293;AF1=0.2581;AC1=2;DP4=77,98,12,40;MQ=18;FQ=119;PV4=0.0092,0.058,1,0.19;AN=8;AC=2 GT:PL:DP:SP:GQ 0/1:58,0,112:49:1:56 0/0:0,102,152:34:0:99 0/1:69,0,180:70:11:67 0/0:0,10,148:74:16:12
1 60332 . T C 97.6 MinMQ DP=239;VDB=0.0192;AF1=0.3089;AC1=2;DP4=77,104,22,32;MQ=17;FQ=97.6;PV4=0.88,1,1,1;AN=8;AC=2 GT:PL:DP:SP:GQ 0/0:0,0,99:61:9:3 0/0:0,111,185:37:0:99 0/1:70,0,145:68:7:70 0/1:33,0,123:69:0:33
1 61219 . T C 42.6 MinMQ DP=180;VDB=0.0243;AF1=0.1351;AC1=1;DP4=41,105,18,15;MQ=24;FQ=42.6;PV4=0.0069,0.27,0.32,1;AN=8;AC=1 GT:PL:DP:SP:GQ 0/0:0,6,128:36:5:11 0/0:0,57,168:31:0:62 0/1:48,0,227:53:21:43 0/0:0,16,255:59:6:21
1 61442 . A G 999 PASS DP=96;VDB=0.0348;AF1=1;AC1=8;DP4=0,0,41,47;MQ=30;FQ=-27;AN=8;AC=8 GT:PL:DP:SP:GQ 1/1:197,42,0:14:0:75 1/1:108,21,0:7:0:54 1/1:255,99,0:33:0:99 1/1:255,102,0:34:0:99
1 61499 . G A 87.1 PASS DP=140;VDB=0.0120;AF1=0.3006;AC1=2;DP4=54,55,18,12;MQ=35;FQ=87.1;PV4=0.41,0.3,2.4e-12,1;AN=8;AC=2 GT:PL:DP:SP:GQ 0/1:61,0,243:24:2:60 0/0:0,1,231:14:7:4 0/1:32,0,255:43:3:31 0/0:0,48,255:58:6:49
1 61579 . G A 107 MinMQ DP=161;VDB=0.0304;AF1=0.25;AC1=2;DP4=88,29,32,10;MQ=20;FQ=107;PV4=1,0.12,0.015,0.43;AN=8;AC=2 GT:PL:DP:SP:GQ 0/0:0,65,171:34:3:67 0/0:0,32,152:18:3:34 0/1:40,0,174:56:0:38 0/1:75,0,207:51:1:73
1 61987 . A G 999 PASS DP=206;VDB=0.0287;AF1=0.5;G3=1.244e-38,1,7.859e-46;HWE=0.0185;AC1=4;DP4=46,65,42,48;MQ=39;FQ=999;PV4=0.48,0.072,0.00033,1;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:255,0,255:58:0:99 0/1:182,0,218:22:2:99 0/1:255,0,255:62:3:99 0/1:220,0,255:59:6:99
1 61989 . G C 999 PASS DP=208;VDB=0.0311;AF1=0.5;G3=3.141e-40,1,8.15e-49;HWE=0.0185;AC1=4;DP4=47,65,42,49;MQ=39;FQ=999;PV4=0.57,0.058,8.7e-05,1;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:255,0,255:59:1:99 0/1:190,0,233:22:2:99 0/1:255,0,255:62:3:99 0/1:216,0,255:60:6:99
1 62203 . T C 999 PASS DP=258;VDB=0.0354;AF1=0.5;G3=3.125e-31,1,5e-52;HWE=0.0185;AC1=4;DP4=77,73,47,52;MQ=40;FQ=999;PV4=0.61,1,2.6e-25,1;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:255,0,255:66:5:99 0/1:145,0,255:38:1:99 0/1:252,0,255:84:1:99 0/1:210,0,255:61:0:99
1 62239 . TACACACACACACACACA TACACACACACACACA 999 PASS INDEL;DP=223;VDB=0.0280;AF1=0.4961;G3=3.069e-08,1,4.923e-103;HWE=0.0192;AC1=4;DP4=83,54,34,25;MQ=41;FQ=999;PV4=0.75,0.056,2.1e-17,1;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:248,0,255:45:1:99 0/1:12,0,255:23:9:15 0/1:183,0,255:68:2:99 0/1:158,0,255:60:6:99
1 62271 . A G 134 PASS DP=187;VDB=0.0101;AF1=0.498;G3=2.233e-09,1,4.959e-103;HWE=0.0189;AC1=4;DP4=92,56,14,20;MQ=41;FQ=134;PV4=0.033,0.013,1.9e-22,0.0028;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:65,0,255:35:1:68 0/1:19,0,255:28:2:22 0/1:17,0,255:60:11:20 0/1:39,0,255:59:14:42
1 62777 . A T 999 MinMQ DP=251;VDB=0.0308;AF1=0.499;G3=3.104e-08,1,1.551e-35;HWE=0.0187;AC1=4;DP4=80,108,35,26;MQ=21;FQ=999;PV4=0.055,1,1,0.39;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:115,0,137:69:1:99 0/1:124,0,169:34:6:99 0/1:87,0,202:76:5:90 0/1:18,0,109:70:3:21
1 63735 . CCTACTA CCTA 999 MinMQ INDEL;DP=141;VDB=0.0354;AF1=0.5;G3=9.836e-28,1,9.22e-15;HWE=0.0185;AC1=4;DP4=36,24,32,41;MQ=20;FQ=216;PV4=0.082,1.2e-09,1,0.023;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:83,0,86:30:6:84 0/1:130,0,56:21:0:59 0/1:184,0,40:43:5:43 0/1:181,0,55:39:3:58
1 64613 . T A 999 MinMQ DP=328;VDB=0.0091;AF1=0.5;G3=7.862e-12,1,3.13e-37;HWE=0.0185;AC1=4;DP4=134,130,23,36;MQ=27;FQ=999;PV4=0.11,1,1,0.36;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:79,0,175:95:0:82 0/1:48,0,208:47:2:51 0/1:193,0,248:95:13:99 0/1:107,0,211:86:0:99
1 66162 . A T 999 PASS DP=215;VDB=0.0231;AF1=0.4998;G3=1.238e-10,1,1.58e-77;HWE=0.0186;AC1=4;DP4=62,67,26,36;MQ=39;FQ=999;PV4=0.44,1,3.3e-21,1;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:170,0,255:45:3:99 0/1:77,0,255:30:0:80 0/1:183,0,255:69:3:99 0/1:26,0,255:47:0:29
1 66442 . TATATAATATA TATATAATATAATATA 132 PASS INDEL;DP=233;VDB=0.0328;AF1=0.3333;AC1=3;DP4=64,69,21,10;MQ=42;FQ=135;PV4=0.071,1,3.5e-27,0.11;AN=8;AC=3 GT:PL:DP:SP:GQ 0/1:62,0,255:30:2:62 0/0:0,50,255:37:0:50 0/1:91,0,255:65:8:91 0/1:3,0,255:32:2:5
1 66507 . T A 999 PASS DP=202;VDB=0.0385;AF1=0.626;AC1=5;DP4=25,14,63,82;MQ=42;FQ=999;PV4=0.03,0.023,1,0.0014;AN=8;AC=5 GT:PL:DP:SP:GQ 0/1:255,0,205:42:7:99 0/1:255,0,20:37:12:21 0/1:255,0,155:57:4:99 1/1:255,72,0:48:0:71
1 66521 . TATATAATATA TATATAATATAATATA 999 PASS INDEL;DP=200;VDB=0.0384;AF1=0.3747;AC1=3;DP4=61,75,25,12;MQ=43;FQ=999;PV4=0.016,1,3.8e-20,0.38;AN=8;AC=3 GT:PL:DP:SP:GQ 0/1:233,0,255:40:7:99 0/1:25,0,255:32:16:26 0/1:178,0,255:56:3:99 0/0:0,75,255:45:3:74
1 69511 . A G 999 MinMQ DP=79;VDB=0.0355;AF1=1;AC1=8;DP4=1,0,44,31;MQ=18;FQ=-30.9;PV4=1,1,1,1;AN=8;AC=8 GT:PL:DP:SP:GQ 1/1:133,42,0:14:0:79 1/1:95,25,0:12:0:62 1/1:170,57,0:19:0:94 1/1:192,93,0:31:0:99
1 70300 . C T 14.2 MinMQ DP=147;VDB=0.0063;AF1=0.1206;AC1=1;DP4=63,68,7,5;MQ=19;FQ=14.2;PV4=0.56,1,0.057,0.22;AN=8;AC=1 GT:PL:DP:SP:GQ 0/1:20,0,181:35:0:15 0/0:0,39,144:25:7:45 0/0:0,122,201:46:0:99 0/0:0,111,207:37:0:99
1 73822 . A G 161 MinMQ DP=175;VDB=0.0290;AF1=0.2494;AC1=2;DP4=67,87,9,10;MQ=25;FQ=161;PV4=0.81,0.5,1,0.45;AN=8;AC=2 GT:PL:DP:SP:GQ 0/1:25,0,203:41:1:23 0/1:144,0,170:29:6:99 0/0:0,178,255:59:0:99 0/0:0,132,255:44:0:99
1 73841 . C T 999 PASS DP=182;VDB=0.0366;AF1=0.3748;AC1=3;DP4=50,64,12,26;MQ=30;FQ=999;PV4=0.25,1.6e-10,0.084,1;AN=8;AC=3 GT:PL:DP:SP:GQ 0/1:95,0,255:33:3:96 0/1:174,0,204:27:9:99 0/1:28,0,255:53:17:29 0/0:0,64,255:39:6:63
1 74092 . G A 26.8 MinMQ DP=158;VDB=0.0267;AF1=0.2721;G3=0.7501,7.846e-07,0.2499;HWE=0.0437;AC1=2;DP4=91,48,4,10;MQ=11;FQ=26.8;PV4=0.0093,0.39,1,1;AN=8;AC=2 GT:PL:DP:SP:GQ 0/0:0,5,79:32:10:7 1/1:40,26,0:33:7:19 0/0:0,105,85:35:0:93 0/0:0,160,38:53:0:46
1 79033 . A G 217 MinMQ DP=19;VDB=0.0139;AF1=1;AC1=8;DP4=0,0,12,7;MQ=18;FQ=-12.7;AN=8;AC=8 GT:PL:DP:SP:GQ 1/1:77,18,0:6:0:36 1/1:56,12,0:4:0:30 1/1:28,18,0:6:0:36 1/1:56,9,0:3:0:27
1 79050 . G T 258 MinMQ DP=29;VDB=0.0203;AF1=1;AC1=8;DP4=0,0,15,11;MQ=16;FQ=-18.3;AN=8;AC=8 GT:PL:DP:SP:GQ 1/1:87,21,0:7:0:45 1/1:52,15,0:5:0:39 1/1:45,24,0:8:0:48 1/1:74,18,0:6:0:42
1 79418 . G C 28.5 PASS DP=99;VDB=0.0139;AF1=0.2016;AC1=2;DP4=31,59,1,5;MQ=39;FQ=28.5;PV4=0.66,0.015,0.045,0.00068;AN=8;AC=2 GT:PL:DP:SP:GQ 0/1:5,0,236:21:3:4 0/1:30,0,229:18:0:27 0/0:0,78,255:26:0:81 0/0:0,93,255:31:0:96
1 79772 . C G 999 PASS DP=138;VDB=0.0342;AF1=0.25;AC1=2;DP4=68,47,11,9;MQ=30;FQ=999;PV4=0.81,1,0.41,1;AN=8;AC=2 GT:PL:DP:SP:GQ 0/1:132,0,127:29:6:99 0/1:84,0,219:27:2:82 0/0:0,117,255:39:0:99 0/0:0,120,255:40:0:99
1 82115 . A G 51.1 MinMQ DP=137;VDB=0.0291;AF1=0.1264;AC1=1;DP4=71,55,3,7;MQ=27;FQ=51.1;PV4=0.19,0.085,1,1;AN=8;AC=1 GT:PL:DP:SP:GQ 0/1:57,0,204:36:6:52 0/0:0,14,206:20:0:19 0/0:0,123,255:41:0:99 0/0:0,117,255:39:0:99
1 82133 . CAAAAAAAAAAAAAAAAAAAAA CAAAAAAAAAAAAAAAAAAA,CAAAAAAAAAAAAAA 81.9 PASS INDEL;DP=107;VDB=0.0354;AF1=1;AC1=8;DP4=0,0,5,17;MQ=37;FQ=-22.7;AN=8;AC=7,1 GT:PL:DP:SP:GQ 1/2:83,73,58,31,0,18:9:0:29 1/1:10,3,0,10,3,10:1:0:17 1/1:32,15,0,32,15,32:7:0:29 1/1:32,15,0,32,15,32:5:0:29
1 82303 . T C 21 PASS DP=111;VDB=0.0241;AF1=0.1243;AC1=1;DP4=47,50,3,8;MQ=38;FQ=21;PV4=0.22,1,1.5e-14,1;AN=8;AC=1 GT:PL:DP:SP:GQ 0/1:27,0,255:20:5:22 0/0:0,24,255:21:0:29 0/0:0,96,255:32:0:99 0/0:0,75,255:35:0:80
1 82456 . A G 49.6 MinMQ DP=151;VDB=0.0367;AF1=0.2495;AC1=2;DP4=77,55,15,0;MQ=20;FQ=49.6;PV4=0.0011,0.29,1,1;AN=8;AC=2 GT:PL:DP:SP:GQ 0/1:31,0,124:31:8:29 0/1:27,0,107:23:9:25 0/0:0,151,251:50:0:99 0/0:0,129,255:43:0:99
1 82676 . T G 999 PASS DP=152;VDB=0.0213;AF1=0.25;AC1=2;DP4=70,59,9,11;MQ=34;FQ=999;PV4=0.48,0.37,0.0004,1;AN=8;AC=2 GT:PL:DP:SP:GQ 0/1:124,0,255:34:0:99 0/1:89,0,255:31:6:87 0/0:0,138,255:46:0:99 0/0:0,114,255:38:0:99
1 83084 . T A 999 PASS DP=84;VDB=0.0128;AF1=1;AC1=8;DP4=0,0,38,37;MQ=37;FQ=-33;AN=8;AC=8 GT:PL:DP:SP:GQ 1/1:255,48,0:16:0:87 1/1:203,27,0:9:0:66 1/1:255,72,0:24:0:99 1/1:255,78,0:26:0:99
1 83514 . C T 88.2 MinMQ DP=139;VDB=0.0336;AF1=0.4441;AC1=4;DP4=54,53,22,8;MQ=14;FQ=87.9;PV4=0.037,0.05,1,0.35;AN=8;AC=4 GT:PL:DP:SP:GQ 1/1:30,5,0:38:2:4 0/0:0,30,24:20:6:24 0/1:17,0,23:36:4:18 0/1:51,0,62:43:8:53
1 83786 . TAAAAAAAA TAAAAAAAAAAA 134 PASS INDEL;DP=144;VDB=0.0396;AF1=0.2505;AC1=2;DP4=43,54,4,9;MQ=40;FQ=137;PV4=0.39,1,2.7e-06,1;AN=8;AC=2 GT:PL:DP:SP:GQ 0/1:47,0,67:25:0:45 0/1:113,0,16:11:0:24 0/0:0,84,98:29:0:86 0/0:0,135,116:45:0:99
1 83829 . GAGAAAGAAAGAAAGAAAGAAAGAAAGAAAGAAAGAAAGAAAGAA GAGAAAGAAAGAAAGAAAGAAAGAAAGAAAGAAAGAAAGAA 999 PASS INDEL;DP=131;VDB=0.0388;AF1=0.5;AC1=4;DP4=20,14,15,17;MQ=39;FQ=999;PV4=0.46,1,7e-05,1;AN=8;AC=4 GT:PL:DP:SP:GQ 1/1:255,39,0:13:0:36 0/1:36,0,124:6:0:39 0/1:255,0,203:27:2:99 0/0:0,60,255:20:0:57
1 83895 . GAGAAAGAAAGAAAGAAAGA GAGAAAGAAAGAAAGA 999 PASS INDEL;DP=151;VDB=0.0325;AF1=0.75;AC1=6;DP4=18,19,38,34;MQ=44;FQ=999;PV4=0.69,0.00045,0.12,1;AN=8;AC=6 GT:PL:DP:SP:GQ 1/1:255,60,0:20:0:62 1/1:255,48,0:16:0:50 0/1:255,0,255:37:1:99 0/1:255,0,255:36:3:99
1 84010 . G A 37 PASS DP=190;VDB=0.0033;AF1=0.125;AC1=1;DP4=85,71,6,3;MQ=38;FQ=37;PV4=0.73,0.2,1,0.14;AN=8;AC=1 GT:PL:DP:SP:GQ 0/0:0,126,255:42:0:99 0/0:0,81,255:27:0:86 0/1:43,0,255:38:2:38 0/0:0,64,255:58:2:69
1 84014 . G A 38 PASS DP=192;VDB=0.0055;AF1=0.125;AC1=1;DP4=89,67,3,4;MQ=39;FQ=38;PV4=0.47,0.21,1,0.021;AN=8;AC=1 GT:PL:DP:SP:GQ 0/1:44,0,255:42:0:39 0/0:0,84,255:28:0:89 0/0:0,72,255:36:0:77 0/0:0,172,255:57:0:99
1 84018 . G A 79.6 PASS DP=188;VDB=0.0107;AF1=0.2497;AC1=2;DP4=77,64,8,3;MQ=41;FQ=79.6;PV4=0.35,0.17,1,0.0032;AN=8;AC=2 GT:PL:DP:SP:GQ 0/0:0,120,255:40:0:99 0/1:28,0,255:24:0:26 0/0:0,108,255:36:0:99 0/1:60,0,255:52:6:58
1 84244 . A C 999 PASS DP=213;VDB=0.0276;AF1=0.25;AC1=2;DP4=83,93,14,22;MQ=41;FQ=999;PV4=0.46,0.29,0.019,1;AN=8;AC=2 GT:PL:DP:SP:GQ 0/1:255,0,255:47:6:99 0/1:208,0,255:49:1:99 0/0:0,178,255:59:0:99 0/0:0,172,255:57:0:99
1 85597 . A C 999 PASS DP=139;VDB=0.0342;AF1=0.25;AC1=2;DP4=47,60,16,14;MQ=30;FQ=999;PV4=0.41,0.057,1,1;AN=8;AC=2 GT:PL:DP:SP:GQ 0/1:198,0,135:26:0:99 0/1:184,0,221:35:5:99 0/0:0,114,255:38:0:99 0/0:0,114,255:38:0:99
1 86018 . C G 999 PASS DP=181;VDB=0.0399;AF1=0.25;AC1=2;DP4=70,69,13,26;MQ=42;FQ=999;PV4=0.07,1,7e-14,0.12;AN=8;AC=2 GT:PL:DP:SP:GQ 0/1:200,0,255:36:11:99 0/1:240,0,255:46:2:99 0/0:0,129,255:43:0:99 0/0:0,160,255:53:0:99
1 86303 . G T 999 PASS DP=182;VDB=0.0329;AF1=0.25;AC1=2;DP4=76,66,17,22;MQ=40;FQ=999;PV4=0.28,1,3.7e-11,0.31;AN=8;AC=2 GT:PL:DP:SP:GQ 0/1:255,0,255:44:3:99 0/1:207,0,255:31:1:99 0/0:0,138,255:46:0:99 0/0:0,181,255:60:0:99
1 86331 . A G 999 PASS DP=187;VDB=0.0331;AF1=0.25;AC1=2;DP4=69,74,18,23;MQ=40;FQ=999;PV4=0.72,1,2.4e-05,1;AN=8;AC=2 GT:PL:DP:SP:GQ 0/1:255,0,255:51:1:99 0/1:216,0,255:34:0:99 0/0:0,120,255:40:0:99 0/0:0,178,255:59:0:99
1 86656 . G T 36 MinMQ DP=148;VDB=0.0132;AF1=0.125;AC1=1;DP4=76,63,0,9;MQ=28;FQ=36;PV4=0.0012,1,0.028,1;AN=8;AC=1 GT:PL:DP:SP:GQ 0/1:42,0,220:30:8:37 0/0:0,39,255:24:6:44 0/0:0,126,255:42:0:99 0/0:0,157,255:52:0:99
1 89677 . A G 195 MinMQ DP=136;VDB=0.0371;AF1=0.25;AC1=2;DP4=46,51,14,17;MQ=23;FQ=195;PV4=0.84,0.28,1,1;AN=8;AC=2 GT:PL:DP:SP:GQ 0/0:0,72,186:24:0:74 0/0:0,63,181:21:0:65 0/1:130,0,145:39:1:99 0/1:73,0,155:44:1:71
1 91072 . A G 200 MinMQ DP=88;VDB=0.0383;AF1=0.25;AC1=2;DP4=44,22,11,10;MQ=25;FQ=200;PV4=0.3,1,1,0.41;AN=8;AC=2 GT:PL:DP:SP:GQ 0/0:0,42,166:14:0:44 0/0:0,33,123:11:0:35 0/1:114,0,203:38:10:99 0/1:94,0,104:24:2:92
1 91075 . T C 187 MinMQ DP=85;VDB=0.0399;AF1=0.25;AC1=2;DP4=41,21,11,9;MQ=26;FQ=187;PV4=0.43,1,1,1;AN=8;AC=2 GT:PL:DP:SP:GQ 0/0:0,42,148:14:0:44 0/0:0,33,136:11:0:35 0/1:100,0,186:36:8:98 0/1:95,0,112:21:5:93
1 91336 . A T 999 MinMQ DP=69;VDB=0.0304;AF1=0.7419;AC1=6;DP4=11,19,1,37;MQ=20;FQ=216;PV4=0.0003,1,1,1;AN=8;AC=6 GT:PL:DP:SP:GQ 1/1:51,10,0:11:3:12 1/1:104,39,0:13:0:41 0/1:54,0,110:17:6:62 0/1:61,0,114:27:13:69
1 98929 . A G 28 MinMQ DP=172;VDB=0.0198;AF1=0.1249;AC1=1;DP4=63,78,12,9;MQ=21;FQ=28;PV4=0.35,0.23,0.15,1;AN=8;AC=1 GT:PL:DP:SP:GQ 0/0:0,90,234:30:0:95 0/1:34,0,191:43:3:29 0/0:0,114,245:38:0:99 0/0:0,28,150:51:6:33
1 98999 . TTTTATTTATTTATTTATTTATTTATTTATTTATTTATTTATTT TTTTATTTATTTATTTATTTATTTATTTATTTATTTATTTATTTATTT 999 PASS INDEL;DP=236;VDB=0.0395;AF1=0.5;G3=1.982e-19,1,1.257e-17;HWE=0.0185;AC1=4;DP4=33,36,32,46;MQ=35;FQ=999;PV4=0.51,1,0.0067,0.38;AN=8;AC=4 GT:PL:DP:SP:GQ 0/1:112,0,192:30:0:99 0/1:255,0,77:31:4:80 0/1:86,0,191:32:2:89 0/1:255,0,100:54:12:99
1 101686 . A G 89 PASS DP=304;VDB=0.0327;AF1=0.125;AC1=1;DP4=99,171,10,15;MQ=30;FQ=89;PV4=0.83,0.3,0.035,1;AN=8;AC=1 GT:PL:DP:SP:GQ 0/1:95,0,255:58:0:90 0/0:0,72,255:24:0:77 0/0:0,101,255:102:0:99 0/0:0,255,255:111:0:99
X 1 . A G 89 PASS DP=304;VDB=0.0327;AF1=0.125;AC1=1;DP4=99,171,10,15;MQ=30;FQ=89;PV4=0.83,0.3,0.035,1;AN=8;AC=1 GT:PL:DP 0/1:95,0,255:11 0/0:0,72,255:11 0/0:0,101,255:11 0:0,255:11
X 2 . A G 89 PASS DP=304;VDB=0.0327;AF1=0.125;AC1=1;DP4=99,171,10,15;MQ=30;FQ=89;PV4=0.83,0.3,0.035,1;AN=8;AC=1 GT:PL:DP 0/1:95,0,255:11 0/0:0,72,255:11 0/0:0,101,255:11 1:255,0:11
X 3 . A G 89 PASS DP=304;VDB=0.0327;AF1=0.125;AC1=1;DP4=99,171,10,15;MQ=30;FQ=89;PV4=0.83,0.3,0.035,1;AN=8;AC=1 GT:PL:DP 0/0:0,95,255:11 0/0:0,72,255:11 0/0:0,101,255:11 1:255,0:11
vcftools-0.1.16/examples/fill-an-ac.out 0000664 0000000 0000000 00000002764 13331057401 0017716 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.0
##INFO=
##FORMAT=
##FORMAT=
##FORMAT=
##FILTER=
##INFO=
##INFO=
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT A
1 100 . GTTT G 1806 q10 AC=1;AN=2;DP=35 GT:GQ:DP 0/1:409:35
1 110 . C T,G 1792 PASS AC=1,0;AN=2;DP=32 GT:GQ:DP 0/1:245:32
1 110 . CAAA C 1792 PASS AC=1;AN=2;DP=32 GT:GQ:DP 0/1:245:32
1 120 . GA G 628 q10 AC=2;AN=2;DP=21 GT:GQ:DP 1/1:21:21
1 130 . G T 1016 PASS AC=1;AN=2;DP=22 GT:GQ:DP 0/1:212:22
1 130 . GAA GG 1016 PASS AC=1;AN=2;DP=22 GT:GQ:DP 0/1:212:22
1 140 . GT G 727 PASS AC=1;AN=2;DP=30 GT:GQ:DP 0/1:150:30
1 150 . TAAAA TA,T 246 PASS AC=1,1;AN=2;DP=10 GT:GQ:DP 1/2:12:10
1 160 . TAAAA TA,T 246 PASS AC=1,1;AN=2;DP=10 GT:GQ:DP 1/2:12:10
2 100 . GTTT G 1806 q10 AC=1;AN=2;DP=35 GT:GQ:DP 0/1:409:35
2 110 . CAAA C 1792 PASS AC=1;AN=2;DP=32 GT:GQ:DP 0/1:245:32
2 120 . GA G 628 q10 AC=2;AN=2;DP=21 GT:GQ:DP 1/1:21:21
2 130 . GAA G 1016 PASS AC=1;AN=2;DP=22 GT:GQ:DP 0/1:212:22
2 140 . GT G 727 PASS AC=1;AN=2;DP=30 GT:GQ:DP 0/1:150:30
2 150 . TAAAA TA,T 246 PASS AC=1,1;AN=2;DP=10 GT:GQ:DP 1/2:12:10
2 160 . TAAAA TA,TC,T 246 PASS AC=0,1,0;AN=2;DP=10 GT:GQ:DP 0/2:12:10
vcftools-0.1.16/examples/filters.txt 0000664 0000000 0000000 00000011445 13331057401 0017467 0 ustar 00root root 0000000 0000000 # Examples of user-defined filters. Edit and run with -f filters.txt.
# The examples below are self-explanatory. Notice the use of the predefined
# variables ($PASS, $FAIL, $MATCH, $RECORD) and methods (error).
# In this example, a minimum value of AF1=0.1 is required
{
tag => 'INFO/AF1', # The VCF tag to apply this filter on
name => 'MinAF', # The filter ID
desc => 'Minimum AF1 [0.01]', # Description for the VCF header
test => sub { return $MATCH < 0.01 ? $FAIL : $PASS },
},
# Filter all indels (presence of INDEL tag is tested)
{
tag => 'INFO/INDEL',
apply_to => 'indels', # Can be one of SNPs, indels, all. Default: [All]
name => 'Indel',
desc => 'INDEL tag present',
test => sub { return $FAIL },
},
# Only loci with enough reads supporting the variant will pass the filter
{
tag => 'INFO/DP4',
name => 'FewAlts',
desc => 'Too few reads supporting the variant',
apply_to => 'SNPs',
test => sub {
if ( !($MATCH =~ /^([^,]+),([^,]+),([^,]+),(.+)$/) )
{
error("Could not parse INFO/DP4: $CHROM:$POS [$MATCH]");
}
if ( 0.1*($1+$2) > $3+$4 ) { return $PASS; }
return $FAIL;
},
},
# Example of filtering based on genotype columns and the QUAL column
{
tag => 'FORMAT/PL',
name => 'NoHets',
desc => 'Inbred homozygous mouse, no hets expected',
apply_to => 'SNPs',
test => sub {
for my $pl (@$MATCH)
{
my @pls = split(/,/,$pl);
if ( $pls[1]<$pls[0] && $pls[1]<$pls[2] ) { return $FAIL; }
}
return $PASS;
},
},
# This example splits the four PV4 values into four tags names PV0, PV1, PV2 and PV3.
# Note the use of the 'header' key, and the $RECORD and $VCF variables.
{
header => [
qq[key=INFO,ID=PV0,Number=1,Type=Float,Description="P-value for strand bias"],
qq[key=INFO,ID=PV1,Number=1,Type=Float,Description="P-value for baseQ bias"],
qq[key=INFO,ID=PV2,Number=1,Type=Float,Description="P-value for mapQ bias"],
qq[key=INFO,ID=PV3,Number=1,Type=Float,Description="P-value for tail distance bias"]
],
tag => 'INFO/PV4',
name => 'SplitPV4',
desc => 'Split PV4',
apply_to => 'all',
test => sub {
my @vals = split(/,/,$MATCH);
$$RECORD[7] = $VCF->add_info_field($$RECORD[7],'PV0'=>$vals[0],'PV1'=>$vals[1],'PV2'=>$vals[2],'PV3'=>$vals[3]);
return $PASS;
},
},
# Do whatever you want with every record and edit it according to your needs. This silly
# example removes the tag SILLY in records where ID is set and depth is bigger than 5.
{
tag => 'Dummy',
test => sub {
if ( $$RECORD[2] eq '.' ) { return $PASS; } # Modify only lines with ID
my $dp = $vcf->get_info_field($$RECORD[7],'DP');
if ( $dp>5 ) { $$RECORD[7] = $VCF->add_info_field($$RECORD[7],'SILLY'=>undef); }
return $PASS;
},
}
# Filter records with the value XY absent or not equal to 42
{
tag => 'Dummy',
header => [
qq[key=FILTER,ID=XY,Description="XY not OK"],
],
test => sub {
my $xy = $VCF->get_info_field($$RECORD[7],'XY');
my $is_bad = ( !defined $xy or $xy!=42 ) ? 1 : 0;
$$RECORD[6] = $VCF->add_filter($$RECORD[6],'XY'=>$is_bad);
return $PASS;
},
},
# Annotate INFO field with SINGLETON flag when one and only one sample is different from the reference
{
header => [
qq[key=INFO,ID=SINGLETON,Number=0,Type=Flag,Description="Only one non-ref sample"],
],
tag => 'FORMAT/GT',
name => 'Dummy',
desc => 'Dummy',
test => sub {
my $nalt = 0;
for my $gt (@$MATCH)
{
my @gt = $VCF->split_gt($gt);
for my $allele (@gt)
{
if ( $allele ne 0 && $allele ne '.' ) { $nalt++; last; }
}
if ( $nalt>1 ) { last; }
}
if ( $nalt==1 ) { $$RECORD[7] = $VCF->add_info_field($$RECORD[7],'SINGLETON'=>''); }
return $PASS;
},
},
# Set genotypes to unknown ("." or "./." depending on ploidy) when coverage is low (by Shane McCarthy).
{
tag => 'FORMAT/DP',
name => 'MinSampleDP',
desc => 'Genotypes set to . for samples with DP < 2',
apply_to => 'all',
test => sub {
my $i = 8;
for my $dp (@$MATCH)
{
$i++;
next unless ($dp<2);
my @format = split(/:/,$$RECORD[$i]);
$format[0] = $format[0] =~ /\// ? "./." : ".";
$$RECORD[$i] = join(":",@format);
}
return $PASS;
},
},
vcftools-0.1.16/examples/fix-ploidy.out 0000664 0000000 0000000 00000003604 13331057401 0020071 0 ustar 00root root 0000000 0000000 61098 M1 0 0,9,72,5,6,7 M2 0 0,15,140,5,6,7 F3 1/1 147,0,5 F4 0/0 0,131,5 M5 0 0,9,83,5,6,7 M6 0 0,6,56,5,6,7
61270 M1 0 8,14,58 M2 0 0,6,52 F3 0/0 0,6,56 F4 0/0 0,15,117 M5 0 0,6,45 M6 0 0,12,87
61275 M1 0 0,3,13 M2 0 0,3,28 F3 0/0 8,0,41 F4 0/0 0,12,97 M5 0 0,6,49 M6 0 0,9,67
61282 M1 0/1 15,3,0 M2 0/0 0,6,51 F3 0 6,0,31 F4 0 0,6,57 M5 0/1 7,0,19 M6 0/1 16,0,20
61795 M1 0/0 0,27,203 M2 0/0 0,21,174 F3 0 0,45,229 F4 0 0,27,199 M5 0/0 0,24,182 M6 0/0 0,9,85
62731 M1 0/0 0,27,194 M2 0/0 0,24,194 F3 0 0,18,141 F4 0 0,30,201 M5 0/0 0,18,153 M6 0/0 0,33,202
63008 M1 0/0 0,42,255 M2 0/0 0,15,128 F3 0 0,15,136 F4 0 0,39,251 M5 0/0 0,15,111 M6 0/0 0,27,200
63231 M1 0/0 0,42,246 M2 0/0 0,18,141 F3 0 0,27,209 F4 0 0,24,186 M5 0/0 0,12,110 M6 0/0 0,21,145
63244 M1 0 0,36,209 M2 0 0,21,174 F3 . 0,27,198 F4 . 0,24,184 M5 0 0,15,132 M6 0 0,21,159
63328 M1 0 0,42,242 M2 0 0,12,110 F3 . 0,36,231 F4 . 0,36,226 M5 0 0,15,135 M6 0 0,18,132
63452 M1 0 0,27,200 M2 0 0,15,123 F3 . 0,33,228 F4 . 0,9,88 M5 0 0,24,171 M6 0 0,15,134
63799 M1 0 0,36,211 M2 0 0,27,205 F3 . 0,18,125 F4 . 0,15,125 M5 0 1,0,150 M6 0 0,12,106
63967 M1 0 0,30,183 M2 0 0,30,206 F3 . 0,30,206 F4 . 0,33,230 M5 0 0,21,160 M6 0 0,12,112
65288 M1 . 0,18,155 M2 . 0,12,113 F3 0 0,18,144 F4 0 0,21,155 M5 . 0,21,176 M6 . 0,6,63
65900 M1 . 162,24,0 M2 . 160,21,0 F3 1 219,30,0 F4 1 213,30,0 M5 . 248,42,0 M6 . 148,21,0
65951 M1 . 0,18,125 M2 . 0,15,132 F3 0 0,24,183 F4 0 0,45,252 M5 . 0,30,217 M6 . 0,15,101
66370 M1 . 255,57,0 M2 . 193,24,0 F3 1 97,12,0 F4 1 208,30,0 M5 . 129,15,0 M6 . 72,9,0
67184 M1 . 0,33,202 M2 . 0,6,57 F3 0 0,42,223 F4 0 0,21,142 M5 . 0,30,181 M6 . 0,33,201
67760 M1 . 0,42,224 M2 . 0,9,77 F3 0 0,45,243 F4 0 0,21,147 M5 . 0,33,205 M6 . 0,12,102
68303 M1 . 0,33,205 M2 . 0,33,236 F3 0 0,33,214 F4 0 0,27,197 M5 . 0,18,149 M6 . 0,9,72
68618 M1 0/0 0,24,176 M2 0/0 0,30,214 F3 0/0 0,21,159 F4 0/0 0,24,191 M5 0/0 0,15,133 M6 0/0 0,18,133
vcftools-0.1.16/examples/fix-ploidy.samples 0000664 0000000 0000000 00000000052 13331057401 0020720 0 ustar 00root root 0000000 0000000 M1 M
M2 M
F3 F
F4 F
M5 M
M6 M
vcftools-0.1.16/examples/fix-ploidy.txt 0000664 0000000 0000000 00000000337 13331057401 0020101 0 ustar 00root root 0000000 0000000 ploidy =>
{
20 =>
[
{ from=>1, to=>61275, M=>1, F=>2 },
{ from=>61282, to=>63231, F=>1 },
{ from=>63244, to=>63967, M=>1, F=>0 },
{ from=>65288, to=>68303, M=>0, F=>1 },
],
}
vcftools-0.1.16/examples/fix-ploidy.vcf 0000664 0000000 0000000 00000021455 13331057401 0020044 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.1
##samtoolsVersion=0.1.16-dev (r969:252)
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##FORMAT=
##FORMAT=
##FORMAT=
##FORMAT=
##FORMAT=
##FORMAT=
##source_20111007.1=/software/vertres/codebase/scripts/vcf-annotate -f +/D=1200
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##FILTER=
##source_20111007.2=/software/vertres/codebase/scripts/vcf-annotate -f +/D=1200
##source_20120109.1=vcf-subset(r660) -c QTL190284,QTL190301,QTL190321,QTL190576,QTL190627,QTL190628
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT M1 M2 F3 F4 M5 M6
20 61098 . C A,T 999 PASS AC1=41;AF1=0.2104;DP4=209,284,67,76;DP=658;FQ=999;MQ=45;PV4=0.39,4.4e-10,0.0034,0.2 GT:PL:DP:SP:GQ 0/1:0,9,72,5,6,7:3:212:12 0/0:0,15,140,5,6,7:5:458752:18 1:147,0,5:7:384:24 0:0,131,5:5:208:18 0/0:0,9,83,5,6,7:3:392:12 0/0:0,6,56,5,6,7:2:204:9
20 61270 . A T 93 PASS AC1=5;AF1=0.02733;DP4=149,185,5,10;DP=398;FQ=93.1;MQ=43;PV4=0.44,3.5e-05,0.028,1 GT:PL:DP:SP:GQ 0/0:8,14,58:3:0:19 0/0:0,6,52:2:0:19 0/0:0,6,56:2:0:19 0/0:0,15,117:5:0:28 0/0:0,6,45:2:0:19 0/0:0,12,87:4:0:25
20 61275 . T G 14.5 PASS AC1=14;AF1=0.07104;DP4=76,141,1,24;DP=345;FQ=14.8;MQ=43;PV4=0.001,4e-16,2.4e-07,0.39 GT:PL:DP:SP:GQ 0/0:0,3,13:1:16908804:12 0/0:0,3,28:1:0:12 0/0:8,0,41:3:201985031:3 0/0:0,12,97:4:134480904:21 0/0:0,6,49:2:117901063:15 0/0:0,9,67:3:33686020:18
20 61282 . T C 120 PASS AC1=35;AF1=0.1794;DP4=50,118,1,30;DP=333;FQ=121;G3=0.6049,0.3951,2.481e-10;HWE=0.0152;MQ=44;PV4=0.0012,1.1e-08,5.7e-11,1 GT:PL:DP:SP:GQ 0/1:15,3,0:1:17302017:6 0/0:0,6,51:2:0:10 0/1:6,0,31:3:437393682:4 0/0:0,6,57:2:134283524:10 0/1:7,0,19:2:404167697:5 0/1:16,0,20:2:67633410:13
20 61795 . G A 999 PASS AC1=63;AF1=0.3239;DP4=193,313,105,127;DP=753;FQ=999;MQ=48;PV4=0.075,4e-23,2.6e-26,1 GT:PL:DP:SP:GQ 0/0:0,27,203:9:-1598497894:27 0/0:0,21,174:7:50331648:21 0/0:0,45,229:15:0:45 0/0:0,27,199:9:1056977028:27 0/0:0,24,182:8:0:24 0/0:0,9,85:3:1678500058:10
20 62731 . C A 999 PASS AC1=24;AF1=0.1207;DP4=326,391,40,53;DP=846;FQ=999;MQ=49;PV4=0.74,6.3e-14,5.2e-07,0.44 GT:PL:DP:SP:GQ 0/0:0,27,194:9:134349316:33 0/0:0,24,194:8:0:30 0/0:0,18,141:6:505290270:24 0/0:0,30,201:10:16908801:36 0/0:0,18,153:6:505290270:24 0/0:0,33,202:11:67175432:39
20 63008 . C A 122 PASS AC1=2;AF1=0.01078;DP4=303,374,3,7;DP=692;FQ=122;MQ=49;PV4=0.52,0.0011,0.093,1 GT:PL:DP:SP:GQ 0/0:0,42,255:14:67371265:59 0/0:0,15,128:5:0:32 0/0:0,15,136:5:505290270:32 0/0:0,39,251:13:67634177:56 0/0:0,15,111:5:505290270:32 0/0:0,27,200:9:33818632:44
20 63231 . T A 999 PASS AC1=5;AF1=0.02753;DP4=289,324,7,18;DP=652;FQ=999;MQ=49;PV4=0.067,8.8e-11,0.004,0.49 GT:PL:DP:SP:GQ 0/0:0,42,246:14:0:54 0/0:0,18,141:6:0:30 0/0:0,27,209:9:0:39 0/0:0,24,186:8:0:36 0/0:0,12,110:4:0:24 0/0:0,21,145:7:0:33
20 63244 . A C 999 PASS AC1=37;AF1=0.1905;DP4=273,269,66,56;DP=670;FQ=999;MQ=49;PV4=0.48,4.1e-21,2.5e-21,1 GT:PL:DP:SP:GQ 0/0:0,36,209:12:858855379:39 0/0:0,21,174:7:7:24 0/0:0,27,198:9:0:30 0/0:0,24,184:8:1055941246:27 0/0:0,15,132:5:0:18 0/0:0,21,159:7:203312940:24
20 63328 . A C 26.4 PASS AC1=1;AF1=0.006786;DP4=439,259,2,3;DP=711;FQ=26.4;MQ=49;PV4=0.37,0.00092,0.2,0.31 GT:PL:DP:SP:GQ 0/0:0,42,242:14:0:61 0/0:0,12,110:4:0:31 0/0:0,36,231:12:0:55 0/0:0,36,226:12:0:55 0/0:0,15,135:5:0:34 0/0:0,18,132:6:0:37
20 63452 . C A 86.5 PASS AC1=4;AF1=0.02128;DP4=399,301,5,6;DP=718;FQ=86.6;MQ=49;PV4=0.54,0.084,0.0093,0.49 GT:PL:DP:SP:GQ 0/0:0,27,200:9:-419321220:41 0/0:0,15,123:5:0:29 0/0:0,33,228:11:0:47 0/0:0,9,88:3:1060858040:23 0/0:0,24,171:8:0:38 0/0:0,15,134:5:1094576049:29
20 63799 . C T 999 PASS AC1=68;AF1=0.347;DP4=215,280,110,139;DP=796;FQ=999;MQ=49;PV4=0.88,7.1e-84,0.16,1 GT:PL:DP:SP:GQ 0/0:0,36,211:12:212:36 0/0:0,27,205:9:50331648:27 0/0:0,18,125:6:384:18 0/0:0,15,125:5:208:15 0/1:1,0,150:8:392:4 0/0:0,12,106:5:204:12
20 63967 . A G 999 PASS AC1=5;AF1=0.02598;DP4=384,427,9,10;DP=833;FQ=999;MQ=49;PV4=1,0.0053,1.9e-05,0.043 GT:PL:DP:SP:GQ 0/0:0,30,183:10:0:43 0/0:0,30,206:10:0:43 0/0:0,30,206:10:0:43 0/0:0,33,230:11:0:46 0/0:0,21,160:7:0:34 0/0:0,12,112:4:0:25
20 65288 . G C 999 PASS AC1=23;AF1=0.1172;DP4=217,304,31,52;DP=612;FQ=999;MQ=49;PV4=0.47,1.3e-40,0.001,1 GT:PL:DP:SP:GQ 0/0:0,18,155:6:212:24 0/0:0,12,113:4:0:18 0/0:0,18,144:6:384:24 0/0:0,21,155:7:208:27 0/0:0,21,176:7:392:27 0/0:0,6,63:2:204:12
20 65900 . G A 999 PASS AC1=156;AF1=0.7977;DP4=98,72,334,335;DP=857;FQ=999;MQ=46;PV4=0.086,1,2.3e-21,1 GT:PL:DP:SP:GQ 1/1:162,24,0:8:-1210247533:27 1/1:160,21,0:7:458754:24 1/1:219,30,0:10:0:33 1/1:213,30,0:10:1056438877:33 1/1:248,42,0:14:0:45 1/1:148,21,0:7:-1125025851:24
20 65951 . T A 142 PASS AC1=4;AF1=0.02096;DP4=349,437,7,6;DP=818;FQ=142;MQ=48;PV4=0.58,0.0014,0.18,1 GT:PL:DP:SP:GQ 0/0:0,18,125:6:1040080185:32 0/0:0,15,132:5:0:29 0/0:0,24,183:8:0:38 0/0:0,45,252:15:1063282696:59 0/0:0,30,217:10:0:44 0/0:0,15,101:5:1099561607:29
20 66370 . G A 999 PASS AC1=160;AF1=0.8125;DP4=76,68,292,297;DP=774;FQ=999;MQ=46;PV4=0.52,1,5.3e-34,1 GT:PL:DP:SP:GQ 1/1:255,57,0:19:0:60 1/1:193,24,0:8:131075:27 1/1:97,12,0:4:0:15 1/1:208,30,0:10:0:33 1/1:129,15,0:5:0:18 1/1:72,9,0:3:0:13
20 67184 . C A 110 PASS AC1=2;AF1=0.01039;DP4=366,477,2,10;DP=863;FQ=111;MQ=48;PV4=0.08,1,0.12,0.24 GT:PL:DP:SP:GQ 0/0:0,33,202:11:-2137326772:50 0/0:0,6,57:2:0:23 0/0:0,42,223:14:0:59 0/0:0,21,142:7:1069487242:38 0/0:0,30,181:10:0:47 0/0:0,33,201:11:1190753285:50
20 67760 . C A 16.9 PASS AC1=2;AF1=0.008593;DP4=412,444,1,5;DP=870;FQ=16.9;MQ=47;PV4=0.22,1,0.0035,1 GT:PL:DP:SP:GQ 0/0:0,42,224:14:-1076013622:60 0/0:0,9,77:3:0:27 0/0:0,45,243:15:0:63 0/0:0,21,147:7:1070943924:39 0/0:0,33,205:11:0:51 0/0:0,12,102:4:783554555:30
20 68303 . A G 65.3 PASS AC1=1;AF1=0.005274;DP4=480,351,4,5;DP=874;FQ=65.3;MQ=49;PV4=0.51,9.5e-09,0.34,1 GT:PL:DP:SP:GQ 0/0:0,33,205:11:-1062059871:53 0/0:0,33,236:11:0:53 0/0:0,33,214:11:0:53 0/0:0,27,197:9:1072683922:47 0/0:0,18,149:6:0:38 0/0:0,9,72:3:1109840439:29
20 68618 . G C 62.6 PASS AC1=2;AF1=0.01065;DP4=409,382,5,5;DP=833;FQ=62.6;MQ=49;PV4=1,1.9e-07,0.00018,0.4 GT:PL:DP:SP:GQ 0/0:0,24,176:9:1561273746:41 0/0:0,30,214:10:0:47 0/0:0,21,159:7:0:38 0/0:0,24,191:8:1056871757:41 0/0:0,15,133:5:0:32 0/0:0,18,133:6:253656329:35
vcftools-0.1.16/examples/floats.vcf 0000664 0000000 0000000 00000001064 13331057401 0017242 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.1
##INFO=
##FORMAT=
##reference=file:/lustre/scratch105/projects/g1k/ref/main_project/human_g1k_v37.fasta
##contig=
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT NA00001
19 14370 . G A 29 PASS FLOATTAG=0.0001 GT 0|0
19 14371 . G A 29 PASS FLOATTAG=1e-4 GT 0|0
19 14372 . G A 29 PASS FLOATTAG=1E-4 GT 0|0
19 14373 . G A 29 PASS FLOATTAG=1e4 GT 0|0
vcftools-0.1.16/examples/indel-stats.out 0000664 0000000 0000000 00000000060 13331057401 0020225 0 ustar 00root root 0000000 0000000 total 20
in-frame 9
frameshift 5
ratio 0.357143
vcftools-0.1.16/examples/indel-stats.tab 0000664 0000000 0000000 00000000030 13331057401 0020161 0 ustar 00root root 0000000 0000000 1 20 30
1 40 50
1 60 80
vcftools-0.1.16/examples/indel-stats.vcf 0000664 0000000 0000000 00000001013 13331057401 0020173 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.1
#CHROM POS ID REF ALT QUAL FILTER INFO
1 15 . ACGT A . PASS .
1 15 . A ACGT . PASS .
1 18 . ACGT A . PASS .
1 18 . A ACGT . PASS .
1 25 . ACGT A . PASS .
1 25 . A ACGT . PASS .
1 27 . ACGTA A . PASS .
1 27 . A ACGT . PASS .
1 29 . ACGT A . PASS .
1 29 . A ACGT . PASS .
1 35 . ACGT A . PASS .
1 35 . A ACGT . PASS .
1 38 . ACGT A . PASS .
1 38 . A ACGT . PASS .
1 45 . ACGT A . PASS .
1 45 . A ACGT . PASS .
1 47 . ACGTA A . PASS .
1 47 . A AACGT . PASS .
1 49 . ACGT A . PASS .
1 49 . A ACGT . PASS .
vcftools-0.1.16/examples/invalid-4.0.vcf 0000664 0000000 0000000 00000004046 13331057401 0017702 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.0
##problem1=The first base of the second ALT allele at 20:1234567 does not match the reference.
##fileDate=20090805
##source=myImputationProgramV3.1
##reference=1000GenomesPilot-NCBI36
##phasing=partial
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##FILTER=
##FILTER=
##FORMAT=
##FORMAT=
##FORMAT=
##FORMAT=
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT NA00001 NA00002 NA00003
19 111 . A C 9.6 . . GT:HQ 0|0:10,10 0|0:10,10 0\1:3,3
19 112 . A G 10 . . GT:HQ 0|0:10,10 0|0:10,10 0\1:3,3
20 14370 rs6054257 G A 29 0 NS=3;DP=14;AF=0.5;DB;H2 GT:GQ:DP:HQ 0|0:48:1:51,51 1|0:48:8:51,51 1/1:43:5:-1,-1
20 17330 . T A 3 q10 NS=3;DP=11;AF=0.017 GT:GQ:DP:HQ 0|0:49:3:58,50 0|1:3:5:65,3 0/0:41:3:-1,-1
20 1110696 rs6040355 A G,T 67 0 NS=2;DP=10;AF=0.333,0.667;AA=T;DB GT:GQ:DP:HQ 1|2:21:6:23,27 2|1:2:0:18,2 2/2:35:4:-1,-1
20 1230237 . T . 47 0 NS=3;DP=13;AA=T GT:GQ:DP:HQ 0|0:54:-1:56,60 0|0:48:4:51,51 0/0:61:2:-1,-1
20 1234567 microsat1 G GA,AAC 50 0 NS=3;DP=9;AA=G;AN=6;AC=3,1 GT:GQ:DP 0/1:-1:4 0/2:17:2 1/1:40:3
20 1235237 . T . -1 . . GT 0\0 0|0 ./.
X 10 rsTest AC A,ATG 10 . . GT 0 0/1 0|2
X 11 rsTest2 T A,G 10 q10;s50 . GT:DP:GQ 0:3:10 .:5:20 0:3:10
vcftools-0.1.16/examples/isec-n2-test.vcf.out 0000664 0000000 0000000 00000002140 13331057401 0020771 0 ustar 00root root 0000000 0000000 Warning: The column names do not match (e.g. B):
B
A
##fileformat=VCFv4.0
##FILTER=
##FORMAT=
##FILTER=
##INFO=
##FORMAT=
##FORMAT=
##FORMAT=
##INFO=
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT A
1 3062915 . GTTT G 1806 q10 DP=35;DP4=1,2,3,4;SF=0f,1f,2 GT:GQ:DP:GL 0/1:409:35:-20,-5,-20
1 3106154 . CAAA C 1792 PASS DP=32;SF=0,1,2 GT:GQ:DP 0/1:245:32
1 3157410 . GA G 628 q10 DP=21;SF=0f,1,2 GT:GQ:DP 1/1:21:21
1 3162006 . GAA G 1016 PASS DP=22;SF=0,1,2 GT:GQ:DP 0/1:212:22
1 3177144 . GT G 727 PASS DP=30;SF=0,1,2 GT:GQ:DP 0/1:150:30
1 3184885 . TAAAA TA,T 246 PASS DP=10;SF=0,1 GT:GQ:DP 1/2:12:10
3 3199815 . G A 353 PASS DP=19;SF=1,2 GT:GQ:DP 0/1:188:19
vcftools-0.1.16/examples/merge-test-a.vcf 0000664 0000000 0000000 00000001705 13331057401 0020246 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.0
##INFO=
##FORMAT=
##FORMAT=
##FORMAT=
##FORMAT=
##FILTER=
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT A
1 3062915 . GTTT G 1806 q10 DP=35;DP4=1,2,3,4 GT:GQ:DP:GL 0/1:409:35:-20,-5,-20
1 3106154 . CAAA C 1792 PASS DP=32 GT:GQ:DP 0/1:245:32
1 3157410 . GA G 628 q10 DP=21 GT:GQ:DP 1/1:21:21
1 3162006 . GAA G 1016 PASS DP=22 GT:GQ:DP 0/1:212:22
1 3177144 . GT G 727 PASS DP=30 GT:GQ:DP 0/1:150:30
1 3184885 . TAAAA TA,T 246 PASS DP=10 GT:GQ:DP 1/2:12:10
2 3199812 . G GTT,GT 481 PASS DP=26 GT:GQ:DP 1/2:322:26
3 3212016 . CTT C,CT 565 PASS DP=26 GT:GQ:DP 1/2:91:26
4 3258448 . TACACACAC T 325 PASS DP=31 GT:GQ:DP 0/1:325:31
vcftools-0.1.16/examples/merge-test-b.vcf 0000664 0000000 0000000 00000001774 13331057401 0020255 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.0
##INFO=
##FORMAT=
##FORMAT=
##FORMAT=
##FORMAT=
##FILTER=
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT B
1 3062915 . GTTT GT 376 q20 DP=14;DP4=1,2,3,4 GT:GQ:DP:GL 0/1:376:14:-10,0,-10
1 3106154 . CAAAA C 677 PASS DP=15 GT:GQ:DP:GL 0/1:277:15:-10,0,-10
1 3157410 . GA G 249 PASS DP=11 GT:GQ:DP 0/1:49:11
1 3162006 . GAA G 663 PASS DP=19 GT:GQ:DP 0/1:589:19
1 3177144 . GT G 460 PASS DP=24 GT:GQ:DP 0/1:236:24
1 3184885 . TAAA T 598 PASS DP=16 GT:GQ:DP 0/1:435:16
2 3188209 . GA G 162 . DP=15 GT:GQ:DP 0/1:162:15
3 3199812 . G GTT,GT 353 PASS DP=19 GT:GQ:DP 1/2:188:19
3 3199815 . G A 353 PASS DP=19 GT:GQ:DP 0/1:188:19
4 3212016 . CTT C 677 q20 DP=15 GT:GQ:DP 0/1:158:15
vcftools-0.1.16/examples/merge-test-c.vcf 0000664 0000000 0000000 00000001436 13331057401 0020251 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.0
##INFO=
##FORMAT=
##FORMAT=
##FORMAT=
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT C
1 3062915 . GTTT G 388 . DP=10;DP4=1,2,3,4 GT:GQ:DP 0/1:149:10
1 3106154 . CAAA C 269 . DP=9 GT:GQ:DP 1/1:25:9
1 3157410 . GA G 212 . DP=10 GT:GQ:DP 0/1:52:10
1 3162006 . GAA G 558 . DP=17 GT:GQ:DP 0/1:163:17
1 3177144 . GT G 211 . DP=14 GT:GQ:DP 0/1:151:14
1 3199812 . G GT . . . GT 1/1
2 3212016 . CTT C 613 . DP=11 GT:GQ:DP 0/1:41:11
3 3199815 . G T 353 PASS DP=19 GT:GQ:DP 0/1:188:19
3 3242491 . TT T . . . GT 1/1
4 3291771 . T TAA,TAAA 336 . DP=12 GT:GQ:DP 1/2:2:12
vcftools-0.1.16/examples/merge-test.vcf.out 0000664 0000000 0000000 00000004300 13331057401 0020630 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.2
##FILTER=
##FORMAT=
##FILTER=
##INFO=
##FORMAT=
##FORMAT=
##FORMAT=
##INFO=
##INFO=
##INFO=
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT A B C
1 3062915 . GTTT G,GT 856.67 q20;q10 AC=2,1;AN=6;DP4=3,6,9,12;DP=59;SF=0f,1f,2 GT:GL:DP:GQ 0/1:-20,-5,-20,.,.,.:35:409 0/2:-10,.,.,0,.,-10:14:376 0/1:.:10:149
1 3106154 . CAAAA CA,C 912.67 PASS AC=3,1;AN=6;DP=56;SF=0,1,2 GT:GL:DP:GQ 0/1:.:32:245 0/2:-10,.,.,0,.,-10:15:277 1/1:.:9:25
1 3157410 . GA G 363.00 q10 AC=4;AN=6;DP=42;SF=0f,1,2 GT:DP:GQ 1/1:21:21 0/1:11:49 0/1:10:52
1 3162006 . GAA G 745.67 PASS AC=3;AN=6;DP=58;SF=0,1,2 GT:DP:GQ 0/1:22:212 0/1:19:589 0/1:17:163
1 3177144 . GT G 466.00 PASS AC=3;AN=6;DP=68;SF=0,1,2 GT:DP:GQ 0/1:30:150 0/1:24:236 0/1:14:151
1 3184885 . TAAAA TA,T 422.00 PASS AC=2,1;AN=4;DP=26;SF=0,1 GT:DP:GQ 1/2:10:12 0/1:16:435 .
1 3199812 . G GT . . AC=2;AN=2;SF=2 GT . . 1/1
2 3188209 . GA G 162.00 . AC=1;AN=2;DP=15;SF=1 GT:DP:GQ . 0/1:15:162 .
2 3199812 . G GTT,GT 481.00 PASS AC=1,1;AN=2;DP=26;SF=0 GT:DP:GQ 1/2:26:322 . .
2 3212016 . CTT C 613.00 . AC=1;AN=2;DP=11;SF=2 GT:DP:GQ . . 0/1:11:41
3 3199812 . G GTT,GT 353.00 PASS AC=1,1;AN=2;DP=19;SF=1 GT:DP:GQ . 1/2:19:188 .
3 3199815 . G T,A 353.00 PASS AC=1,1;AN=4;DP=38;SF=1,2 GT:DP:GQ . 0/2:19:188 0/1:19:188
3 3212016 . CTT C,CT 565.00 PASS AC=1,1;AN=2;DP=26;SF=0 GT:DP:GQ 1/2:26:91 . .
3 3242491 . TT T . . AC=2;AN=2;SF=2 GT . . 1/1
4 3212016 . CTT C 677.00 q20 AC=1;AN=2;DP=15;SF=1f GT:DP:GQ . 0/1:15:158 .
4 3258448 . TACACACAC T 325.00 PASS AC=1;AN=2;DP=31;SF=0 GT:DP:GQ 0/1:31:325 . .
4 3291771 . T TAA,TAAA 336.00 . AC=1,1;AN=2;DP=12;SF=2 GT:DP:GQ . . 1/2:12:2
vcftools-0.1.16/examples/parse-test.vcf 0000664 0000000 0000000 00000000346 13331057401 0020043 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.1
#CHROM POS ID REF ALT QUAL FILTER INFO
1 100 . GTTT g 1806 PASS DP=35
1 104 . C g,,t, 1792 PASS DP=32
1 104 . C ,g,,t 1792 PASS DP=32
1 104 . C ,g,,t, 1792 PASS DP=32
vcftools-0.1.16/examples/perl-api-1.pl 0000664 0000000 0000000 00000002637 13331057401 0017465 0 ustar 00root root 0000000 0000000 #!/usr/bin/env perl
#
# Example code for generating a minimal VCF file using the perl API
#
# Author: pd3@sanger
#
use strict;
use warnings;
use Carp;
use Vcf;
my $sample = 'Sample1';
my $vcf_out = Vcf->new();
$vcf_out->add_columns($sample);
$vcf_out->add_header_line({key=>'FORMAT',ID=>'GT',Number=>'1',Type=>'String',Description=>"Genotype"});
$vcf_out->add_header_line({key=>'ALT',ID=>'DEL',Description=>"Deletion"});
$vcf_out->add_header_line({key=>'ALT',ID=>'DEL:ME:ALU',Description=>"Deletion of ALU element"});
$vcf_out->add_header_line({key=>'ALT',ID=>'DEL:ME:L1',Description=>"Deletion of L1 element"});
$vcf_out->add_header_line({key=>'ALT',ID=>'DUP',Description=>"Duplication"});
$vcf_out->add_header_line({key=>'INFO',ID=>'DP',Number=>1,Type=>'Integer',Description=>"Total Depth"});
$vcf_out->add_header_line({key=>'INFO',ID=>'H2',Number=>0,Type=>'Flag',Description=>"HapMap2 membership"});
print $vcf_out->format_header();
my $pos = 1;
for my $gt qw(A/A C/C /C / / /)
{
$pos++;
my %out;
$out{CHROM} = '1';
$out{POS} = $pos;
$out{ID} = '.';
$out{ALT} = [];
$out{REF} = 'C';
$out{QUAL} = '.';
$out{FILTER} = ['.'];
$out{INFO} = { DP=>3, H2=>undef };
$out{FORMAT} = ['GT'];
$out{gtypes}{$sample}{GT} = $gt;
$vcf_out->format_genotype_strings(\%out);
print $vcf_out->format_line(\%out);
}
vcftools-0.1.16/examples/query-test.out 0000664 0000000 0000000 00000000377 13331057401 0020133 0 ustar 00root root 0000000 0000000 1:100100 ref=G alt=C qual=0 1 A=G|C B=G/C
1:100200 ref=G alt=C qual=0 1 A=G|C B=G/G
1:100300 ref=G alt=C qual=0 1 A=C/C B=./.
1:100400 ref=C alt=G,T qual=35 1 A=G/G B=C/T
1:100500 ref=A alt=G qual=0 1 A=G/G B=A/A
1:100600 ref=C alt=G qual=0 1 A=G/G B=C/C
vcftools-0.1.16/examples/shuffle-test.vcf 0000664 0000000 0000000 00000001176 13331057401 0020367 0 ustar 00root root 0000000 0000000 ##fileformat=VCFv4.0
##INFO=