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
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
Comment