generating a genome reference fasta file
Entering edit mode
8 days ago
erl1977 ▴ 60

Hi all,

I have some files (CRAM) that were given to me that are ONT reads aligned to a genome. I would like to generate a 'personal' genome reference for each file. So, essentially, I would like to take the mappings and extract out a consensus fasta sequence for each sample, sort of like generating individual genomes in the pangenome projects. I'm playing around with using "samtools consensus", but not sure if this is the best approach, or if there is another better method out there.

I should mention, that in the output, I would like the new fasta file to incorporate any potential SNPs that may be present in the sample.

(I did try searching for approaches, but haven't been too successful in finding many solutions). Just interested if anyone has done anything similar and has any thoughts.

Thanks for any help.

FASTA BAM reference • 118 views

Login before adding your answer.

Traffic: 3635 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6