This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Using RSem of Trinity, align_and_estimate_abundance.pl script, all folders for sample/condition is not being created

Hi,

Folders for all the samples are not being created, only 1/4 is created using align_and_estimate_abundance.pl script provided by Trinity, is giving output.

What to do? What must be the error? Can I run the script for each sample/ condition separately?

Thank you. Hope to get a response soon

Anupama

error

rsem quantification trinity
• 1,197 views • link
updated by Ram Ram • written by Roy.anupama11 • Roy.anupama11

Error, fewer reads in file specified with -2 than in file specified with -1

That's pretty clear is it? Some of your fastq file pairs are out of sync. Did you do any manipulation like trimming of paired-end data with a non paired-end aware tool? Please describe how data were obtained and processed before running that analysis here.

• 0 views • link
written by ATpoint ATpoint

Firstly thank you for responding. I did trimming for adapter sequences using Trimmomatic. In the script for quantification, I am using only trimmed paired end data, leaving the upaired data formed after trimming. I have used HISAT2, then STRINGTIE and TRINITY for preparing reference.

• 0 views • link
written by Roy.anupama11 • Roy.anupama11

In the script for quantification, I am using only trimmed paired end data, leaving the upaired data formed after trimming.

What does that mean? Please show code. In any case, unpaired data should be gone after the trimming as this is exactly what causes this issue here. The tool expects paired data, no singleton reads.

• 0 views • link
written by ATpoint ATpoint

0 answers

No answers yet.

Log in to answer this question.

More posts like this