Skip to content

Commit

Permalink
updated foldernames inputmodes <4 in readcounts
Browse files Browse the repository at this point in the history
  • Loading branch information
merlyescalona committed Feb 15, 2018
1 parent 15a6551 commit 39ef9c5
Showing 1 changed file with 20 additions and 10 deletions.
30 changes: 20 additions & 10 deletions ngsphy/readcounts.py
Original file line number Diff line number Diff line change
Expand Up @@ -1424,16 +1424,26 @@ def launchCommand(self, referenceForCurrST, indexREP,indexLOC, individuals,cover
nLoci=self.numLociPerReplicate[indexREP-1]
self.appLogger.debug("Iterating over nLoci: {0}/{1}".format(indexLOC,nLoci))
numGeneTreeDigits=len(str(nLoci))
filepathLoc=os.path.join(\
self.settings.basepath,\
"{0:0{1}d}".format(indexREP,self.numReplicatesDigits),\
"{0}_{1:0{2}d}_TRUE.fasta".format(\
self.settings.simphyDataPrefix,\
indexLOC,\
numGeneTreeDigits\
)\
)

if self.settings.inputmode<4:
filepathLoc=os.path.join(\
self.settings.basepath,\
"REPLICATE_{0:0{1}d}".format(indexREP,self.numReplicatesDigits),\
"{0}_{1:0{2}d}_TRUE.fasta".format(\
self.settings.simphyDataPrefix,\
indexLOC,\
numGeneTreeDigits\
)\
)
else:
filepathLoc=os.path.join(\
self.settings.basepath,\
"{0:0{1}d}".format(indexREP,self.numReplicatesDigits),\
"{0}_{1:0{2}d}_TRUE.fasta".format(\
self.settings.simphyDataPrefix,\
indexLOC,\
numGeneTreeDigits\
)\
)
# dictionary. indices=variable sites, content=variable nucleotide set
self.appLogger.debug(\
"Extracting variable positions from Replicate - {0:0{1}d}/ Locus - {2:0{3}d}".format(\
Expand Down

0 comments on commit 39ef9c5

Please sign in to comment.