This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Do You Guys Have Any Recommendation On Variant Calling Using Long Contigs Or Scaffold

I have assembled a human genome and used bwa mem aligning the contigs back to the reference genome. Do you guys have any recommendation of variant caller?

variant-calling snp sv indel

2 answers

GATK
samtools
VarScan2
SNVmix

I wonder if the regular tools would be appropriate in this scenario since these are not raw reads mapped to reference but assembled contigs and as such the error models would likely be inaccurate.

Not to mention some of the variants here might actually be mis-assembly events.

For the SV/CNV calling, I've used CLEVER and CNVnator.

Log in to answer this question.