mskilab-org / fragCounter

GC and mappability corrected fragment coverage for paired end whole genome sequencing
MIT License
7 stars 11 forks source link

Correcting coverage for both hg19 and GRCh37 #5

Open kcygan opened 4 years ago

kcygan commented 4 years ago

Just a suggestion. For people who mapped reads to hg19 (i.e. "chr" included in seqnames), it would be sufficient to add one line: cov = gr.sub(cov)

right after line 627 https://github.com/mskilab/fragCounter/blob/575af9926e5177a39b45a31ad37048953a680ca4/R/fragCounter.R#L627