Working With Fasta And Bam
1
0
Entering edit mode
9.4 years ago
Pellix • 0

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

c fasta bam • 2.1k views
ADD COMMENT
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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 REPLY
0
Entering edit mode
9.4 years ago
Pellix • 0

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

Thanks

ADD COMMENT
2
Entering edit mode

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

ADD REPLY

Login before adding your answer.

Traffic: 2131 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6