Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Bowtie and ERANGE

    Hi,

    We're planning on using the Mortazavi ERANGE package for quantifying gene expression in our rna-seq data. An aligner must be used for this process, and we are going to use Bowtie.

    The ERANGE documentation describes how Bowtie should be used to create output that is compatible with the ERANGE process, but part of it is confusing me:
    We use the following settings to emulate Eland output:

    $BOWTIEDIR/bowtie zzz -v 2 -k 11 -m 10 -t --best -f s1.32mer.query.txt --unfa s1.unmapped.fa --maxfa s1.repeat.fa s1.zzz.bowtie.txt

    where zzz is the genome prefix that you gave when building the
    genome. In particular, we ask bowtie to map all multireads up
    to 11 ("-k") with up to 2 mismatches ("-v" and "--best"), however
    we will only import all multireads up to 10x multiplicity ("-m").

    My understanding is that the "-m 10" will make the "-k 11" non-sensical, since any reads with more than 10 matches will be discarded, therefore making it impossible to report up to 11.

    Have I misunderstood?

    Thanks!

  • #2
    Specifying "-m 10" without "-k 11" would cause (a) all alignments for reads with more than 10 alignments to be suppressed, and (b) one random alignment to be reported for reads with 1 to 10 alignments. Specifying both (equivalent to "-m 10 -a") causes (a) all alignments for reads with more than 10 alignments to be suppressed, and (b) *all* alignments to be reported for reads with 1 to 10 alignments

    Hope that helps,
    Ben

    Comment


    • #3
      Hi Ben,

      Thanks for the explanation.

      Does this mean that "-k 10" or "-a" would have been more appropriate here, since "-k 11" would be as redundant as "-k 1234" ?

      I suppose I'm trying to determine if I'm missing something from the ERANGE documentation. Does their explanation in blue make sense?

      Thanks!

      Comment


      • #4
        I think you should also try out TopHat and Cufflinks as it does much the same thing as Erange.

        In addition, TopHat and Cufflinks were (at least for me) much easier to set up, and the latest versions seem to have sorted out the most significant bugs. Never having succesfully run Erange I can't tell you if there are speed differences, but I think cufflinks is probably much faster than Erange (I believe much of the code is in Python, but there might be C/C++ code for the speed sensitive parts).

        Looking at the list of Cufflinks contributors I see Ali Mortazavi and Barbara Wold, so Cufflinks might even have design elements in common with Erange. Ben, you must know a lot more about this?

        Comment


        • #5
          @Bio.X2Y:

          Yes, their explanation makes sense, (incidentally, it's similar to what TopHat does). I would have said "-a -m 10", but, as I say, the effect is exactly the same as "-k 11 -m 10" or "-k 1234 -m 10".

          @Thomas Doktor

          Yes, TopHat and Cufflinks are highly recommended. Ali and Barbara and Cole (and others) collaborate very tightly on them. Please check them out! If you have questions about those tools, feel free to post them here. Cole reads this forum with some regularity.

          Thanks,
          Ben

          Comment


          • #6
            Hi,
            Could anyone explain to me where these input files come from
            s1.32mer.query.txt
            s1.unmapped.fa
            s1.repeat.fa

            I only have the fastq files downloaded from GEO-database and the bowtie d_melanogaster index files.
            Thanks

            Comment

            Latest Articles

            Collapse

            • 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
            • seqadmin
              Multiomics Techniques Advancing Disease Research
              by seqadmin


              New and advanced multiomics tools and technologies have opened new avenues of research and markedly enhanced various disciplines such as disease research and precision medicine1. The practice of merging diverse data from various ‘omes increasingly provides a more holistic understanding of biological systems. As Maddison Masaeli, Co-Founder and CEO at Deepcell, aptly noted, “You can't explain biology in its complex form with one modality.”

              A major leap in the field has
              ...
              02-08-2024, 06:33 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, 02-28-2024, 06:12 AM
            0 responses
            30 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 02-23-2024, 04:11 PM
            0 responses
            74 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 02-21-2024, 08:52 AM
            0 responses
            87 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 02-20-2024, 08:57 AM
            0 responses
            71 views
            0 likes
            Last Post seqadmin  
            Working...
            X