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

            • kinpict.com
              Reply to 16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
              by kinpict.com
              The part about this topic feels surprisingly grounded, and the details on the page is the bit that made me slow down. I’d keep <a href="https://kinpict.com/">ai family photo generator</a> bookmarked for later.
              Today, 06:37 PM
            • 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

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by SEQadmin2, 11-12-2026, 07:59 AM
            0 responses
            27 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 11-05-2026, 09:25 AM
            0 responses
            35 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 11-03-2026, 11:54 AM
            2 responses
            56 views
            0 reactions
            Last Post kinpict.com  
            Started by SEQadmin2, 10-28-2026, 09:57 AM
            0 responses
            175 views
            0 reactions
            Last Post SEQadmin2  
            Working...