Just started learning about NGS, and using samtools to understand the data provided.
Please let me know how can i write a perl script to call somatic variants from a sample and my final output file contains only somatic calls.
•
written
by
shilpy •
1 answer
is not possible to detect somatic variants from an unique list only, because some comparison pattern must be established in order to remove the germline background. a good idea to find somatic variants out is to use trios, where you could easily set filters among them (in perl if you are fond of it, although simple bash scripting would do) to find which variants of the child were not inherited from the parents. have you consider using existing tools to help you on this task, such as GATK's PhaseByTransmission?
Log in to answer this question.
-
How to find nu. of variants present in a gene for my Exome seq. sample
written by bioinf_data-analyst9449 •
-
Java Samtools Api: Write Bam To Stdout
written by dariober
-
Variant Calling With Bcftools Call
written by tayebwajb
-
Somaticsniper And Filtering With Included Perl Scripts
written by DoubleD
-
Samtools Filter Reads Cigar Field
written by dfernan
-
Are My Samtools Commands Correct
written by shilpy •
-
Haloplex & Allele Calling
written by Pierre Lindenbaum
-
Reference Genome Path In Samtools
written by win
-
How To Predict Which Next Gen Variants Would Not Validate By Sanger
written by Biomed
-
Samtools + Picard Markduplicates
written by Nery •
Hi, you should provide a little more information if you want us to help you. Can you paste an example of the input and of the output that you want to get?
Moreover, please note that some people may get annoyed if you ask for help, without demonstrating that you have at least made an effort. Why don't you post at least some code that you wrote for this task? Even if the code is not working, you should show that you have at least tried to do something, before asking for help here.
Writing a variant caller in perl is a BAD IDEA.