In hg38 genome builds contig name for mitochondria is named chrM/M whereas in hg19 it's named chrMT/MT.
Currently scout does not support the M notation and as such we cannot filter against the contig.
Would it be possible to support both notations?
I wasn't aware of that! Of course it should be fixed!
Seems we are also using the CHROMSOME ones (i.e. MT) for fetching from biomart. There are not that many genes on M/MT, but are they properly fetched for hg38?
@dnil It appears that the genes are loaded properly from biomart. Ensembl is still calling the chromosome MT.
I'm quite confused by this. Maybe we're using some strange version of the genome build (but I don't think so)?
@ViktorHy Would renaming M to MT in our fasta solve all these issues or are there other things to consider?
I think VEP feelt it would be to challenging to change this in all their databases. However, you can use the --synonyms to make VEP use a map of chromosome names instead. This is what we do in MIP.
Most helpful comment
I think VEP feelt it would be to challenging to change this in all their databases. However, you can use the
--synonymsto make VEP use a map of chromosome names instead. This is what we do in MIP.