Hi,
I have a reference genome and reads from a cancer cell line aligned to that reference in a BAM and SAM file. I need to convert the BAM file into a FASTA file where the cancer consensus sequence is aligned to the reference genome so I can identify . Other threads on BioStars addressing this question explain how to generate FASTQ and FASTA files based off the BAM file, but don't explain how to align those sequences to the reference. Thanks so much in advance!
•
written
by
ryan.englander •
1 answer
To align your consensus sequence to your reference sequence you will need to use an alignment tool to create a multiple sequence alignment, such as:
MAFFT, MEGA, CLUSTAL, and loads of others
These tools don't care that one is a consensus, and one is reference, it just aligns all the sequences you give it together
Log in to answer this question.
-
Question about `vg giraffe`
written by zhengluo •
-
Create a reference genome from aligned bam file
written by Vincent •
-
How to align and visualize data with .fasta and .gff3 files in IGV?
written by Phuc •
-
Aligning PacBio hifiasm CCS to reference genome using pbmm2
written by K •
-
Generate consensus from contig alignment to reference sequence
written by daewowo •
-
filtering alignned seqeunces
written by priscilladraj •
-
How to remove bam files that don't contain matching reads
written by samlambrechts299
-
Consensus sequence with characters apart from ATGC
written by sambioinfo2018 •
-
Generation of Consenses Fasta from SAM/BAM
written by mail2steff
-
Convert BAM format file to FASTA format
written by rodbari.zahra •
