Hi all,
I have RNA-seq paired-end data from Illumina and I am using BWA-mem for mapping and samtools for SNP calling. I wanted to understand and get your opinion regarding following.
1. What is the criteria for calling the variant? Specifically, at a particular position, how many reads should contain variation for it to be called as a variant.
I want to call variant at a position if it supported by more than 20% of reads at that position.
What parameter can I use to achieve this.
Any help is appreciated!!
I have RNA-seq paired-end data from Illumina and I am using BWA-mem for mapping and samtools for SNP calling. I wanted to understand and get your opinion regarding following.
1. What is the criteria for calling the variant? Specifically, at a particular position, how many reads should contain variation for it to be called as a variant.
I want to call variant at a position if it supported by more than 20% of reads at that position.
What parameter can I use to achieve this.
Any help is appreciated!!
Comment