Tweet
Anybody know a best fit method for determining mask/indexes for the mouse genome? I saw on another post they recommended the 25 bp/50 bp sets from the manual supplementation for a 35 bp plant sample. I figured I could just apply the same masks to my own. I'm running both the 25 and 50 specified index sets right now. I imagine I might have to use 'btestindex' Any help is appreciated and thanks!
-
-
Since the genome size of mouse and human is close enough in size, you can use the indexes for human recommended in the manual. Of course, the polymorphism rate is ~10 times greater in mouse but that may be why you are using BFAST (?).Originally posted by JohnK View Post -
Affirmative- I'm using BFAST primarily for concordance with SOLiD data analysis. It'd be nice to figure out btestindex... I'll keep looking at the manual and maybe I'll get it down.Originally posted by nilshomer View PostComment
-
Interpreting Btestindex Output?
Nils,
Below is the command and sample output of Btest index, I am guessing that the two way table has the sensitivity numbers, given the error model defined by row and column.
Can you expand on what BP:0,1:del,2i1,2i2 etc mean?
Is it one insertion of size 2 or 2 insertions of size 1 etc?
Also is there a way to get the cumulative answer (i.e after being filtered by all masks?)
$ btestindexes -a 1 -r 80 -A 0 -S 10000 -I 2 -M 8 -f maskset1
...MASKS IN BFAST FORMAT
1111111111111111
1011001101010101100111011001
11010101010011010010100101000111
1110101001101111101111
N Masks = 1
MM BP:0 1=del 2i1 2i2
0 1.000 1.000 1.000 1.000
1 1.000 1.000 1.000 1.000
2 1.000 1.000 1.000 1.000
3 1.000 1.000 1.000 1.000
4 1.000 1.000 0.999 1.000
5 0.999 0.995 0.994 0.992
6 0.994 0.978 0.975 0.975
7 0.976 0.942 0.937 0.933
8 0.936 0.890 0.882 0.875
:::::::::::::::::::::::::::::::::::::::::
:::::::::::::::::::::::::::::::::::::::::Comment
-
Also do you have any mask suggestions vis a vis ion torrent ? (It may be moot at this point given the number of reads, but may be some day in the future.)Comment
-
There is no manual for this and I will not support it, but read the supplemental materials in the original BFAST paper for more info. Use the recommended masks in the manual's appendix if you are unsure.Originally posted by cdry7ue View Post
For BFAST, I would recommend using a smaller key size and width, but this is at the cost of run time. If you are unsure of what this means or how to tune this, use the recommended masks. Two other mappers that may be useful for Ion Torrent are SHRiMP2 or TMAP (available and documented on the Ion Torrent website).Originally posted by cdry7ue View PostComment
-
Love this - good data definitely starts from good input, and poor input can only give relatively poor data. I particularly like the mention of Nanodrop/absorbance based methods for quantification. It's such a toss up if you'll get an accurate reading or what amounts to a randomly generated number, and a lot of library/sequencing related issues can be traced back to poor quant.07-01-2026, 11:43 AM
-
I’m not a sequencing expert. I’m a purification scientist who uses NGS to evaluate workflows my group develops. With this perspective, we think about the sample first and the NGS workflow second. The sequencer is an exceptionally honest reporter, but it can only report on what you give it, so whether you get clean, interpretable data from an NGS workflow is largely determined before you begin.
Here are nine questions we think about, in roughly the order they matter, before...
Comment