GC-BIAS correction for BAM file
1
2
Entering edit mode
5.5 years ago
amitgsir ▴ 50

Hi,

I am looking for any software which can correct the GC-Bias in the BAM file and provides the output as BAM.

I know that deeptools (correctGCBias) does the exactly same job, but the output of the deeptools is variable in each run.

So, I am looking for more specific tool.

Basically, I want to count the sequence reads for each chromosome in NIPT (low depth, WGS sample) BAM file after the GC-Bias correction.

Many LOESS based or R-based GC-bias correction methods are available, but they provides the output as normalized data, not as total number of reads in each chromosome/ specified region.

Any suggestions?

genome RNA-Seq R deepTools • 2.9k views
ADD COMMENT
4
Entering edit mode
5.5 years ago

The reason correctGCBias can produce slightly different results each time is due to using a random number at one point. If you really need to produce the exact same results each time then you could hardcode a seed in. I expect that adding the following near the top of deeptools/correctGCBias.py will do that:

np.random.seed(1234)
ADD COMMENT
2
Entering edit mode

Yeah,

Fixing the seed using np.random.seed(1234) is giving the exact same results each time.

Thanks

ADD REPLY
1
Entering edit mode

Great, glad that worked!

ADD REPLY
1
Entering edit mode

Thanks Ryan, I will try this and will inform you if this works.

ADD REPLY
0
Entering edit mode

@ Devon This post of yours and the deeptools paper are both from 2016.

Fast forward to the present day (Jun 2020), and the advice I seek from you is whether correctGCBias from deeptools is still good enough a tool for RNA-Seq GC bias detection and correction?

If not, which software / tool would you recommend? ALPINE or something else?

Thanks, in advance.

ADD REPLY

Login before adding your answer.

Traffic: 2116 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6