Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • bio_informatics
    Senior Member
    • Nov 2013
    • 182

    Error while running clustalw2 in biopython

    Code:
    import re, sys, argparse, os, subprocess
    from subprocess import call
    from Bio import SeqIO
    from Bio import AlignIO
    from Bio.Seq import Seq
    from Bio.Align.Applications import ClustalwCommandline
    
    for key in store_seq:
        fileN=str(key)+".txt" #"1.txt"
        call(["touch",fileN])
        outFile=str(key)+".fasta" #"1.fasta"
        
        with open(fileN,'w') as f:
            f.write(">"+str(store_seq[key])+"\n") # print seq id
            f.write(str(chemo_seq[store_seq[key]] )+"\n")#print sequence
            f.write( ">"+str(brig_orth[store_seq[key]])+ "\n")
            f.write ( str(chemo_seq[brig_orth[store_seq[key]]]) +"\n")
            f.write( ">"+str(reman_orth[store_seq[key]]) +"\n" )
            f.write( str(chemo_seq[reman_orth[store_seq[key]]]) + "\n")
             #run clustalw2 for the file..--->      
    cline=ClustalwCommandline("clustalw2",infile=fileN,type="PROTEIN",output="FASTA",outfile=outFile,quiet)
        cline()
    
        #one key ends..
    Error:cline=ClustalwCommandline("clustalw2",infile=fileN,type="PROTEIN",output="FASTA",outfile=outFile,quiet)
    SyntaxError: non-keyword arg after keyword arg
    In the above code, cline=Clus.. is not indenting with with open, kindly pardon me.

    However, when I run on python shell it runs smoothly. I have clustaw2 in my path. Instead of fileN and outFile, I give file name in double quotes.

    Editor using: emacs
    What is the cause for this?
    How do I correct it?
    Python 2.7
    Bioinformaticscally calm
  • maubp
    Peter (Biopython etc)
    • Jul 2009
    • 1544

    #2
    The error is in your Python syntax, probably you need quiet=True at the end of the line defining cline.

    Right now you have several keyword arguments (key=value) followed by a non-keyword argument (value). In Python any non-keyword arguments (also called positional arguments) must come before any keyword arguments.
    Last edited by maubp; 07-28-2014, 06:49 AM. Reason: Correct type (quite vs quiet)

    Comment

    • maubp
      Peter (Biopython etc)
      • Jul 2009
      • 1544

      #3
      This is a duplicate of https://www.biostars.org/p/107528/ - please don't double post like this without cross linking. It saves duplication of effort.

      Comment

      • bio_informatics
        Senior Member
        • Nov 2013
        • 182

        #4
        Hi maubp,
        Thank you for pointing that.
        I didn't understand this point:

        Right now you have several keyword arguments (key=value) followed by a non-keyword argument (value). In Python any non-keyword arguments (also called positional arguments) must come before any keyword arguments.

        I apologize for multi posts.
        Bioinformaticscally calm

        Comment

        • maubp
          Peter (Biopython etc)
          • Jul 2009
          • 1544

          #5
          Try this:

          cline=ClustalwCommandline("clustalw2", infile=fileN, type="PROTEIN", output="FASTA", outfile=outFile, quiet=True)

          Comment

          • bio_informatics
            Senior Member
            • Nov 2013
            • 182

            #6
            Thank you.
            Originally posted by maubp View Post
            Try this:

            cline=ClustalwCommandline("clustalw2", infile=fileN, type="PROTEIN", output="FASTA", outfile=outFile, quiet=True)
            Bioinformaticscally calm

            Comment

            Latest Articles

            Collapse

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by SEQadmin2, Today, 06:09 AM
            0 responses
            11 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 06-09-2026, 11:58 AM
            0 responses
            33 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 06-05-2026, 10:09 AM
            0 responses
            38 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 06-04-2026, 08:59 AM
            0 responses
            43 views
            0 reactions
            Last Post SEQadmin2  
            Working...