SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Bowtie. Different outputs from equivalent(?) inputs. AEB Bioinformatics 1 08-18-2011 12:34 AM
question about bwa sam file peachgil Bioinformatics 1 04-06-2011 08:53 AM
the meaning of CIGAR column in SAM file outputs by BWA holywoool Bioinformatics 2 01-04-2011 04:34 AM
gsAssembly (Newbler) de novo behaviour, inputs and outputs nicolallias 454 Pyrosequencing 6 10-29-2010 12:16 AM
Aligner that outputs H2 flag in SAM file phalaenopsis Bioinformatics 0 08-04-2010 05:14 AM

Reply
 
Thread Tools
Old 01-06-2010, 08:07 AM   #1
ikrier
Member
 
Location: Lausanne

Join Date: Dec 2009
Posts: 19
Default bwa question on file outputs/inputs

Reposted from an old thread in the SOLiD forum :

Hi, I did the indexing with bwtsw and no -p and I got the following files :
Mouse_genome.fa.amb
Mouse_genome.fa.ann
Mouse_genome.fa.bwt
Mouse_genome.fa.pac
Mouse_genome.fa.rbwt
Mouse_genome.fa.rpac
Mouse_genome.fa.rsa
Mouse_genome.fa.sa

I managed to get the .sai file from the aln command, but now I'm stuck because the samse command gives me the error:
fail to open file '../Mouse_genome.fa.nt.ann'. Abort!

But I never get the .nt.ann file with indexing. I'm confused.
ikrier is offline   Reply With Quote
Old 01-19-2010, 01:01 AM   #2
dnusol
Senior Member
 
Location: Spain

Join Date: Jul 2009
Posts: 133
Default

Hi Ikrier,
did you manage to solve this problem? I am interested in how you got round it, because I am having problems in the previous step

I am trying to use the bwa aln command with default paramenters, I start getting lots of weird characters on the screen and finally I get no out.sai output. I get no warning that I can see.

Can anyone drop some hints?

Cheers
dnusol is offline   Reply With Quote
Old 01-19-2010, 01:26 AM   #3
ikrier
Member
 
Location: Lausanne

Join Date: Dec 2009
Posts: 19
Default

i'm still confused because the problem solved itself. I don't even know why, but after trying several times it didn't ask for the same files anymore.
ikrier is offline   Reply With Quote
Old 01-19-2010, 02:16 AM   #4
dnusol
Senior Member
 
Location: Spain

Join Date: Jul 2009
Posts: 133
Default

Did you actually get the out.sai file in the bwa aln step?

my command is

$ ./bwa aln myref.fasta /path_to/s_1_sequence.fastq out.sai

I cannot get that out.sai file
dnusol is offline   Reply With Quote
Old 01-19-2010, 02:18 AM   #5
ikrier
Member
 
Location: Lausanne

Join Date: Dec 2009
Posts: 19
Default

Did you try
./bwa aln myref.fasta /path_to/s_1_sequence.fastq > out.sai
ikrier is offline   Reply With Quote
Old 01-19-2010, 02:36 AM   #6
dnusol
Senior Member
 
Location: Spain

Join Date: Jul 2009
Posts: 133
Default

Shoot! I missed the >

Thanks, I'll try again
dnusol is offline   Reply With Quote
Old 01-10-2011, 03:13 AM   #7
huma Asif
Member
 
Location: Japan

Join Date: Oct 2010
Posts: 53
Default error in alignment

Dear All
while running BWA aln i got this error can anybody help in fixing it
dr@dr-Revision-A:~/Desktop/bwa-0.5.7$ ./bwa aln myref.fasta /path_to/s_1_sequence.fastq > out.sai
[bwa_aln] 17bp reads: max_diff = 2
[bwa_aln] 38bp reads: max_diff = 3
[bwa_aln] 64bp reads: max_diff = 4
[bwa_aln] 93bp reads: max_diff = 5
[bwa_aln] 124bp reads: max_diff = 6
[bwa_aln] 157bp reads: max_diff = 7
[bwa_aln] 190bp reads: max_diff = 8
[bwa_aln] 225bp reads: max_diff = 9
[bwa_seq_open] fail to open file '/path_to/s_1_sequence.fastq'. Abort!
Aborted

Regards
Huma
huma Asif is offline   Reply With Quote
Old 01-10-2011, 06:29 AM   #8
dnusol
Senior Member
 
Location: Spain

Join Date: Jul 2009
Posts: 133
Default

Hi Huma,

you actually have to write the full path to your fastq file.

So if for example you have the fastq file in

/home/huma/data/s_1_sequence.fastq

then your command should say so

./bwa aln myref.fasta /home/huma/data/s_1_sequence.fastq > out.sai

If your fastq file is in your current drectory then you don't need to specify it.

also, if you do not specify the full path where you want to redirect the output, this will be saved in the current directory

A useful trick to use is to complete the name of the files requested using the tab key. If you try to do the following from the command line

./bwa aln myref.fasta /path_to/s_1_sequence.fas [press tab to complete]

then the shell will not complete the name automatically

HTH.
dnusol is offline   Reply With Quote
Old 01-10-2011, 08:34 PM   #9
huma Asif
Member
 
Location: Japan

Join Date: Oct 2010
Posts: 53
Default use of index file in BWA

Thank you so much for your kind reply .Actually i am confused about the index files that i generated in my command
./bwa aln ref/putidaKT2440.fasta reads/pseudo.fq out.sai
I am not using index file anywhere so I am confused what do i need to do to make use of it
Regards

Last edited by huma Asif; 01-11-2011 at 12:08 AM.
huma Asif is offline   Reply With Quote
Old 01-11-2011, 12:05 AM   #10
huma Asif
Member
 
Location: Japan

Join Date: Oct 2010
Posts: 53
Default BWA output

finally i have run BWA aligner successfully
i have got an output file in SAM format
now i opened that sam file in samtool pseudo.bam,
pseudo.sorted.bam
pseudo.sorted.bam.bai
raw.pileup

please guide me which file shal i use to view alignment

Regards
huma Asif is offline   Reply With Quote
Old 09-11-2012, 02:39 AM   #11
Orr Shomroni
Member
 
Location: Netherlands

Join Date: Oct 2011
Posts: 26
Smile Problem solved

ikier, I had the same problem as you with the core dumped and the fa.nt.ann failure to open. I found the solution. When you run the bwa index, it should be run with the following code:
Code:
bwa index -a bwtsw database.fasta
For some reason, bwa index by default uses algorithm IS, which doesn't work as well as the BWTSW algorithm. See the explanation on the website
HTML Code:
http://bio-bwa.sourceforge.net/bwa.shtml
__________________
"Though it may seem that all's been said and done, originality still lives on" - some unoriginal guy who had nothing better to write as his signature
Orr Shomroni is offline   Reply With Quote
Old 10-12-2012, 04:40 AM   #12
dtc
Junior Member
 
Location: ECNU

Join Date: Oct 2012
Posts: 3
Default

Hi Ikrier,
do you know how it fixed itself now,I have the same problem,but I tried many times,it still can't work.
Cheers
dtc is offline   Reply With Quote
Old 10-12-2012, 08:22 AM   #13
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

Quote:
Originally Posted by Orr Shomroni View Post
ikier, I had the same problem as you with the core dumped and the fa.nt.ann failure to open. I found the solution. When you run the bwa index, it should be run with the following code:
Code:
bwa index -a bwtsw database.fasta
For some reason, bwa index by default uses algorithm IS, which doesn't work as well as the BWTSW algorithm. See the explanation on the website
HTML Code:
http://bio-bwa.sourceforge.net/bwa.shtml
The default index on bwa index works fine for small genomes. I use all the time with no problems on bacterial genomes.

But for larger genomes, like mammals, yes, you have to use bwtsw.
swbarnes2 is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 10:25 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO