def updateGeneName(dFN,
fFN,
wigDir,
chrom,
strand,
prefix,
switchStrand=False):
NX = Nexus(dFN, fFN)
NX.load(['geneNames', 'tcc'])
if switchStrand:
strand = -strand
strand = str(strand)
coord_gName = cgWig.loadSingleWigTranscript(wigDir, chrom, strand, prefix)
while NX.nextID():
chrom, strand, start, end = bioLibCG.tccSplit(NX.tcc)
overlappingGenes = coord_gName.get(start, ".")
if overlappingGenes == "NONE":
NX.geneNames = []
else:
NX.geneNames = overlappingGenes.split(',')
NX.save()
def updateGeneName(dFN, fFN, wigDir, chrom, strand, prefix, switchStrand = False):
NX = Nexus(dFN, fFN)
NX.load(['geneNames', 'tcc'])
if switchStrand:
strand = -strand
strand = str(strand)
coord_gName = cgWig.loadSingleWigTranscript(wigDir, chrom, strand, prefix)
while NX.nextID():
chrom, strand, start, end = bioLibCG.tccSplit(NX.tcc)
overlappingGenes = coord_gName.get(start, ".")
if overlappingGenes == "NONE":
NX.geneNames = []
else:
NX.geneNames = overlappingGenes.split(',')
NX.save()
