python - python，按坐标从DAS提取序列

ucsc DAS server，通过坐标获得DNA序列。

网址：http://genome.ucsc.edu/cgi-bin/das/hg19/dna?segment=chr20:30037432,30038060

样本文件：

<DASDNA>
  <SEQUENCE id="chr20" start="30037832" stop="30038060" version="1.00">
  <DNA length="229">
    gtggcacccaaagatgctggaatctttatggcaaatgccgttacagatgc
    tccaagaaggaaagagtctatgtttactgcataaataataaaatgtgctg
    cgtgaagcccaagtaccagccaaaagaaaggtggtggccattttaactgc
    tttgaagcctgaagccatgaaaatgcagatgaagctcccagtggattccc
    acactctatcaataaacacctctggctga
  </DNA>
  </SEQUENCE>
</DASDNA>

我想要的是这部分：

    gtggcacccaaagatgctggaatctttatggcaaatgccgttacagatgc
    tccaagaaggaaagagtctatgtttactgcataaataataaaatgtgctg
    cgtgaagcccaagtaccagccaaaagaaaggtggtggccattttaactgc
    tttgaagcctgaagccatgaaaatgcagatgaagctcccagtggattccc
    acactctatcaataaacacctctggctga

我想从成千上万个此类网址中获取序列部分，该怎么办？
我试图将数据写入文件并解析文件，它工作正常，但是有什么方法可以直接解析类似xml的字符串吗？我尝试了其他帖子中的一些示例，但是它们没有用。

在这里，我添加了解决方案。感谢下面的2个答案。

解决方案1：

def getSequence2(chromosome, start, end):
    base = 'http://genome.ucsc.edu/cgi-bin/das/hg19/dna?segment='
    url = base + chromosome + ':' + str(start) + ',' + str(end)
    doc = etree.parse(url,parser=etree.XMLParser())
    if doc != '':
        sequence = doc.xpath('SEQUENCE/DNA/text()')[0].replace('\n','')
    else:
        sequence = 'THE SEQUENCE DOES NOT EXIST FOR GIVEN COORDINATES'
    return sequence

解决方案2：

def getSequence1(chromosome, start, end):
    base = 'http://genome.ucsc.edu/cgi-bin/das/hg19/dna?segment='
    url = base + chromosome + ':' + str(start) + ',' + str(end)
    xml = urllib2.urlopen(url).read()
    if xml != '':
        w = open('temp.xml', 'w')
        w.write(xml)
        w.close()
        dom = parse('temp.xml')
        data = dom.getElementsByTagName('DNA')
        sequence = data[0].firstChild.nodeValue.replace('\n','')
    else:
        sequence = 'THE SEQUENCE DOES NOT EXIST FOR GIVEN COORDINATES'
    return sequence

当然，他们将需要导入一些必要的库。

最佳答案

>>> from lxml import etree
>>> doc = etree.parse("http://genome.ucsc.edu/cgi-bin/das/hg19/dna?segment=chr20:30037432,30038060",parser=etree.XMLParser())
>>> doc.xpath('SEQUENCE/DNA/text()')
['\natagtggcacatgtctgttgtcctagctcctcggggaaactcaggtggga\ngagtcccttgaactgggaggaggaggtttgcagtgagccagaatcattcc\nactgtactccagcctaggtgacagagcaagactcatctcaaaaaaaaaaa\naaaaaaaaaaaaaagacaatccgcacacataaaggctttattcagctgat\ngtaccaaggtcactctctcagtcaaaggtgggaagcaaaaaaacagagta\naaggaaaaacagtgatagatgaaaagagtcaaaggcaagggaaacaaggg\naccttctatctcatctgtttccattcttttacagacctttcaaatccgga\ngcctacttgttaggactgatactgtctcccttctttctgctttgtgtcag\ngtggcacccaaagatgctggaatctttatggcaaatgccgttacagatgc\ntccaagaaggaaagagtctatgtttactgcataaataataaaatgtgctg\ncgtgaagcccaagtaccagccaaaagaaaggtggtggccattttaactgc\ntttgaagcctgaagccatgaaaatgcagatgaagctcccagtggattccc\nacactctatcaataaacacctctggctga\n']

关于python - python，按坐标从DAS提取序列，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/22328319/