Question: Red error from RSEM Trinity -- Solution: input matched paired end reads
0
gravatar for j_pasook
11 weeks ago by
j_pasook40
j_pasook40 wrote:

Hi,

I've constructed a trinity assembly using ~30M PE read. Then I was trying to estimate abundance using "Align reads and estimate abundance on a de novo assembly of RNA-Seq data" tool in usegalaxy Europe server. I performed this exact tool before using 1M PE read sample before, and it was perfectly fine When it comes to my real 30M PE read data, the tool was run for 8 hours and it showed red error with this messages;

Fatal error: Exit code 2 ()
CMD: touch /data/dnb01/galaxy_db/job_working_directory/004/118/4118029/working/input.fa.bowtie.started
CMD: bowtie-build /data/dnb01/galaxy_db/job_working_directory/004/118/4118029/working/input.fa /data/dnb01/galaxy_db/job_working_directory/004/118/4118029/working/input.fa.bowtie
CMD: touch /data/dnb01/galaxy_db/job_working_directory/004/118/4118029/working/input.fa.RSEM.rsem.prepped.started
CMD: rsem-prepare-reference  --transcript-to-gene-map /data/dnb01/galaxy_db/job_working_directory/004/118/4118029/working/gene_to_trans.map /data/dnb01/galaxy_db/job_working_directory/004/118/4118029/working/input.fa /data/dnb01/galaxy_db/job_working_directory/004/118/4118029/working/input.fa.RSEM
$VAR1 = [
          {
            'left' => '/data/dnb01/galaxy_db/job_working_directory/004/118/4118029/working/paired_left.fq',
            'output_dir' => 'output',
            'right' => '/data/dnb01/galaxy_db/job_working_directory/004/118/4118029/working/paired_right.fq'
          }
        ];
CMD: set -o pipefail && bowtie -q --all --best --strata -m 300 --chunkmbs 512 -X 800 -S -p 1 /data/dnb01/galaxy_db/job_working_directory/004/118/4118029/working/input.fa.bowtie -1 /data/dnb01/galaxy_db/job_working_directory/004/118/4118029/working/paired_left.fq -2 /data/dnb01/galaxy_db/job_working_directory/004/118/4118029/working/paired_right.fq | samtools view -F 4 -S -b | samtools sort -n -o bowtie.bam 
reads processed: 54708014
reads with at least one reported alignment: 52205 (0.10%)
reads that failed to align: 54655809 (99.90%)
Reported 135391 paired-end alignments
CMD: touch bowtie.bam.ok
CMD: convert-sam-for-rsem bowtie.bam bowtie.bam.for_rsem
Number of first and second mates in read 5175620's full alignments (both mates are aligned) are not matched!
Error, cmd: convert-sam-for-rsem bowtie.bam bowtie.bam.for_rsem died with ret: 65280 at /usr/local/tools/_conda/envs/mulled-v1-eea4a55dc36ed6afc8679898c026ba6276e343800300cbf167af08a314bb9f65/bin/align_and_estimate_abundance.pl line 790.

I think it has to be with my data but I don't know why. Can anyone enlighten me about this ? Thanks, j_pasook

ADD COMMENTlink modified 11 weeks ago by Jennifer Hillman Jackson25k • written 11 weeks ago by j_pasook40
1
gravatar for Jennifer Hillman Jackson
11 weeks ago by
United States
Jennifer Hillman Jackson25k wrote:

Hello,

The same number of reads must be entered for both forward and reverse inputs. The sequence identifiers should also be a match and be in the same order.

If you performed some QA and reads were lost, cycling the data through these two NGS: QC and manipulation tools can help you to get matched paired input:

  • FASTQ interlacer paired end reads
  • FASTQ de-interlacer on paired end reads

That said, I am not exactly certain if that is your problem or not. I asked the developers at the IUC Gitter chat (where many of the .eu server admins check-in). Please feel free to join in or just follow the conversation: https://gitter.im/galaxy-iuc/iuc?at=5b5f502b17c942036b839c85

Thanks! Jen, Galaxy team

ADD COMMENTlink modified 11 weeks ago • written 11 weeks ago by Jennifer Hillman Jackson25k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 134 users visited in the last hour