Hi everyone,
I am a newbie to the NGS studies. I have PE data (fastq files), and after performing mapping with bowtie2 I got bam files. Then I extracted mapped reads using samtools command:
samtools view -b -F 4 file.bam > mapped.bam
So now I am having bam file of mapped reads. Are there any tools or programs available to extract information which I need?
I need two things in a single text file:
- First Column: a list (a text file) of mapped reads(IDs) and
- Second Column: a list of regions on the reference genome where the reads are mapped.
For example, "this" is the read or read id which is mapped to "this region" on the reference genome.
read id1 , region 1
read id2 , region 2
So how can I achieve the same? Any help is highly appreciated.
Thank you!
•
written
by
SV •
2 answers
remove the '-b' option from samtools. It will produces a text/sam file rather than a binary/bam file. use cut to extract the columns read-name,read-chr,read-pos
or use http://bedtools.readthedocs.io/en/latest/content/tools/bamtobed.html
Pysam (http://pysam.readthedocs.io/en/latest/api.html) lets you freely parse the BAM file directly and extract the information you need. Downside is you will have to know some python to use it.
Log in to answer this question.
-
How to Extract Specific Region on Bam file
written by santos48 •
-
Extract reads mapped within chosen intronic region from BAM file (no junctions)
written by gspirito •
-
How to join two txt files in unix
written by xiaoyonf •
-
meta file table
written by figueroacg •
-
Samtools view Exec format error Help
written by aman-coasab •
-
extracting orphan reads from bam file
written by Assa Yeroslaviz
-
Extract parts of reads from BAM file that overlap a specific region of genome
written by Chadi Saad
-
Multiple sam to bam followed by raw read count
written by Bioinfonext
-
Filtering nth length from mapped SAM/BAM file.
written by unique379
-
Can I Merge The Good Quality (Unpaired Pe ) Bam File And Paired Pe Bam File With Samtools?
written by waswangs •
Why do you want that, what's your end goal?
further I have to do statistical analysis to compare data from different samples.
That is probably the vaguest possible reply. What I'm getting at is that there's likely a much much more efficient way to go about things, but you would have to provide enough details about your goals for us to help you find out.