コード例 #1
0
def illumina_files(runobj):
    start = time()
    #    if os.uname()[1] == 'ashipunova.mbl.edu':
    #        import shutil
    #        shutil.rmtree('/Users/ashipunova/BPC/py_mbl_sequencing_pipeline/test/data/fastq/illumina_files_test/output/analysis/')
    illumina_files = IlluminaFiles(runobj)
    illumina_files.split_files(compressed=runobj.compressed)
    illumina_files.perfect_reads()
    illumina_files.uniq_fa()

    elapsed = time() - start
    print "illumina_files time = %s" % str(elapsed)
コード例 #2
0
def illumina_files(runobj):  
    start = time()
    illumina_files = IlluminaFiles(runobj)
    illumina_files.open_dataset_files()
    illumina_files.split_files(compressed = runobj.compressed)
    if runobj.do_perfect: 
        illumina_files.perfect_reads()
    else:
        illumina_files.partial_overlap_reads()
        illumina_files.filter_mismatches()
    illumina_files.uniq_fa()
    
    elapsed = (time() - start)
    print "illumina_files time = %s" % str(elapsed)