Question: remove low number of reads from BAM file
0
gravatar for Mehmet
3.7 years ago by
Mehmet460
Japan
Mehmet460 wrote:

Dear All,

I want to remove low number of reads from alignment BAM file. How can I do that?

rna-seq alignment gene genome • 1.2k views
ADD COMMENTlink modified 3.7 years ago by Pierre Lindenbaum116k • written 3.7 years ago by Mehmet460

Please elaborate what you mean by "low number of reads". Are these reads with low mapping quality ?

ADD REPLYlink written 3.7 years ago by Ashutosh Pandey11k

Hi, I am looking my alignment file on Artemis. I have seen low mapping quality (low number of reads per an  exon). I don't know how to explain this. For example, There is only one read per an exon, which affects my gene statistics obtained by using RNAseq hints. I used Augustus for gene prediction.

ADD REPLYlink written 3.7 years ago by Mehmet460

Hi again, I want to remove low coverage maps from BAM file. How can I do that?

ADD REPLYlink written 3.7 years ago by Mehmet460
2
gravatar for Pierre Lindenbaum
3.7 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum116k wrote:

* use GATK depthofcoverage to get the coverage of the bam

* get a BED file of the high-cov regions

* use samtools view to extract the reads overlaping the BED file.

ADD COMMENTlink written 3.7 years ago by Pierre Lindenbaum116k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 845 users visited in the last hour