SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Shortread Package INSTALL Wei-HD Bioinformatics 21 04-30-2014 12:33 AM
QA with ShortRead dawe Bioinformatics 5 02-26-2014 04:44 AM
ShortRead R package Aligned Read Output?? nazanin Bioinformatics 0 01-31-2011 12:34 AM
Problem with ShortRead (R, Bioconductor) function readIntensities marlei Bioinformatics 1 11-26-2009 11:23 AM
ShortRead: a bioconductor package for input, quality assessment and exploration Chien-Yuan Chen Literature Watch 0 11-09-2009 02:36 PM

Reply
 
Thread Tools
Old 08-27-2008, 07:52 AM   #1
bioinfosm
Senior Member
 
Location: USA

Join Date: Jan 2008
Posts: 482
Default ShortRead package

Hi,

Any experiences with the ShortRead R package?
http://bioconductor.org/packages/2.3...c/Overview.pdf

I can finally make it to run, but people's experience using it, tweaking it, and QC.ing data from the functions will be helpful..

sm
bioinfosm is offline   Reply With Quote
Old 11-28-2008, 12:30 PM   #2
bioinfosm
Senior Member
 
Location: USA

Join Date: Jan 2008
Posts: 482
Default

<bump>

Also more stuff coming out in R for NGS data
http://bioconductor.org/workshops/2008/SeattleNov08/
bioinfosm is offline   Reply With Quote
Old 04-30-2014, 12:37 AM   #3
archi
Member
 
Location: India

Join Date: Apr 2014
Posts: 16
Default

I'm getting the following errors in using the readFastq function
readFastq("home/linux/sratoolkit.2.34-2-centos_linux64/bin",pattern="SRR039194_1.fastq")
Error: Input/Output
no input files found
dirPath: home/linux/sratoolkit.2.34-2-centos_linux64/bin
pattern: SRR039194_1.fastq

any suggestions on how to fix this?
please reply asap
thanks!
archi is offline   Reply With Quote
Old 04-30-2014, 04:52 AM   #4
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,800
Default

You may be missing a leading "/" in your path home/linux/sratoolkit.2.34-2-centos_linux64/bin. Are the sequence files in this path?
GenoMax is offline   Reply With Quote
Old 05-23-2014, 02:36 AM   #5
archi
Member
 
Location: India

Join Date: Apr 2014
Posts: 16
Default

I'm able to generate the fastq version of the sra file in the terminal.But how do I save its copy in my computer??

please reply asap
archi is offline   Reply With Quote
Old 05-23-2014, 03:52 AM   #6
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,800
Default

Have you looked in the same directory where your SRA file was?
GenoMax is offline   Reply With Quote
Old 05-24-2014, 02:46 AM   #7
archi
Member
 
Location: India

Join Date: Apr 2014
Posts: 16
Default

okay,got it.
thanks
archi is offline   Reply With Quote
Old 05-24-2014, 02:49 AM   #8
archi
Member
 
Location: India

Join Date: Apr 2014
Posts: 16
Default

While using ShortRead package to read a fastq file what are the exact arguments to be put in the readFastq function?
and how do we segregate from so many reads only those which have a particular sequence?
archi is offline   Reply With Quote
Old 05-24-2014, 03:30 AM   #9
mastal
Senior Member
 
Location: uk

Join Date: Mar 2009
Posts: 667
Default

If you type '?readFastq' at the R prompt on the commandline, you should get a help page that explains the usage for the function.

Download the manuals for the ShortRead package from:

http://www.bioconductor.org/packages...ShortRead.html

I think that readFastq() can only be used to read in complete fastq files, but then you can apply various filters to the data to get only the reads you want. See the documentation.

Last edited by mastal; 05-24-2014 at 04:07 AM.
mastal is offline   Reply With Quote
Old 06-01-2014, 10:25 AM   #10
archi
Member
 
Location: India

Join Date: Apr 2014
Posts: 16
Default

I have around 4 Lac lines in my fastq file.So while using readFastq I'm constantly getting an error message of unexpected number of lines.
Can someone please suggest how do I manipulate my reads now,my final aim is to get a set of reads which have a particular sequence at the end.

please reply asap
archi is offline   Reply With Quote
Reply

Tags
bioconductor, quality, solexa

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 06:05 AM.


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