Mapping of CDS sequence of genes using minimap2

Question

[harisanjanas]

Hi, I have downloaded the CDS sequence of some Bos taurus genes from NCBI in fasta format. My goal is to view the different exons in IGV after mapping it to the reference genome. The problem now is that minimap2 (even after using the --splice option) is only mapping the first exon sequence. Does anyone know a way to map these reads which includes the splice junctions and hence all the exons in the gene. ADD INFO: when mapping and viewing the complete sequence it is mapping fine but the splice junctions are not properly marked hence it becomes hard to find the exact position of the exons.

Answer 1

[lieven.sterck]

I suggest to have a look at the more 'old-school' mappers: GenomeThreader, GMAP, EST2genome, ... they are definitely slower than minimap but are very capable to identify (or correctly splice-map) CDS/cDNA sequences

[harisanjanas]

Thank you very much. Will do.