DoaneAS / rseqc

Automatically exported from code.google.com/p/rseqc
0 stars 0 forks source link

same number reads for both read 1 and read2 #25

Open GoogleCodeExporter opened 9 years ago

GoogleCodeExporter commented 9 years ago
What steps will reproduce the problem?
1. bam_stat.py
2.
3.

What is the expected output? What do you see instead?
Different reads for read 1 and read 2

What version of the product are you using? On what operating system?

Rseqc 2.3.7
Please provide any additional information below.
What does non-primary hits mean?

Why do I have same number of reads on  minus and plus strands?

#==================================================
Total Records:                38534316

QC failed:                    0
Optical/PCR duplicate:        0
Non Primary Hits              6308104
Unmapped reads:               0
Multiple mapped reads:        3067306

Uniquely mapped:              29158906
Read-1:                       14579453
Read-2:                       14579453
Reads map to '+':             14579453
Reads map to '-':             14579453
Non-splice reads:             25139832
Splice reads:                 4019074
Reads mapped in proper pairs: 29158906

Original issue reported on code.google.com by biotechd...@gmail.com on 3 Jul 2013 at 10:25

GoogleCodeExporter commented 9 years ago
I also want to know.
I just find the explanation as below, but I can't understand the exact 
biological meaning of non primary hits:
non primary hits:a read having split hits may have multiple primary alignment 
records

Original comment by liucuihu...@gmail.com on 13 Nov 2013 at 8:30