calculate ti/tv ratio in WGS data
2
0
Entering edit mode
3.1 years ago
chengyi31000 ▴ 10

I am dealing with WGS data with vcf files. Wondering is it necessary to calculate ti/tv ratio to do QC? Can I calculate ti/tv ratio using vcftools? I have only see the function in Variants tools, but it is not that convenience to upload the big data.

sequencing ti/tv ratio • 3.1k views
ADD COMMENT
2
Entering edit mode

ts/tv ratio is important for judging the overall quality of your calls. It should not be too different from known values for your species. I usually employ it as one of the gold standards to reach while filtering.

ADD REPLY
0
Entering edit mode

how can i determine the standards of ts/tv ratio? How can I filtering?

ADD REPLY
1
Entering edit mode

The expected ratio of transitions to transversions in the whole genome is near 2.0 for humans.

ADD REPLY
0
Entering edit mode

Is it meaningful to calculate ti/tv ratio by each chromosome ? or it must calculate ti/tv ratio by all chromosome

ADD REPLY
1
Entering edit mode
3.1 years ago
trausch ★ 1.6k
bcftools stats snv.vcf.gz | grep "TSTV"
ADD COMMENT
0
Entering edit mode
13 days ago
davidc • 0

snpSift tstv my.vcf.gz

gives individual tstv values for each samples in a multi sample vcf and takes into account the genotype when calculating

I believe bvcftools counts the chr/pos/ref/alt but if the genotype is 0/0 still counts it.

ADD COMMENT

Login before adding your answer.

Traffic: 2428 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6