Question: Pangenome from bam files
0
gravatar for ssj262
6 months ago by
ssj2620
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 • 215 views
ADD COMMENTlink written 6 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 6 months ago by jrj.healey11k
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: 760 users visited in the last hour