Unconfigured Ad

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

    #1

    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

            • SEQadmin2
              Beyond CRISPR/Cas9: Understand, Choose, and Use the Right Genome Editing Tool
              by SEQadmin2



              CRISPR/Cas9 sparked the gene editing revolution for both research and therapeutics.1 But this system still showed severe issues that limited its applications. The most prominent were the heavy reliance on PAM sequences, delivery limitations, double-stranded breaks that prompt unintended edits and cell death, and editing inefficiency (both in targeting and in knock-in reliability).

              Despite this, “CRISPR helped turn genome editing from a specialized technique into
              ...
              07-31-2026, 11:01 AM
            • SEQadmin2
              Proteomic Platforms: How to Choose the Right Analytical Strategy to Improve Detection and Clinical Applications
              by SEQadmin2


              Proteomics platforms are evolving rapidly, with advances in mass spectrometry and affinity-based approaches expanding what researchers can detect and at what scale. As the field moves toward deeper proteome coverage and clinical applications, scientists face an increasingly complex landscape of tools. This article will explore how researchers are navigating these choices to find the right platform for their work.

              The systematic characterization of the human proteome has
              ...
              07-20-2026, 11:48 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by SEQadmin2, Yesterday, 12:22 PM
            0 responses
            12 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 08-11-2026, 10:35 AM
            0 responses
            14 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 08-06-2026, 07:41 AM
            0 responses
            31 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 08-03-2026, 10:13 AM
            0 responses
            49 views
            0 reactions
            Last Post SEQadmin2  
            Working...