Difference between revisions of "JP Feb 16 16"

From GcatWiki
Jump to: navigation, search
Line 4: Line 4:
  
 
Ran Blast2Go using NCBI database and the following parameters (default):
 
Ran Blast2Go using NCBI database and the following parameters (default):
−
 
−
 
−
 
−
 
−
 
−
 
−
 
−
 
−
 
−
<nowiki>
 
−
#file_name = raw_input('Enter file name: ').strip()
 
−
file_name = 'first80.txt'
 
−
 
−
Seqlist= []
 
−
f = open(file_name)
 
−
 
−
for pos, line in enumerate(f):
 
−
    if pos % 4 == 1:
 
−
        Seqlist.append(line)       
 
−
 
−
f.close()
 
−
NewSeqlist = []
 
−
 
−
ID = 1
 
−
 
−
for seq in Seqlist:
 
−
    seq = ">" + str(ID) + "\n" + seq
 
−
    NewSeqlist.append(seq)
 
−
    ID += 1
 
−
   
 
−
print NewSeqlist
 
−
 
−
newfilename = "renamed_" + file_name
 
−
 
−
with open (newfilename, 'w') as newfile:
 
−
    for seq in NewSeqlist:
 
−
        newfile.write(seq)
 
−
</nowiki>
 

Revision as of 19:26, 16 February 2016

Julia Preziosi

Wrote a python code to rename the raw fastq files. (at bottom of page)

Ran Blast2Go using NCBI database and the following parameters (default):