Question: RSEM on single replicate with multiple fasta files
0
gravatar for gozrom
21 months ago by
gozrom50
gozrom50 wrote:

Hi, there,

I'm running RSEM on 12 replicates each replicate has 8 fasta files, paired-end reads.

This is inside 2 For loops with i from 1 to 4 and j from 1 to 4

rsem-calculate-expression -p 8 --bowtie2 --paired-end --strandedness forward --append-names --estimate-rspd --output-genome-bam ~/RNAseq/RSEM_fasta/control-"$i"*_L00"$j"_R1_001.fastq.gz ~/RNAseq/RSEM_fasta/control-"$i"*_L00"$j"_R2_001.fastq.gz  ~/RNAseq/ref/mouse_ensemble ~/RNAseq/RSEM/Control"$i"/sample_"$j"

This command will generate 4 fasta files of single replicate,

is there a way to create a combined one file that includes results from all 4 fasta files for a single replicate

I used this previously For loop with i running from 1 to 12 and j running from 1 to 4 but then I realized that each time it creates a new file with the same name and doesn't really combine with the previous one. So at the end I have files for the last two fasta files for each sample (that's my guess).

rsem-calculate-expression -p 8 --bowtie2 --paired-end --strandedness forward --append-names --estimate-rspd --output-genome-bam ~/RNAseq/RSEM_fasta/*S"$i"*_L00"$j"_R1_001.fastq.gz ~/RNAseq/RSEM_fasta/*-"$i"*_L00"$j"_R2_001.fastq.gz  ~/RNAseq/ref/mouse_ensemble ~/RNAseq/RSEM/sample_"$i"

I still want to combine all fasta files for each replicate into one file and then all those files for 12 replicates into one file

Is there a nice way to do it with RSEM?

Thank you.

rna-seq • 818 views
ADD COMMENTlink modified 21 months ago • written 21 months ago by gozrom50
1
gravatar for swbarnes2
21 months ago by
swbarnes27.2k
United States
swbarnes27.2k wrote:

Why can't you cat the files together before putting them in RSEM? Catting gzipped fastq files works fine.

ADD COMMENTlink written 21 months ago by swbarnes27.2k

Thanks, that seems to work...

ADD REPLYlink modified 21 months ago • written 21 months ago by gozrom50
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: 953 users visited in the last hour