INNER CODE UNIT · Python

report

BirolLab/straglr · extract_repeats.py:54

def report(seqs, out_fa):
    with open(out_fa, 'w') as out:
        for locus in sorted(seqs.keys()):
            for read_name, repeat_size, seq in seqs[locus]:
                out.write('>{} {} {}\n{}\n'.format(read_name, locus, repeat_size, seq))

def parse_args():
    parser = argparse.ArgumentParser()
    parser.add_argument("tsv", type=str, help="Straglr tsv")
    parser.add_argument("bam", type=str, help="bam file") 
    parser.add_argument("out", type=str, help="output fasta")
    parser.add_argument("--locus", type=str, help="UCSC-format coordinate")
    parser.add_argument("--flank", type=int, default=100, help="flank size. Default:100")
    args = parser.parse_args()
    return args
    
def main():
    args = parse_args()

View source record →