SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > SOLiD

Similar Threads
Thread Thread Starter Forum Replies Last Post
GTF to BED converter? mmuratet Bioinformatics 7 02-20-2012 10:19 AM
GFF 2 genbank converter deMan Bioinformatics 3 02-16-2012 01:33 PM
XSQ format smarkel SOLiD 4 02-21-2011 05:00 AM
pileup2vcf converter? a11msp Bioinformatics 2 11-21-2010 04:29 AM
soap2bowtie converter? bioenvisage Bioinformatics 1 12-03-2009 07:10 AM

Reply
 
Thread Tools
Old 12-29-2011, 04:58 AM   #1
paolo.kunder
Member
 
Location: Milano, Italy

Join Date: Aug 2011
Posts: 73
Default XSQ converter

I have a PE run, 6 lane, 12 barcoded samples.
ICS SOLiD software generated 6 XSQ file ( one per lane).

Now I will use XSQ_Tools to split each XSQ file in 12 indexed XSQ.

I will end up with 72 XSQ indexed files (6 lanes x 12 samples).

Is there a way to combine the 6 XSQ files prior to split them in order to generate at the end only 12 XSQ indexed files?

May I use merge Linux command to combine the 6 XSQ files or I may lose information?

Thanks,
Paolo

Last edited by paolo.kunder; 12-29-2011 at 05:22 AM.
paolo.kunder is offline   Reply With Quote
Old 01-02-2012, 06:01 PM   #2
genehunter
Junior Member
 
Location: California

Join Date: Aug 2009
Posts: 3
Default Xsq

No, you cannot combine them using linux commands. The files are binary.
You may want to use LifeScope or HDF (format) APIs.
genehunter is offline   Reply With Quote
Old 01-03-2012, 05:17 AM   #3
colindaven
Senior Member
 
Location: Germany

Join Date: Oct 2008
Posts: 300
Default

We have been doing alignments of all and recombining - merging - the SAM files with Picards MergeSamFiles.

This is admittedly not a great solution !

Let us know if you have any better ideas.
colindaven is offline   Reply With Quote
Old 01-03-2012, 11:30 PM   #4
paolo.kunder
Member
 
Location: Milano, Italy

Join Date: Aug 2011
Posts: 73
Default

Finally I converted my indexed XSQ files in csfasta (and QV.qual) files with XSQconverter and merged each individual csfasta (and QV.qual) with cat function,
information seems to be maintained,
paolo
paolo.kunder is offline   Reply With Quote
Old 01-04-2012, 04:40 AM   #5
colindaven
Senior Member
 
Location: Germany

Join Date: Oct 2008
Posts: 300
Default

Yeah, did that too. Bioscope and Lifescope wouldn't use all the reads for alignment, but just the first set.
Not sure why that was.

NovoalignCS doesn't seem to pick up on the different read sets and uses all reads.

Colin
colindaven is offline   Reply With Quote
Old 01-05-2012, 01:46 AM   #6
kevleb
Member
 
Location: Sophia-Antipolis

Join Date: Jun 2009
Posts: 10
Default

Quote:
Originally Posted by paolo.kunder View Post
Finally I converted my indexed XSQ files in csfasta (and QV.qual) files with XSQconverter and merged each individual csfasta (and QV.qual) with cat function,
information seems to be maintained,
paolo
When you cat files, you should obtain a final file with 6 times the same bead_id for different sequence comming from each of your 6 lanes. I'm not sure it won't be a problem then in the future bam file.
For external pipelines using .csfasta, instead of cat file, i open each .csfasta (and QV.qual) and update the panel number:
lane 1 : 1 to 708
lane 2 : 709 to 1416
lane 3 ...

In case of you use lifescope you should not do anything, import .xsq in lifescope and create reads set directly in lifescope. You can merge the .bam files from each lane with samtools after mapping. There you will have the bead_id issues in your final bam file. It seems not to be problem in case of only visualization in genome browser but it certainly depends on what you plan to do with the merge .bam file.

kevin.

Last edited by kevleb; 01-05-2012 at 02:00 AM.
kevleb 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 11:04 PM.


Powered by vBulletin® Version 3.8.6
Copyright ©2000 - 2014, Jelsoft Enterprises Ltd.