SEQanswers

Go Back   SEQanswers > General



Similar Threads
Thread Thread Starter Forum Replies Last Post
Contig vs contig or map against contig lib? JackieBadger Bioinformatics 1 05-30-2016 05:34 AM
Finding contig repeat counts by mapping contigs to the reference genome misagh Bioinformatics 1 01-09-2014 05:41 AM
contig mapping reference blat psl result handle Jane Lee Bioinformatics 0 06-19-2013 07:31 PM
reference-based contig assembly? lukas1848 Bioinformatics 5 02-03-2013 06:55 AM
SRMA Problem SAMRecord contig does not match the current reference sequence contig gavin.oliver Bioinformatics 5 07-05-2011 05:28 AM

Reply
 
Thread Tools
Old 04-23-2015, 01:03 PM   #1
gnarls
Junior Member
 
Location: Pennsylvania

Join Date: Apr 2015
Posts: 1
Default BWA-SW contig to reference -- what was unaligned?

I used BWA-SW to align 60 contigs to a reference genome. The resulting SAM file has 339 entries (I did `samtools view -c myfile.sam`). All of these are "mapped" (`samtools view -c -F 4 myfile.sam`).

So BWA-SW, while aligning, broke the contigs into on average 5 or 6 pieces, is what I understand from that. I can visualize the alignments in IGV and they look good. What I'm trying to find out is: are there any sub-segments of those 60 contigs that were not successfully aligned to the reference genome I passed to BWA? In other words I have 339 sub-segments and I'm not sure if they fully constitute the 60 contigs I originally passed to BWA-SW.

I've searched for how to figure this out as best I can, but it seems like the majority of use cases for BWA are for aligning short reads, which is not what I'm doing at the moment. So I can find info on how to count unmapped reads but it doesn't seem to apply to my situation.

Thanks....
gnarls is offline   Reply With Quote
Reply

Tags
bwa, sam

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 09:17 AM.


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