view gc_skew.py @ 9:df7356989ac1 draft

"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/circos commit 041a4d1618298f540e3917f00781e09ec7d71b75"

import sys

import pyBigWig
from Bio import SeqIO
from Bio import SeqUtils


span = int(sys.argv[2])
bw = pyBigWig.open(sys.argv[3], "w")

# Prepare header separately because ugh
data = []
for rec in SeqIO.parse(sys.argv[1], "fasta"):
    data.append((rec.id, len(rec)))
bw.addHeader(data)

for rec in SeqIO.parse(sys.argv[1], "fasta"):
    gc = SeqUtils.GC_skew(rec.seq, span)

    bw.addEntries(rec.id, 0, values=list(gc), span=span, step=span)

bw.close()