Question: Interpreting bam-files from Prokka annotation
0
gravatar for Hansen_869
6 weeks ago by
Hansen_86910
Hansen_86910 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 • 129 views
ADD COMMENTlink modified 6 weeks ago • written 6 weeks ago by Hansen_86910

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 6 weeks ago by ATpoint26k

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 6 weeks ago by Hansen_86910

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 6 weeks ago by Asaf6.5k

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 6 weeks ago by Hansen_86910
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: 1322 users visited in the last hour