How to make consensus sequence from aligned reads ?
0
0
Entering edit mode
23 days ago
Info.shi ▴ 10

We already have sequence reads aligned to a reference genome. How can we extract the final reads that map to a specific gene? To make it clear I want to make a consensus sequence from the reads aligned to each gene. SAMTOOLS suggests the command:

samtools mpileup -uf ref.fa aln.bam | bcftools view -cg - | vcfutils.pl vcf2fq

But it gives many wearied characters after + sign of sequence.

I will thankful if anyone kindly guides me.

sequence Transcriptome aligned consensus reads • 159 views
ADD COMMENT
2
Entering edit mode

Generating consensus sequence from bam file

How to get consensus sequences in FASTA format from BAM

create consensus sequence from bam file

How To Generate A Consensus Fasta Sequence From Sam Tools Pileup?

etc.

I think the "wearied characters" you are talking about are fastq quality score values. Note that vcf2fq outputs fastq, not fasta..

ADD REPLY

Login before adding your answer.

Traffic: 1163 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6