Varscan Filtering scripts: very little overlap between output of varscan filtering scripts
0
1
Entering edit mode
9.2 years ago
Chirag Nepal ★ 2.4k

Hi there,

I am using varscan2 (for matched N/T pairs) to call somatic SNVs. For filtering, I used script from varscan to filter high confidence somatic call. I used two scripts from varscan, and I get completely different results with very little overlap.

Filter script-1

java -jar VarScan.v2.3.7.jar processSomatic $name

This will be two output files, namely $name.Somatic and $name.Somatic.hc (high confidence output)

Filter script-2

java -jar VarScan.jar somaticFilter $name.Somatic --min-coverage 15 --min-reads2 5 --min-strands2 2 --p-value 0.001 --output-file $name.Somatic.hc"Criteria2"

Overlap between two high-confidence somatic calls,

overlap is about 10%. Why two filtering criteria are so discordant.

Any insight will be helpful.

cheers
Chirag

snv exome varscan filter-varscan • 3.1k views
ADD COMMENT

Login before adding your answer.

Traffic: 2013 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