Variant Calling for Organism aligned to a reference genome of not the same species
0
0
Entering edit mode
7 months ago

I aligned exon capture data from multiple individuals to a closely related reference genome of my study organism. I was wondering if there is a best pipeline for variant calling under these circumstances. I am looking at GATK, but am wondering if generating gvcf files for all of my samples would be poor practice because they are not the same species as the reference I am using. Any feedback would be great!

genome alignment variant calling • 226 views
ADD COMMENT
1
Entering edit mode

best pipeline for variant calling under these circumstances

None, you cannot distinguish between intra e inter-species variants.

Maybe it is better to try to assemble your sequenced regions, create a consensus sequence and do variant calling in such areas

ADD REPLY

Login before adding your answer.

Traffic: 1454 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6