ó QŽfXc@s²ddlZddlZd„Zd„Zd„Zddd„ƒYZddd„ƒYZd dd „ƒYZd dd „ƒYZd dd„ƒYZ d„Z ddd„ƒYZ dS(iÿÿÿÿNcCs|t|ƒ |kS(N(tlen(tstringtprefix((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pytStringStartsFromPrefix sc Cs“ddddddddd g }tƒ}tƒ}d }xM|D]E}||kr{|j|ƒ|jt|ƒƒd }q@||7}q@W||gS( NtMtItDtNtStHtPt=tXt(tlisttappendtint(tcigartdelimst cigar_listtnum_listtcur_numts((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pyt CIGAR_to_List s!      cCsJ|dks$|dks$|dkr1||7}n|dkrF|}n|S(NRRRR ((t align_lent cigar_chart cigar_lentseq_len((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pytUpdateAlignmentLengths $   t SAM_ConfigcBszeZdZdZdZdZdZdZdZdZ dZ dZ d Z d Z d Zd Zd ZdZdZdZdZRS(s@SQs@HDs@RGs@PGs@COisSN:i iiiiiiiii i (t__name__t __module__t sq_prefixt hd_prefixt rg_prefixt pg_prefixt oc_prefixtsq_tname_indextsq_tname_prefixtnum_mand_fieldst query_indext flag_indext target_indext pos_indext mapq_indext cigar_indextmate_target_indextmate_pos_indext tlen_indext seq_indext qual_index(((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pyR$s&t SAM_entrycBskeZdZdZdZdZdZdZdZdZ dZ dZ dZ dZ eƒZd„Zd„ZRS(R icCsst|jƒ}|d}|d}xItdt|ƒƒD]2}t|j||||t|jƒƒ|_q9WdS(Nii(RRtrangeRRtalentseq(tselftlistst char_listt lens_listti((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pytComputeAlignmentLengthRs   cCsA|jƒ}t|ƒ|jjkrDtjjdƒtjdƒn||jj|_ t ||jj ƒ|_ ||jj |_t ||jjƒ|_t ||jjƒ|_||jj|_||jj|_t ||jjƒ|_t ||jjƒ|_||jj|_||jj|_|jƒdS(Ns9ERROR: Mandatory fields of alignment were not specified! i( tsplitRt sam_configR'tsyststderrtwritetexitR(t query_nameRR)tflagR*t target_nameR+tposR,tmapping_qualityR-RR.tmate_target_nameR/tmate_posR0ttlenR1R6R2tqualR<(R7talignment_stringtsplits((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pyt__init__Ys  (RRRCRDRERFtmapping_qialityRRHRIRJR6RKR5RR>R<RN(((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pyR3Bs  t SAMEntryInfocBseZd„Zd„ZRS(c Csš||_||_||_||_||_|d?d@dk|_||_|d?d@dk|_|d?d@dk|_||_ ||_ | |_ dS(Niii( ttidttnameRFR5R6t is_unmappedtnametproper_alignmentt secondaryRDRKR( R7RQRRRFR5R6RDRTRKR((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pyRNts        cCs’tjj|jdt|jƒdt|jƒdt|jƒdt|jƒdt|j ƒdt|j ƒdt|j ƒdƒdS(Nt s ( R?tstdoutRARTtstrRQRFR5RSRURDRV(R7((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pytPrintƒs(RRRNRZ(((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pyRPss tSamItercBs#eZd„Zd„Zd„ZRS(cCsd|_||_dS(Ni(t iter_indextsam(R7R]((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pyRN‡s cCs |jƒS(N(tnext(R7((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pyt__next__Œsc Csªxš|j|jjƒkrœ|jj|j}|jj|j}t||j|j|j|j |j |j |j |j ƒ }|jd7_|js|SqWtƒ‚dS(Ni(R\R]t NumEntriestentriest target_mapRERPRFR5R6RDRCRKRRVt StopIteration(R7tentryRQt sam_iterator((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pyR^s< (RRRNR_R^(((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pyR[†s  tSamfilecBsÅeZeƒZeƒZeƒZeƒZeƒZeƒZe ƒZ e ƒZ d„Z d„Zd„Zd„Zd„Zd„Zd„Zd„Zd„Zd „Zd „Zd „Zd „ZRS( cCst||jjƒS(N(RR>R (R7tline((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pytIsLineReferenceDescr§scCst||jjƒS(N(RR>R!(R7Rg((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pytIsLineHeaderDescrªscCst||jjƒS(N(RR>R"(R7Rg((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pytIsLineReadGroupDescr­scCst||jjƒS(N(RR>R#(R7Rg((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pytIsLineProgramDescr°scCst||jjƒS(N(RR>R$(R7Rg((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pyt IsLineComment³scCst|ƒ}|S(N(R3(R7Rgt sam_entry((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pyt GetSAMEntry¶s cCs`|jj|ƒ|jƒ}||jj}|t|jjƒ}t|jƒd|j|R%RR&Rb(R7RgRMRE((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pytUpdateTargetFieldsºs  cCsXtƒ|_tƒ|_tƒ|_tƒ|_tƒ|_tƒ|_tƒ|_dS(N( RRotheaderstqueriestprogramstcommentsRatdictRb(R7((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pyt InitFieldsÅs      cCsîtjjdtt|jƒƒdƒtjjdtt|jƒƒdƒtjjdtt|jƒƒdƒtjjdtt|jƒƒdƒtjjdtt|j ƒƒdƒtjjdtt|j ƒƒdƒdS(Ns # targets: s s # headers: s # queries: s # programs: s # comments: s # entries: ( R?RXRARYRRoRqRrRsRtRa(R7((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pyt PrintStatsÎs '''''cCsWtjj|ƒs:tjjd|dƒtjdƒn|jƒt|dƒ}|j ƒ}d|j dRhRiRjRkRlRnRpRvRwRNR`R„R…(((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pyRf›s*                "  ccs+x$|D]}x|D] }|VqWqWdS(N((t iteratorstittelement((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pyt chain_iters  tSamChaincBs#eZd„Zd„Zd„ZRS(cCs ||_dS(N(t sam_files(R7R‹((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pyRNscCs t|jƒS(N(R‰R‹(R7((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pyR„ scCs:x3|jD](}|j|ƒdkr |j|ƒSq WdS(N(R‹R…tNone(R7RRR]((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pyR…s(RRRNR„R…(((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pyRŠs  ((((((( RzR?RRRRR3RPR[RfR‰RŠ(((sg/home/psgendb/BIRCHDEV/pkg/SPAdes-3.13.0/linux-x86_64/share/spades/spades_pipeline/common/sam_parser.pyts     1g