Hi all,
I found the genotype is also available for the sites which without any reads
passed my quality filter criterion. for example, the genotype is 1/1 but
without any reads 1/1:0,0,0:0:0:0:0,0,0,0:0,0. So why does samtools report
results like that, whether these sites should be marked as ./. .
My command is as follow (version 1.31)
samtools mpileup -q 20 -t DP,DPR,DV,DP4,INFO/DPR,SP -g -f ref.fa -b bam.list |
bcftools call -V indels -v -c -O z -o raw.vcf.gz
Can you give me some suggestions on how to deal with these sites?
best
Zheng zhuqing
------------------------------------------------------------------------------
What NetFlow Analyzer can do for you? Monitors network bandwidth and traffic
patterns at an interface-level. Reveals which users, apps, and protocols are
consuming the most bandwidth. Provides multi-vendor support for NetFlow,
J-Flow, sFlow and other flows. Make informed decisions using capacity planning
reports.http://sdm.link/zohodev2dev
_______________________________________________
Samtools-help mailing list
Samtools-help@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/samtools-help
