
Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Building bowtie index with mirBase hairpin.fa file  Gators  RNA Sequencing  6  05072015 12:43 PM 
Building bfast index with btestindexes  bre  Bioinformatics  35  09292011 05:34 AM 
strange bowtie index building and mapping problem  Gangcai  Bioinformatics  0  08042010 06:02 PM 
tophatbowtie building index  repinementer  Bioinformatics  1  07172010 11:53 PM 
BWA building index of full human (ensembl) fails  inijman  Bioinformatics  4  12232009 06:00 AM 

Thread Tools 
09072011, 12:39 PM  #1 
Member
Location: USA Join Date: Sep 2011
Posts: 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 bfastbook, 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. 
09072011, 04:02 PM  #2 
Member
Location: USA Join Date: Sep 2011
Posts: 11

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.

Tags 
bfast, btestindexes, homer, index 
Thread Tools  

