I have a set of whole genome sequencing samples, some of which have matched normal while some not. I want to call somatic mutations on non-coding regions.

I was looking for GATK best practices (like this) but it seems it is not yet ready for production. Which alternative tools do you suggest?

  • $\begingroup$ Have you looked at this question? bioinformatics.stackexchange.com/questions/2744/… $\endgroup$
    – Bioathlete
    Dec 13, 2017 at 17:59
  • $\begingroup$ @Bioathlete thanks, I actually didn't find that question before. However, I found that question too specific on tumor-only samples. I also have normals so maybe there are better tools? $\endgroup$
    – gc5
    Dec 13, 2017 at 18:09

1 Answer 1


Those will have to be treated as separate analysis. You can use Mutect or VarScan for the matched samples and the tools listed in the other question for the unmatched samples however that means that you won't be able to compare the results easily. If the inter-sample comparison is important you will have to run all samples as unmatched. Personally I have used LoFreq for somatic variant detection with some success. It is very sensitive so the false positive rate can be high depending on how noisy your data is.


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.