If you happen to see a question you know the answer to, please do chime in and help your fellow community members. We encourage our fourm members to be more involved, jump in and help out your fellow researchers with their questions. GATK forum is a community forum and helping each other with using GATK tools and research is the cornerstone of our success as a genomics research community.We appreciate your help!
Test-drive the GATK tools and Best Practices pipelines on Terra
Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.
We will be out of the office on November 11th and 13th 2019, due to the U.S. holiday(Veteran's day) and due to a team event(Nov 13th). We will return to monitoring the GATK forum on November 12th and 14th respectively. Thank you for your patience.
Comparison of pairs of tumor samples
Dear GATK, I would like to analyse data from pairs of samples. Each pair is a sample from a tumor and from a cell line derived from that tumor. We expect the tumor samples to be quite heterogenous, and the derived cell lines be less heterogenous. I want to find what differs between the tumor and the derived cell line. In general we think most of the variants present in the derived cell lines will also present in the original tumor (although possibly at a different frequency). Of primary interest would be variants present in tumor but absent ( or present at a reduced frequency ) in the cell lines.
I thought I could either:
1) Run mutect2 in matched pair mode, using the cell lines as the normal samples. I would have to set ploidy suitably high to account for low frequency variants. I think I could omit the PON as I don't think it is very important whether or not the the variant is likely to be germline, and anyway I don't have normal samples. It is possible that the "normal" sample in this scenario will contain large deletions not found in in the "tumour" samples, will this cause problems?
2) Run mutect2 in tumor only mode, then compare the output using another tool such as "GenotypeConcordance"
Do either of these options seem sensible? Also is it possible to get mutect2 to output all sites (not just variant sites).
Thank you and sorry for the long question