Analysis of Next Generation Sequencing Data

[pochi04]

baminfo: you can get the length of each chromosome.
bamread(File,RefSeq,Range): File is 'Sample.bam', RefSeq is chromosome 21, Range is from 1 to the length of chromosome 21. The output of bamread gives the information of how the short reads mapped to chromosome 21 regions.

[nilshomer]

What is your question?

[pochi04]

How to generate the number of reads mapped to each position of chromosome 21 in sample.bam file using matlab buit in function bamread.
eg:- bam = bamread('sample.bam', chr21, [1 21])