I have a Python script which takes as input a list of SNPs, then calls
tabix to retrieve .VCF file and does some stuff with it.
The problem is that i need to manually add (using UCSC genome browser web*) the coordinates of the SNP for tabix to work properly.
Does any one know any simple python API / function / module to easy add this information without having to perform the search manually?
Hi Peixe, you should accept one of these excellent answers! :)