Question: Getting consensus fasta sequence for a chromosome from a bam file
gravatar for nagarsaggi
22 months ago by
nagarsaggi10 wrote:

I have a bam file per sample which I have generated by mapping Illumina short reads on ~1100 single copy orthologues genes (busco). I want to extract consense sequence from bam file for a given single copy orthologue, in another way for a given chromosome. The aim is to generate a multiple sequence alignment out of the consensus sequence of a particular single-copy ortholog gene.

Following command have been recommended for extracting consense from a bam file but I am not able to use it for a particular chromosome

samtools mpileup -uf reference.fa aligment.bam | bcftools view -cg - | vcfutils vcf2fq >out.fasta
alignment • 669 views
ADD COMMENTlink modified 22 months ago by finswimmer13k • written 22 months ago by nagarsaggi10

So essentially what you need is to get a single chromosome from the bam file? Have you looked at the manual of samtools?

ADD REPLYlink modified 22 months ago • written 22 months ago by WouterDeCoster44k
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: 743 users visited in the last hour