Mercurial > repos > jay > pdaug_peptide_length_distribution

comparison PDAUG_TSVtoFASTA/PDAUG_TSVtoFASTA.py @ 2:5b1950c46dd5 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload for repository https://github.com/jaidevjoshi83/pdaug commit 8b18552f6d2b2261efebe1075ff4c18a295b94dd"
author jay
date Tue, 29 Dec 2020 18:29:31 +0000
parents 5be6abf137ba
children b5545b41ab70
comparison
equal deleted inserted replaced
1:5be6abf137ba 2:5b1950c46dd5
19 19
20 for line in lines[1:]: 20 for line in lines[1:]:
21 l.append(line.split('\t')[1].strip('\n').strip('\r')) 21 l.append(line.split('\t')[1].strip('\n').strip('\r'))
22 l = list(set(l)) 22 l = list(set(l))
23 23
24 print(l)
25
24 for line in lines: 26 for line in lines:
25 27
26 if l[0] in line.split('\t')[1].strip('\n').strip('\r'): 28 if l[1] in line.split('\t')[1].strip('\n').strip('\r'):
27 n= n+1 29 n= n+1
28 of1.write('>peptide_'+str(n)+'_'+str(l[0])+'\n') 30 of1.write('>peptide_'+str(n)+'_'+str(l[0])+'\n')
29 of1.write(line.split('\t')[0]+'\n') 31 of1.write(line.split('\t')[0]+'\n')
30 32
31 if l[1] in line.split('\t')[1].strip('\n').strip('\r'): 33 if l[0] in line.split('\t')[1].strip('\n').strip('\r'):
32 m= m+1 34 m= m+1
33 of2.write('>peptide_'+str(m)+'_'+str(l[1])+'\n') 35 of2.write('>peptide_'+str(m)+'_'+str(l[1])+'\n')
34 of2.write(line.split('\t')[0]+'\n') 36 of2.write(line.split('\t')[0]+'\n')
35 37
36 elif Method == 'NoClassLabel': 38 elif Method == 'NoClassLabel':