SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
How to extract uniquely mapped paired end reads from bam file serenaliao Bioinformatics 1 08-21-2014 02:02 PM
How to extract uniquely mapped paired end reads from bam file serenaliao RNA Sequencing 0 08-20-2014 10:42 AM
Convert Paired End Reads BAM file to GRanges ayushraman Bioinformatics 1 10-10-2012 03:06 AM
Filter paired & mapped reads from SAM/BAM file pirates.genome Bioinformatics 2 06-19-2012 01:21 AM
Filter paired end BAM file based on iSize Leif Bergsagel Bioinformatics 2 12-16-2010 11:50 AM

Reply
 
Thread Tools
Old 09-16-2017, 02:06 AM   #1
Vca80553
Junior Member
 
Location: Stockholm

Join Date: Sep 2017
Posts: 5
Default Filter paired end reads by ID in bam file

Hello everyone,

I am new in bioinformatics, so I might have to say sorry if the question is a bit odd. I have looked everywhere and didn't find a solution for it.

I am mapping my paired end reads to a viral genome(Nextgenmap). In the bam file, I extracted those that were both mapped and were paired end (I discarded those mapped with mate not unmapped, so I only have paired end data).

With this viral_paired_end_mapped.bam, what I want to do is filter out possible human reads.

What I did is to convert the viral_paired_end_mapped.bam into fastq file (R1 and R2) and then map it again to hg19. And in this final bam file, I filter only the unmapped paired end (Viralmapped_Humanunmapped paired end). Those are the sequences that I want (mapped to virus and not to human, paired end)

I have to continue now to do coverage analysis,etc. So I need to merge several bam files, according to sample id. My question is then:

Is it possible to merge the "Viralmapped_Humanunmapped paired end.bam"and then do a coverage analysis with the viral reference?
Or should I somehow go back to my first bam file (viral_paired_end_mapped.bam) and somehow filter out the sequence IDs from "Viralmapped_Humanunmapped paired end.bam"?

If so, does anyone know how to do it?
Vca80553 is offline   Reply With Quote
Reply

Tags
bam, filter reads, read id

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 04:43 PM.


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