- VariantFiltration
Filter variant calls based on INFO and/or FORMAT annotaitions.
gatk VariantFiltration \
-R reference.fasta\
-V input.vcf.gz\
-O output.vcf.gz\
--fitler-expression "AB<0.2 || MQ0 >50" \
--fitler-name "my_filters"
不能使用VQSR，但需要对HaplotypeCaller的结果进行handfilter时可以用这个工具。
- FixMateInformation(Picard)
Verify mate-pair information between mates and fix if needed.
java -jar --Xmx4g picard.jar FixMateInformation \
I=input.bam \
O=fixed_mate.bam \
ADD_MATE_CIGAR=true AS=true \
SO=coordinate
这一步是放在MarkDuplicates后面，查到有别人说如果MarkDuplicates 把重复去掉了,会对mate信息产生影响。如果MarkDuplicates只是标记而没有去除重复，则不会对mate信息产生影响，理论上可以不用做FixMate.
我比较了输入文件和输出文件的大小，fix之后的文件要大一些。
3.SelectVariants
Select a subset of variants from a VCF file.
gatk SelectVariants \
-R reference.fasta \
-V input.vcf \
-selectType SNP \
-O output.vcf
-selectType ：INDEL,SNP,MIXED,MNP,SYMBOLIC,NO_VARIATION.（can be specified multiple times.)
常用的是从vcf中把snp和indel分别输出到两个文件。
