Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • #16
    It appears that your index files are not matching the reference files. That is likely the source of all the problems you have posted across several threads.

    I would suggest getting the complete set of index, GTF and sequence files from one source so that they all match: http://cufflinks.cbcb.umd.edu/igenomes.html

    Comment


    • #17
      I have used all ensembl .. I downloaded the whole pack of tar.gz files from ensembl cos last time you recommended me... something is wrong... it is only at cuffmerge and cuffdiff.... I have tried both UCSC and Ensembl.
      Well, I can it one more time and get back to you.
      Thank you...

      Comment


      • #18
        If you downloaded the "Ensembl" archive file from iGenomes then the genome file is located at /Homo_sapiens/Ensembl/GRCh37/Sequence/WholeGenomeFasta/genome.fa.

        This genome file contains no "GL000*" fasta entries. So I am not sure where you are getting that error from if you are using this Ensembl genome file for cuffdiff analysis.

        Can you post the command line you are using?

        Comment


        • #19
          Dear GenoMax
          I thank you very much for your kind reply and willingness to help me out of this problem.

          My commands and error are attached in the file named 'Errors_RNA-seq_Cuffmerge'. I also think that I didnt sorted .bam files, could that be a reason gettin errors ?
          Please have a look and suggest me something...
          Thank you..


          Originally posted by GenoMax View Post
          If you downloaded the "Ensembl" archive file from iGenomes then the genome file is located at /Homo_sapiens/Ensembl/GRCh37/Sequence/WholeGenomeFasta/genome.fa.

          This genome file contains no "GL000*" fasta entries. So I am not sure where you are getting that error from if you are using this Ensembl genome file for cuffdiff analysis.

          Can you post the command line you are using?
          Attached Files
          Last edited by Charitra; 04-23-2013, 12:33 AM.

          Comment


          • #20
            Thanks for providing the detailed log file. You appear to be on right track.

            In spite of those "Warnings" (they are not errors) are you getting output from the cuffmerge step?

            I checked a few of the ID's in the "Warning" messages and most of those look to be pseudogenes or lincRNA's etc. Though they are present in the gtf file no fasta records are created for them when you use the "gffread" utility (http://cufflinks.cbcb.umd.edu/gff.html)

            How much memory do you have on this machine? Are you using 32-bit or 64-bit ubuntu?

            Comment


            • #21
              I am using Server 64 bit ubuntu has enough memory for multiple tasks ~50GB..
              your mail is a relief for me. I have visualised the outputs of cufflinks using cummeRbund... but I was not sure about my cuffmerge and cuffdiff...

              Would you like to recommend me something please ?




              Originally posted by GenoMax View Post
              Thanks for providing the detailed log file. You appear to be on right track.

              In spite of those "Warnings" (they are not errors) are you getting output from the cuffmerge step?

              I checked a few of the ID's in the "Warning" messages and most of those look to be pseudogenes or lincRNA's etc. Though they are present in the gtf file no fasta records are created for them when you use the "gffread" utility (http://cufflinks.cbcb.umd.edu/gff.html)

              How much memory do you have on this machine? Are you using 32-bit or 64-bit ubuntu?

              Comment


              • #22
                Are you getting output from Cuffmerge/Cuffdiff or there is nothing being produced?

                Comment


                • #23
                  Results are there... I have not yet analysed properly... I just got the output from R but didnt read the results from R... results are not empty... there are genes but there is not significant differential gene expression but still not sure...
                  do you think.. i can continue with these warnings and analyse the R results ??

                  Comment


                  • #24
                    Dear GenoMax

                    I got one doubt, if there are lincRNA's and I just ignore them, then I am trying to differentiate the expression level in my samples. Sometimes, these are needed for differentiation ? is that okay if I continue... or is there anything where I can modify my commands to get more accurate results.
                    Thanks




                    Originally posted by GenoMax View Post
                    Thanks for providing the detailed log file. You appear to be on right track.

                    In spite of those "Warnings" (they are not errors) are you getting output from the cuffmerge step?

                    I checked a few of the ID's in the "Warning" messages and most of those look to be pseudogenes or lincRNA's etc. Though they are present in the gtf file no fasta records are created for them when you use the "gffread" utility (http://cufflinks.cbcb.umd.edu/gff.html)

                    How much memory do you have on this machine? Are you using 32-bit or 64-bit ubuntu?

                    Comment


                    • #25
                      Originally posted by Charitra View Post
                      Dear GenoMax

                      I got one doubt, if there are lincRNA's and I just ignore them, then I am trying to differentiate the expression level in my samples. Sometimes, these are needed for differentiation ? is that okay if I continue... or is there anything where I can modify my commands to get more accurate results.
                      Thanks
                      You should create a separate posting with this question since I can't answer your question (and biologists on the forum may not read this thread). With a new thread you will have a better chance of getting an answer.

                      Comment


                      • #26
                        Originally posted by GenoMax View Post
                        You should create a separate posting with this question since I can't answer your question (and biologists on the forum may not read this thread). With a new thread you will have a better chance of getting an answer.
                        Dear GenoMax,
                        I got cuffdiff results and , i think, gene name or somthing went wrong which was giving warnings because i can not fined my gene in cummeRbund but gene is there in database when i checked opening in excel. is this problem with cuffdiff ?

                        kind help appreciated..

                        Comment

                        Latest Articles

                        Collapse

                        • seqadmin
                          Strategies for Sequencing Challenging Samples
                          by seqadmin


                          Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                          03-22-2024, 06:39 AM
                        • 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

                        ad_right_rmr

                        Collapse

                        News

                        Collapse

                        Topics Statistics Last Post
                        Started by seqadmin, 03-27-2024, 06:37 PM
                        0 responses
                        13 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 03-27-2024, 06:07 PM
                        0 responses
                        11 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 03-22-2024, 10:03 AM
                        0 responses
                        53 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 03-21-2024, 07:32 AM
                        0 responses
                        69 views
                        0 likes
                        Last Post seqadmin  
                        Working...
                        X