ó êOü_c@sðddlZddlZddlZddlZdddddddd gZd d gZd „Zd dd„ƒYZddd„ƒYZd„Z d„Z d„Z d„Z d„Z d„Zd„Zd„Zd„Zd„Zd„ZdS(iÿÿÿÿNs.fas.fass.fastas.seqs.fsas.fnas.ffns.frns.fqs.fastqcCsN|jƒjdƒr4tjdƒtj||ƒƒStj||ddƒSdS(Ns.gzs.gzipsUTF-8tencodingsutf-8(s.gzs.gzip(tlowertendswithtcodecst getreadertgziptopen(tftmode((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pytOpenstReadercBsPeZd„Zd„Zd„Zd„Zd„Zd„Zd„Zd„Z RS(cCs||_d|_dS(N(thandlertNonetcash(tselfR ((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pyt__init__s cCs(|jdkr$|jjƒ|_ndS(N(R R R treadline(R((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pytFillCashscCs d|_dS(N(R R (R((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pyt TrashCash scCs|jƒ|jS(N(RR (R((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pytTop#s cCs$|jƒ|jƒ}|jƒ|S(N(RRR(Rtresult((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pyR's   cCsng}xXtr`|jƒ}|dks3||ƒr@dj|ƒS|jƒ|j|jƒƒq Wdj|ƒS(Nt(tTrueRtjoinRtappendtstrip(RRRtnext((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pyt ReadUntill-s    cCs»d}g}x„|jƒ r’|jƒdkr’|t|jƒjƒƒ|kr’|j|jƒjƒƒ|t|jƒjƒƒ7}|jƒqW||kr®tdƒ‚ndj|ƒS(NiRs€The sequence and quality strings for one of reads have different length in file {FILE}. Please check the correctness of the file(tEOFRtlenRRRt ExceptionR(Rtbuf_sizetcntR((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pytReadUntillFill7sD cCs|jƒdkS(NR(R(R((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pyRDs( t__name__t __module__RRRRRRR!R(((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pyR s      t SeqRecordcBs8eZdd„Zd„Zd„Zd„Zd„ZRS(cCsˆ|dkr?t|ƒt|ƒkr?tjjd|dƒn|dksit|ƒt|ƒksit‚||_||_||_dS(Ntoppa( R RtsyststdouttwritetAssertionErrortidtseqtqual(RR+R*R,((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pyRIs $*  cCs t|jƒS(N(RR+(R((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pyt__len__QscCs |j|S(N(R+(Rtkey((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pyt __getitem__TscCs!|jdkr|j||!SdS(N(R,R (Rtltr((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pyt QualSubseqWscCsv|dks!|t|jƒkrnt|j||!|jdt|dƒdt|ƒd|j||ƒƒS|SdS(Nit(it-t)(RR+R$R*tstrR2(RR0R1((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pytsubseq\s!MN(R"R#R RR-R/R2R7(((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pyR$Hs     cCsB|dkst‚|dkr(t|ƒS|dkr>t|ƒSdS(Ntfastatfastq(R8R9(R)t parse_fastat parse_fastq(R t file_type((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pytparsecs    ccsŠt|ƒ}xw|jƒs…|jƒjƒ}t|ƒdksO|ddkr^tdƒ‚n|jd„ƒ}t||dƒVqWdS(Niit>s {FILE} is not a valid FASTA filecSs |jdƒS(NR>(t startswith(ts((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pytqs(R RRRRRRR$(R treadertrec_idtrec_seq((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pyR:ks "ccsßt|ƒ}xÌ|jƒsÚ|jƒjƒ}t|ƒdksO|ddkr^tdƒ‚n|jd„ƒ}|jƒ}t|ƒdksž|ddkr­tdƒ‚n|jt|ƒƒ}t||d|ƒVqWdS(Niit@s {FILE} is not a valid FASTQ filecSs |jdƒS(Nt+(R?(R@((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pyRA{sRF( R RRRRRRR!R$(R RBRCRDttmptrec_qual((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pyR;us " "cCs›|dkr;|jd|jdƒ|j|jdƒn\|dkr—|jd|jdƒ|j|jdƒ|jdƒ|j|jdƒndS(NR8R>s R9RERFs+ (R(R*R+R,(trecR R<((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pyR(ƒs   cCs=x6t||ƒD]%}||ƒrt|||ƒqqWdS(N(R=R((t input_handlertoutput_handlerRR<tcontig((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pyt FilterContigsŽs cCsÃx¼t|dƒD]«}d}x0|t|ƒkrN||dkrN|d7}qWt|ƒ}x.||kr‹||ddkr‹|d8}q^W||krtt|j||!|jƒƒqqWdS(NR8itNi(R=RR(R$R+R*(RJRKRLR0R1((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pytRemoveNs”s% # cCsb|jƒ}tjj|ƒ\}}|dkr@|d}}ntjj|ƒ\}}||kS(Ns.gzs.gzipR(s.gzs.gzip(Rtostpathtsplitext(tfnametextension_listtbasename_plus_innerextt outer_exttbasenamet inner_ext((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pytcheck_extension s   cCs t|tƒS(N(RYt fasta_ext(t file_name((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pytis_fasta¬scCs t|tƒS(N(RYt fastq_ext(R[((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pytis_fastq°scCst|dgƒS(Ns.bam(RY(R[((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pytis_bam´scCs8t|ƒrdSt|ƒr dSt|ƒr0dSdSdS(NR9R8tbam(R^R\R_R (R[((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pytget_read_file_type¸s   (((RRR&RPRZR]R R R$R=R:R;R(RMRORYR\R^R_Ra(((sp/oak/stanford/groups/akundaje/marinovg/programs/SPAdes-3.15.0-Linux/share/spades/spades_pipeline/common/SeqIO.pyts&      1