ó gbc@sEdZddlmZddlmZd„Zd„Zd„ZdS(s, Created on 23 nov. 2011 @author: merciece iÿÿÿÿ(tsequence(ttaxonomycCs%t|ƒ}t||ƒ}||fS(N(tecoPCRDatabaseReadertecoPCRFileReader(tentriestoptionstfilteredDataFromDBt filteredData((sAlib/python2.7/site-packages/obitools/barcodecoverage/readFiles.pytmain s cs¶tjˆjƒ}tjˆjd|ƒ}|jdƒ}|jdƒ}|jdƒ}|jdƒ}t|||gƒ}i‰i‰i}xä|D]Ü} d| kr”| d} tƒ} x´|j| ƒD] } | d|kró| j| dƒnt|| @ƒdkrÉx(| D] } ˆj | d ƒdˆ| |]4}tˆ|ƒttˆƒƒˆjkr|VqdS(N(tfloattlent rankthresold(t.0tr(tingroupRt usedrankid(sAlib/python2.7/site-packages/obitools/barcodecoverage/readFiles.pys 3sRtoutgrouptranks( Rt EcoTaxonomyDBRtEcoPCRDBSequenceIteratortfindRankByNametsettparentalTreeIteratortaddRtgett isAncestorRtid(Rttaxtseqstnorankidt speciesidtgenusidtfamilyidt minrankseqRtsR tallranktpRt keptranks((RRRsAlib/python2.7/site-packages/obitools/barcodecoverage/readFiles.pyRs8     " cCsSg}xF|D]>}d|kr |j}||dkrK|j|ƒqKq q W|S(NR R(Rtappend(RRRR&tseqId((sAlib/python2.7/site-packages/obitools/barcodecoverage/readFiles.pyR=s   N(t__doc__tobitools.ecopcrRRRRR(((sAlib/python2.7/site-packages/obitools/barcodecoverage/readFiles.pyts   +