SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Extracting Reads from a BAM/SAM file that map back to a subset of contigs ediane Bioinformatics 2 03-01-2017 01:24 PM
Counting Reads in BAM file per genomic position C9r1y Bioinformatics 4 11-24-2015 02:11 PM
Lack of paired mate in BAM file PJohn Bioinformatics 0 07-02-2015 02:36 PM
Weird PNEXT and RNEXT values in Bam from BWA sudders Bioinformatics 1 10-29-2013 12:41 PM
Can you read a bam file by mate pairs? kgulukota Bioinformatics 5 08-29-2013 10:12 AM

Reply
 
Thread Tools
Old 07-01-2019, 04:58 AM   #1
Notorious
Junior Member
 
Location: Italy

Join Date: Jun 2019
Posts: 1
Default BAM file: RNEXT read mate map to another position

Dear All,

I am analyzing a PE RNAseq experiment . Looking at the resulting sorted bams, I have noticed that some of the aligned reads map on another reference (e.g. in the attached figure, the other read map on chrM).

Any idea on what can cause this mismatch? Could it be raise problems during reads counting?


Bam files are obtained as follows:

PE fastq are checked with FASTQC (good overall quality)
PE fastq are aligned to mm10 genome index using hisat2 (mean overall alignment rate >97%)
samtools is used to convert sams to bams, sort and index them.

Thank you very much for your help.
Attached Images
File Type: jpg Example.jpg (89.7 KB, 6 views)
Notorious 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 05:09 PM.


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