SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   How to map SOLiD paired end reads by Bfast (http://seqanswers.com/forums/showthread.php?t=8616)

beliefbio 12-28-2010 10:10 PM

How to map SOLiD paired end reads by Bfast
 
1. I found that bfast localalign can't support -1 and -2 options. So I'm confused how it works in
http://seqanswers.com/forums/showthread.php?t=7100
My bfast version is 0.6.4e. Is -1 and -2 work in other versions?

2. Then I tried to merge 2 paired end fastq files into 1 file and run bfast match. But in the bfast postprocess step how to set -r ReadGroup.txt? Could you please show an simple example?

3. The MAPQ of sam files generated by bfast is also Phred scaled, right? If so I can set an cutoff to filter them.
Thanks a lot for you help and time!

nilshomer 12-29-2010 01:55 AM

Quote:

Originally Posted by beliefbio (Post 32077)
1. I found that bfast localalign can't support -1 and -2 options. So I'm confused how it works in
http://seqanswers.com/forums/showthread.php?t=7100
My bfast version is 0.6.4e. Is -1 and -2 work in other versions?

Check out the bfast+bwa version to get that functionality (I assumed you used the bwa functionality instead of bfast match).

For clarity, if you wish to use BWA's "bwa aln" functionality for paired end reads, then you need to align each end independently, merging them together in "bfast localalign" with the "-1/-2" options. Otherwise, the paired end reads should be merged (automatically in solid2fastq) into one fastq file to be processed by "bfast match".

Quote:

Originally Posted by beliefbio (Post 32077)
2. Then I tried to merge 2 paired end fastq files into 1 file and run bfast match. But in the bfast postprocess step how to set -r ReadGroup.txt? Could you please show an simple example?

See the SAM spec on how to create an RG line.

Quote:

Originally Posted by beliefbio (Post 32077)
3. The MAPQ of sam files generated by bfast is also Phred scaled, right? If so I can set an cutoff to filter them.
Thanks a lot for you help and time!

You can set a cutoff to filter them yes. Better (higher) mapq should correlate with higher stringency.


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

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.