Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • priyanka_15
    Member
    • May 2015
    • 17

    #16
    Tried to export PATH but still not worrking

    Hi

    I tried exporting the PATH But it is still not working ..PIctures attched
    Attached Files

    Comment

    • GenoMax
      Senior Member
      • Feb 2008
      • 7142

      #17
      If trf is the file the program is trying to execute then $PATH directory declaration needs to end at the level of the directory which contains the trf file (e.g. /path_to/repeatscout/RepeatScout-1). Open a new terminal, re-export PATH and try again.
      Last edited by GenoMax; 10-02-2015, 11:28 AM.

      Comment

      • priyanka_15
        Member
        • May 2015
        • 17

        #18
        ok , one more question should this be permanent ? because every time I close terminal it goes back to previous PATH..

        I am new to this all ,perl,linux and repeatscout ; may be these are very easy questions

        Comment

        • priyanka_15
          Member
          • May 2015
          • 17

          #19
          Hello

          I tried to do what you said ,is this the way ? It still is not reading the trf file
          Attached Files

          Comment

          • GenoMax
            Senior Member
            • Feb 2008
            • 7142

            #20
            Two questions.

            1. Have you compiled the software by doing "make" in this directory?
            2. I see the following help when I try to run RepeatScout. Are you trying to follow directions given by someone/some website?

            EDIT: I see that you are on step three of this process. You are trying to process the output of RepeatScount like this?

            $ cat repeats.fa | filter-stage-1.prl > repeats-filtered.prl
            Last edited by GenoMax; 10-02-2015, 11:57 AM.

            Comment

            • GenoMax
              Senior Member
              • Feb 2008
              • 7142

              #21
              Looking at the help further these are the steps involved. What steps have you completed?

              First, build_lmer_table
              creates a file that tabulates the frequency of all l-mers in the
              sequence to be analyzed. Second, RepeatScout takes this table and
              the sequence and produces a fasta file that contains all the repetitive
              elements that it could find. Third, the "filter-stage-1.prl" script
              is run on the output of RepeatScout to remove low-complexity and
              tandem elements; RepeatMasker is run on the sequence of interest using
              this filtered RepeatScout library. The program "filter-stage-2.prl"
              then filters out any repeat element that does not appear a certain number
              of times (by default, 10). Finally, the locations of the repeats found
              by RepeatMasker are used, in conjuction with GFF files that describe
              segmental duplications or exons or other such "uninteresting" regions
              to remove sequences from the library that are likely to not be mobile
              elements; the program "compare-out-to-gff.prl" does exactly this.

              Comment

              • priyanka_15
                Member
                • May 2015
                • 17

                #22
                Hi, yes I did compile it by make ..
                I am trying to learn while doing this so I am using commands from the following link :
                http://http://seqanswers.com/forums/showthread.php?t=5448
                Yes, I am trying to fliter low complexity reads on stage 1.so this is my command:
                cat kib_repeats.fa| ./filter-stage-1.prl >kib_repeats_filtered1.fasta

                Comment

                • GenoMax
                  Senior Member
                  • Feb 2008
                  • 7142

                  #23
                  Have you completed the first 2 steps? Did those run without errors? Can you post the commands you used for those steps? This post has a nice recap of the commands: http://seqanswers.com/forums/showpos...0&postcount=10

                  Comment

                  • priyanka_15
                    Member
                    • May 2015
                    • 17

                    #24
                    Hello

                    yes ,first two ran without error :
                    these are my commands for first two [./build_lmer_table -l 14 -sequence kib.fa -freq ~/Desktop/Kib15.freq
                    ]
                    [/Storage_1/repeatscout/RepeatScout-1$ ./RepeatScout -sequence kib.fa -output kib_repeats.fa -freq Kib15.freq -l 14]
                    I did looked for the commands you have shown ,but I cannot solve this trf problem .

                    Comment

                    • priyanka_15
                      Member
                      • May 2015
                      • 17

                      #25
                      Hi

                      I could run it this time but the fasta file is empty ..

                      Comment

                      • GenoMax
                        Senior Member
                        • Feb 2008
                        • 7142

                        #26
                        You are going to have to install the TRF program from here: http://tandem.bu.edu/trf/trf.download.html Once that is done make sure the trf program is in your path or export TRF_COMMAND='/path_to/trf' is set.

                        You will need nseg as well. I think it is here: ftp://ftp.ncbi.nih.gov/pub/seg/
                        Last edited by GenoMax; 10-02-2015, 01:05 PM.

                        Comment

                        • priyanka_15
                          Member
                          • May 2015
                          • 17

                          #27
                          I did downloaded the trf from here , and also have nseg there..

                          as i think trf is in path right now but empty file again and ahain

                          Comment

                          • GenoMax
                            Senior Member
                            • Feb 2008
                            • 7142

                            #28
                            No other error messages? I tried a small file and step 3 worked for me. Did you rename (or softlink) the executable from trf407b.linux64 to just trf?

                            Comment

                            • priyanka_15
                              Member
                              • May 2015
                              • 17

                              #29
                              yes ,no other error message and yes I did change the name to trf.

                              Comment

                              • priyanka_15
                                Member
                                • May 2015
                                • 17

                                #30
                                did u had nseg in your path too ? I just have it installed the same directory as trf , is that a difference ? should I try to export nseg path too ?

                                Comment

                                Latest Articles

                                Collapse

                                • SEQadmin2
                                  From Collection to Sequencing: Why Sample Preparation and Preservation Define Sequencing Data
                                  by SEQadmin2


                                  Data variability is still an issue in sequencing technologies despite the advances in reproducibility and accuracy of these platforms. But the problem does not originate in the sequencing itself, but in the previous steps, before the sample reaches the sequencer.


                                  The first step is collection, followed by preservation and sample preparation for analysis. Most scientists overlook those steps, but not being careful might just be skewing the experiment’s results.
                                  ...
                                  Yesterday, 10:05 AM
                                • SEQadmin2
                                  Single-Cell Sequencing at an Inflection Point: Early Impacts of New Platforms and Emerging Trends
                                  by SEQadmin2


                                  With the launch of new single-cell sequencing platforms in 2026, the field stands at an exciting inflection point. This article surveys the most impactful advances in the field and discusses how they’re reshaping research in cancer, immunology, and beyond.


                                  Introduction

                                  Single-cell sequencing technologies have undergone remarkable advances over the past decade, transitioning from low-throughput experimental approaches to highly scalable platforms capable of...
                                  05-22-2026, 06:42 AM
                                • SEQadmin2
                                  Environmental Genomics in the Age of NGS: From Microbes to Conservation Strategies
                                  by SEQadmin2

                                  Studying ecosystems means dealing with complex, multi-species communities that are hard to observe at scale. This complexity, however, hides many important questions to be answered, from how biogeochemical cycles work and how climate change can affect species distribution to how conservation strategies can work best.


                                  Genomics, particularly since the expansion of NGS, has transformed ecosystem ecology. By sequencing environmental DNA, we can now assess biodiversity without direct...
                                  05-06-2026, 09:04 AM

                                ad_right_rmr

                                Collapse

                                News

                                Collapse

                                Topics Statistics Last Post
                                Started by SEQadmin2, Yesterday, 12:03 PM
                                0 responses
                                19 views
                                0 reactions
                                Last Post SEQadmin2  
                                Started by SEQadmin2, Yesterday, 11:40 AM
                                0 responses
                                14 views
                                0 reactions
                                Last Post SEQadmin2  
                                Started by SEQadmin2, 05-28-2026, 11:40 AM
                                0 responses
                                29 views
                                0 reactions
                                Last Post SEQadmin2  
                                Started by SEQadmin2, 05-26-2026, 10:12 AM
                                0 responses
                                31 views
                                0 reactions
                                Last Post SEQadmin2  
                                Working...