Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

outputing statistics different with the qcreport #79

Open
Wenbin-Wang opened this issue Nov 7, 2019 · 2 comments
Open

outputing statistics different with the qcreport #79

Wenbin-Wang opened this issue Nov 7, 2019 · 2 comments

Comments

@Wenbin-Wang
Copy link

I am using hichipper to analyze the hichip data of GM12878 which is provided in the qc reports,but I have some questions about the result.I got a different result than the qc report except the Total_PETs.I used the peaks' option : EACH,SELF .

@caleblareau
Copy link
Contributor

Could you provide the report and indicate which one you are comparing to?

@Wenbin-Wang
Copy link
Author

Wenbin-Wang commented Nov 8, 2019

My report is as follows:
------------------------------------GM12878_Rep1---GM12878_Rep2
Total_PETs--------------------------332,279,257--------299,478,736
Mapped_unique_quality-----------285,550,786--------262,042,752
Mapped_unique_quality_valid-----107,434,251--------103,031,346
Valid_Intra---------------------------95,935,175---------92,121,438
Valid_Intra_<5Kb--------------------30,036,124---------28,320,268
Valid_Intra_5Kb-2Mb----------------54,808,566---------53,090,007
Valid_Intra_>2Mb--------------------11,090,485---------10,711,163
Valid_Intra_anchor_5Kb-2Mb--------9,737,321-----------9,457,383
Number of Anchors------------------66,819--------------66,819

The report which is provided in qcreport is as follows:
------------------------------------GM12878_Rep1---GM12878_Rep2
Total_PETs--------------------------332,279,257--------299,478,736
Mapped_unique_quality-----------286,612,971--------262,546,720
Mapped_unique_quality_valid-----107,833,087--------102,943,458
Valid_Intra---------------------------98,170,528---------93,725,004
Valid_Intra_<5Kb--------------------40,389,599---------37,473,831
Valid_Intra_5Kb-2Mb----------------48,060,076---------46,810,968
Valid_Intra_>2Mb--------------------9,720,853----------9,440,205
Valid_Intra_anchor_5Kb-2Mb--------2,937,509----------2,783,661
Number of Anchors------------------84,598-------------84,598

Is it because I use hichipper 0.7.7 but the report which is provided uses hichipper 0.7.1?

Thanks!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants