No announcement yet.

how to create an output file using perl

  • Filter
  • Time
  • Show
Clear All
new posts

  • how to create an output file using perl

    I have two files.One file contains blast result of contigs and the next file contains the contigs that are used for blast.I want an output text file of sequences in a separate folder which gives no hits in blast result from original fasta file with all contigs.I have no idea about this.Can you please help me to create perl script for this problem.

    Help me pleaseee...
    Attached Files
    Last edited by vineetha; 06-30-2015, 09:22 PM.

  • #2
    I want to ask this question too.


    • #3
      I don't think anyone wants to write a blast parser for you. If you are using blast from the command line, you can tell it to format the output a number of ways, one of those ways might be easier to parse for what you want.


      • #4
        If you want an output file with the "no hits", you could just grep the blast file for "no hits", with -B to grab the query, then grep that to get all the querys only, then grep the blast file with the query file using -f and with grep -A to get the sequence itself. Need to convert the fasta file to single line fasta, though.
        Providing nextRAD genotyping and PacBio sequencing services.