Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • BLAST parameters for large datasets

    I am mapping contig sequences with plant database (size 2 GB). I am using tblastx and blastx search tool. Can anyone suggest what parameter for blastx/tblastx i should use for mapping with database?. I am having about 1 million contig sequences.
    parameter in case for e value, length to be covered by query sequence etc.

  • #2
    For e-value it depends on how many false positives you can stand. E=1 means that you would expect 1 of your contigs to map to the database by random chance. And many others to be weakly associated to the database. You may wish to cast such a wide net and then take the time to track down the false positives. Personally I use E of 10 to the minus 6 as a first pass. I may miss some interesting correlations but at the same time I can focus in on what I do obtain since those results are going to be strong.

    Comment


    • #3
      The culling_limit and max_target_sequences should also affect the speed. culling_limit discards hits that have n number of reads that are better than it. max_target_sequences are the number of hits to retain. I generally set the max_target_sequences > culling_limit but you can tune it to what you want.

      Hop this helps!

      Comment


      • #4
        why max_target_seqs > culling_limit?

        It isn't immediately obvious to me how these two parameters are related. The reason I believe Twaddlac suggested these parameters is because for any given query, you'd like to retrieve at least as many hits as the number of better matches than the alignment in question. This seems like a recursive definition, because when looking across many species most of the hits for a query will not be the best hit, so you'll just end up printing all the hits no matter what...

        I'm pretty sure I'm wrong about something here, but cannot rectify my current understanding of how these parameters function:

        -max_target_seqs <Integer, >=1>
        Maximum number of aligned sequences to keep

        Example 1. max_target_sequences = 5
        If target1 is the best match for query1, but targets 2,3,4,5 are also good matches - this will retrieve the 5 alignments in a target-centric manner.

        -culling_limit <Integer, >=0>
        If the query range of a hit is enveloped by that of at least this many
        higher-scoring hits, delete the hit

        Example 2. culling_limit = 5
        Queries (q1, q2, q3, q4, q5) are all good matches for target 1. However, if any of these queries has higher scoring matches to 5+ other targets, then the query->target1 alignment hit will be deleted.

        If the query matches X genes better than a given target, therefore, the query-target relationship is deleted, it seems to me this has nothing to do with number of hits you want to retrieve at the end.

        Can anyone explain if/how these two parameters are related? Thanks in advance.
        Last edited by flow_science; 05-06-2015, 04:19 PM.

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