vcf file summary

Question

[rheab1230]

Hello everyone, I have vcf files for chr 1 to 22. I want to do a summary statistics for these vcf files. Like find out the sample size, no. of SNP,mean and standard deviation. Is there any software or tools which can plot graph for these calculations? Thank You.

[Ram]

You should mention that you've tried the vcfstats software and are facing issues; and also point to your post: error related to vcfstats

[rheab1230]

Yes, I tried vcfstats but its showing error. So now I am trying bcftools stats and then plot-vcfstats to analyse the output produce by bcftools stats. But its showing error the error:

Parsing bcftools stats output: test.vchk
Plotting graphs: python plot.py
Traceback (most recent call last):
  File "plot.py", line 53, in <module>
    import matplotlib as mpl
ModuleNotFoundError: No module named 'matplotlib'
The command exited with non-zero status 256:
        python plot.py

 at /home/kxj190026/anaconda3/envs/chip-seq/bin/plot-vcfstats line 99.
        main::error("The command exited with non-zero status 256:\x{a}\x{9}python plot.py\x{a}\x{a}") called at /home/anaconda3/envs/chip-seq/bin/plot-vcfstats line 287
        main::plot(HASH(0x55a9a5d69e80)) called at /home/anaconda3/envs/chip-seq/bin/plot-vcfstats line 69

[GenoMax]

ModuleNotFoundError: No module named 'matplotlib'

Install matplotlib: https://matplotlib.org/stable/users/installing.html

[rheab1230]

Thank you so much. I am able to do it now. I installed matplotlib

[rheab1230]

I got this as summary:

[rheab1230]

I am not really able to understand how to get information from this summary? Like what does n or ts/tv means. Like what does this summary overall predict?

[DBScan]

ts/tv is the ratio between transition/transversion, see for instance here: Ti Tv ratio and their usefulness in exome sequencing n is the number of SNPs (2893523) or indels (109293) I think.