How to apply bam file into assembly sequence
2
0
Entering edit mode
22 months ago
nonaddldy ▴ 10

http://p2.so.qhimgs1.com/bdr/300_115_/t02d3afc4c46284a492.jpg

I used the bam file into UCSC Genome Browser, then I want to assembly these reads into the sequence at chromosome level, what should I do?

I mean how to extract the reads according to the nucleotide order of GRCh38 reference sequence

sequence assembly 1000 human genomes • 758 views
ADD COMMENT
0
Entering edit mode

I want to assembly these reads into the sequence at chromosome level,

it's not clear to me.

ADD REPLY
0
Entering edit mode

A BAM file typically contains the genome location where individual (short) reads have aligned to. This seems to have been successful in your case as the snapshot shows.

What exactly is it that you're after? Do you simply want a text file with the ACGTs that make up, say, chromosome 1? That would be the reference sequence, i.e. the FASTA file that you used to generate the BAM file in the first place.

ADD REPLY
0
Entering edit mode

I want to assembly an individual human chromosome Y rather than reference sequence, so I should sort the reads and output the sequence according the order of reference human Y, while the reads in the bam file are in disorder, thus, my question is how sort the reads aligned to the chromosome Y by the nucleotide order.

ADD REPLY
0
Entering edit mode

You can use samtools sort -o sorted.bam input.bam on your aligned BAM. Samtools does co-ordinate sort by default.

ADD REPLY
0
Entering edit mode

Thanks a lot

ADD REPLY

Login before adding your answer.

Traffic: 2344 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