Hi everyone,
I have aligned bam files with the reference genome. All reads are single reads.
I want to analyze fragment size using bamPEFragmentSize.
My questions are as follows:
How does it work for single reads?. Because while running the program I saw a message says No pairs were found. Is the data from a paired-end sequencing experiment?
- I find that bamPEFragmentSize read the sorted bam in my case. Do I need to pass the unsorted bam file instead of the sorted bam file? regards adr
1 answer
It is not a tool question. Think of it that way, you are trying to hammer a nail into the wall. You have 10 different hammers, but no nail. You cannot nail something without nails. Same here, you want insert size stats and have no insert sizes because it is single-end. At best you can ask the people who made the library for the details on the fragment size assessed by bioanalyzer or gel electrophoresis. You do not have the data to analyse what you want to analyze.
Log in to answer this question.
-
Does the bigwig file of a single-end chip-seq shows 2 peaks for a fragment whereas the paired-end s…
written by Rajendra KC •
-
Warning: Mate records missing HTSEQ
written by san96 •
-
unsorting a bam file
written by wiscoyogi •
-
how to do basic statistics for bam files
written by adR
-
Demultiplexing before merging paired-end reads?
written by asidhu •
-
CNV detection using BAM file
written by reza
-
samtools - Finding Average Fragment size for pair?
written by oars
-
Difference in htseq-count values for sorted bam and unsorted bam by name
written by nalandaatmi •
-
finding copy number variation from bam file.
written by Varun Gupta
-
Why is the output of Artificial FASTQ Generator not aligning properly with Bowtie 2?
written by John Smith
