Hi all,
Sorry to resurrect this topic, but my workflow depends on Bowtie (not Bowtie2) so I'd like to ask for an advice on the following issue:
In a couple of published paired-end datasets, I encountered a problem that long reads are being mapped OK, but reads shorter than 100 bp in one dataset (or shorter than 150 bp in another dataset) don't get mapped. I am using Bowtie with the following parameters:
-t -v 2 -m 1 --solexa-quals hg19 -1 [reads file 1] -2 [reads file 2]
Could you please suggest what's going wrong?
Thanks!
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
I'm obviously not Ben, but "--un unmapped.fastq" or "--un sample.unmapped.fastq" or something along those lines would be common. Pick a name that makes sense to you, it doesn't matter what it is.
Leave a comment:
-
Hi Ben,
really happy that i can talk with you here because at first when i was working with bowtie2 i asked myself how much you can be clever that created bowtie and how much i am not who cant run bowtie properly...
anyway i have a question about --un option:
if i want to separate mapped and unmapped reads when aligning, which code i should type???
bowtie2 -x [name of the bowtie2-build indicized file containing the rRNA sequence] --un [name of the fastq file which will contain the UNMAPPED reads] -U [name of the fastq file containing the reads] -S [name of the .sam file that will contain the MAPPED and UNMAPPED reads]
I could not understand about --un option because i don't know which i should type instead of [name of the fastq file which will contain the UNMAPPED reads]
Leave a comment:
-
Guest repliedHi everybody!
I starting using bowtie today, i wanted to align csfasta + qual file width the bowtie.
I build the reference fasta file width the bowtie-build, after that i try to align the csfasta+qual file to the reference file(s), but i have error massege.
The bowtie-build command:
bowtie-build -C reference_genom.fa ref/reference_genom
bowtie -C ref/reference_genom -f read.csfasta -Q quality.qual -S align.sam
bowtie -C ref/reference_genom -f read.csfasta -Q quality.qual -S align.sam
/usr/include/seqan/sequence/string_base.h:237 Assertion failed : static_cast<TStringPos>(pos) < static_cast<TStringPos>(length(me)) was: 48 >= 48 (Trying to access an element behind the last one!)
Aborted
My question is that what is the error mean? I try to search this error message but don't found anything.
I installed the bowtie width the following way:
sudo apt-get install bowtie
Thank you!Last edited by Guest; 04-03-2014, 08:34 AM.
Leave a comment:
-
Originally posted by sparks View PostAgree you can merge Bowtie and calculate a vague alignment quality. Perhaps you can give Angie the formulae for it.
Oh, the config.mode just denotes --end-to-end or --local. You'd need to change that to be a function input rather than relying on a global struct I think the remainder should work, though!
Code:/****************************************************************************** * * Calculate a MAPQ, given AS, XS, and the minimum score (ala bowtie2) * *******************************************************************************/ int calc_MAPQ_BT2(int AS, int XS, int scMin) { int diff, bestOver, bestdiff; diff = abs(scMin); //Range of possible alignment scores bestOver = AS-scMin; //Shift alignment score range, so worst score is 0 //The method depends on config.mode bestdiff = (int) abs(abs((float) AS)-abs((float) XS)); //Absolute distance between alignment scores if(config.mode == 0) { //--end-to-end (default) if(XS < scMin) { if(bestOver >= diff * (double) 0.8f) return 42; else if(bestOver >= diff * (double) 0.7f) return 40; else if(bestOver >= diff * (double) 0.6f) return 24; else if(bestOver >= diff * (double) 0.5f) return 23; else if(bestOver >= diff * (double) 0.4f) return 8; else if(bestOver >= diff * (double) 0.3f) return 3; else return 0; } else { if(bestdiff >= diff * (double) 0.9f) { if(bestOver == diff) { return 39; } else { return 33; } } else if(bestdiff >= diff * (double) 0.8f) { if(bestOver == diff) { return 38; } else { return 27; } } else if(bestdiff >= diff * (double) 0.7f) { if(bestOver == diff) { return 37; } else { return 26; } } else if(bestdiff >= diff * (double) 0.6f) { if(bestOver == diff) { return 36; } else { return 22; } } else if(bestdiff >= diff * (double) 0.5f) { if(bestOver == diff) { return 35; } else if(bestOver >= diff * (double) 0.84f) { return 25; } else if(bestOver >= diff * (double) 0.68f) { return 16; } else { return 5; } } else if(bestdiff >= diff * (double) 0.4f) { if(bestOver == diff) { return 34; } else if(bestOver >= diff * (double) 0.84f) { return 21; } else if(bestOver >= diff * (double) 0.68f) { return 14; } else { return 4; } } else if(bestdiff >= diff * (double) 0.3f) { if(bestOver == diff) { return 32; } else if(bestOver >= diff * (double) 0.88f) { return 18; } else if(bestOver >= diff * (double) 0.67f) { return 15; } else { return 3; } } else if(bestdiff >= diff * (double) 0.2f) { if(bestOver == diff) { return 31; } else if(bestOver >= diff * (double) 0.88f) { return 17; } else if(bestOver >= diff * (double) 0.67f) { return 11; } else { return 0; } } else if(bestdiff >= diff * (double) 0.1f) { if(bestOver == diff) { return 30; } else if(bestOver >= diff * (double) 0.88f) { return 12; } else if(bestOver >= diff * (double) 0.67f) { return 7; } else { return 0; } } else if(bestdiff > 0) { if(bestOver >= diff * (double) 0.67f) { return 6; } else { return 2; } } else { if(bestOver >= diff * (double) 0.67f) { return 1; } else { return 0; } } } } else { //--local if(XS < scMin) { if(bestOver >= diff * (double) 0.8f) return 44; else if(bestOver >= diff * (double) 0.7f) return 42; else if(bestOver >= diff * (double) 0.6f) return 41; else if(bestOver >= diff * (double) 0.5f) return 36; else if(bestOver >= diff * (double) 0.4f) return 28; else if(bestOver >= diff * (double) 0.3f) return 24; else return 22; } else { if(bestdiff >= diff * (double) 0.9f) return 40; else if(bestdiff >= diff * (double) 0.8f) return 39; else if(bestdiff >= diff * (double) 0.7f) return 38; else if(bestdiff >= diff * (double) 0.6f) return 37; else if(bestdiff >= diff * (double) 0.5f) { if (bestOver == diff) return 35; else if(bestOver >= diff * (double) 0.5f) return 25; else return 20; } else if(bestdiff >= diff * (double) 0.4f) { if (bestOver == diff) return 34; else if(bestOver >= diff * (double) 0.5f) return 21; else return 19; } else if(bestdiff >= diff * (double) 0.3f) { if (bestOver == diff) return 33; else if(bestOver >= diff * (double) 0.5f) return 18; else return 16; } else if(bestdiff >= diff * (double) 0.2f) { if (bestOver == diff) return 32; else if(bestOver >= diff * (double) 0.5f) return 17; else return 12; } else if(bestdiff >= diff * (double) 0.1f) { if (bestOver == diff) return 31; else if(bestOver >= diff * (double) 0.5f) return 14; else return 9; } else if(bestdiff > 0) { if(bestOver >= diff * (double) 0.5f) return 11; else return 2; } else { if(bestOver >= diff * (double) 0.5f) return 1; else return 0; } } } }
Leave a comment:
-
Agree you can merge Bowtie and calculate a vague alignment quality. Perhaps you can give Angie the formulae for it.
Originally posted by dpryan View PostWhile using BWA or Novoalign are certainly the better solutions, one can relatively simply recalculate MAPQs from multiple alignment files to different references with bowtie. The bowtie MAPQ score is dependent primarily on the AS:i: and XS:i: score of each read, so you can just rerun the algorithm on that (bowtie MAPQs are more of a vague approximation than you may think). This is the approach I took in bison, where there are multiple parallel alignments of each read to different bisulfite converted genomes.
Leave a comment:
-
Originally posted by sparks View PostHi Angie,
I think if you use split reference you'll have issues calculating the alignment quality of the best alignment during the merge, it's a bit more complicated than just selecting the best alignment.
You will likely get more accurate alignment qualities if you don't split the reference and instead use an aligner like BWA or Novoalign that can handle genomes >4Gbp.
KR, Colin
Leave a comment:
-
Thanks for the reply Colin and rshina. As suggested I have indexed the reference without splitting it with BWA so I will proceed with this approach
Cheers
Angela
Leave a comment:
-
Reply to thread 'Bowtie, an ultrafast, memory-efficient, open source short read align
Hi Angie,
I think if you use split reference you'll have issues calculating the alignment quality of the best alignment during the merge, it's a bit more complicated than just selecting the best alignment.
You will likely get more accurate alignment qualities if you don't split the reference and instead use an aligner like BWA or Novoalign that can handle genomes >4Gbp.
KR, Colin
Originally posted by angie_red View PostHi Ben,
Sorry to resurrect an old post. I am getting the error Error: Reference sequence has more than 2^32-1 characters!. I know this means I need to split my reference in order to use bowtie2-build but I am wondering about mapping my reads to this reference which has been split. Is it possible to concatenate the split-indexed files and map the reads to this concatenated file or will I have to map the reads to each indexed files separately and write scripts to find which has the best hit.
Thank you
Angela
Leave a comment:
-
Bowtie2
I think you need to map your reads to divided indexes and then write script to bring them together.
Leave a comment:
-
Hi Ben,
Sorry to resurrect an old post. I am getting the error Error: Reference sequence has more than 2^32-1 characters!. I know this means I need to split my reference in order to use bowtie2-build but I am wondering about mapping my reads to this reference which has been split. Is it possible to concatenate the split-indexed files and map the reads to this concatenated file or will I have to map the reads to each indexed files separately and write scripts to find which has the best hit.
Thank you
Angela
Leave a comment:
-
Dear Ben,
Why is the last version of Bowtie using the mm9 rather than mm10?
What is better Bowtie or Bowtie2 for alighment of 50 nt HiSeq Illumina ChIP-Seq redas?
I have read that Bowtie is good for short reads up to 100 nt, but Bowtie2 from 50 nt and higher. Still 50 nt reads are on the border for the programms.
If Bowtie2 is used, how to get rid of ununique reads?
Many thanks in advance
Leave a comment:
-
Originally posted by sahiilseth View PostFrom Bowtie website:
They also say:
'If your computer has more than 3-4 GB of memory and you would like to exploit that fact to make index building faster, use a 64-bit version of the bowtie2-build binary. The 32-bit version of the binary is restricted to using less than 4 GB of memory. If a 64-bit pre-built binary does not yet exist for your platform on the sourceforge download site, you will need to build one from source.'
I thought 64 bit binary, should be able to handle more characters as well; not true?
Leave a comment:
Latest Articles
Collapse
-
by seqadmin
Innovations in next-generation sequencing technologies and techniques are driving more precise and comprehensive exploration of complex biological systems. Current advancements include improved accessibility for long-read sequencing and significant progress in single-cell and 3D genomics. This article explores some of the most impactful developments in the field over the past year.
Long-Read Sequencing
Long-read sequencing has...-
Channel: Articles
12-02-2024, 01:49 PM -
-
by seqadmin
The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben MartÃnez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...-
Channel: Articles
11-06-2024, 07:24 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 12-02-2024, 09:29 AM
|
0 responses
140 views
0 likes
|
Last Post
by seqadmin
12-02-2024, 09:29 AM
|
||
Started by seqadmin, 12-02-2024, 09:06 AM
|
0 responses
50 views
0 likes
|
Last Post
by seqadmin
12-02-2024, 09:06 AM
|
||
Started by seqadmin, 12-02-2024, 08:03 AM
|
0 responses
38 views
0 likes
|
Last Post
by seqadmin
12-02-2024, 08:03 AM
|
||
Started by seqadmin, 11-22-2024, 07:36 AM
|
0 responses
70 views
0 likes
|
Last Post
by seqadmin
11-22-2024, 07:36 AM
|
Leave a comment: