bVf3ddlZddlZddlZddlZddlZddlZGddZGddZGddZ Gdd Z y) Nc^eZdZdZdZdZdZdZdZdZ dZ d Z d Z d Z d Zd ZdZy) TestRemoteNhttp://raw.githubusercontent.com/dpryan79/pyBigWig/master/pyBigWigTest/test.bwcLtj|j}|J|S)N)pyBigWigopenfnameselfbws Z/var/www/html/software/conda/envs/catlas/lib/python3.12/site-packages/pyBigWigTest/test.pydoOpenzTestRemote.doOpen s# ]]4:: &~~ ctj|j5}|jdddk(sJ dddy#1swYyxYw)Nc Q?110)rrr chromsr s r doOpenWithzTestRemote.doOpenWithsF ]]4:: & E"99; "CC DC E E Es AA c|jdddk(sJ|jddk(sJ|jdJy)Nrrrrc)rr s r doChromszTestRemote.doChromssEyy{IY??@?yy~*+*yy~%&%rc >|jddddddddk(sJy) Nir)maxValsumDataminValversion sumSquarednLevels nBasesCovered)headerr s r doHeaderzTestRemote.doHeaders9yy{saTUehuvJMNN ONrcZ|jddddgk(sJ|jdddddgk(sJ|jddd dd d d gk(sJ|jdtjdtjd dd d d gk(sJ|jddgk(sJy)NrrgUUU?max)type@33?crr-nBins`ff??gq\?)statsnpint64r s r doStatszTestRemote.doStatssxxQ"'9&::;:xxQx.3F2GGHGxxBsax8=NPS|jdddddddgk(sJy) NrrBrDsumrr1gA@g@B@)r5r s r doSumzTestRemote.doSum-s'xxS#Ex;d|KLKrctjd}|j}|jt j |d}|Jd|j dfd|j dfg}t|j dk(sJ|j|d|D]}|j|d }g}g} g} g} |D]R} |j|d | j| d | j| d| j| dT|j|| | | |jt j |}|j|jk(sJ|j |j k(sJ|D]1}|j|d } |j|d }| |k(r1J|j|jtj|y) NFdeletewrrrrmaxZoomsrendsr<)tempfileNamedTemporaryFilenamecloserrrlen addHeaderrEappend addEntriesr(osremove)r r ofileonamebw2rrintschroms2startsrQr<entryints1ints2s r doWritezTestRemote.doWrite0s++59  mmE3' #'$ $)@A299;1$%$ fq ) FA<<!%DGFDF (qt$ eAh' E!H% eAh'  ( NN7FfN E F mmE"yy{cjjl*+*yy{cjjl*+* #ALL1&EMM!A$'EE> "> #    %rcdgdz}gd}gd}gd}tjd}|j}|jt j |d}|j d d g|j|d d |d d |d d |d d  |j|d d|d d|d d|d d |jd|dd|ddd|jd|dd|ddd|jdd|dddd|jdd|dddd |j|d d |d d |d d |d d  dd k(sJ |jd|dd|ddddd k(sJ |jd|dd|ddddd k(sJ |jdd|dddddd k(sJ |jdgd z|d d |d d |d d  |jtjt |djj}|dk(sJtj|y#t$rYwxYw#t$rYwxYw#t$rYwxYw#t$rYwxYw)z Test all three modes of storing entries. Also test to ensure that we get error messages when doing something silly This is a modified version of the writing example from libBigWig r) rrB}r0riXiqii iR)x~)?i@b@9@rqrrrsrtrurv4rvrwrxrvFrKrMr@B2`rr+rP )r<span i)r<rstepir3r|rb ef104f198c6ce8310acc149d0377fc16N)rRrSrTrUrrrWrY RuntimeErrorhashlibmd5read hexdigestrZr[) r rrarQr<r\r]r hs r doWrite2zTestRemote.doWrite2Us qK+++59  ]]5# & nn56 fQqk6!A;T!AYvaPQ{ S fQqk6!A;T!AYvaPQ{ S c6!A;va{ D c6!B<q 2 F c3vb}2B G c3vb}2B G  MM&1+va{a6RSTU;M Wa4L4  MM#va{6!A;RM Ha4L4  MM#va{6!A;RM Ha4L4  MM#s6"R=rM Ka4L4 seAgva{a6!A; O  KKUD)..0 1 ; ; =4454 %3               sH(H6"I'"I I#6 II II I I # I/.I/ctjd}|j}|jt j |d}|j ddg|jtjt |djj}|dk(sJt j |}|jdd d k(sJ|jd dk(sJ|jd d dgk(sJ|jd d dd ddgk(sJ|jtj |y)NFrKrMryr{r 361c600e5badf0b45d819552a7822937rzr})rr|rrr)r2)rRrSrTrUrrrWrrrrrrEr<r5rZr[)r r\r]r rs r doWriteEmptyzTestRemote.doWriteEmptys++59  ]]5# & nn56   KKUD)..0 1 ; ; =4454]]5 !yy{G'::;:||C D()(yya)R/0/xxQqx1dD\ABA   %rctjd}|j}|jt j |d}|j gdtjdgdzdgdzzdgdzzd gdzz}tjgd tj }tjgd tj }tjtjjd tj }|j|||||jt|||D cgc] \}} } || | f } } }} t j |}|jdddddk(sJt!gdD]q\} } t!|j#| D]O\}}|d|d| z|zk(sJ|d|d| z|zk(sJtj$|d|d| z|zrOJs|jt'j(|ycc} } }w)NFrKrM))chr1rB)chr2rC)chr3r0)chr4rrrrr)r (2<FPZdtype)rkr-7AKU_rPrBrCr0r)rrrrrr)rRrSrTrUrrrWr6arrayr7random random_samplefloat64rYzipr enumeraterEiscloserZr[)r r\r]r rrarQvalues0xyzvalsidx1chromidx2tups r doWriteNumpyzTestRemote.doWriteNumpys++59  ]]5# & QR6(Q,&A51 DxRS|ST9Jxx7rxxH((2992215RZZH ff4 @  ),VT7)CDDgaAAq DD ]]5 !yy{sCcRRSR$%EF EKD%&r||E':; E c1vD4!88981va$ho!6676zz#a&'!d(T/*BCDC E E   %EsH:c|j}|j||jjds|j ||j ||j ||j||j||j||j|j|j|j|jy)Nhttp)rrr startswithr)r8rIr=rFrerrrrrUr s r testAllzTestRemote.testAlls [[] bzz$$V, MM"  LL  JJrN MM"    R LL  OO  MMO          rN)__name__ __module__ __qualname__r rrrr)r8r=rFrIrerrrrrr rrsO \E E' O6| m M#J9v.4rrceZdZdZy) TestLocalct}tjjtj dz|_|jy)Nz/pyBigWigTest/test.bw)rrZpathdirnamer__file__r r)r blahs r testFoozTestLocal.testFoos2|WW__X%6%67:QQ  rN)rrrrrrr rrsrrceZdZdZy) TestBigBedc:tjjtjdz}tj |}|J|j dk(sJ|jdk(sJd}|j}t|tr|jd}||k(sJ|jddd}gd }||k(sJ|jdtjdtjd}||k(sJ|jy) Nz/pyBigWigTest/test.bigBedrratable RnaElements "BED6 + 3 scores for RNA Elements data " ( string chrom; "Reference sequence chromosome or scaffold" uint chromStart; "Start position in chromosome" uint chromEnd; "End position in chromosome" string name; "Name of item" uint score; "Normalized score from 0-1000" char[1] strand; "+ or - or . for unknown" float level; "Expression level such as RPKM or FPKM. Set to -1 for no data." float signif; "Statistical significance such as IDR. Set to -1 for no data." uint score2; "Additional measurement/count e.g. number of reads. Set to 0 for no data." ) ASCIIrii))ii(z61035 130 - 0.026 0.42 404)i7͘iQΘz61047 136 - 0.029 0.42 404)iΘisz61048 630 - 5.420 0.00 2672399)rZrrrrrisBigWigisBigBedSQL isinstancebytesdecodeentriesr6r7rU)r r bbroutputoexpecteds r testBigBedzTestBigBed.testBigBeds 1 125PP ]]5 !~~{{}!"!{{}!"!  fe $]]7+F}} JJvhx 0CH}} JJvbhhx0(1C DH}}  rN)rrrrrrr rrsrrceZdZdZdZy) TestNumpyc ddl}tjdk(ryddltjdd j dgdj dgdz}j gdj }j gd j }j jjdj } j|||| j gd j }j gd j }j jjdj } j|||| j gdj }j jjdj } jd|j d|j gdj }j jjdj } jd|j d|j jjdj } jdj dj dj d| j jjdj } jdj dj dj d|  jtjd J fd} | d|| d|| d|| d|| d| | d|  jdddd} t| tsJ jdddd} t| jsJ j|j dy)Nrz/tmp/delete.bwrM)rirNrr) rrrrrrrrrrr) rkr#rrrrrrrP) rBnrlrC) isrh) r0rirjriiii"r)rr<) ,i6i@iJiTi^ihiri|iir)rrr<rcg}jd||dzD]}|j|dj|}jj||z dksJy)NrrBrgh㈵>)rErXrallabs)startv2vtr r6s r compyz"TestNumpy.testNumpy..compy!smA\\#ueck: 1  A266&"&&R.4/0 10rrBr0rFnumpyT)rZrrrrWrr7rrrrYrUr<rlistndarrayr[)r rZrrarQrvalues1values2values3values4values5rfoor r6s @@r testNumpyzTestNumpy.testNumpysf >>Q  ]]+S 1 k]Q /3%"*%ARrxx?rxxP"((2992226bjjI ff4 @LTVT\T\]rxxJRTRZRZ["((2992226bjjI ff4 @LTVT\T\]"((2992226bjjI c6 G DLTVT\T\]"((2992226bjjI c6 G D"((2992226bjjI c8288C=xrxx{!U\ ]"((2992226bjjI c8288C=xrxx{!U\ ]  ]]+ ,~~ 2 a c7 c7 c7 c7 c7iiQ5i1#t$%$iiQ4i0#rzz*+*   "#rc tjdk(ryddl}d}tj|d}|j|j dddd|j |j ddd|j d sJ|j|jdddd d d|j |jdddd d |jd sJy) NrrrrrTrr) equal_nanmeanrk) rrrallcloser<rfloat32r5r)r r6r r s r testNumpyValueszTestNumpy.testNumpyValues8s >>Q ` ]]5# &r{{ IIc1bI - BHHRYYsAr*"** =   r{{ HHS!R$H 7 BHHRXXc1b&!4BJJ G   rN)rrrrr rrr rrs >$@ rr) rrRrZsysrrr6rrrrrrr r sE GGR @S S r