def caf_to_fastq(infile, outfile, min_length=0, trim=False):
    '''Convert a CAF file to fastq. Reads shorter than min_length are not output. If clipping information is in the CAF file (with a line Clipping QUAL ...) and trim=True, then trim the reads'''
    caf_reader = caf.file_reader(infile)
    fout = utils.open_file_write(outfile)

    for c in caf_reader:
        if trim:
            if c.clip_start is not None and c.clip_end is not None:
                c.seq.seq = c.seq.seq[c.clip_start:c.clip_end + 1]
                c.seq.qual = c.seq.qual[c.clip_start:c.clip_end + 1]
            else:
                print('Warning: no clipping info for sequence', c.id, file=sys.stderr)


        if len(c.seq) >= min_length:
            print(c.seq, file=fout)

    utils.close(fout)
def caf_to_fastq(infile, outfile, min_length=0, trim=False):
    '''Convert a CAF file to fastq. Reads shorter than min_length are not output. If clipping information is in the CAF file (with a line Clipping QUAL ...) and trim=True, then trim the reads'''
    caf_reader = caf.file_reader(infile)
    fout = utils.open_file_write(outfile)

    for c in caf_reader:
        if trim:
            if c.clip_start is not None and c.clip_end is not None:
                c.seq.seq = c.seq.seq[c.clip_start:c.clip_end + 1]
                c.seq.qual = c.seq.qual[c.clip_start:c.clip_end + 1]
            else:
                print('Warning: no clipping info for sequence', c.id, file=sys.stderr)


        if len(c.seq) >= min_length:
            print(c.seq, file=fout)

    utils.close(fout)
def caf_to_fastq ( infile , outfile , min_length = 0 , trim = False ) : caf_reader = caf . file_reader ( infile ) fout = utils . open_file_write ( outfile ) for c in caf_reader : if trim : if c . clip_start is not None and c . clip_end is not None : c . seq . seq = c . seq . seq [ c . clip_start : c . clip_end + 1 ] c . seq . qual = c . seq . qual [ c . clip_start : c . clip_end + 1 ] else : print ( 'Warning: no clipping info for sequence' , c . id , file = sys . stderr ) if len ( c . seq ) >= min_length : print ( c . seq , file = fout ) utils . close ( fout )
Convert a CAF file to fastq. Reads shorter than min_length are not output. If clipping information is in the CAF file (with a line Clipping QUAL ...) and trim=True, then trim the reads
Convert a CAF file to fastq . Reads shorter than min_length are not output . If clipping information is in the CAF file ( with a line Clipping QUAL ... ) and trim = True then trim the reads