# Input for Gene Body Coverage (BAM) **URL:** \ **Category:** usegalaxy.org support\ **Created:** [March 8, 2021, 10:36am UTC](https://help.galaxyproject.org/t/input-for-gene-body-coverage-bam/5549 "2021-03-08T10:36:34Z")\ **Posts on this page:** 5\ **Page:** 1
**Author:** ![stephen\_li](https://avatars.discourse-cdn.com/v4/letter/s/7c8e57/32.png) [@stephen\_li](https://help.galaxyproject.org/u/stephen_li)\ **Post date:** [March 8, 2021, 10:36am UTC](https://help.galaxyproject.org/t/input-for-gene-body-coverage-bam/5549/1 "2021-03-08T10:36:34Z")
Hi, I have been trying to generate a gene coverage plot from my bam file using **Gene Body Coverage (BAM)** from [Galaxy.org](http://Galaxy.org), but keep getting results that I am not expecting (see the attached heatmap png). While the png file suggested all samples have a shart peak near the 40th percentile, it is very different from what I see when I visualise my data on IGV. So I am wondering if my input bed12 file is in the wrong format. I have attached a screenshot of my bed12 file (see the screenshot png below). Would you mind checking if the bed12 file is in the right format? ![output.geneBodyCoverage.heatMap](https://us1.discourse-cdn.com/flex020/uploads/galaxy/original/2X/1/1fbff50b33a55038abd465bd9a703c5086fbea00.png) ![Screen Shot 2021-03-08 at 10.35.00](https://us1.discourse-cdn.com/flex020/uploads/galaxy/original/2X/b/bb314cdb19e11fbdcee152ea9be98eec85b77238.png) best wishes, Stephen ---
**Author:** ![gallardoalba](https://sea2.discourse-cdn.com/flex020/user_avatar/help.galaxyproject.org/gallardoalba/32/1903_2.png) [@gallardoalba](https://help.galaxyproject.org/u/gallardoalba)\ **Post date:** [March 8, 2021, 6:35pm UTC](https://help.galaxyproject.org/t/input-for-gene-body-coverage-bam/5549/2 "2021-03-08T18:35:47Z")
Hi @stephen_li, the format is correct, but that signal seems to be the result of some kind of artifact since it is too narrow and homogeneous to be an over-represented transcript. Could you analyze [the read distribution across features](https://training.galaxyproject.org/training-material/topics/transcriptomics/tutorials/rna-seq-reads-to-counts/tutorial.html#read-distribution-across-features-exons-introns-intergenic)? I could provide some additional information. Regards ---
**Author:** ![stephen\_li](https://avatars.discourse-cdn.com/v4/letter/s/7c8e57/32.png) [@stephen\_li](https://help.galaxyproject.org/u/stephen_li)\ **Post date:** [March 9, 2021, 6:12am UTC](https://help.galaxyproject.org/t/input-for-gene-body-coverage-bam/5549/3 "2021-03-09T06:12:44Z")
Hi, Sure. Heres the link to the read distribution result : [https://usegalaxy.org/datasets/bbd44e69cb8906b59133118c8ae76acd/display?to\_ext=html](https://usegalaxy.org/datasets/bbd44e69cb8906b59133118c8ae76acd/display?to_ext=html) I am puzzled by the fact that none of the reads are classified as “CDS\_Exons”. Wonder if that may provide some hint to the strange result from the previous post? The full input bed12 file is attached in the screenshot below ![Screen Shot 2021-03-09 at 06.11.49](https://us1.discourse-cdn.com/flex020/uploads/galaxy/original/2X/b/b5f95ffc97147ca5341be26df7dbe01cb47e72d1.png) best wishes, Stephen ---
**Author:** ![stephen\_li](https://avatars.discourse-cdn.com/v4/letter/s/7c8e57/32.png) [@stephen\_li](https://help.galaxyproject.org/u/stephen_li)\ **Post date:** [March 10, 2021, 6:18am UTC](https://help.galaxyproject.org/t/input-for-gene-body-coverage-bam/5549/4 "2021-03-10T06:18:11Z")
Hi, I have just done a quick test using just one RNASeq sample for one gene. The IGV screenshot of the selected gene (rnpB) is shown as following. It shows the data should be skewed towards the 3’end ![rnpB_IGV](https://us1.discourse-cdn.com/flex020/uploads/galaxy/original/2X/4/4b33b254c374d18341cad33825d686ec365412a1.png) . However, the output for gene body coverage showed one sharp peak in about the 50th percentile ![rnpB_(Curves_pdf)](https://us1.discourse-cdn.com/flex020/uploads/galaxy/original/2X/1/1ce4836bf0ef23839566aaf3e4b8ba4064d406d4.png) Is that suggesting error from my input files? best wishes, Stephen ---
**Author:** ![gallardoalba](https://sea2.discourse-cdn.com/flex020/user_avatar/help.galaxyproject.org/gallardoalba/32/1903_2.png) [@gallardoalba](https://help.galaxyproject.org/u/gallardoalba)\ **Post date:** [March 10, 2021, 12:36pm UTC](https://help.galaxyproject.org/t/input-for-gene-body-coverage-bam/5549/5 "2021-03-10T12:36:14Z")
According to the features of your BED12 file, it seems that the transcripts correspond to non-coding RNAs (tracrRNA, crRNA, 6S RNA). It can explain the fact that none of them are classified as CDs. Which kind of RNA-seq are you using?