GC & AT content calculation
1
0
Entering edit mode
5.1 years ago
anu014 ▴ 190

Hello Biostars,

Can anyone tell me how can I get GC & AT percentage in chip-seq peak bed file?
I have seen that GC content can be calculated using UCSC browser & EMBOSS suit itself. But I am not getting any tool for AT calculation.

Please suggest me some methods. Thank you :)

ChIP-Seq next-gen R genome • 2.9k views
ADD COMMENT
1
Entering edit mode

If you're able to get GC content, then AT content is just 1 minus that.

ADD REPLY
0
Entering edit mode

Thanks Devon. I realized it now..

ADD REPLY
2
Entering edit mode
5.1 years ago

You can use bedtools.

bedtools nuc -fi /data/genomes/hg38/hg38.fa -bed input.bed | cut -f1-5

#1_usercol  2_usercol   3_usercol   4_pct_at    5_pct_gc
chr1    28216   31238   0.465255    0.534745
chr1    51379   52104   0.478621    0.521379
chr1    136251  145899  0.473984    0.526016
ADD COMMENT
0
Entering edit mode

Thank you Simon. It's perfect!

ADD REPLY

Login before adding your answer.

Traffic: 1799 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