Try using the "-U" option.
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
with the -U option added, it finished within 1 min!
but it means BFAST will not consider the results for "mated-pair". How much it will affect the final result (compared with not using -U option)?
If I use the -U option, may I use the SAMTools or other software to further filter the result before SNP discovery?
Thanks!
Comment
-
Originally posted by NanYu View Postwith the -U option added, it finished within 1 min!
but it means BFAST will not consider the results for "mated-pair".
Originally posted by NanYu View PostHow much it will affect the final result (compared with not using -U option)?
Originally posted by NanYu View PostIf I use the -U option, may I use the SAMTools or other software to further filter the result before SNP discovery?
Thanks!
Comment
-
Originally posted by nilshomer View PostI am not sure why you are concerned. Did you look at the data and see something odd?
I'm wondering in the postprocess step, if BFAST will do something like: choose the best for F and second best for R within a pair because that is within the expected distance, yet the best for R is in another chromsome.
I know it is a rare scenario and maybe I should not worry about it.
Comment
-
Originally posted by Chipper View PostThe problem seems to be the estmation of distances and standar deviations. If you set these parameters it will be much faster.
you mean I should provide a value for "-S"?
We are expecting ~1.5kb, so should I set "-S 3000" (allow the distance to be a bit bigger).
Thanks!
Comment
-
Originally posted by nilshomer View PostYes, that is a possibility with any pairing step.
As mentioned above, if we are expecting ~1.5kb (insertion size, from the lab protocol) , do you think if I set "-S 3000" (very loose criteria) is an OK value for "-S"?
I saw the "-v 160 -s 20 " in #3 but could not find them in the BFAST manual (postprocess section). Could you please let me know what they are?
I really appreciate your help in this forum.
Comment
-
Originally posted by nilshomer View PostIt may take a very long time to run. Note that "-S" is the standard deviation, not the mean (etc.).
The best way, as a scientist, to find out these answers is to try it out for yourself! I cannot predict the behavior
It did takes much longer time to run (compared to use "-U") option.
There is no option for me to tell BFAST at the postprocess step the expected insertion size (or range of the insertion size) for pair-end sequences.
Maybe what I should do is to use "-U" and use a perl script to look for column #8 in the sam file and separate the results according to the value of #8 (and other information). By doing so, it means this information (expected insertion size) is wasted when using BFAST to map the paired-end sequences.
Comment
-
Originally posted by nilshomer View PostSee "-v" and "-s". Make sure you have latest BFAST.
Nils
The version I have is 0.6.5a and the bfast-book.pdf I am reading also shows "version 0.6.5a" but I could not find -v option and -s is common parameter and has nothing to do with insertion size.
"-s INTEGER, --startReadNum=INTEGER
Specifies the first read in which to process. This may be useful when distributing a large data set across a cluster."
Maybe I didn't get things right?
Comment
-
Originally posted by nilshomer View PostGood point, the manual and binaries are not in sync. Try "bfast postprocess -h".
Thanks for your help! I will try them now.
-v insertSizeAvg Specifies the mean insert size to use when rescuing
-s insertSizeStdDev Specifies the standard deviation of the insert size to use when rescuing
Comment
-
Hi. everyone.
Should we transform the fastq file from solid2fastq (Bfast) before mapping with bwa?
I mean from "0123123" to "ACTGACT", something likewise plus some other format issues?
Before, I proced without any transformation, then got malformed .sai file and some "segmentation fault" which I haven't specified yet.
Any suggestions? Thank you
Comment
Latest Articles
Collapse
-
by seqadmin
The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...-
Channel: Articles
05-06-2024, 07:48 AM -
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...-
Channel: Articles
04-22-2024, 07:01 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 06:35 AM
|
0 responses
14 views
0 likes
|
Last Post
by seqadmin
Yesterday, 06:35 AM
|
||
Started by seqadmin, 05-09-2024, 02:46 PM
|
0 responses
19 views
0 likes
|
Last Post
by seqadmin
05-09-2024, 02:46 PM
|
||
Started by seqadmin, 05-07-2024, 06:57 AM
|
0 responses
17 views
0 likes
|
Last Post
by seqadmin
05-07-2024, 06:57 AM
|
||
Started by seqadmin, 05-06-2024, 07:17 AM
|
0 responses
19 views
0 likes
|
Last Post
by seqadmin
05-06-2024, 07:17 AM
|
Comment