Question: Trinity for single end library distributed among different lanes.
0
gravatar for waqaskhokhar999
7 months ago by
waqaskhokhar99940 wrote:

RNA-seq data is single ended and distributed among different lanes. I am trying to assembel the transcriptome using this command:

trinity/trinityrnaseq-Trinity-v2.8.3/Trinity --seqType fq --max_memory 10G --single SRR3460466.fastq,SRR3460467.fastq,SRR3460468.fastq,SRR3460469.fastq,SRR3460470.fastq,SRR3460471.fastq,SRR3460472.fastq,SRR3460473.fastq,SRR3460474.fastq,SRR3460475.fastq,SRR3460476.fastq,SRR3460477.fastq

But I am getting this error

Error, not recognizing read name formatting: [SRR3460466.1]

Command line summary:

waqas@waqas-Hybris:/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq$ trinity/trinityrnaseq-Trinity-v2.8.3/Trinity --seqType fq --max_memory 10G --single SRR3460466.fastq,SRR3460467.fastq,SRR3460468.fastq,SRR3460469.fastq,SRR3460470.fastq,SRR3460471.fastq,SRR3460472.fastq,SRR3460473.fastq,SRR3460474.fastq,SRR3460475.fastq,SRR3460476.fastq,SRR3460477.fastq


     ______  ____   ____  ____   ____  ______  __ __
    |      ||    \ |    ||    \ |    ||      ||  |  |
    |      ||  D  ) |  | |  _  | |  | |      ||  |  |
    |_|  |_||    /  |  | |  |  | |  | |_|  |_||  ~  |
      |  |  |    \  |  | |  |  | |  |   |  |  |___, |
      |  |  |  .  \ |  | |  |  | |  |   |  |  |     |
      |__|  |__|\_||____||__|__||____|  |__|  |____/

    Trinity-v2.8.3



Single read files: $VAR1 = [
          'SRR3460466.fastq',
          'SRR3460467.fastq',
          'SRR3460468.fastq',
          'SRR3460469.fastq',
          'SRR3460470.fastq',
          'SRR3460471.fastq',
          'SRR3460472.fastq',
          'SRR3460473.fastq',
          'SRR3460474.fastq',
          'SRR3460475.fastq',
          'SRR3460476.fastq',
          'SRR3460477.fastq'
        ];
Trinity version: Trinity-v2.8.3
** NOTE: Latest version of Trinity is Trinity-v2.8.4, and can be obtained at:
    https://github.com/trinityrnaseq/trinityrnaseq/releases

Wednesday, November 7, 2018: 16:00:20   CMD: java -Xmx64m -XX:ParallelGCThreads=2  -jar /media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/trinity/trinityrnaseq-Trinity-v2.8.3/util/support_scripts/ExitTester.jar 0
Wednesday, November 7, 2018: 16:00:20   CMD: java -Xmx4g -XX:ParallelGCThreads=2  -jar /media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/trinity/trinityrnaseq-Trinity-v2.8.3/util/support_scripts/ExitTester.jar 1


----------------------------------------------------------------------------------
-------------- Trinity Phase 1: Clustering of RNA-Seq Reads  ---------------------
----------------------------------------------------------------------------------

---------------------------------------------------------------
------------ In silico Read Normalization ---------------------
-- (Removing Excess Reads Beyond 200 Coverage --
---------------------------------------------------------------

# running normalization on reads: $VAR1 = [
          [
            '/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460466.fastq',
            '/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460467.fastq',
            '/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460468.fastq',
            '/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460469.fastq',
            '/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460470.fastq',
            '/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460471.fastq',
            '/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460472.fastq',
            '/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460473.fastq',
            '/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460474.fastq',
            '/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460475.fastq',
            '/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460476.fastq',
            '/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460477.fastq'
          ]
        ];


Wednesday, November 7, 2018: 16:00:20   CMD: /media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/trinity/trinityrnaseq-Trinity-v2.8.3/util/insilico_read_normalization.pl --seqType fq --JM 10G  --max_cov 200 --min_cov 1 --CPU 2 --output /media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/trinity_out_dir/insilico_read_normalization   --max_pct_stdev 10000  --single /media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460466.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460467.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460468.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460469.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460470.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460471.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460472.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460473.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460474.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460475.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460476.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460477.fastq
CMD: seqtk-trinity seq -A /media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460466.fastq >> single.fa
Error, not recognizing read name formatting: [SRR3460466.1]

If your data come from SRA, be sure to dump the fastq file like so:

    SRA_TOOLKIT/fastq-dump --defline-seq '@$sn[_$rn]/$ri' --split-files file.sra 

Error, cmd: seqtk-trinity seq -A /media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460466.fastq >> single.fa died with ret 512 at /media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/trinity/trinityrnaseq-Trinity-v2.8.3/util/insilico_read_normalization.pl line 762.
Error, cmd: /media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/trinity/trinityrnaseq-Trinity-v2.8.3/util/insilico_read_normalization.pl --seqType fq --JM 10G  --max_cov 200 --min_cov 1 --CPU 2 --output /media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/trinity_out_dir/insilico_read_normalization   --max_pct_stdev 10000  --single /media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460466.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460467.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460468.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460469.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460470.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460471.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460472.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460473.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460474.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460475.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460476.fastq,/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/SRR3460477.fastq died with ret 512 at trinity/trinityrnaseq-Trinity-v2.8.3/Trinity line 2684.
    main::process_cmd("/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/trinity"...) called at trinity/trinityrnaseq-Trinity-v2.8.3/Trinity line 3230
    main::normalize("/media/waqas/waqas_data/19_Accessions_w/new/BUR/fastq/trinity"..., 200, ARRAY(0x56443710b4b0)) called at trinity/trinityrnaseq-Trinity-v2.8.3/Trinity line 3177
    main::run_normalization(200, ARRAY(0x56443710b4b0)) called at trinity/trinityrnaseq-Trinity-v2.8.3/Trinity line 1319
rna-seq assembly • 286 views
ADD COMMENTlink modified 7 months ago by h.mon26k • written 7 months ago by waqaskhokhar99940
0
gravatar for h.mon
7 months ago by
h.mon26k
Brazil
h.mon26k wrote:

Reading the error message carefully, you will see:

If your data come from SRA, be sure to dump the fastq file like so:

SRA_TOOLKIT/fastq-dump --defline-seq '@$sn[_$rn]/$ri' --split-files file.sra

See also the Trinity FAQ.

ADD COMMENTlink written 7 months ago by h.mon26k

I have used SRA_TOOLKIT/fastq-dump to dump the fastq files, since my data is single ended, do i still need to use --split-files as this option is used for paired ended data

ADD REPLYlink written 7 months ago by waqaskhokhar99940
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 2.3.0
Traffic: 1020 users visited in the last hour