ó £žbc@s[dd„Zdd„Zdd„Zdd„Zded„Zdd„Zdd„ZdS( gI@cCst|||ƒ\}}|S(s8 Abstract: Returns the NG50 value of the passed list of numbers. Comments: Works for any percentage (e.g. NG60, NG70) with optional argument Usage: NG50(numlist, reference_length) Based on the definition from this SEQanswers post http://seqanswers.com/forums/showpost.php?p=7496&postcount=4 (modified Broad Institute's definition https://www.broad.harvard.edu/crd/wiki/index.php/N50) See SEQanswers threads for details: http://seqanswers.com/forums/showthread.php?t=2857 http://seqanswers.com/forums/showthread.php?t=2332 (t NG50_and_LG50(tnumlisttreference_lengtht percentagetng50tlg50((s</home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/N50.pytNG50scCst|||ƒ\}}|S(sÄ Abstract: Returns the LG50 value of the passed list of numbers. Comments: Works for any percentage (e.g. LG60, LG70) with optional argument Usage: LG50(numlist, reference_length) (R(RRRRR((s</home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/N50.pytLG50scCst|t|ƒ|ƒS(s® Abstract: Returns the N50 value of the passed list of numbers. Comments: Works for any percentage (e.g. N60, N70) with optional argument Usage: N50(numlist) (Rtsum(RR((s</home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/N50.pytN50%scCst|t|ƒ|ƒS(s® Abstract: Returns the L50 value of the passed list of numbers. Comments: Works for any percentage (e.g. L60, L70) with optional argument Usage: L50(numlist) (RR(RR((s</home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/N50.pytL50.sc Cs |dkst‚|dks$t‚|r=|jdtƒn|}|d|d}d}x>|D]6}||8}|d7}||krb|}||fSqbWdS(NggY@treverseii(NN(tAssertionErrortsorttTruetNone( RRRt need_sorttstlimitRtlR((s</home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/N50.pyR7s    cCst|t|ƒ|ƒS(N(RR(RR((s</home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/N50.pyt N50_and_L50IscCsâyÉt|ƒdks!tdƒ‚tg|D]}t|tttfƒ^q+ƒs^tdƒ‚|rst|ƒ}ntt|ƒƒ}|dkstdƒ‚ttg|D]}|d^qªƒƒ|SWntk rÝdSXdS(s€ numlist - length of contigs metric explanation http://lh3.github.io/2020/04/08/a-new-metric-on-assembly-contiguity isEmpty list as inputsNon-numerical inputgsall contigs are 0iN( tlenR tallt isinstancetinttfloattcomplexRR(RRtitdenumtn((s</home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/N50.pyt au_metricMs=/ N( RRR R tFalseRRRR(((s</home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/N50.pyts