QuasR error message
1
0
Entering edit mode
ta_awwad ▴ 10
@ta_awwad-11382
Last seen 15 days ago
Frankfurt am Main

Hi everybody, I am trying to align my paired end reads in R using QuasR. everything was easy and smooth until the alignment step. I am using the following line for alignment:

> proj <- qAlign(sampleFile, genome= genomeFile, splicedAlignment= TRUE, alignmentsDir= "bam", clObj = clObj, paired = 'fr')

however there is an error message in the log file says: Error in checkForRemoteErrors(val) : one node produced an error: Error on Galaxy2 processing sample /media//KO_day3_1_1.fq : failed while aligning 25mers

this is happening only when I am applying splicedAlignment=T.

any ideas about what is going on?

thanks much Tamer

QuasR RNA_Seq • 540 views
ADD COMMENT
1
Entering edit mode
@hotz-hans-rudolf-3951
Last seen 15 months ago
Switzerland

Hi Tamer
What is your read length? SpliceMap (the aligner used, when 'splicedAlignment= TRUE') struggles with reads shorter than 50nts, see also: SpliceMap in QuasR package failing when sorting 25mer-alignments


As a workaround, I recommend to use STAR to create the BAM files outside of QuasR and then create the qproject with a sample file listing the BAM files (see QuasR vignete: 5.1 Create a sample file)


Regards, Hans-Rudolf

 

ADD COMMENT

Login before adding your answer.

Traffic: 223 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6