Extract base from bam file

[empyrean]

I have a bam file and i wanted to extract the following information.

For a read in a bamfile, can i extract the base based on the position? For example, if i say i need the base at 100th position and it should output whether its A/T/G/C. Can i get this using samtools or bamtools??

[pbluescript]

Try this:

samtools view file.bam | awk '{print substr($10, 100, 1)}'

[empyrean]

thank you.. its doing exactly as i asked .. but looking at the output i realized that i asked my question wrongly.

I guess i should have made my question clear..sorry about that.. So once the read maps to the reference, i wanted to give input a position based on the reference say 100th position on reference, i need all bases from the reads which are mapping at that particular position including the readname. if i use the solution below, its printing the 100th base on the reads which is not equal to what i need because if there are indels in reference or reads, then the position of bases changes completely. So is there anyway to get what i need? i tried with mpileup. Its printing out all the bases at that positions but not the readnames.

[SeekAnswers]

Something like this?

samtools view <in.bam> 'chr1:start_pos-end_pos' | awk '{print $1"\t"$10}'

[pbluescript]

thank you.. its doing exactly as i asked .. but looking at the output i realized that i asked my question wrongly.

I guess i should have made my question clear..sorry about that.. So once the read maps to the reference, i wanted to give input a position based on the reference say 100th position on reference, i need all bases from the reads which are mapping at that particular position including the readname. if i use the solution below, its printing the 100th base on the reads which is not equal to what i need because if there are indels in reference or reads, then the position of bases changes completely. So is there anyway to get what i need? i tried with mpileup. Its printing out all the bases at that positions but not the readnames.

You should look into bedtools. You can use intersectBed to get the reads that overlap with any defined interval, including single bases in a vcf file.