Hi,
I have two samples mapped to reference genome and would like to detect difference between two samples (not SNPs vs. the reference sequences).
I'm trying to use the "bcftools -bvcgT pair sample1.bvf sample2.bvf -> s1_vs_s2.bcf" option to do it. Is the "-T pair" the right option to use?
Also, the version I have does not have this option. Does it mean that I need a newer version or this option is not valid anymore?
If -T is not the right one, is there a way to do so using samtools? Or I will have to write some scripts by myself?
Thanks!
I have two samples mapped to reference genome and would like to detect difference between two samples (not SNPs vs. the reference sequences).
I'm trying to use the "bcftools -bvcgT pair sample1.bvf sample2.bvf -> s1_vs_s2.bcf" option to do it. Is the "-T pair" the right option to use?
Also, the version I have does not have this option. Does it mean that I need a newer version or this option is not valid anymore?
If -T is not the right one, is there a way to do so using samtools? Or I will have to write some scripts by myself?
Thanks!