Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • VarScan2. Strange behaviour of somatic-p-value

    Hi community of SEQanswers,

    I'm running varscan2 for somatic mutation calling following the commands of its documentation: mpileup + varscan somatic

    In order to test the effect of --somatic-p-value parameter in the number of variants called by varscan somatic, I have tried different p-values and logically I obtain different number of variants in the final output. However, checking the differents outputs I see that the value of the field somatic p-value (col. 15 in the tabulated output or value SPV of col. INFO in the VCF format) is greater than the specified threshold. I expected that all the somatic variants of the output has a somatic-p-value below my threshold.

    Is there any reason to explain this behaviour? Am I doing something wrong? I would like to know whether --somatic-p-value parameter is filtering variants according its value or not.

    Thank you in advance for you answers,
    Last edited by slois; 04-03-2014, 08:01 AM.

  • #2
    Hi,

    According the documentation:
    --p-value - P-value threshold to call a heterozygote [0.99]
    --somatic-p-value - P-value threshold to call a somatic site [0.05]

    Are these parameters filtering the variants identified according its corresponding values?

    Different tests applied using different somatic p-values: (outputs are modified to fit in this post). Last two columns are the variant-p-value and somatic-p-value:

    Results obtained for --somatic-p-value = 1
    result1.snp:chr1 1019753 A G 10 0 0% A 8 2 20% R Somatic 1.0 0.2368
    result1.snp:chr1 1111270 C T 36 2 5,26% C 52 17 24,64% Y Somatic 1.0 0.0090
    result1.snp:chr1 1111272 C T 35 2 5,41% C 55 16 22,54% Y Somatic 1.0 0.0182

    Results obtained for --somatic-p-value = 0.05 (Default)
    result2.snp:chr1 1019753 A G 10 0 0% A 8 2 20% R Somatic 1.0 0.2368
    result2.snp:chr1 1111270 C T 36 2 5,26% C 52 17 24,64% Y Somatic 1.0 0.0090
    result2.snp:chr1 1111272 C T 35 2 5,41% C 55 16 22,54% Y Somatic 1.0 0.0182

    Results obtained for --somatic-p-value = 0.01
    result3.snp:chr1 1019753 A G 10 0 0% A 8 2 20% R Somatic 1.0 0.2368
    result3.snp:chr1 1111270 C T 36 2 5,26% C 52 17 24,64% Y Somatic 1.0 0.0090

    Results obtained for --somatic-p-value = 0.001
    result4.snp:chr1 1019753 A G 10 0 0% A 8 2 20% R Somatic 1.0 0.2368

    Thank you

    Comment


    • #3
      I also have a same question,

      I have also ran the VarScan2 on my normal/tumor pair with p-value as 0.1 somatic p-value 0.001 and then ran the somatic filter with parameters -min-strands2 2 -min-avg-qual 30. I get very less variants that are significant. I get around 250 variants. Can anyone tell me what is the standard threshold for the p-value and somatic pvalue? Which is the paramter you have been supplying?

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Strategies for Sequencing Challenging Samples
        by seqadmin


        Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
        03-22-2024, 06:39 AM
      • seqadmin
        Techniques and Challenges in Conservation Genomics
        by seqadmin



        The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

        Avian Conservation
        Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
        03-08-2024, 10:41 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, 03-27-2024, 06:37 PM
      0 responses
      13 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-27-2024, 06:07 PM
      0 responses
      11 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-22-2024, 10:03 AM
      0 responses
      53 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-21-2024, 07:32 AM
      0 responses
      69 views
      0 likes
      Last Post seqadmin  
      Working...
      X