Hi all,
as asked above, do you think would it be possible or would there be any significant reason for the inclusion of singletons in the assembly?
Literally, zero read contigs does not make sense.
I assembled 25M, 90bp fastq reads using trinity assembler and obtained 0.13 Million contigs/transcripts.
I have mapped fastq reads to those transcripts i got by assembling the reads using trinity.
Now when i found the number of reads mapped to each contig/transcript using bowtie2 and samtools, i am confused to see few contigs also had 1 read as well as 0 read.
I also observed that these 1 read, 0 read contig had length of 100bp-500bp length. If the reads are not overlapped into contigs why would the length be >90bp. But, my question is why did bowtie2 showed 0 reads for few of the contig still?
All ideas are welcome and help is appreciated.
Thanks in advance.
•
written
by
Prakki Rama
3 answers
Realigning is not really the litmus test you think it is. You should try to see what read-tracking export formats are in Trinity.
Graph-based assemblers assemble kmers, not reads, and will willingly break reads up at the points of ambiguity (repeat junctions). So it is possible some of your reads have their feet planted in two or more contigs.
Thank you Jeremy,lh3. Yes i agree that kmers are used in the assembly. But could i please know why i am getting transcripts without any kmer? I used a kmer of 43 in velvet and used inbuilt (.afg) file to get the following results. The numbers below show Contig_ID, Reads( must be kmers from different reads),length of the contig.
93764 0 87 95957 0 88 105528 0 89 509 1 85 2015 1 85 6601 1 85 . . . 897 2 85 2728 2 87 4386 2 85 . .
How is it possible that the length of contig becomes 87 without any overlap with other kmers? Please clarify and help me if i am misunderstanding.
did you actually turn on read tracking?
what was your velvetg command?
Hi Jeremy. I used the following commands to assemble: ./velveth dir 39,65,4 -fastq -short s_read_1.fq s_read_2.fq ./velvetg dir_43 -read_trkg yes -amos_file yes -unused_reads yes
we would have to see the afg file. AFG is a pretty weird format in that scaffolds and contigs are treated similarly. Make sure you review the spec carefully. http://biostar.stackexchange.com/questions/10745/questions-about-the-asm-file-produced-by-the-velvet-program
•
written
by
Jeremy Leipzig
My understanding is that unused reads are those reads which are not input to the assembly algorithm, normally because they have failed some quality control test. Singletons, on the other hand, are "good" reads which are put into the assembly but do not overlap any other reads.
As for "contigs with zero reads", I think you must have misunderstood something that you have read. Contigs are, by definition, made up of overlapping reads. There's no such thing as a contig with zero reads.
hello neilfws, thank you for your answer. The reason why i raised the issue of zero read contig is, when i mapped the fastq reads to assembly transcripts using bowtie2 and tried to find out the number of reads falling in each contig, i found a significant amount of contigs with no reads mapped. Should'nt i consider them zero read contigs? would be thankful to your reply.
sorry for my typo. I meant those Transcripts obtained after assembling the fastq reads.
Sorry for the typo. Yes. You seem to be correct neilfws. Literally, zero read contigs does not make sense.
I assembled 25M, 90bp fastq reads using trinity assembler and obtained 0.13 Million contigs/transcripts.
I have mapped fastq reads to those transcripts i got by assembling the reads using trinity.
Now when i found the number of reads mapped to each contig/transcript using bowtie2 and samtools, i am confused to see few contigs also had 1 read as well as 0 read.
I also observed that these 1 read, 0 read contig had length of 100bp-500bp length. As you said, If they were not overlapped into contigs why would the length be >90bp. But, my question is why did bowtie2 showed 0 reads for few of the contig still?
•
written
by
Prakki Rama
Log in to answer this question.
-
Total trinity genes, Total trinity transcripts ad unigene
written by sansan •
-
BUSCO duplicated
written by san96 •
-
RNA de novo assembly - blasts - KEGG - GO
written by dimitrischat
-
How to split paired end read as length distribution like (0-50,50-150, >150)?
written by siqizhao2018 •
-
Trinity transcriptome quality assessment
written by cwbenson1993 •
-
htseq read count zero ?
written by cvu
-
Comparing Trinity Assembly Results
written by tanz.renner •
-
assessing de novo Trinity assembly quality
written by trish0401 •
-
14 contigs recieved from SMRTportal using HGAP.3
written by komwit.s •
-
Cleaning up trinity assemblies
written by dpearton •
