Dear friends,
i was doing the trinity assembly, but i made a small mistake that for the --library type parameter, i used the FR, but actually it was RF, till i did the blast that i found the known gene matched to the reverse strand, cause the trinity take times and CPU, can i just reverse the Trinity.fasta file, and no need redo it again? does it works???
thank you so much for your help.
•
written
by
biologo •
1 answer
It might be low quality data, but please tell us step by step what was done exactly. Also, you can run it through a trinity workflow for example on Galaxy main.
How many reads do you have ?
Another assembler - CLC is easy and quite good, or soapdenovo-trans has a decent reputation - might help too.
Also, just map your transcripts to the genome (I like Gmap) then visualize with Jbrowse or IGV etc. Are multiple redundant transcript fragments present ?
Log in to answer this question.
-
What is correct percentage of sense-stranded transcripts in directional de novo assemblies ?
written by Lada •
-
converting RF strand-specific mRNA sequencing data to FR stand-specific
written by upadhyay.maulik •
-
why forward reads were much more than reverse reads In strand specific RNAseq data?
written by yuaner •
-
trinity assembly error
written by biologo •
-
Strand specific pair end sequencing
written by Bioinfonext
-
Change paired end library from FR to RF
written by Jake Warner
-
Trinity.fasta file missing "does not exist in the out directory:
written by bilal.sarwar •
-
Trinity Error with the the flag of --run_as_paired
written by seta
-
de novo assembly v.s. genome guided assembly (Trinity) -- better N50 from de novo than genome guided
written by Yang Li
-
Futher analysis on Trinity fasta output
written by Whoknows
No, you'll need to redo it, but you should get a much better assembly.
Colin
i got it, thank you, but your extral advice remind me that actually the quality is not good, even i isolate the llongest unigene, i got over 430,000 transcriptome.
Thats a lot of transcripts, and they are probably highly fragmented. I'm sure you did read trimming before, right ? Check your data with FASTQC before and after trimming too.
yes, i did. trim the adaptor and filter some short or low quality reads, but the trinity result seems no better options, thank you for your patience and your kind reply.
You might try out Bridger
yes, i also considered that, but my colleague suggested me use ingap-cdg, have you ever use that, he told me it really works, and i also did the test, and only left 7,000,000 transcripts, seems good, but i am still doubt about the result.