Mercurial > repos > yating-l > jbrowse_hub
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 ''' |