SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Merge sam/bam file genelab Bioinformatics 11 06-04-2013 04:33 PM
a question about merge bam files camelbbs Bioinformatics 2 10-24-2011 09:00 AM
perl -reading from 2 files and merge them robelb4 General 5 10-07-2011 04:51 AM
Merge individual vcf files francy Bioinformatics 5 06-21-2011 02:10 AM
Can we merge 2 csfasta files ? tdm SOLiD 9 12-10-2010 09:10 AM

Reply
 
Thread Tools
Old 10-19-2011, 10:19 PM   #1
pengchy
Senior Member
 
Location: China

Join Date: Feb 2009
Posts: 116
Default how to merge sam files from mapping the same fq to different reference

Dear all,

I want to know if there some utilities to merge sam files produced from mapping the same reads file to the different reference sequences. The necessity occurred when the reference sequence so large that the mapper can not process.

Thank you in advance.

pengchy
pengchy is offline   Reply With Quote
Old 10-20-2011, 12:19 AM   #2
maubp
Peter (Biopython etc)
 
Location: Dundee, Scotland, UK

Join Date: Jul 2009
Posts: 1,540
Default

That is tricky.

What would you expect to happen when a read is mapped nicely to both reference files?

Are you using paired end data? Consider the case where each read has mapped nicely to a different reference.

P.S. How long are your reference sequences? I'm trying to encourage the SAM/BAM community to think about supporting very large chromosomes as in plants.
maubp is offline   Reply With Quote
Old 10-20-2011, 02:07 AM   #3
pengchy
Senior Member
 
Location: China

Join Date: Feb 2009
Posts: 116
Default

thank maubp for your reply.

Because the reads can only be produced from one place, so it is expected giving only one randomly. We care about the unique mapped reads or the randomly one place for the multimapped reads. So, if the both ends uniquely mapped to one place, it will be not a matter. For the multi hit pe-reads, the proper paired position, or randomly one if there are many, will be expected.

This is direction will be a question for the future large genome that current aligner can not tackle.
pengchy is offline   Reply With Quote
Old 10-20-2011, 02:25 AM   #4
maubp
Peter (Biopython etc)
 
Location: Dundee, Scotland, UK

Join Date: Jul 2009
Posts: 1,540
Default

Would you also expect some read should map across the break point? Or can you choose where to break the long reference sequence (e.g. a large region NNNN from scaffolding, or an area annotated with no genes).

[I am assuming you have had to break up a long chromosome - things are easier if you just have divided up the chromosomes into one FASTA reference file per chromosome]
maubp 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:49 PM.


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