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
During the COVID-19 pandemic, scientists observed that while some individuals experienced severe illness when infected with SARS-CoV-2, others were barely affected. These disparities left researchers and clinicians wondering what causes the wide variations in response to viral infections and what role genetics plays.
Jean-Laurent Casanova, M.D., Ph.D., Professor at Rockefeller University, is a leading expert in this crossover between genetics and infectious...-
Channel: Articles
09-09-2024, 10:59 AM -
-
by seqadmin
The first FDA-approved CRISPR-based therapy marked the transition of therapeutic gene editing from a dream to reality1. CRISPR technologies have streamlined gene editing, and CRISPR screens have become an important approach for identifying genes involved in disease processes2. This technique introduces targeted mutations across numerous genes, enabling large-scale identification of gene functions, interactions, and pathways3. Identifying the full range...-
Channel: Articles
08-27-2024, 04:44 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Today, 06:25 AM
|
0 responses
13 views
0 likes
|
Last Post
by seqadmin
Today, 06:25 AM
|
||
Started by seqadmin, Yesterday, 01:02 PM
|
0 responses
12 views
0 likes
|
Last Post
by seqadmin
Yesterday, 01:02 PM
|
||
Started by seqadmin, 09-18-2024, 06:39 AM
|
0 responses
14 views
0 likes
|
Last Post
by seqadmin
09-18-2024, 06:39 AM
|
||
Started by seqadmin, 09-11-2024, 02:44 PM
|
0 responses
14 views
0 likes
|
Last Post
by seqadmin
09-11-2024, 02:44 PM
|
Comment