aFN = sys.argv[2] 

seqFN = sys.argv[3] 
tFN = '/home/chrisgre/scripts/alignSeqs/dRNA.results.updated'

timer = bioLibCG.cgTimer()
timer.start()


#initialize oRNA
blankIDs.blankIDs(seqFN, oFN)

#initializing alignments
print '...appending T Info'
parRun(30, 3, '/home/chrisgre/scripts/alignSeqs/cgAlignmentFlat.py', 'appendTInfo', aFN, tFN)
parClean(aFN, 30)
print '...appending Tran Info'
parRun(30, 3, '/home/chrisgre/scripts/alignSeqs/cgAlignmentFlat.py', 'appendTranInfo', aFN, tFN)
parClean(aFN, 30)
print timer.split()


print '...updating paired interactions: centered mismatches and center expression'
parRun(30, 5, '/home/chrisgre/scripts/alignSeqs/updateMismatchAndMiddleFlat.py', 'markCenterExpression', aFN, '/home/chrisgre/smallLibs/siRNA/degradome/wigsk50')
parClean(aFN, 30)
print timer.split()
parRun(30, 3, '/home/chrisgre/scripts/alignSeqs/updateMismatchAndMiddleFlat.py', 'markMismatchedPairs', aFN)
parClean(aFN, 30)
print timer.split()

print 'Linking oRNA to Targets'
import bioLibCG
import sys
from parRun import parRun
from checkExit import parClean
import subprocess

#fN and Init
aFN = sys.argv[1] 
aFilteredFN = sys.argv[2]
cLevel = sys.argv[3]

timer = bioLibCG.cgTimer()
timer.start()

#filter targets
print 'filtering targets'
parRun(5, 3, '/home/chrisgre/scripts/endoClip/filteringFlat.py', 'filterTargetsInPlace', aFN, 'True', '1', '1', '.%s' % cLevel)
parClean(aFN, 5)
print timer.split()

#make db smaller
print 'truncating db'
subprocess.Popen(['/home/chrisgre/scripts/endoClip/truncate.filterAlignments.sh', aFN, aFilteredFN]).wait()
print timer.split()




aFN = runName + '/all.aligned.0.5.tabs'

seqFN = runName + '/oRNA.sequences'
tFN = runName + '/dRNA.results'

timer = bioLibCG.cgTimer()
timer.start()

#initialize oRNA
blankIDs.blankIDs(seqFN, oFN)

#initializing alignments
print '...appending T Info'
parRun(40, 3, '/home/chrisgre/scripts/alignSeqs/cgAlignmentFlat.py',
       'appendTInfo', aFN, tFN)
parClean(aFN, 40)
print '...appending Tran Info'
parRun(40, 3, '/home/chrisgre/scripts/alignSeqs/cgAlignmentFlat.py',
       'appendTranInfo', aFN, tFN)
parClean(aFN, 40)
print timer.split()

print '...updating paired interactions: centered mismatches and center expression'
parRun(30, 4,
       '/home/chrisgre/scripts/alignSeqs/updateMismatchAndMiddleFlat.py',
       'markCenterExpression', aFN,
       '/home/chrisgre/smallLibs/siRNA/degradome/wigsk50')
parClean(aFN, 30)
print timer.split()
parRun(50, 3,
       '/home/chrisgre/scripts/alignSeqs/updateMismatchAndMiddleFlat.py',
import blankIDs
from parRun import parRun
from checkExit import parClean

#files and init
aFN = sys.argv[1]
dFN = sys.argv[2]

timer = bioLibCG.cgTimer()
timer.start()

#initializing alignments
'''
print '...appending T Info'
parRun(50, 3, '/home/chrisgre/myLibs/cgAlignmentFlat.py', 'appendTInfoFlat', aFN, dFN)
parClean(aFN, 50)
print timer.split()
'''

print '...updating paired interactions: centered mismatches and center expression'
parRun(40, 5, '/home/chrisgre/scripts/endoClip/updateMismatchAndMiddleFlat.py',
       'markCenterExpression', aFN,
       '/home/chrisgre/smallLibs/siRNA/degradome/wigsk50')
parClean(aFN, 40)
print timer.split()
'''
parRun(50, 3, '/home/chrisgre/scripts/endoClip/updateMismatchAndMiddleFlat.py', 'markMismatchedPairs', aFN)
parClean(aFN, 50)
print timer.split()
'''
from parRun import parRun
from checkExit import parClean


#files and init
aFN = sys.argv[1] 
dFN = sys.argv[2]

timer = bioLibCG.cgTimer()
timer.start()


#initializing alignments
'''
print '...appending T Info'
parRun(50, 3, '/home/chrisgre/myLibs/cgAlignmentFlat.py', 'appendTInfoFlat', aFN, dFN)
parClean(aFN, 50)
print timer.split()
'''

print '...updating paired interactions: centered mismatches and center expression'
parRun(40, 5, '/home/chrisgre/scripts/endoClip/updateMismatchAndMiddleFlat.py', 'markCenterExpression', aFN, '/home/chrisgre/smallLibs/siRNA/degradome/wigsk50')
parClean(aFN, 40)
print timer.split()

'''
parRun(50, 3, '/home/chrisgre/scripts/endoClip/updateMismatchAndMiddleFlat.py', 'markMismatchedPairs', aFN)
parClean(aFN, 50)
print timer.split()
'''
#files and init
aFN = sys.argv[1]
dFN = sys.argv[2]

timer = bioLibCG.cgTimer()
timer.start()

#initializing alignments
print '...appending T Info'
parRun(5, 3, '/home/chrisgre/myLibs/cgAlignmentFlat.py', 'appendTInfoFlat',
       aFN, dFN)
parCleanSplit(aFN, 5)
print timer.split()

print '...updating paired interactions:  center expression'
#splitRun(aFN, 5, '/home/chrisgre/scripts/endoClip/updateMismatchAndMiddleFlat.py', 'markCenterExpression', 'splitFN', '/home/chrisgre/smallLibs/siRNA/degradome/wigsk50')
#splitRun(aFN, 5, '/home/chrisgre/scripts/endoClip/updateMismatchAndMiddleFlat.py', 'markCenterExpression', 'splitFN', '/home/chrisgre/smallLibs/siRNA/degradome/HeLa/wigs.3t7n0m20k20b.contigFiltered')
splitRun(aFN, 5,
         '/home/chrisgre/scripts/endoClip/updateMismatchAndMiddleFlat.py',
         'markCenterExpression', 'splitFN',
         '/home/chrisgre/data/Lab_Data/U87/result_primary/wigs.n0m20k20b')
parCleanSplit(aFN, 5)
print timer.split()

print '...updating paired interactions:  center mismatch'
splitRun(aFN, 3,
         '/home/chrisgre/scripts/endoClip/updateMismatchAndMiddleFlat.py',
         'markMismatchedPairs', 'splitFN')
parClean(aFN, 5)
print timer.split()
Exemple #7
0
import updateDegPeaks as degPeaks
import sys
import blankIDs
from parRun import parRun
from checkExit import parClean


dRNAFN = sys.argv[1]
peakFN = sys.argv[2]
seqFN = sys.argv[3]

timer = bioLibCG.cgTimer()
timer.start()

#initialize oRNA database
print 'make blank IDs'
blankIDs.blankIDs(seqFN, dRNAFN)
print 'sequence'
degPeaks.updateSequence(dRNAFN, seqFN)
print 'tcc'
degPeaks.updateTcc(dRNAFN, peakFN)
print 'eLevel'
degPeaks.updateELevel(dRNAFN, '/home/chrisgre/smallLibs/siRNA/degradome/wigsk50')
print 'gSequence'
degPeaks.updateGSequence(dRNAFN)

print 'gScore'
parRun(50, 3, '/home/chrisgre/scripts/endoClip/updateDegPeaks.py', 'updateGScore', dRNAFN)
parClean(dRNAFN, 50)

oFN = runName + "/oRNA.data"
aFN = runName + "/all.aligned.0.5.tabs"

seqFN = runName + "/oRNA.sequences"
tFN = runName + "/dRNA.results"

timer = bioLibCG.cgTimer()
timer.start()

# initialize oRNA
blankIDs.blankIDs(seqFN, oFN)

# initializing alignments
print "...appending T Info"
parRun(40, 3, "/home/chrisgre/scripts/alignSeqs/cgAlignmentFlat.py", "appendTInfo", aFN, tFN)
parClean(aFN, 40)
print "...appending Tran Info"
parRun(40, 3, "/home/chrisgre/scripts/alignSeqs/cgAlignmentFlat.py", "appendTranInfo", aFN, tFN)
parClean(aFN, 40)
print timer.split()


print "...updating paired interactions: centered mismatches and center expression"
parRun(
    30,
    4,
    "/home/chrisgre/scripts/alignSeqs/updateMismatchAndMiddleFlat.py",
    "markCenterExpression",
    aFN,
    "/home/chrisgre/smallLibs/siRNA/degradome/wigsk50",
)