Go Back   SEQanswers > Bioinformatics > Bioinformatics
Similar Threads
Thread Thread Starter Forum Replies Last Post
extracting reads from sam/bam file gene wise chris_bioinfo Bioinformatics 5 09-20-2013 02:09 AM
PubMed: Ultra-high resolution HLA genotyping and allele discovery by highly multiplex Newsbot! Literature Watch 0 08-08-2012 03:20 AM

Thread Tools
Old 03-14-2017, 08:11 AM   #1
Location: FL

Join Date: Sep 2012
Posts: 13
Question Extracting HLA-E allele sequences from BAM for gene editing

I have WGS data (2x100bp) that has been aligned to hg19. Using HLA-VBSeq, I HLA typed the data, and now the researcher I am working with is interested in two allele sequences for HLA-E for gene editing purposes. I would really appreciate some feedback on my workflow since I have never done this type of analysis before.

Using samtools, I subset the list of reads covering the genomic region of HLA-E, then used picard FilterSamReads to output the filtered bam file. I converted the bam to paired fastq and then aligned the fastq to the HLA-E genomic DNA sequences for alleles E*01:03:01:01 and E*01:01:01:01 from the IMGT/HLA database using BWA-MEM. I then sorted and indexed my bam files for each allele for viewing in IGV. Now, I have a couple questions:

1. Should I have aligned the filtered bam to the nucleotide coding sequences for both alleles instead of the genomic sequences?
2. The researcher wants all sequences aligned to E*01:03:01:01 and E*01:01:01:01. Would it be appropriate to convert the bam files to fasta files using seqtk to deliver the sequences?
samtools bam2fq input.bam | seqtk seq -A > output.fa
I'm not sure if my workflow is correct, or if delivering the fasta for each allele is appropriate for gene editing downstream. Help?
emily202 is offline   Reply With Quote

allele database, bam, hla typing, wgs

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 10:41 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2022, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO