Is there any way to exclude the read lengths below certain length say 100bp for denovo assembly using Velvet
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Originally posted by nxtgenkid10 View PostIs there any way to exclude the read lengths below certain length say 100bp for denovo assembly using Velvet
Out of curiosity, why are you wanting to do this? If your data is 454 then Newbler will probably do a better job, in my experience, and if your data is from an Illumina 100bp library then you are not really going to have any variation in read length.
-
Originally posted by SES View PostI don't think there is a way to specify the reads that go into the contigs with velvet (like you can with Newbler). You can specify the minimum contig length in Velvet, but that is not exactly what you wanted. The best thing to do would be to filter the reads before assembly.
Out of curiosity, why are you wanting to do this? If your data is 454 then Newbler will probably do a better job, in my experience, and if your data is from an Illumina 100bp library then you are not really going to have any variation in read length.
Comment
-
Originally posted by nangillala View PostSo why don't you trust your reads below 100 bp? (Maybe there's a reason for this).
Velvet uses k-mers, so the reads will be "splitted" into shorter fragments anyways.
Comment
-
Originally posted by nxtgenkid10 View PostThe thing is that I'm lokking for the reads with Q30+ and feels it may have some adapter contamination; aS i told before I'm new to this and looking for options in the tools to make an clear interpretation from the data in hand
I'm not sure what kind of reads you have: Illumina? 454? SFF format? Fastq?
Assuming you have files in fastq format you could use the fastx toolkit
to get rid of adapters (Fastq_Clipper). You can also use it to discard sequences shorter than a threshold you choose.
You may also like the Fastq Quality Filter. You can choose the minimum quality score to keep and the minimum percent of bases that must have this quality.
Hope that helps!?
Comment
-
Originally posted by nangillala View PostHi,
I'm not sure what kind of reads you have: Illumina? 454? SFF format? Fastq?
Assuming you have files in fastq format you could use the fastx toolkit
to get rid of adapters (Fastq_Clipper). You can also use it to discard sequences shorter than a threshold you choose.
You may also like the Fastq Quality Filter. You can choose the minimum quality score to keep and the minimum percent of bases that must have this quality.
Hope that helps!?
Comment
Latest Articles
Collapse
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...-
Channel: Articles
04-22-2024, 07:01 AM -
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Today, 08:47 AM
|
0 responses
12 views
0 likes
|
Last Post
by seqadmin
Today, 08:47 AM
|
||
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
60 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
59 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
54 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
Comment