INNER CODE UNIT · Python

_check_fasta

Gabaldonlab/redundans · redundans.py:474

def _check_fasta(lastOutFn, minSize=1000, log=sys.stderr):
    """Exit if empty FastA file"""
    faidx = FastaIndex(open(lastOutFn))
    if not len(faidx) or faidx.genomeSize < minSize:
        log.write("[ERROR] Empty FastA file encountered: %s !\n" % lastOutFn)
        sys.exit(1)


def redundans(
    fastq,
    longreads,
    fasta,
    reference,
    outdir,
    mapq,
    threads,
    mem,
    resume,

View source record →