Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • How to extract the longest sequence from a fasta file

    I have a fasta file with many sequences, an I would like to get just the longest one (does not matter the id). How can I do that?

    I would like to have as output a fasta file with only the longest sequence.

    thanks!

  • #2
    You can make a read length histogram with BBMap's readlength.sh (or reformat.sh lhist), then run:

    Code:
    reformat.sh in=file.fa out=longest.fa minlength=X

    Comment


    • #3
      Must. Not. Write. Perl. One. Liner.

      --
      Phillip

      Comment


      • #4
        Originally posted by clarissaboschi View Post
        I have a fasta file with many sequences, an I would like to get just the longest one (does not matter the id). How can I do that?
        What about the case where there is more than one sequence which match the longest length? Save all? Only the first or last?

        Comment


        • #5
          Originally posted by pmiguel View Post
          Must. Not. Write. Perl. One. Liner.

          --
          Phillip
          Go right ahead

          Originally posted by kmcarr
          What about the case where there is more than one sequence which match the longest length? Save all? Only the first or last?
          These are interesting questions that make me wonder whether the original problem is worth solving. I'm guessing the OP wants the main chromosome from a perfect bacterial assembly, but perhaps the Clarissa should chime in on what she's trying to do.

          Comment


          • #6
            This sort of task is nearly trivial in Python, it is well worth your while to learn it and the important packages such as BioPython which make jokes like seen in this thread all the more funny.

            So here is a quick and dirty way. A little thought could make it more elegant and/or generalizable. This would just parse through your multifasta file and spit back the longest sequence.

            Code:
            from Bio import SeqIO
            
            myList = []
            
            for seq_record in SeqIO.parse("test.fa", "fasta"):
                myList.append([seq_record.id, str(seq_record.seq), len(seq_record)])
            
            myList.sort(key=lambda x: x[2])
            
            print("the longest sequence is:")
            print(">", myList[-1][0], sep='')
            print(myList[-1][1])

            Comment


            • #7
              Thanks for the suggestions.
              I have a list of protein ids, and I need to retrieve the respective nucleotide sequence. I am using eutils to retrieve them based on these ids, but for some ids, I have many sequences for one id. So for this specific case, I am selecting the longest sequence.

              Comment


              • #8
                /*
                how to compile: gcc -Wall -O2 -o longestseq longestseq.c
                example: ./longestseq example.fa

                */

                #include <stdlib.h>
                #include <string.h>
                #include <stdio.h>

                char s[100000]; // input buffer
                int main(int argc, char *argv[])
                {
                FILE *fp;
                int i;
                size_t cur_spot;
                size_t last_spot;
                size_t best_spot;
                size_t cursize;
                size_t bestsofarsize;

                if (argc != 2) { fprintf(stderr,"Error: need inputfilefasta. This program finds first longest fasta sequence(just one, no same length "ties"). Usage: ./longestseq example.fa\n"); exit(0); }
                fp = fopen(argv[1],"r");
                best_spot = cur_spot = last_spot = ftell(fp);
                cursize = bestsofarsize = 0;
                while (fgets(s,99999,fp))
                {
                for (i=0;s[i];i++) { if ((s[i]=='\r')||(s[i]=='\n')) {s[i] = (char)0;break;}} //handle dos or unix
                if (s[0]=='>')
                {
                last_spot = cur_spot;
                cursize = 0;
                continue;
                }
                cursize += strlen(s);
                if (cursize > bestsofarsize)
                {
                bestsofarsize = cursize;
                best_spot = last_spot;
                }
                cur_spot = ftell(fp);
                }

                fseek(fp,best_spot,SEEK_SET);
                fgets(s,99999,fp);
                printf("%s",s);
                while (fgets(s,99999,fp))
                {
                if (s[0] == '>') break;
                printf("%s",s);
                }
                fclose(fp);
                return 0;
                }


                [FONT="Courier New"][SIZE="1"]____

                Comment

                Latest Articles

                Collapse

                • seqadmin
                  Techniques and Challenges in Conservation Genomics
                  by seqadmin



                  The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                  Avian Conservation
                  Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                  03-08-2024, 10:41 AM
                • seqadmin
                  The Impact of AI in Genomic Medicine
                  by seqadmin



                  Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
                  02-26-2024, 02:07 PM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, 03-14-2024, 06:13 AM
                0 responses
                32 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-08-2024, 08:03 AM
                0 responses
                71 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-07-2024, 08:13 AM
                0 responses
                80 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-06-2024, 09:51 AM
                0 responses
                68 views
                0 likes
                Last Post seqadmin  
                Working...
                X