SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
SRMA error Ori Bioinformatics 6 08-19-2011 08:09 AM
SRMA - Error - cs TAG mrxcm3 Bioinformatics 3 06-20-2011 10:12 AM
SRMA memory leak Berlinq Bioinformatics 0 03-02-2011 02:46 AM
another question about SRMA rlh Bioinformatics 2 10-23-2010 11:51 PM
SRMA Error dmurdock Bioinformatics 14 10-21-2010 06:02 PM

Reply
 
Thread Tools
Old 10-21-2010, 06:33 PM   #1
rlh
Member
 
Location: China

Join Date: Oct 2010
Posts: 12
Post question about SRMA

i am a freshman in this field, i am trying to use SRMA to realign my data, but i don't know how to generate a ranges file for SRMA, it is a txt file, the content is like:
chr1 1889866 1890066
chr1 12561395 12561595
chr1 34999494 34999694
chr1 43681831 43682031
chr1 74810345 74810545
chr1 74810352 74810552
chr1 89245929 89246129
chr1 143585219 143585419
chr1 144037256 144037456
chr1 150462252 150462452
chr1 156418029 156418229
chr1 169823406 169823606
chr1 227839825 227840025

any body any suggest?i will appreciate your help!
rlh is offline   Reply With Quote
Old 10-21-2010, 07:18 PM   #2
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

Quote:
Originally Posted by rlh View Post
i am a freshman in this field, i am trying to use SRMA to realign my data, but i don't know how to generate a ranges file for SRMA, it is a txt file, the content is like:
chr1 1889866 1890066
chr1 12561395 12561595
chr1 34999494 34999694
chr1 43681831 43682031
chr1 74810345 74810545
chr1 74810352 74810552
chr1 89245929 89246129
chr1 143585219 143585419
chr1 144037256 144037456
chr1 150462252 150462452
chr1 156418029 156418229
chr1 169823406 169823606
chr1 227839825 227840025

any body any suggest?i will appreciate your help!
That looks fine, what was the error?
nilshomer is offline   Reply With Quote
Old 10-21-2010, 07:39 PM   #3
rlh
Member
 
Location: China

Join Date: Oct 2010
Posts: 12
Post

Quote:
Originally Posted by nilshomer View Post
That looks fine, what was the error?
perhaps i didn't say clearly, i just wanna know how to generate that RANGES file. can you offer some advices?
rlh is offline   Reply With Quote
Old 10-21-2010, 07:49 PM   #4
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

Oh, well if you have regions that you are interested in, use a ranges file. Otherwise, you can run it without the ranges file, which will re-align your entire dataset.
nilshomer is offline   Reply With Quote
Old 10-21-2010, 08:07 PM   #5
rlh
Member
 
Location: China

Join Date: Oct 2010
Posts: 12
Post

Quote:
Originally Posted by nilshomer View Post
Oh, well if you have regions that you are interested in, use a ranges file. Otherwise, you can run it without the ranges file, which will re-align your entire dataset.
thankyou for your reply! my database is 56G, if i run re-align in the entire database, the memory it needs is extreme, so i wanna run it with RANGES file.
is there any tools can help me generate RANGES file from my database, if i write it by myself will cost such a long time.
rlh is offline   Reply With Quote
Old 10-21-2010, 08:09 PM   #6
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

Quote:
Originally Posted by rlh View Post
thankyou for your reply! my database is 56G, if i run re-align in the entire database, the memory it needs is extreme, so i wanna run it with RANGES file.
is there any tools can help me generate RANGES file from my database, if i write it by myself will cost such a long time.
It doesn't load the entire reference and reads into memory, it loads them in as needed. I assume by database, you mean reference.
nilshomer is offline   Reply With Quote
Old 10-21-2010, 08:30 PM   #7
rlh
Member
 
Location: China

Join Date: Oct 2010
Posts: 12
Post

Quote:
Originally Posted by nilshomer View Post
It doesn't load the entire reference and reads into memory, it loads them in as needed. I assume by database, you mean reference.
i run a 5G testing data, 64* coverage, -Xmx8G, it print "Error: java heap space", and the records processing stop at about 11000000, someone suggest me to do it with RANGES file.
rlh is offline   Reply With Quote
Old 10-22-2010, 06:05 AM   #8
krobison
Senior Member
 
Location: Boston area

Join Date: Nov 2007
Posts: 747
Default

A Nils (who wrote the program, so he might understand it :-) said, if you want to speed things up you can use the ranges file to specify specific regions of interest. You, and only you, can decide what you find interesting.

Probably you want to pick a bunch of genes of interest & then extract the columns for them from the annotation file off UCSC for the genome you are interested in. That's pretty trivial with Perl/Python/Java/whathaveyou (or even with grep & cut in UNIX) & the time you invest now to learn how to do it will pay off many times over.
krobison is offline   Reply With Quote
Old 10-22-2010, 07:06 AM   #9
rlh
Member
 
Location: China

Join Date: Oct 2010
Posts: 12
Default

Quote:
Originally Posted by krobison View Post
A Nils (who wrote the program, so he might understand it :-) said, if you want to speed things up you can use the ranges file to specify specific regions of interest. You, and only you, can decide what you find interesting.

Probably you want to pick a bunch of genes of interest & then extract the columns for them from the annotation file off UCSC for the genome you are interested in. That's pretty trivial with Perl/Python/Java/whathaveyou (or even with grep & cut in UNIX) & the time you invest now to learn how to do it will pay off many times over.
thank you! i understand it now, i am really appreciate!
rlh is offline   Reply With Quote
Old 10-22-2010, 07:08 AM   #10
rlh
Member
 
Location: China

Join Date: Oct 2010
Posts: 12
Default

i understand what you mean now. thanks very much!
rlh is offline   Reply With Quote
Reply

Tags
srma

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 07:51 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO