Entering edit mode
7.7 years ago
bioguy24
▴
230
Is there a way to query the usuc mysql database using a file and return the fasta sequence of each line in the file and get the same output as in the table browser? This is DNA, CDS exons, one fasta record per region with 10 bases 5' and 3', exons in upper case, everything else in lower case. Thank you :).
file
A2M
output
>hg19_refGene_NM_000014_0 range=chr12:9268350-9268455 5'pad=10 3'pad=10 strand=- repeatMasking=none
tttctgcaacATGGGGAAGAACAAACTCCTTCATCCAAGTCTGGTTCTTC
TCCTCTTGGTCCTCCTGCCCACAGACGCCTCAGTCTCTGGAAAACCgtga
gttcca
>hg19_refGene_NM_000014_1 range=chr12:9265946-9266149 5'pad=10 3'pad=10 strand=- repeatMasking=none
tcattcccagGCAGTATATGGTTCTGGTCCCCTCCCTGCTCCACACTGAG
ACCACTGAGAAGGGCTGTGTCCTTCTGAGCTACCTGAATGAGACAGTGAC
TGTAAGTGCTTCCTTGGAGTCTGTCAGGGGAAACAGGAGCCTCTTCACTG
ACCTGGAGGCGGAGAATGACGTACTCCACTGTGTCGCCTTCGCTgtgagt
gtgg
Maybe:
mysql --user=genome --host=genome-mysql.cse.ucsc.edu -A -D hg37 -e 'select /home/user/file()'
using mysql only ? no
What do you suggest? Thanks.