How parse collectReadCounts output in r?
0
0
Entering edit mode
19 months ago
Nemo • 0

I need to get read counts per gene. I am using GATK collectReadCounts method via following piece of code:

"gatk",
  "CollectReadCounts",
  "-I", paste0("./sorted-bams/D10.bam"),
  "-L", "chr1:11869-14409",
  "-imr", "OVERLAPPING_ONLY",
  "--format", "TSV",
  "-O", paste0("./collectReadCounts/DDX11L1.tsv")

It is running without any problem. I would like to parse the output in R (rstudio). How I can do it? I used read.table but gives me error and I think its not an appropriate method to do it. I would appreciate any idea.

CollectReadCounts r gatk rstudio • 426 views
ADD COMMENT
0
Entering edit mode

add a few lines from the output to your post, plus be specific as to the error that you are getting

ADD REPLY

Login before adding your answer.

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