SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
using bwa to map illumina paired end reads mikeworth Bioinformatics 6 08-13-2013 06:02 PM
SOLiD reads map using BFAST and Samtools jsun529 SOLiD 17 02-24-2011 12:55 PM
PE SOLiD reads alignment by bwa m_elena_bioinfo Bioinformatics 8 02-21-2011 08:09 AM
How to map SOLiD paired end reads by Bfast beliefbio Bioinformatics 1 12-29-2010 12:55 AM
map solid reads to ncRNA database crh SOLiD 2 07-19-2010 11:22 AM

Reply
 
Thread Tools
Old 12-03-2010, 10:53 PM   #1
Chandana
Junior Member
 
Location: Singapore

Join Date: Dec 2010
Posts: 8
Default How to map NCBI SOLiD reads with BWA

Hi,
I have downloaded some SOLiD reads from NCBI and needs to map them using BWA. Should I simply convert 0->A, 1->C etc. or do I need to chop the first two bases in addition?
Chandana is offline   Reply With Quote
Old 12-04-2010, 12:34 AM   #2
Chipper
Senior Member
 
Location: Sweden

Join Date: Mar 2008
Posts: 324
Default

If you really need to use BWA you have to use the solid2fastq.pl script that comes with it. You also need to change the parameters to increase the number of mismatches allowed.
Chipper is offline   Reply With Quote
Old 12-04-2010, 01:27 AM   #3
Chandana
Junior Member
 
Location: Singapore

Join Date: Dec 2010
Posts: 8
Default

Thanks Chipper for the quick reply, but the script does not seem to work on my case. the reads are like


@ERR012630.16 VAB_solid0029_20100115_QR1000044_50X2_QR10000441_14_61/1
T0010..330..00013.00.0.20...00.0..02..0000..00.....
+
!/95%!!<).!!:%,/)!17!9!54!!!45!-!!!!68!!9;!!!!!

but the scripts expects the quality values and the reads separated. Can you help me in this case?
Chandana is offline   Reply With Quote
Old 12-04-2010, 05:57 AM   #4
Chipper
Senior Member
 
Location: Sweden

Join Date: Mar 2008
Posts: 324
Default

Ok, since it is FastQ already you could just change 0123 to ACGT and remove the first two bases as you suggested. You could also try bowtie, it's faster and I think it will accept the file as is.
Chipper is offline   Reply With Quote
Reply

Tags
bwa ncbi solid

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 03:22 AM.


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