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
Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...-
Channel: Articles
10-18-2024, 07:11 AM -
-
by seqadmin
Non-coding RNAs (ncRNAs) do not code for proteins but play important roles in numerous cellular processes including gene silencing, developmental pathways, and more. There are numerous types including microRNA (miRNA), long ncRNA (lncRNA), circular RNA (circRNA), and more. In this article, we discuss innovative ncRNA research and explore recent technological advancements that improve the study of ncRNAs.
Nobel Prize for MicroRNA Discovery
This week,...-
Channel: Articles
10-07-2024, 08:07 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 06:09 AM
|
0 responses
10 views
0 likes
|
Last Post
by seqadmin
Yesterday, 06:09 AM
|
||
New Model Aims to Explain Polygenic Diseases by Connecting Genomic Mutations and Regulatory Networks
by seqadmin
Started by seqadmin, 10-30-2024, 05:31 AM
|
0 responses
13 views
0 likes
|
Last Post
by seqadmin
10-30-2024, 05:31 AM
|
||
Started by seqadmin, 10-24-2024, 06:58 AM
|
0 responses
22 views
0 likes
|
Last Post
by seqadmin
10-24-2024, 06:58 AM
|
||
New AI Model Designs Synthetic DNA Switches for Targeted Gene Expression in Specific Cell Types
by seqadmin
Started by seqadmin, 10-23-2024, 08:43 AM
|
0 responses
52 views
0 likes
|
Last Post
by seqadmin
10-23-2024, 08:43 AM
|
Comment