consensus fasta and multifasta file of human genomes

[huma Asif]

Dear All,
i desperately need help in converting a thought to work.

i have sequence six human samples and have created six alignment file (BAM) using UCSC reference.
i need to do multiple alignment of all these samples how can i do that

how i can convert every bam to consensus fasta
how i can merge six consensus fasta into multifasta
any software or reference will be appreciated

[GenoMax]

From Samtools manual page:

Generate the consensus sequence for one diploid individual:

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

Use seqtk to convert from fastq to fasta.

Concatenate (end to end, you are not "merging" the files this way) the fasta files with "cat" to generate multifasta.

If you truly want to do whole genome alignments: http://genomewiki.ucsc.edu/index.php...lignment_howto

[huma Asif]

thank you so much will try now

[huma Asif]

well i am done with first step but look at my fastq file with many NNNNNNNNNN
is this because my data is targeted exome sequencing
how can i get rid of it
if i use this fastq in seqtk it gives me empty output

[vivek_]

If your data is exome sequencing use the -L option in samtools and specify your capture regions file in the mpileup step.

[huma Asif]

give the same result with -L

see my code
samtools mpileup -uf human_hg19.fa reads.bam -L targetSeq_exome_target_regions_hg19.bed |bcftools view -cg - |/usr/local/genome/samtools-0.1.18/bcftools/vcfutils.pl vcf2fq > cns.fastq



fq with many n and empty fasta with
seqtk seq –a cns.fastq > cns.fa

[GenoMax]

That should have been a lower case "l".

-l FILE list of positions (chr pos) or regions (BED) [null]

-L means something else for samtools mpileup.

[huma Asif]

thank you so much i am going to try this now

[huma Asif]

tried with small l still n in final output please suggest where i can get wrong

see my code
samtools mpileup -uf human_hg19.fa reads.bam -l targetSeq_exome_target_regions_hg19.bed |bcftools view -cg - |/usr/local/genome/samtools-0.1.18/bcftools/vcfutils.pl vcf2fq > cns.fastq