Skip to content

Commit d9ae6b1

Browse files
authored
Merge pull request #21 from Quarkins/master
Stranded two contig bug
2 parents e0c309b + 725b58d commit d9ae6b1

File tree

2 files changed

+10
-1
lines changed

2 files changed

+10
-1
lines changed

BuildSuperTranscript.py

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -233,7 +233,7 @@ def BuildGraph(fname,transcripts,verbose=False):
233233
#Filter psl table where two transcripts can have multiple rows (usually because blat does T1 vs T2 then T2 vs T1 later)
234234
bData, trandir = filt_dir(bData)
235235

236-
if(len(bData['strand'].unique()) > 1): #That is we have both pos and neg strands
236+
if(len(bData['strand'].unique()) > 1 or bData['strand'].unique() == '-'): #That is we have both pos and neg strands or potentially two transcripts with a negative strand
237237
print("Double Stranded Contigs\n")
238238

239239
#Re-correct the transcripts to be the reverse compliments if one of the transcripts has a negative directionality

Lace.py

Lines changed: 9 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -165,6 +165,15 @@ def Split(genome,corsetfile,ncore,maxTran,outdir):
165165

166166

167167
if __name__ == '__main__':
168+
169+
#Print Lace Version
170+
print(" __ __ ____ ____ ")
171+
print("( ) / _\ / )( __)")
172+
print("/ (_/\/ \( (__ ) _) ")
173+
print("\_____/\_/\_/\_____)(____)")
174+
print("Lace Version: 0.80")
175+
print("Last Editted: 07/11/16")
176+
168177

169178
#Make argument parser
170179
parser = argparse.ArgumentParser()

0 commit comments

Comments
 (0)