Question: Interpreting bam-files from Prokka annotation
0
gravatar for Hansen_869
10 months ago by
Hansen_86920
Hansen_86920 wrote:

Hi!

I have used Prokka on some bins, which I found with Maxbin2. On the annotated genes, I then used BWA together with the raw FASTQ files, to obtain a BAM file. Now, how do I interpret the bam file? The main goal is to find the depth for each of the genes that were annotated with Prokka.

Thanks!

bwa prokka samtools bam • 288 views
ADD COMMENTlink modified 10 months ago • written 10 months ago by Hansen_86920

Instead of posting basically the same question twice you might consider investing more effort into one of them, adding details and context. People seem to lack the necessary information to help you. Finding depth of genes

ADD REPLYlink written 10 months ago by ATpoint38k

Thanks for your response. In the first post I was asking for the approach of finding the depth, which I have now decided. Now I just need help interpreting the BAM-file. So i find the two posts quite different. But people are free to ask for further information if in doubt.

ADD REPLYlink written 10 months ago by Hansen_86920

It's basically the same question. I started to reply to the other one but wasn't sure if my advice was the best. There are several tools to count how many reads span each gene, from the simple htseq-count to RSEM and stringTie.

ADD REPLYlink written 10 months ago by Asaf8.4k

Thanks for your reply. I will take down the other post then, to avoid confusion. Do these tools just require a BAM-file as input?

ADD REPLYlink written 10 months ago by Hansen_86920
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1599 users visited in the last hour