Actually, while it's not documented, the flag "minlength" also works with BBMap. Reads shorter than that will be discarded completely (they won't be output as unmapped).
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Originally posted by Brian Bushnell View PostActually, while it's not documented, the flag "minlength" also works with BBMap. Reads shorter than that will be discarded completely (they won't be output as unmapped).
Is there a way to make BBMap require the smallest read to be some fraction of the longest read length? I know that's a niche use but BBMap always suprises me with it's built in functions.
Comment
-
Originally posted by Shini Sunagawa View PostDear Brian,
I have been looking for a tool that would quickly dereplicate (100% containments) nucleotide sequences and track for each unique sequence the identifiers of the removed duplicates.
Something like:
dedupe.sh in=in.fa out=out.fa outd=outd.fa mid=100 mop=100
where:
in.fa:
seq1
seq2 (contained in seq1)
seq3 (contained in seq1)
seq4
out.fa:
seq1
seq4
outd.fa:
seq2
seq3
I am interested in:
seq1<tab>seq2,seq3
seq4
dedupe.sh does a fantastic job in returning out and outd, but I cannot find any option that would return the information I am interested in. Is this something that I am missing? Otherwise, I believe this could be a great feature, since compared to other tools that return this information, dedupe is so much faster.
Best,
Shini
Comment
-
Hmmm... no, not yet, though I did add into Clumpify's dereplication step the ability to count the number of duplicate reads and add "count=3", for example, to the name of a read representing 3 total reads (itself and 2 duplicates). It would not be difficult to modify that to report read identifiers. I'll add it to my todo list.
Comment
-
Originally posted by Brian Bushnell View PostHmmm... no, not yet, though I did add into Clumpify's dereplication step the ability to count the number of duplicate reads and add "count=3", for example, to the name of a read representing 3 total reads (itself and 2 duplicates). It would not be difficult to modify that to report read identifiers. I'll add it to my todo list.
Comment
-
Clumpify can consider sequences as duplicates if they have at most X substitutions, but it's not as flexible as Dedupe. For example, Clumpify requires duplicates to overlap 100% with neither overhanging, while Dedupe allows containments (this only matters when using variable-length sequences) and also allows indels. What I actually added to my todo list was to update both of them with that capability, since it seems useful.
Comment
-
Originally posted by Brian Bushnell View PostClumpify can consider sequences as duplicates if they have at most X substitutions, but it's not as flexible as Dedupe. For example, Clumpify requires duplicates to overlap 100% with neither overhanging, while Dedupe allows containments (this only matters when using variable-length sequences) and also allows indels. What I actually added to my todo list was to update both of them with that capability, since it seems useful.
Comment
-
Hello Brian,
I would like to ask you a suggestion for bbmap.
I am trying to reassemble a bin from a metagenomic data, hoping that I will get better assembly if I use just the mapped reads.
I tried bbmap.sh on normal parameters and used the outm to take all the aligned reads, then I normalized the reads with bbnorm.sh, and reassembled with SPAdes. I want to note that the initial metagenomic assembly was done on normalized reads, and SPAdes does error correction, but I did not use these libraries, I used the adapter and quality trimmed libraries but not normalized and error corrected (this is why I do normalization after mapping).
I got better assembly, (some longer scaffolds, and slightly larger N50) but checking briefly the SSU I noticed that some "contaminants" were present. Also, the amount of SSU sequences was much higher than expected. (I expect 4, 3 complete ones and one near complete).
In the metagemic data (assembled using all the data) I have 10 SSUs, but here I got a lot of them (15+) and most of them are really partial.
What I am thinking is that bbmap, includes in the output some(not all) reads from other bacterial SSU which map to a certain degree to my reference (since it can have very conserved regions) then SPAdes is somehow confused by these, and fragments my SSU sequences in multiple places due to these reads. Sorry if this sounds strange, but I am just speculating, I am not sure if this is the case, and unfortunately I am not an expert bioinformatician.
I was thinking to add to the command line the parameters minidentity=0.98 idfilter=0.98 hoping that I will somehow avoid the mapping of "non-specific" reads.
I avoid using the perfectmode=t due to the fact that SPAdes does error correction, and this would somehow cause some small mismatches, and thus I would lose some reads.
Would you have any suggestion for better parameters of bbmap?
Thank you!
PS: I am using PE 2X250bp and PE2x300bp libraries for the mapping.
Comment
-
Non-Deterministic BBmap results... how ensure deterministic?
Hi Brian,
It is important that my pipeline's results can be perfectly reproduced. I notice that the non-deterministic behavior is coming from the human read removal from my datasets... here are the parameters I have specified in this call:
bbmap.sh\
-Xmx23g\
minid=0.9\
idfilter=0.9\
maxindel=3\
bwr=0.16\
bw=12\
minhits=2\
printunmappedcount=t\
The numbers are close (ie. 672510 vs 672492)
I run the PE reads through BBduk to remove low-quality pairs first. It looks like that output is sorted the same and deterministic upon re-runs.
Thanks for your thoughts, KateLast edited by sk8bro; 05-03-2017, 08:32 AM.
Comment
-
Originally posted by Brian Bushnell View PostActually, while it's not documented, the flag "minlength" also works with BBMap. Reads shorter than that will be discarded completely (they won't be output as unmapped).
I get this error: Read of length 36 outside of range 50--1. Paired input is incompatible with 'breaklength'
Comment
-
Originally posted by darthsequencer View PostDoes this work with paired end reads?
I get this error: Read of length 36 outside of range 50--1. Paired input is incompatible with 'breaklength'
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 | ||
---|---|---|---|---|
New Model Aims to Explain Polygenic Diseases by Connecting Genomic Mutations and Regulatory Networks
by seqadmin
Started by seqadmin, Yesterday, 05:31 AM
|
0 responses
10 views
0 likes
|
Last Post
by seqadmin
Yesterday, 05:31 AM
|
||
Started by seqadmin, 10-24-2024, 06:58 AM
|
0 responses
20 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
48 views
0 likes
|
Last Post
by seqadmin
10-23-2024, 08:43 AM
|
||
Started by seqadmin, 10-17-2024, 07:29 AM
|
0 responses
58 views
0 likes
|
Last Post
by seqadmin
10-17-2024, 07:29 AM
|
Comment