View Single Post
Old 01-26-2017, 09:56 AM   #1
clarissaboschi
Member
 
Location: US

Join Date: Apr 2010
Posts: 63
Default CrossMap analysis

I am trying to convert some genome coordinates of SNPs in one vcf file using the CrossMap program (http://crossmap.sourceforge.net/).

I am having difficulties with the conversion. My vcf file does not have all chromosomes present in the chain file (for example chromosome random, etc), so the conversion is not performed.

Do I need to remove all chromosomes not present in my vcf file from the chain file? (and maybe from the fasta file as well?)

My error message was:
KeyError: "sequence 'chr10_NT_461738v1_random' not present"

This sequence is present only in the chain file, but I am not sure if I should edit the chain file.

Also, Is the chromosome format need to be the same in the 3 files (input, chain and fasta file): Chr1, chr1 or 1?

thanks
clarissaboschi is offline   Reply With Quote