BLAT

Perform a BLAT search using the UCSC Genome Browser.

Arguments

NameShortLongDescription
Sequence Type-s--seqtypeSpecify the structure format [default: dna] [possible values: dna, protein, translated-rna, translated-dna]
Database Name-d--db-nameSpecifies the database name to query [default: hg38]
Output-o--outputoptional filepath to write output to [default=stdout]

Command Line Interface

# query UCSC genome browser for the first 121 bp of AP2S1
ggetrs ucsc blat GGGCCCTACAACTGCACCCTGAGCCGGAGCTGCCCAGTCGCCGCGGGACCGGGGCCGCTGGGGTCTGGACGGGGGTCGCCATGGTAACGGGGGAGCGCTACGCCGGGGACTGGCGGAGGG

Python

import ggetrs

# query UCSC genome browser for the first 121 bp of AP2S1
ggetrs.ucsc.blat(
  "GGGCCCTACAACTGCACCCTGAGCCGGAGCTGCCCAGTCGCCGCGGGACCGGGGCCGCTGGGGTCTGGACGGGGGTCGCCATGGTAACGGGGGAGCGCTACGCCGGGGACTGGCGGAGGG"
)

# query UCSC genome browser with amino acid sequence
ggetrs.ucsc.blat(
  "MIRFILIQNRAGKTRLAKWYMQFDDDEKQKLIEEVHAVVTVRDAKHTNFVEFRNFKIIYRRYAGLYFCICVDVNDNNLAYLEAIHNFVEVLNEYFHNVCELDLVFNFYKVYTVVDEMFLAGEIRETSQTKVLKQLLMLQSLE",
  seqtype="protein"
)