Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • varscan2: copynumber detection

    Hi all,

    I would like to know if anybody has experience using VarScan2 to call copynumber (assuming we have Tumor/Normal matched sample). Specifically, I want to know how you calculate the parameter "data ratio".

    Any suggestion will be appreciated.

  • #2
    Thank you for your message. We developed VarScan and use it for exome-based copy number analysis. The data-ratio value could simply be calculated by looking at the SAMtools "flagstat" for tumor and normal BAM files to calculate the number of reads that are mapped and unique (not marked duplicate). The ratio should be normal OVER tumor:

    r = (normal_unique_reads / tumor_unique_reads)

    Of course, this assumes that reads are the same length in normal and tumor datasets. If they're not, one would have to adjust for uniquely mapped base pairs.

    You may notice, after calling copy number and performing CBS, some samples with a systematic bias in copy number calls where neutral regions seem shifted away from the zero-axis. We believe these to be instances of different exome capture efficiency between samples.

    You can normalize raw VarScan copy number values up or down to address these situations; see the VarScan copyCaller command usage for details.

    Please feel free to e-mail me directly if you have further questions:
    dkoboldt (at) genome [dot] wustl (dot) edu.

    Comment


    • #3
      I have a question regarding the mergeSegments.pl that is a part of the copy number workflow from Vrascan.

      After a lot of searching, I generated the file with the following columns (Please note, nowhere, on the website is it mentioned that we need the segments.p function from the DNAcopy package to get these columns ):
      ID, chrom, loc.start, loc.end, num.mark, seg.mean, bstat, pval, lcl, ucl.

      and when I try to run this file with mergeSegments.pl ( perl ../../mergeSegments.pl out.file --output-basename basenametest --verbose 1) it gives me an error:

      Use of uninitialized value $input in <HANDLE> at ../../mergeSegments.pl line 446.
      readline() on unopened filehandle at ../../mergeSegments.pl line 446.
      Can't use an undefined value as a symbol reference at ../../mergeSegments.pl line 456.

      What could I change here?

      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, Yesterday, 06:37 PM
      0 responses
      8 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, Yesterday, 06:07 PM
      0 responses
      8 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-22-2024, 10:03 AM
      0 responses
      49 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-21-2024, 07:32 AM
      0 responses
      66 views
      0 likes
      Last Post seqadmin  
      Working...
      X