Consensus sequence from SAM or BAM file?
I am trying to perform reference-based assembly. Most of the tutorials teach how to create a bam file and view alignemnts in IGV or Tablet. But, I want a assembled genome sequence in fasta format. How can I do that ? Is there any tool or linux syntax available for this purpose ?
• 2,463 views
•
link
updated
by
Biostar
•
written
by
lokraj2003
•
written
by
lokraj2003
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
How to subset large BAM files specifically/ extract specific subsets?
written by ella •
-
Converting .delta file from NUCMER to format accepted by IGV? (.bam/.vcf)
written by Joanna •
-
IGV Fasta Error: Uneven Line Lengths
written by eds35016 •
-
Creating SingleCellExperiment class from scratch !
written by lokraj2003
-
Seeing only "N's" in tablet genome viewer?
written by tuj92050 •
-
FASTA to SAM for Tablet Assembly Visulisations
written by peflanag •
-
How can I create LASTdb for nr.gz data?
written by MAPK
-
BWA align contigs with Clustal-like output
written by biotech
-
View annotations in igv of bam file with multiple chromosomes
written by yellowbike •
-
How To Use Tablet To Visualization Bam File?
written by Haiping •
Generate consensus from BAM file
Generate Consensus sequence from BAM file ;
Hello lokraj2003!
It appears that your post has been cross-posted to another site: https://bioinformatics.stackexchange.com/questions/7441/consensus-sequence-from-sam-or-bam-file
This is typically not recommended as it runs the risk of annoying people in both communities.
Assembling a genome is not what you get by aligning to a reference. For that, you need an assembler.
Creating a consensus from an alignment will also give you a fasta consensus, but it is not the same as a real de novo assembly.