Coverage For Pair-End Rna-Seq, Extend Reads Or Not?
HI, all I have a question when computing region coverage of pair-end RNASeq data. As showed by the sketch map, when computing region coverage, whether I should use actually mapped reads or extended...
View ArticleGsnap: How To Output Only Alignments With Minimal Number Of Mismatches?
I am aligning paired-end ChIP-Seq data of drosophila transcription factors.I would like to have a .sam file with the following alignments reported: - A maximum of three mismatches, i.e. three wrong...
View ArticleAnalyzing Older Illumina Paired-End Data
Hi all, I've got some reads from mid-2010 that I'd like to re-align. I'm not sure how best to proceed. These are from a Illumina Genome Analyzer IIx. They're 40-bp, paired-end. Here's what the text...
View ArticleTophat - Understated Number Of Reads In The "Align_Summary.Txt" File
Hi all. I'm working with paired-end rna-seq data to assemble transcriptome of my species of interest. I've just realized that Tophat is understating the number of reads that I actually have and...
View ArticleAlign Paired End Reads Using Blast
Hi all, Has anyone align illumina paired end reads using BLAST, I used gsnap to do the alignment first, then use BLAST to align the reads which were not mapped by gsnap. It seems that BLAST can only...
View ArticleMapping Applications On Hadoop Cluster
Hello,I would like to know if there are any new mapping applications for paired end data designed to work on hadoop cluster. The mapping applications that i know are CloudBurst (single ended), Crossbow...
View ArticleTool: Trim Adapters Of Paired-End Reads (Fastq)
Trimming adapter sequences of paired-end experiments is sometimes a problem. If you clip the mates in two steps, it migh happen that you loose one mate, but not the corresponding one, resulting in two...
View ArticleForum: Attempting to understand pair end sequencing, mainly illumina
Hi I am currently attempted to understand how pair end sequencing works. I understand the basics that it sequences from both ends and if there is an overlapping region they can be joined to create...
View ArticleEfficiently join paired-end read coordinates in the same line?
Dear community, I have a huge paired-end HiC dataset (BAM format) which I want to format like this way: HWI-D00283:117:C5KKJANXX:2:1101:1139:77789 chr6 153338506 153338556 37 + chr6 153338031 153338081...
View ArticleDna Sequencing Using Abyss-Pe Of Paired End Reads
Hello! I've got to run abyss-pe on 2 files I've got and find the best parameters (k and so on) that would create the best contig coverage (sorry if I'm messing things up, I'm a programmer that had a...
View ArticleIs There Any Advantage Of Paired End Sequencing For Chip-Seq?
Hi,I hope that some of you may have a recommendation regarding the use of paired-end sequencing over single-end for ChIP-seq. In principle we expect only a minor improvement of using paired-end...
View ArticleHigh level of duplicate in one reads of paired-end data
Hi,We are doing some transcriptomic analysis on bovine immune blood cells and we seem's to have some problem with high levels of duplicate in our data. Our library were prepared with the illumina...
View ArticleAligning Reads To Specific Chromosome Using Bwa
Hello Everyone,I have whole genome illumina paired end reads and I want to align my reads to specific chromosome (chr 21) using BWA. I was thinking of aligning the entire reads to fasta file of the...
View ArticleCounting sense reads in bacterial paired-end RNA-seq data
Hi, I'm trying to count reads mapping to sense strand. I have doubts which counts file I should chose from this pipeline. I think is "plate_R.counts" because has more reads counted in total. Am I...
View ArticleFilter Paired-End Sam File For Xt:A:U
Dear all,I have a sam file (BWA output, paired-end reads). I would like to retain only reads which are "properly paired". This I would do by:samtools view -f 0x002 file.sam >...
View ArticleHow To Extract Information From Fastq Pair End Files
dear BioStars users,I would like to extract from my pair-end fastq files information how many times my read is occurring in my fastq file. So output could look - my read (sequence) - how many times I...
View ArticleSorting Fastq Files After Trimming (Orphans And Pe)
I have a bunch of Illumina PE data that has been run through fastx trimmer and clipper. I am ready to map these reads, but am needing to create 2 files for paired end reads (the left and right hand...
View ArticleDoes The Mapping Of A Read'S Pair Affect Its Own Mapping Score In Bwa?
I would like to know more about how BWA treats the mapping score of paired end reads. I am familiar with the process used to assign mapping scores by BWA (thanks to this post). I see that the mapping...
View ArticleChIPseq paired-end bowtie2 concern regarding biological replicates
I am working with ChIP-seq paired-end data where there is concern that one or more of the biological replicates may not be very good, but it is unknown which replicate may have a problem (I suspect...
View ArticleHow To Interpret Crossbow Data
Hello,May I know how to interpret Crossbow output. Is there any chance of building a SAM or BAM file from Crossbow output. I am looking into genome mapping. I know there are many tools available for...
View Article