Question: Pangenome from bam files
gravatar for ssj262
23 months ago by
ssj2620 wrote:

I have whole genome sequences of 300 bacterial isolates and i want to do the gene presence absence analysis. After quality and adapter trimming, I have been able to generate alignment bam files from the reference genome for all isolates. One of the options for looking at gene presence absence is to do pangenome analysis, but i believe for that one needs whole genome fasta files. I wanted to know how a single fasta file with the consensus sequence can be generated from sam/bam files so that it can be used in a pangenome analysis tool like ROARY. I am new to NGS analysis and I will be very happy to try other approaches or options out there. Thank you for reading.

bam pangenome wgs sequence • 493 views
ADD COMMENTlink written 23 months ago by ssj2620

This doesn’t seem like the most sensible way to do this.

I would do the following:

  1. Convert your BAMs to fastq files.
  2. Assemble all the FASTQs for each bacterial genome
  3. Annotate your assemblies with Prokka
  4. Give those GFFs or genbanks directly to Roary.
ADD REPLYlink written 23 months ago by Joe17k
Please log in to add an answer.


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