tFN = 'dRNA.results.e.t' timer = bioLibCG.cgTimer() timer.start() #initialize oRNA database print 'initializing oRNA' si.updateID(oRNADir, peakFN) print timer.split() print 'truncating alignments' #Truncate alignment file for specified mismatches #si.truncateAlignments(alignFN, minMismatches, maxMismatches, alignTrunFN) print 'loading alignments' #initialize alignments #cgAlignment.loadAlignments2(aDir, alignTrunFN) #cgAlignment.appendTInfo(aDir, tFN) #cgAlignment.appendTranInfo(aDir, tFN) print timer.split() print 'updating initial targets' #update the initial targets for small rnas based off of only local alignment #si.updateTargetIDs(oRNADir, aDir) print timer.split() print 'updating paired interactions: centered mismatches and center expression' mm.markCenterExpression(aDir, 'siDegradome.conf') mm.markMismatchedPairs(aDir) print timer.split()
alignTrunFN = runName + '.all.aligned.nodups' + '.' + str(minMismatches) + '.' + str(maxMismatches) sFN = runName + '.small.results' tFN = runName + '.degradome.results' #Get expression of the small RNAs and the targets print 'updating expression levels' #si.updateSmallExpression(sFN, 'siPeaks.conf', 1, 3, sFN) #si.updateSmallExpression(tFN, 'siDegradome.conf', 1, 2, tFN) print 'truncating alignments' #Truncate alignment file for specified mismatches si.truncateAlignments(alignFN, minMismatches, maxMismatches, alignTrunFN) print 'updating initial targets' #update the initial targets for small rnas based off of only local alignment si.updateTargetIDs(sFN, alignTrunFN, 0, 4, sFN) print 'update microrna/transcript overlaps' #update microRNA biogenesis for small and transcript #si.updateMicroRNAOverlap(sFN, 'mirBaseHumanNew.gff.Tcc', 1, 5, sFN) #si.transcriptSetOverlap(sFN, False, 1, 6, sFN) #si.transcriptSetOverlap(tFN, True, 1, 3, tFN) print 'updating paired interactions: centered mismatches and center expression' mm.markCenterExpression(sFN, tFN, alignTrunFN, 'siDegradome.conf', runName + '.pair.center.data') mm.markMismatchedPairs(sFN, alignTrunFN, runName + '.pair.mismatch.data')
tFN = "/home/chrisgre/scripts/alignSeqs/dRNA.results.updated" timer = bioLibCG.cgTimer() timer.start() # initialize oRNA database print "initializing oRNA" si.updateIDFromQuery(oRNADir, seqFN) print timer.split() # initializing alignments print "Initializing alignments" print "...loading alignments" cgAlignment.loadAlignments2(aDir, alignTrunFN) print "...appending T Info" cgAlignment.appendTInfo(aDir, tFN) print "...appending Tran Info" cgAlignment.appendTranInfo(aDir, tFN) print timer.split() print "...updating paired interactions: centered mismatches and center expression" mm.markCenterExpression(aDir, "siDegradome.conf") mm.markMismatchedPairs(aDir) print timer.split() print "Linking oRNA to Targets" # update the initial targets for small rnas based off of only local alignment si.updateTargetIDs(oRNADir, aDir) print timer.split(), "done", runName
import updateMismatchAndMiddle as mm mm.markCenterExpression('tester.tester','degradome.results', 'all.aligned.nodups.0.5', 'siDegradome.conf', 'test.pair.center.data')
import updateMismatchAndMiddle as mm mm.markCenterExpression('tester.tester', 'degradome.results', 'all.aligned.nodups.0.5', 'siDegradome.conf', 'test.pair.center.data')