text: user4790: I found a good answer related to my question here:
  How to get a list of genes corresponding to the list of SNPs (rs ids)?
  But it says about small number of SNPs. I want to get gene names based on approximately 12000 to 19000 SNPs. NCBI says about 10 read per second with API key. Is there any efficient way to get gene names for thousands of SNPs? Is that database available to store locally? Any advice please?
  
  user776: Modify SNPs and gene annotations for your genome assembly of interest. This example shows how one might do this for hg19, for instance:
  
  Download SNPs and convert to BED:
  $ wget -qO- ftp://ftp.ncbi.nih.gov/snp/organisms/human_9606_b151_GRCh37p13/VCF/common_all_20180423.vcf.gz \
      | gunzip -c \
      | convert2bed --input=vcf --output=bed --sort-tmpdir=${PWD} - \
      > hg19.snp151.bed
  
  Filter this BED file for entries of interest:
  $ grep -wFf snps-of-interest.txt hg19.snp151.bed | cut -f1-6 > hg19.snp151.filtered.bed6
  
  Download gene annotations and convert to BED:
  $ wget -O - ftp://ftp.ebi.ac.uk/pub/databases/gencode/Gencode_human/release_19/gencode.v19.annotation.gtf.gz \
      | gunzip -c \
      | awk '($3 == "gene")' \
      | gtf2bed \
      > gencode.v19.genes.bed
  
  Note: In the future, newer versions of Gencode annotations may be available for your reference assembly. Always check their site to confirm what is most recent.
  Map genes to SNPs:
  $ bedmap --echo --echo-map-id hg19.snp151.filtered.bed6 gencode.v19.genes.bed > snps_with_associated_gene_names.bed
  
  
  A list of genes that overlap a SNP of interest will be in the seventh column of the output.
  If this must be done within Python, you could use something like subprocess.call(...), where ... includes commands wrapped in bash -c. 
  I don't recommend this because of the need to escape quote marks and deal with Python's quirks in running command-line tasks, which makes a Python-based approach very difficult to set up and maintain. I'd suggest learning some basic bash scripting to make this easy and fast.
  
original_text: 4790: <p>I found a good answer related to my question here:
  <a href="https://bioinformatics.stackexchange.com/questions/2503/how-to-get-a-list-of-genes-corresponding-to-the-list-of-snps-rs-ids">How to get a list of genes corresponding to the list of SNPs (rs ids)?</a></p>
  
  <p>But it says about small number of SNPs. I want to get gene names based on approximately 12000 to 19000 SNPs. NCBI says about 10 read per second with API key. Is there any efficient way to get gene names for thousands of SNPs? Is that database available to store locally? Any advice please?</p>
   776: <p>Modify SNPs and gene annotations for your genome assembly of interest. This example shows how one might do this for <code>hg19</code>, for instance:</p>
  
  <ol>
  <li><p>Download SNPs and convert to BED:</p>
  
  <pre><code><span class="math-container">$ wget -qO- ftp://ftp.ncbi.nih.gov/snp/organisms/human_9606_b151_GRCh37p13/VCF/common_all_20180423.vcf.gz \
      | gunzip -c \
      | convert2bed --input=vcf --output=bed --sort-tmpdir=$</span>{PWD} - \
      &gt; hg19.snp151.bed
  </code></pre></li>
  <li><p>Filter this BED file for entries of interest:</p>
  
  <pre><code>$ grep -wFf snps-of-interest.txt hg19.snp151.bed | cut -f1-6 &gt; hg19.snp151.filtered.bed6
  </code></pre></li>
  <li><p>Download gene annotations and convert to BED:</p>
  
  <pre><code><span class="math-container">$ wget -O - ftp://ftp.ebi.ac.uk/pub/databases/gencode/Gencode_human/release_19/gencode.v19.annotation.gtf.gz \
      | gunzip -c \
      | awk '($</span>3 == "gene")' \
      | gtf2bed \
      &gt; gencode.v19.genes.bed
  </code></pre>
  
  <p><strong><em>Note</strong>: In the future, newer versions of Gencode annotations may be available for your reference assembly. Always check their site to confirm what is most recent.</em></p></li>
  <li><p>Map genes to SNPs:</p>
  
  <pre><code>$ bedmap --echo --echo-map-id hg19.snp151.filtered.bed6 gencode.v19.genes.bed &gt; snps_with_associated_gene_names.bed
  </code></pre></li>
  </ol>
  
  <p>A list of genes that overlap a SNP of interest will be in the seventh column of the output.</p>
  
  <p>If this must be done within Python, you could use something like <code>subprocess.call(...)</code>, where <code>...</code> includes commands wrapped in <code>bash -c</code>. </p>
  
  <p>I don't recommend this because of the need to escape quote marks and deal with Python's quirks in running command-line tasks, which makes a Python-based approach very difficult to set up and maintain. I'd suggest learning some basic <code>bash</code> scripting to make this easy and fast.</p>
  