Question: Working With Fasta And Bam
0
gravatar for Pellix
7.9 years ago by
Pellix0
Pellix0 wrote:

I have to open a fasta file and match it with a bam file in c++. Unfortunately i open the bam file correctly, but i have no idea how to open the fasta file and match it with the bam. Before beginning the matching, i havo to take care of some informations of the bam file: the quality of the bases, the read quality, the orientation of the read, the initial and final position of the read and the coverage of the reads alignments. After the matching, i would be able to extract the positions where the bam differs from the reference fasta, with the bases differences in those positions. I know i'ts a lot of stuff, but i come from perl and i'm not as good in c++. Thanks in advance

fasta bam C • 1.8k views
ADD COMMENTlink modified 7.9 years ago • written 7.9 years ago by Pellix0

If you still need help, write me examples of files and of desired output.

ADD REPLYlink written 7.9 years ago by Biomonika (Noolean)3.1k

the fasta file is taken form here ftp://hgdownload.cse.ucsc.edu/goldenPath/hg19/chromosomes/ . it's the hg19 human genome. with a bam sorted file in input i want to match it with the fasta and the program will create a text file taking the positions where the a base of the bam file differs from the fasta one. and of course it will sign what is the difference, and if the mutation is heterozygosis or homozygosis. thanks in advance, all help will be very welcome

ADD REPLYlink written 7.9 years ago by Pellix0
0
gravatar for Pellix
7.9 years ago by
Pellix0
Pellix0 wrote:

Thanks, i didn't see it. But how to open the FASTA and work with it?

Thanks

ADD COMMENTlink written 7.9 years ago by Pellix0
2

this should be a comment, not an answer. See google "parsing fasta C"

ADD REPLYlink written 7.9 years ago by Pierre Lindenbaum127k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1922 users visited in the last hour