:>"^ã@s¢dZddlmZddlmZddlmZdd„Zdd„ZeƒZ d d „Z d d „Z d d„Z e dkržddlmZeƒdS)z3Functions to calculate assorted sequence checksums.é)Úprint_function)Úcrc32)Ú _as_bytesc CsCyttt|ƒƒƒSWn"tk r>tt|ƒƒSYnXdS)zÅReturn the crc32 checksum for a sequence (string or Seq object). Note that the case is important: >>> crc32("ACGTACGTACGT") 20049947 >>> crc32("acgtACGTacgt") 1688586483 N)Ú_crc32rÚstrÚAttributeError)Úseq©r ú:/tmp/pip-build-ww9dw3qa/biopython/Bio/SeqUtils/CheckSum.pyrs rcCs“g}x†tdƒD]x}|}d}xVtdƒD]H}|d@}|dL}|d@r`|dO}|dL}|r2|dN}q2W|j|ƒqW|S)Nérééél0l)ÚrangeÚappend)Ú_table_hÚiZpart_lZpart_hÚjZrflagr r r Ú _init_table_h)s     rcCsyd}d}x\|D]T}|d@d>}|d?}|d?|B}|t|ƒAd@}|t|A}|}qWd||fS)zßReturn the crc64 checksum for a sequence (string or Seq object). Note that the case is important: >>> crc64("ACGTACGTACGT") 'CRC-C4FBB762C4A87EBD' >>> crc64("acgtACGTacgt") 'CRC-DA4509DC64A87EBD' réÿér z CRC-%08X%08X)Úordr)ÚsZcrclZcrchÚcZshrZtemp1hZtemp1lÚidxr r r Úcrc64>s    rc Cs~yt|ƒ}Wntk r$YnXd}}xD|D]<}|d7}||t|jƒƒ7}|dkr6d}q6W|dS)aáReturn the GCG checksum (int) for a sequence (string or Seq object). Given a nucleotide or amino-acid secuence (or any string), returns the GCG checksum (int). Checksum used by GCG program. seq type = str. Based on BioPerl GCG_checksum. Adapted by Sebastian Bassi with the help of John Lenton, Pablo Ziliani, and Gabriel Genellina. All sequences are converted to uppercase. >>> gcg("ACGTACGTACGT") 5688 >>> gcg("acgtACGTacgt") 5688 rr é9i')rrrÚupper)rÚindexZchecksumÚcharr r r ÚgcgVs      r cCsËddl}ddl}|jƒ}yt|ƒ}Wntk rHYnX|jt|jƒƒƒy8|j|j ƒƒ}|j ƒj ddƒj dƒSWntk r®YnX|j |j ƒƒj dƒS)aReturn the SEGUID (string) for a sequence (string or Seq object). Given a nucleotide or amino-acid secuence (or any string), returns the SEGUID string (A SEquence Globally Unique IDentifier). seq type = str. Note that the case is not important: >>> seguid("ACGTACGTACGT") 'If6HIvcnRSQDVNiAoefAzySc6i4' >>> seguid("acgtACGTacgt") 'If6HIvcnRSQDVNiAoefAzySc6i4' For more information about SEGUID, see: http://bioinformatics.anl.gov/seguid/ https://doi.org/10.1002/pmic.200600032 rNÚ Úú=)ÚhashlibÚbase64Úsha1rrÚupdaterrÚ encodebytesÚdigestÚdecodeÚreplaceÚrstripÚ b64encode)rr$r%ÚmÚtmpr r r Úseguidws    # r0Ú__main__)Ú run_doctestN)Ú__doc__Ú __future__rÚbinasciirrZ Bio._py3krrrrr r0Ú__name__Z Bio._utilsr2r r r r Ús     ! '