Dear all,
I have RNASeq data. I have to estimate insert size for multiple samples. For this, I am using the CollectInsertSizeMetrics picard application. Here I have to provide input of BAM file (reads mapped to reference). I am confused whether I should do preprocessing of data and then estimate insert size or I should directly move with raw data to estimate insert size. please look into following points/conditions....
- consider all the reads for mapping and then calculate insert size ?
- perform quality filtration,adapter removal, mapping and then estimate the insert size ?
- only perform the adapter removal and then estimate the insert size ?
Every step will make the difference. I just want to follow the accurate way to estimate insert size.
Waiting for reply
Thank you in advance
•
written
by
Omics data mining
2 answers
If you're going to filter the data anyway then filter it before performing the estimation. You care about the insert size of the data you'll actually use. Do track how much data is filtered, though, since that will also indicate whether there was an issue during library prep.
BTW, you might need to align against the transcriptome to get an accurate number.
If all you are interested in is finding insert size then use BBMap suite like this. There are two ways you can estimate the insert size.
Log in to answer this question.
-
Illumina reads preprocessing best practice for snp calling applications
written by Denis
-
inner-distance plot created by inner_distance.py is always right-skewed
written by Ati •
-
TCGA array expression and rnaseq expression batch correction
written by Omics data mining
-
filters need to apply for Exome bam files after mapping
written by Omics data mining
-
Markduplicates in RNASEQ
written by Omics data mining
-
Flexibility of inner distance parameter for Tophat
written by acorella •
-
the defination of insert size in RNA-seq, including the length of read, or not?
written by zju.whw •
-
Splice-aware insert size for RNA-seq
written by shuelga •
-
long insert size estimation
written by Damian Kao
-
Insert Size For Each Bam Entry
written by Nicolas Rosewick
You might want to take a look at this blog: https://mikelove.wordpress.com/2016/09/26/rna-seq-fragment-sequence-bias/ and at the
alpinebioconductor package.