Can I extract a complete sequence of gene from mapped reads bam file

0

Entering edit mode

9 months ago

singhankit973 • 0

Can I extract gene sequences from mapped reads?

I have mapped Illumina paired WGS against a reference genome. By using locus ID of gene I extracted a bam file (gene.bam) of that gene from my mapped reads. Now what should I do next to extract complete sequence of gene from gene.bam file.

I have Locus ID, protein ID, Chromosomal location of gene in reference genome and reference genome annotation file, and the size of reads is 90 bp.

bam illumina mapping • 726 views

ADD COMMENT • link updated 9 months ago by rfran010 ▴ 900 • written 9 months ago by singhankit973 • 0

0

Entering edit mode

see FastaAlternateReferenceMaker https://gatk.broadinstitute.org/hc/en-us/articles/360037594571-FastaAlternateReferenceMaker

ADD REPLY • link 9 months ago by Pierre Lindenbaum 161k

0

Entering edit mode

I'm sorry, I couldn't understand, i'm new to this, can you elaborate what you want to say?

ADD REPLY • link 9 months ago by singhankit973 • 0

0

Entering edit mode

What have you tried so far?

ADD REPLY • link 9 months ago by swbarnes2 14k

0

Entering edit mode

It sounds like you want to get the consensus sequence for a particular locus? Does this thread help you? Generating consensus sequence from bam file or maybe http://www.htslib.org/doc/samtools-consensus.html ?

ADD REPLY • link 9 months ago by rfran010 ▴ 900

Login before adding your answer.