Question: How to get the gene symbols for my bam file or bed file
0
gravatar for vinaysbharadhwaj
3 months ago by
vinaysbharadhwaj10 wrote:

Hello, I am new to NGS and wanted to know which genes are present in my aligned file. I have performed BWA MEM to obtain a sam file and used picard to get a sorted.bam file. I also used bedtools to get a BED file and used its annotate utility to find the gene id or gene symbol for my reads. But it was only giving my the coverage. How do I get my gene symbols or gene ids using bedtools or any other tool?

Please Help

ADD COMMENTlink written 3 months ago by vinaysbharadhwaj10
1

Can you explain what your final goal is?

ADD REPLYlink written 3 months ago by WouterDeCoster29k

I want to get the genes that are present in my bacterial isolate and check with a database if that gene is virulent or rather the bacteria is virulent

ADD REPLYlink written 3 months ago by vinaysbharadhwaj10

Do you want the number of reads/pairs aligning to each gene or do you want to annotate each read with the gene it aligns to?

ADD REPLYlink written 3 months ago by Devon Ryan81k

I want to annotate each read with the gene it aligns to

ADD REPLYlink written 3 months ago by vinaysbharadhwaj10

just curious: what would you with that information ?

ADD REPLYlink written 3 months ago by Pierre Lindenbaum108k

I want to check if there are any virulent genes present in my bacterial sample

ADD REPLYlink written 3 months ago by vinaysbharadhwaj10
1

I strongly suspect that featureCounts would be easier then.

ADD REPLYlink written 3 months ago by Devon Ryan81k

Thank you so much it worked very well

ADD REPLYlink written 3 months ago by vinaysbharadhwaj10

so as said Devon , what about just counting " the number of reads/pairs aligning to each gene" ?

ADD REPLYlink written 3 months ago by Pierre Lindenbaum108k
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: 1590 users visited in the last hour