Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • feederbing
    Member
    • Sep 2011
    • 11

    Building bfast/btestindexes index for 15% divergence

    I am trying to determine whether BFAST is appropriate for mapping 101 bp unpaired illumina reads to a reference with expected divergence of 15% (http://seqanswers.com/forums/showthread.php?t=13871 has related info on my problem). Primarily, I'm hoping Nils can advise me on whether or not BFAST is suited for this type of job. Assuming that it is, below I describe what I have tried and where I am stuck.

    I am at the stage where I am trying to build an index for BFAST, using btestindexes, but I am stuck at a lack of understanding of the output of btestindexes in "evaluate" mode. I understand the concept of spaced seeds, what I don't understand is how to interpret the output of btestindexes. I have looked at the other four threads here that mention btestindexes, and I have read the supplementary info from the BFAST paper.

    Following the advice in section 6.1 of the bfast-book, I get k+2=21 for a genome of size 2.4G.

    I then ran an index search:
    btestindexes -A 0 -a 0 -S 10000 -s 10 -r 101 -M 20 -n 10 -l 21 -w 31
    I used -M 20 because I think my data will contain some unique matches out to 20% divergence. I used -n 10 to get 10 masks, expecting that the evaluation run of btestindexes will indicate how many I need.

    The resulting masks are
    111111111111111111111
    11111011001101101111111111
    1001111111101110011101001010111
    11001111111011101101101101011
    110011111001111110111010101011
    10011110101111101111101101101
    11111101011111101101011111
    11111010101101110111111111
    1111111101111111111111
    11110001011111011111011110011

    I then ran the evaluation:
    btestindexes -A 1 -a 1 -S 10000 -r 101 -M 10 -f filename

    Looking at the output of the evaluation is where I am stuck. Clearly it is a table with one row per mask and one column per mismatch count deom 0 to 20 (plus a column for a deletion but let's ignore that). There is also a column labeled "CE" which is always zero (perhaps "cumulative error"?). The values are undoubtedly probabilities, but probabilities of what? My initial assumption was that row m gave the probability, for the combination of masks 1 thru m, that a homologous read would be discovered using that set of seeds. This assumption is apparently wrong, because if I shuffle the list of masks, I don't get the same results in the final row.

    Looking back at section 6.1, it advises that I "seelect the minimum number of masks sufficient to tolerate" my desired accuracy. But it gives no advice on how to interpret the output so as to make this decision.

    I have also hunted through the supplement and the distribution to see if there is are any masks recommended for this type of divergence. The supplement states that there distribution includes mask sets for reads up to L=100. I have distribution 0.6.5a from sourceforge and I haven't been able to find them. (http://sourceforge.net/projects/bfas...-0.6.5a.tar.gz).

    At this point, I'm just hoping to get some reassurance that BFAST will be useful for this problem.
  • feederbing
    Member
    • Sep 2011
    • 11

    #2
    Originally posted by feederbing View Post
    I then ran the evaluation:
    btestindexes -A 1 -a 1 -S 10000 -r 101 -M 10 -f filename
    Retracing my steps, I see that should be -A 0 (nt space instead of color space). I've rerun the same masks now. Output is in a different format. Am trying to see if it makes more sense now.

    Comment

    Latest Articles

    Collapse

    • seqadmin
      New Genomics Tools and Methods Shared at AGBT 2025
      by seqadmin


      This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

      The Headliner
      The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
      03-03-2025, 01:39 PM
    • seqadmin
      Investigating the Gut Microbiome Through Diet and Spatial Biology
      by seqadmin




      The human gut contains trillions of microorganisms that impact digestion, immune functions, and overall health1. Despite major breakthroughs, we’re only beginning to understand the full extent of the microbiome’s influence on health and disease. Advances in next-generation sequencing and spatial biology have opened new windows into this complex environment, yet many questions remain. This article highlights two recent studies exploring how diet influences microbial...
      02-24-2025, 06:31 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 03-20-2025, 05:03 AM
    0 responses
    18 views
    0 reactions
    Last Post seqadmin  
    Started by seqadmin, 03-19-2025, 07:27 AM
    0 responses
    20 views
    0 reactions
    Last Post seqadmin  
    Started by seqadmin, 03-18-2025, 12:50 PM
    0 responses
    19 views
    0 reactions
    Last Post seqadmin  
    Started by seqadmin, 03-03-2025, 01:15 PM
    0 responses
    186 views
    0 reactions
    Last Post seqadmin  
    Working...