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
              16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
              by SEQadmin2





              NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification,...
              10-30-2026, 06:46 AM
            • SEQadmin2
              Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
              by SEQadmin2



              Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

              Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
              ...
              10-14-2026, 09:55 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by SEQadmin2, 11-05-2026, 09:25 AM
            0 responses
            20 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 11-03-2026, 11:54 AM
            0 responses
            34 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 10-28-2026, 09:57 AM
            0 responses
            82 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 10-28-2026, 09:51 AM
            0 responses
            50 views
            0 reactions
            Last Post SEQadmin2  
            Working...