Identify point mutations from each read in bam files and change to N
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
To extract uniques SNPs froma. multisample VCF file
written by sam •Hi, I have a multi-sample VCF file. This VCF contains around 20 samples (Sample1 - Sample20). I intend to extract unique SNPs of Sample1 compared …
-
What is the common practice to view WGS data on IGV?
written by chictu •Hi, I would like to manually examine mutations on specific chromosomes in around 50 whole-genome sequencing (WGS) data samples. Each BAM file is larger than …
-
Cap base quality in bam file
written by junnna •Hello, I would like to change the base quality scores in PacBio HiFi bam files, such that the base qualities are capped to 33. Is …
-
Is it possible to remove reads from BAM which match the reference at a position
written by yryan •I'm interested in looking for modified bases and for the downstream analysis I'm doing I'd like to remove the reads which match the reference base …
-
How to introduce artificial mutation in bam
written by bassanio •Hi, Is there a way I can introduce mutation in bam file at specific position of a chromosome without running the alignment. For example I …
-
Edit bam files using Pysam
written by zaidnab •Hi there! I am trying to edit values in a bam file, using Pysam, or any other interface. I want to change the base quality, …
-
Base call from BAM file
written by Simo •Hi, I have a BAM file and a reference (let's say 1000 genomes data), I would like to extract from the BAM file all the …
-
bam-readcount skips genomic positions contained in the bedfile
written by jml96 •Hi, I am using bam-readcount to calculate the read counts for each base in a list of genomic positions. However, in most of the samples …
-
find similar sequences between different sources
written by sam •Hello, I have different sets of sequences from different sources (e.g. I have around 20 fasta files (each fasta file correspond to one source) where …
-
Bam File: Change Chromosome Notation
written by Lisanne •<p>Dear All,</p> <p>I am using <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> files for chip-seq analysis. The chromosome notation in a usual <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> file is like: chr1. In my …
why ? what are you trying to do ?
I have short-read data for a cultivar and want to generate assembly for the genic region based on a reference genome. Apart from the uniquely mapped reads, I am also using filtered multi-mapped reads (to improve coverage) with a certain number of mutations in these reads. However, I want to substitute the mismatched base in the muti-mapped read against reference to 'N'.