Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • #16
    Originally posted by ychen View Post
    Hi Folks,

    I feel lucky to find this thread because I have been struggling with the same problems. After splitting the unusual FASTQ files, my TopHat results are still quite different from what reported in the recent published paper. Can you tell me where to find the provided SAM file? I want to try the the reported command line.

    Thanks a lot,

    Yi-Shiou
    Hi,

    you could find the SAM file in their suplemental material online.

    Comment


    • #17
      Hi syslm01,

      I look into the supplementary information page many times but didn't find any SAM file. Did I miss something very obvious or just look into a wrong place?

      Thanks again,

      Yi-Shiou

      Comment


      • #18
        hi ychen,

        here is the link http://www.ncbi.nlm.nih.gov/geo/quer...i?acc=GSE20846
        GSE20846_RAW.tar contains the SAM and gtf files

        Comment


        • #19
          Hi syslm01,

          Thanks so much for your help, I really appreciate it.


          Yi-Shiou

          Comment


          • #20
            Hi,

            I splitted the fastq files with fastx into the original paired end files (positions 1-75 and 77-151) and ran TopHat with default settings. In the paper they did use:

            tophat -p 8 -F 0.0 -r 50 -m 1 --no-novel-juncs -j ../../pooled_200bp_frags.juncs -o pooled_tophat2 -a 8 /fs/szasmg3/cole/ebwts/m_musculus/mm9/fast_mm9/mm9+wold_spikes s1_1.query75.txt,s2_1.query75.txt s1_2.query75.txt,s2_2.query75.txt
            So the difference to my run is, that I left out "--no-novel-juncs" and "-j ../../pooled_200bp_frags.junc". Further I used the default bowtie mm9 reference.

            I loaded the results into IGV and took a look at the results. The result is very different in many ways.
            1. The coverage seems to be lower for almost every positions. accepted_hits.sam contains about half as much entries as the result from the paper.
            2. My base phred quality varies while the quality from the paper is everywhere 40.
            3. So many reads could be paired mapped (in IGV: "Pair is mapped = No")

            I did run Tophat with different parameters, but the result should not be that different? Maybe I missed something when I separated fastq file into two? Does anyone know how to generate such a splice junction file?

            Many thanks for any advise.

            Comment


            • #21
              Pascal,

              Have you been able to improve the matches and got similar results as in the paper ?

              Comment


              • #22
                Sorry, I have stopped working on this data set. If I continue I will let you know.

                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, Yesterday, 06:37 PM
                0 responses
                8 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, Yesterday, 06:07 PM
                0 responses
                8 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-22-2024, 10:03 AM
                0 responses
                49 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-21-2024, 07:32 AM
                0 responses
                66 views
                0 likes
                Last Post seqadmin  
                Working...
                X