comparison trackObject.py @ 12:e1f188b43750 draft

planemo upload for repository https://github.com/Yating-L/jbrowse_hub commit f18ea51d27ec7addfa6413716391cfefebc8acbc-dirty
author yating-l
date Tue, 14 Mar 2017 14:31:36 -0400
parents e7c80e9b70ae
children a4a54b925c73
comparison
equal deleted inserted replaced
11:22efe9f6e8a8 12:e1f188b43750
28 Convert gff3, BED, blastxml and gtf files into gff3 files 28 Convert gff3, BED, blastxml and gtf files into gff3 files
29 and store converted files in folder 'raw' 29 and store converted files in folder 'raw'
30 ''' 30 '''
31 fileName = os.path.basename(dataFile) 31 fileName = os.path.basename(dataFile)
32 des_path = os.path.join(self.raw_folder, fileName) 32 des_path = os.path.join(self.raw_folder, fileName)
33 if dataType == 'gff3' or dataType == 'fasta' or dataType == 'bam' or dataType == 'bigwig' or dataType == 'bai': 33 if dataType == 'gff3_mrna' or dataType == 'gff3_transcript' or dataType == 'fasta' or dataType == 'bam' or dataType == 'bigwig' or dataType == 'bai':
34 if dataType == 'gff3':
35 self.checkGff3(dataFile)
36 try: 34 try:
37 shutil.copyfile(dataFile, des_path) 35 shutil.copyfile(dataFile, des_path)
38 except shutil.Error as err1: 36 except shutil.Error as err1:
39 print "Cannot move file, error({0}: {1})".format(err1.errno, err1.strerror) 37 print "Cannot move file, error({0}: {1})".format(err1.errno, err1.strerror)
40 except IOError as err2: 38 except IOError as err2:
51 'fileName': fileName, 49 'fileName': fileName,
52 'dataType': dataType 50 'dataType': dataType
53 } 51 }
54 self.tracks.append(track) 52 self.tracks.append(track)
55 53
56 def checkGff3(self, dataFile): 54
55
56 '''
57 def checkGff3(self, dataFile, dataType):
57 with open(dataFile, 'r') as f: 58 with open(dataFile, 'r') as f:
58 for line in f: 59 for line in f:
59 if not line.startswith('#'): 60 if not line.startswith('#'):
60 seq_type = line.rstrip().split('\t')[2] 61 seq_type = line.rstrip().split('\t')[2]
61 if seq_type == 'transcript': 62 if seq_type == 'transcript':
62 dataType = 'gff3-transcript' 63 return 'gff3-transcript'
63 break
64 if seq_type == 'mRNA': 64 if seq_type == 'mRNA':
65 break 65 return 'gff3'
66 '''