Is there a command to output the kmers of each sequence in a multifasta file?
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Trouble parsing header
Dear BBMap team:
I tried to use filterbytile.sh to remove the reads with low quality, but I encountered an error message saying that there was a trouble parsing the header. I've read the description of the script and Brian Bushnell said that was possible when the reads were renamed (such as in SRA) and to contact him if such error happened.
I downloaded the sequencing data (SRA) from ncbi and used fastq-dump to get the fastq files. I wonder if there is a solution to this?
Thank you very much!
Rose
Comment
-
BBsketch alltoall is incomplete
Can I ask a question about bbsketch?
I want to compare the ANI between many genomes (1000+) to each other.
I did
Code:bbsketch.sh perfile genome_folder/*.fasta out=sketch.gz k=31,24 threads=16 comparesketch.sh alltoall sketch.gz k=31,24 prealloc=0.75 format=3 threads=16 out=table.tsv
Code:Set threads to 16 Loading sketches. Loaded 1157 sketches in 59.541 seconds. Total Time: 59.784 seconds.
Code:Set threads to 16 Loading sketches. Executing kmer.KmerTableSet [ways=31, tabletype=10, prealloc=0.75] Initial size set to 45218398 Initial: Ways=31, initialSize=45218398, prefilter=f, prealloc=0.75 Memory: max=91268m, total=91268m, free=90848m, used=420m 3.713 seconds. Indexed 2880884 unique and 10513099 total hashcodes. Loaded 1157 sketches in 8.457 seconds. Ran 1225005 comparisons in 9.344 seconds. Total Time: 17.801 seconds.
- Genomes are highly similar.
#Query Ref ANI QSize RefSize QBases RBases QTaxID RTaxID KID WKID SSU
genome1.fasta genome2.fasta 94.223 1984118 1796930 1987598 1797650 -1 -1 24.952 27.523 .
- It is not simply due to the naming: I neither find "genome1 vs genome2" nor "genome 2 vs genome1"
Any idea?
Comment
-
I'm trying to use BBmap to find all perfect hits or hits with an indel length 1.
Code:bbmapskinner.sh in=kmer.fasta out=result.sam ambiguous=all strictmaxindel=1
Is there something that I am doing wrong?
Comment
Latest Articles
Collapse
-
by seqadmin
While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...-
Channel: Articles
Yesterday, 07:15 AM -
-
by seqadmin
Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.
Somatic Genomics
“We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...-
Channel: Articles
05-24-2024, 01:16 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Today, 06:58 AM
|
0 responses
13 views
0 likes
|
Last Post
by seqadmin
Today, 06:58 AM
|
||
Started by seqadmin, Yesterday, 08:18 AM
|
0 responses
19 views
0 likes
|
Last Post
by seqadmin
Yesterday, 08:18 AM
|
||
Started by seqadmin, Yesterday, 08:04 AM
|
0 responses
18 views
0 likes
|
Last Post
by seqadmin
Yesterday, 08:04 AM
|
||
Started by seqadmin, 06-03-2024, 06:55 AM
|
0 responses
13 views
0 likes
|
Last Post
by seqadmin
06-03-2024, 06:55 AM
|
Comment