Question: Extract gene sequence from BAM.
0
gravatar for Pedro Morell
2.9 years ago by
Elche
Pedro Morell • 0 wrote:

Hi, I have several BAM files with different samples of a chromosome and I want to end up with one file for each sample (two if they are paired) including only a certain region of that chromosome in order to study the variance. I've managed to convert the original BAM I'm using as a test into FASTQ, perform the quality asses, that I already knew was correct, and trim the low quality areas on the edges of each sequence, but I don't know how to cut from the whole chromosome to my target gene sequence. 

Any advise?

fastq sequence gene bam • 1.2k views
ADD COMMENT • link • modified 2.9 years ago by Jennifer Hillman Jackson ♦ 25k • written 2.9 years ago by Pedro Morell • 0
2
gravatar for Jennifer Hillman Jackson
2.9 years ago by
United States
Jennifer Hillman Jackson ♦ 25k wrote:

Hello,

Tools in the groups NGS: SAMtools, NGS: Picard, and BEDtools are the most commonly used to manipulate BAM content (before extracting sequence or other operations). Try Slice BAM first for your case. 

Thanks, Jen, Galaxy team

ADD COMMENT • link written 2.9 years ago by Jennifer Hillman Jackson ♦ 25k

I've been working with NGS:SAMtools and NGS:Picard. My first step was to pas from BAM to SAM, and then to FASTQ, a format I'm used to use, but I don't know how to get just the region I want to study. I'll try Slice BAM and see.

Thanks.

ADD REPLY • link written 2.9 years ago by Pedro Morell • 0
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 16.09
Traffic: 168 users visited in the last hour