This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to get the gene symbols for my bam file or bed file

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

genome alignment sequencing next-gen
• 3,232 views • link
written by vinaysbharadhwaj • vinaysbharadhwaj

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?

• 1 views • link
written by Devon Ryan Devon Ryan

just curious: what would you with that information ?

• 1 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

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

• 1 views • link
written by vinaysbharadhwaj • vinaysbharadhwaj

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

• 1 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

I strongly suspect that featureCounts would be easier then.

• 1 views • link
written by Devon Ryan Devon Ryan

Can you explain what your final goal is?

• 1 views • link
written by WouterDeCoster WouterDeCoster

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

• 1 views • link
written by vinaysbharadhwaj • vinaysbharadhwaj

0 answers

No answers yet.

Log in to answer this question.

More posts like this