Adrian,
You can try running repair.sh to split the file into paired and unpaired reads, and then map twice, once for the paired and once for the unpaired, and then merge the bam files. That will allow maximal use of the available information.
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Yes thats a disadvantage I agree.
Unfortunately, the bam file does not have enough PE reads.
When I used bamtofastq for PE fastq files, interestingly I obtained 0 fastq reads.
Leave a comment:
-
Yes, you'll be expected to decrease your mapping efficiency a bit, since one mate can act as an anchor to rescue the other. Further, it's much easier to use paired-end reads to find isoforms, since you're then not relying solely on alignments over a splice junction.
Leave a comment:
-
converting paired-end (PE) bam file to single-end (SE) fastq
Hi:
while working with COAD TCGA BAM files, I find the very annoying to find PE reads. These files are mashed up and not consistent.
for example:
1. read lengths are not consistent. Some are 34 some 76 reads.
2. Many reads miss mate or pair.
I want to identify novel splicing differences however TCGA BAM files are mapped to known transcripts (known exon pairing from known isoforms gtf) thus limiting the discovery of novel isoforms.
I decided convert BAM to fastq and realign to full genome.
While doing this, because of loss of many pair and mates in bam, I converted them to single end fastq.
Any ideas if converting a paired-end bam to single end fastq pose any problem in philosophical ways.
thanksTags: None
Latest Articles
Collapse
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...-
Channel: Articles
04-22-2024, 07:01 AM -
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 08:47 AM
|
0 responses
14 views
0 likes
|
Last Post
by seqadmin
Yesterday, 08:47 AM
|
||
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
60 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
60 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
54 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
Leave a comment: