biopython - Convert FASTA to GenBank -

is there way use biopython convert fasta files genbank format? there many answers on how convert genbank fasta, not other way around.

before convert, must asign alphabet sequence (dna or protein)

from bio import seqio bio.alphabet import generic_dna, generic_protein  input_handle = open("test.fasta", "ru") output_handle = open("test.gb", "w")  sequences = list(seqio.parse(input_handle, "fasta"))  #asign generic_dna or generic_protein seq in sequences:   seq.seq.alphabet = generic_dna  count = seqio.write(sequences, output_handle, "genbank")  output_handle.close() input_handle.close() print "coverted %i records" % count 

for input:

 >i28q9a102fii8j rank=0668881 x=2144.0 y=1105.0 length=418 acgtcatgagagtttgatcatggctcaggacgaacgctggcggcgtgcttaacacatgcaagtcgaacgatgaa gctccagcttgctggggtggattagtggcgaacgggtgagtaacacgtgagtaacctgcccttgactctgggat aagcgttggaaacgacgtctaataccggatatgacgaccgatggcatcatctggttgtggaaagaattttggtc aaggatggactcgcggcctatcaggtagttggtgaggtaatggctcaccaagcctacgacgggtagccggcctg agagggtgaccggccacactgggactgagacacggcccagactcctacgggaggcagcagtggggaatattgca caatgggcgaaagcctgatgcagcaacgccgcgtgagggatgacggcc >i28q9a102jmh72 rank=0320459 x=3829.0 y=3120.0 length=512 acgtcatgagagtttgatcctggctcaggatgaacgctagcggcaggcttaacacatgcaagtcgagggtagaa atagcttgctattttgagaccggcgcacgggtgcgtaacgcgtatgcaatctgccttttacaggggaatagccc agagaaatttggattaatgccccatagcgctgcagggcggcatcgccgagcagctaaagtcacaacggtaaaga tgagcatgcgtcccattagctagttggtaaggtaacggcttaccaaggcgatgatgggtagggtcctgagaggg agatcccccacactggtactgagacacggaccagactcctacgggaggcagcagtgaggaatattggtcaatgg gcgcaagcctgaaccagccatgccgcgtgcaggatgaaggccttcgggttgtaaactgcttttgacggaacgaa aaagct 

you get:

 locus       i28q9a102fii8j           418 bp    dna              unk 01-jan-1980 definition  i28q9a102fii8j rank=0668881 x=2144.0 y=1105.0 length=418 accession   i28q9a102fii8j version     i28q9a102fii8j keywords    . source      .   organism  .             . features             location/qualifiers origin         1 acgtcatgag agtttgatca tggctcagga cgaacgctgg cggcgtgctt aacacatgca        61 agtcgaacga tgaagctcca gcttgctggg gtggattagt ggcgaacggg tgagtaacac       121 gtgagtaacc tgcccttgac tctgggataa gcgttggaaa cgacgtctaa taccggatat       181 gacgaccgat ggcatcatct ggttgtggaa agaattttgg tcaaggatgg actcgcggcc       241 tatcaggtag ttggtgaggt aatggctcac caagcctacg acgggtagcc ggcctgagag       301 ggtgaccggc cacactggga ctgagacacg gcccagactc ctacgggagg cagcagtggg       361 gaatattgca caatgggcga aagcctgatg cagcaacgcc gcgtgaggga tgacggcc // locus       i28q9a102jmh72           450 bp    dna              unk 01-jan-1980 definition  i28q9a102jmh72 rank=0320459 x=3829.0 y=3120.0 length=512 accession   i28q9a102jmh72 version     i28q9a102jmh72 keywords    . source      .   organism  .             . features             location/qualifiers origin         1 acgtcatgag agtttgatcc tggctcagga tgaacgctag cggcaggctt aacacatgca        61 agtcgagggt agaaatagct tgctattttg agaccggcgc acgggtgcgt aacgcgtatg       121 caatctgcct tttacagggg aatagcccag agaaatttgg attaatgccc catagcgctg       181 cagggcggca tcgccgagca gctaaagtca caacggtaaa gatgagcatg cgtcccatta       241 gctagttggt aaggtaacgg cttaccaagg cgatgatggg tagggtcctg agagggagat       301 cccccacact ggtactgaga cacggaccag actcctacgg gaggcagcag tgaggaatat       361 tggtcaatgg gcgcaagcct gaaccagcca tgccgcgtgc aggatgaagg ccttcgggtt       421 gtaaactgct tttgacggaa cgaaaaagct //