Question: Bowtie2 aligning output
0
gravatar for huiwenteh
5 months ago by
huiwenteh0
huiwenteh0 wrote:

Hi, I am currently using bowtie2 to align my sequence. I am using a assembled sequence to align with my reference sequence.Unfortunately, the output after aligning is stated as killed and error. I am currently stuck in this situation. Below is my command line which is inserted and the output.

$ bowtie2 -f -x Mycobacterium -U R1_final.contigs.fa -S R1.sam Killed (ERR): bowtie2-align exited with value 137

I am not sure what is the problem and how to solve it. I hope I can get some help here :)

p/s: I am still abit new to bioinformatics and so sorry for the inconvenience caused :)

Cheers, HuiWen

alignment • 156 views
ADD COMMENTlink modified 5 months ago by Asaf8.4k • written 5 months ago by huiwenteh0

I think you are looking for multiple sequence alignment. Bowtie2 supports reads up to 1000 bp. What's the length of of your assembled sequence?

ADD REPLYlink written 5 months ago by Arup Ghosh2.7k

Hi, my sequence length is around 4MB. Any suggestion for multiple sequence alignment?

ADD REPLYlink written 5 months ago by huiwenteh0

You don't need Multiple Sequence Alignment, just pairwise alignment.

ADD REPLYlink written 5 months ago by Asaf8.4k

If I would be using pairwise alignment would I be able to preform variant calling?

ADD REPLYlink written 5 months ago by huiwenteh0

You better use the raw reads, not the assembly in that case

ADD REPLYlink written 5 months ago by Asaf8.4k

means that I run my raw reads by using pairwise alignment and then straight variant calling?

ADD REPLYlink written 5 months ago by huiwenteh0
1

I strongly recommend properly training on how to do that, it's a complex topic. You might want to collaborate on this aspect.

ADD REPLYlink written 5 months ago by Asaf8.4k
1
gravatar for Asaf
5 months ago by
Asaf8.4k
Israel
Asaf8.4k wrote:

As Arup mentioned, Bowtie2 is meant for short reads mapping, not assembly contigs. You can try other software like minimap2 or mummer.

ADD COMMENTlink written 5 months ago by Asaf8.4k

Hi, I tried mummer before to align and SNP calling but I got stuck when it comes to visualize the SNPs. Due to mummer produce delta file, I wasnt able to convert delta file into sam/bam file to visualize it in IGV/bamview. Thank you for the suggestion on minimap2. I give this software a try :)

ADD REPLYlink written 5 months ago by huiwenteh0

What's your primary goal, variant calling or (pairwise/multiple) sequence alignment?

ADD REPLYlink modified 5 months ago • written 5 months ago by Arup Ghosh2.7k

my primary goal is variant calling. But if I continue my work with pairwise alignment is it possible to do variant calling after the pairwise alignment? Because to my best knowledge for variant calling we need to perform multiple sequence alignment. Correct me if I am wrong.

ADD REPLYlink modified 5 months ago • written 5 months ago by huiwenteh0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1310 users visited in the last hour