Commit c2485eb7 authored by Kai's avatar Kai

basic structure

parent 8011d64e
import csv
from Bio.Seq import Seq
from Bio.Alphabet import generic_dna
import sys
if __name__ == "__main__":
reader = csv.reader(open("./Pseudomonas_aeruginosa_UCBPP-PA14_109.csv","r"))
web_gene = list(open("./Pseudomonas_aeruginosa_UCBPP-PA14_109.fna","r"))[1:]
web_gene = "".join([e.strip() for e in web_gene])
data = list(reader)
for e in data[3:]:
if len(e)>=5:
_from = int(e[3]) - 1
_to = int(e[4])
_name = e[1][:-1]
_ba = e[-2]
_aa = e[-1]
seq = Seq(_ba,generic_dna)
hybrid = Seq(web_gene[_from-180:_to],generic_dna)
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment