Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • #31
    @Ryan: I totally agree with you and that was my fear before I run alignments with bowtie2. However, i reasoned that since I'm dealing with bacterial transcriptome, the complexity of the reference genome is way smaller then that of human RNA-seq. So, I didn't expect a huge number of multimappers as defined by the --score-min function (i.e. reads that have more then one valid alignment).

    I turns out that I was probably right: I lose about 3% of my total reads by filtering for multimappers with XS: tag which is not dramatic. However, one library was more of a mess... It had a deepest sequencing depth (61 M reads) but I lose 14% on "multimappers".

    In my opinion, bowtie1 had somewhat more flexible options. Reporting could be controlled with "-a --best --strata" or "--best -k 1".

    TP

    Comment


    • #32
      @ThePresident: That makes more sense then. Yeah, bowtie1 did have its advantages.

      Comment


      • #33
        Originally posted by ThePresident View Post
        I calculated RPKM values manually, still didn't compare it with Cufflinks RPKMs, but visually I have a lognormal distribution. Sounds logical for me...

        TP
        If/when you do compare to cufflinks RPKM, I would be interested to hear
        the results.

        I also see a (skewed) lognormal computing RPKMs via htseq, but I
        see a bimodal lognormal from cufflinks FPKMs. [Yes, I mean RPKM
        in the first and FPKM in the latter--an oversimplification on my part
        using these paired-end data.]

        Comment


        • #34
          Can someone please explain where the "10^9" in the calculation comes from? Why isn't it 10^6?

          Thanks in advance!

          Comment


          • #35
            Originally posted by raymb View Post
            Can someone please explain where the "10^9" in the calculation comes from? Why isn't it 10^6?

            Thanks in advance!
            L is usually measured in base pairs but the RPKM/FPKM metric requires things in KB, so the 10^9 comes from that conversion. If you input a length measure in KB then just use 10^6.

            Comment


            • #36
              I'm curious about calculating RPKM / FPKM from counts on the gene level. What would you use as "gene length" for a gene with multiple isoforms...? The longest exon chain (i.e. longest transcript isoform)? Or a median transcript length based on all the annotated isoforms of a given gene...?

              Comment


              • #37
                Usually one uses the union of all exons and gets the length of that. At least unless you want to determine an expected gene length given alignment data (some tools can do that).

                Comment


                • #38
                  Can we get the FPKM value only without through a long process of Cuffdiff? I do have 100 samples and each sample has 3 biological replicates. I am only interested in FPKM value which will later use for network analysis. Thanks

                  Comment


                  • #39
                    Yes, but there's more than one way to get and even define FPKMs. You could just get the counts with something like featureCounts and then convert them to FPKM (either before or after normalization). That's the fastest method. You could also get expected counts with something like eXpress and convert those to FPKMs. That would take longer and given you different (through pretty similar) results. I can conceive of still other routes.

                    Comment

                    Latest Articles

                    Collapse

                    • seqadmin
                      Recent Advances in Sequencing Technologies
                      by seqadmin







                      Innovations in next-generation sequencing technologies and techniques are driving more precise and comprehensive exploration of complex biological systems. Current advancements include improved accessibility for long-read sequencing and significant progress in single-cell and 3D genomics. This article explores some of the most impactful developments in the field over the past year.

                      Long-Read Sequencing
                      Long-read sequencing has...
                      12-02-2024, 01:49 PM
                    • seqadmin
                      Genetic Variation in Immunogenetics and Antibody Diversity
                      by seqadmin



                      The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...
                      11-06-2024, 07:24 PM

                    ad_right_rmr

                    Collapse

                    News

                    Collapse

                    Topics Statistics Last Post
                    Started by seqadmin, 12-02-2024, 09:29 AM
                    0 responses
                    150 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 12-02-2024, 09:06 AM
                    0 responses
                    51 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 12-02-2024, 08:03 AM
                    0 responses
                    42 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 11-22-2024, 07:36 AM
                    0 responses
                    74 views
                    0 likes
                    Last Post seqadmin  
                    Working...
                    X