samtools : print single columns at position

How can I print only 1 column from a BAM file at specific coordinate.

For example:

ACGTGACGGACGGGGCCA
   TGACGGACGGGGCCA
       GGACGGGGCCA
         ^

For example print the column at coord 'x' will return AAA (see example above)

print position samtools

• 157 views • link

updated 4.1 years ago by Ram

• written 11.3 years ago by sacha

0 answers

No answers yet.

More posts like this

Extracting specific cancer types from MAF (Mutation Annotation Format) files from TCGA(The Cancer G…
written 5.4 years ago by peter •

I am new to bioinformatics therefore, I apologise I have two basic questions. But I have a MAF file. It is mutational information available for …
Get the top X number of lines per unique value in one column, once you've sorted a text file using …
written 5.6 years ago by niccolo.alfano •

Hi everybody, I have a text file with 19 columns (divided by tab) which I have sorted using a command such as: sort -t$'\t' -k1,1 …
Extracting only reads mapped to specific genomic region
written 5.7 years ago by agata •

I have a specific genomic region saved as `fasta`. I want to check how many reads mapped to this specific region and extract them, as …
How to merge chr+strand+coordinate+score from multiple files?
written 8.2 years ago by Rae •

I have many files (~50) with the same format of four columns: chr, strand, coordinate, score Each file has about 1 Million lines. Now I …
reads y-position in a bam viewer
written 8.4 years ago by sacha

Hi, I m doing a BAM viewer. I wonder how to set read position in vertical axis ? I mean, given N reads , How …
Tutorial: How to display Agilent CGH array in IGV
written 8.7 years ago by sacha

CGH array output file from Agilent technology contains several lines per probe position. From line 11, each line contains all information available for one probe. …
how to construct CDS by using annotation file (coordinates) properly ???
written 8.9 years ago by adeena_hassan •

Assalam o aliakum everyone, I have a BAM file of dog genome and I have generated consensus FASTA from it. BAM is aligned against Canfam3.1 …
remove a percent of reads from a BAM file
written 10.0 years ago by sacha

I would like to remove x% of reads from a BAM file at a specific region . For instance : remove 50% of reads in …
Reads in BAM file overlapping coordinates
written 10.3 years ago by madhavkurup1991 •

Hi guys, I'm trying to create a program that can count the number of reads in a BAM file that overlap a specific genomic coordinate(for …
samtools : print single columns at position
written 11.3 years ago by sacha

How can I print only 1 column from a BAM file at specific coordinate. For example: ``` ACGTGACGGACGGGGCCA TGACGGACGGGGCCA GGACGGGGCCA ^ ``` For example print …