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):