I have the CDS sequences for my genome. I want to align Illumina DNA seq reads from related fungal strains to these. This is to identify SNPs in just the coding sequences. I am not sure if I should use a DNA seq aligner like BWA and Bowtie, or to use a split aware aligner like STAR.
I have the full genome. But I wanted the SNPs for each of the CDS sequences, so I can create consensus sequences for these CDS for the related strains.
Can you provide examples or links to get alignment from CDS region or using VCFtools to get SNPs and INDELs from the CDS region? I have tried looking but haven't found anything useful. Thank you.