Is 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 ArticleBest way to map paired-end, uniquely mapped reads with Bowtie2?
This has come up before, but I have never been satisfied with, or completly understood, the answers. My issue is how can Bowtie2 be used to map uniquely mapping paired-end reads (Illumina) to a genome...
View ArticleSequence Variation Detection For Prokaryotic Genomes
I have resequenced a bacterial genome using SOLiD paired end reads. After looking at SNPs and small Indels, I'm now interested in large Insertions/Deletions (Size of genes and gene clusters).Which...
View ArticleTophat 2 - Both Pairs Map Concordantly
Hi,I have just run the following command tophat2 --solexa1.3-quals -p 12 -r 80 --max-multihits 1 --no-mixed --no-discordant /home/Turkey/Index/turkeyindex /home/Turkey/WTCHG24920061sequence1.fa...
View ArticleHow To Read Maq Paired-End Alignemnt Data Using Bioconductor Packages
Hi all I have maq paired-end alignment files that I want to read into R. I have tried to browse several packages and they all seem to depend on ShortRead package of bioconductor which does not...
View Articlejoin non-overlapping paired-end reads
I have a batch of illumina paired end reads with a shorter than anticipated reverse sequence for a ~550bp amplicon.The joining tools which I have used before (i.e. USEARCH, ea-utils) talk about merging...
View ArticleHow To Count Stand-Specific Paired-End Rna-Seq Reads Overlapping Known...
Dear BiostarsDoes any one how to overlap stand-specific paired-end RNA-Seq reads (BAM) with known protein coding genes (BED) ?I tried the following but I think it is not the correct way ? Would...
View ArticleUsing Paired End And Orphaned Singles For De Novo Assembly
I have been using FastX to process reads prior to de novo assembly and mapping. What I have discovered and few have pointed out is the FastX will delete reads leaving reads unpaired which changes the...
View ArticleHow To Determine If Library Is Strand-Specific
Hello all, this might be a very basic question but I gave it some thought and don't see a satisfactory answer.Let's say I have a FASTQ file from a sequencing experiment. How can one quickly determine...
View ArticleExisting Tools For Post-Processing After Aligning Paired-End Reads With Blat?
Hi,I'm wondering if there are existing tools that can do post-processing on paired-end reads' BLAT output.More specifically, I'd like the tool to "merge" the alignments from the alignments of each...
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 ArticleWhat Better Way To Get The Paired Reads Aligned Against The Reference Genome?
Hi everybody,I have a group of paired reads sequenced using Solid 4 (50bp each mate). I discovered that reads are contaminated by E.coli. My strategy is to align the reads against the reference genome...
View ArticleShould We Dump Illumina Pair-End Mapping Results In Sam With Mapq=0, But Good...
hi, everyone! I am working on illumina pair-end sequencingAfter mapping by bwa, I got a pair of reads with MAPQ=0, with both reads mapped to more than one place. But the Template Length is OK, and I...
View ArticleCoverage 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 ArticleCollect Read Pairs Where At Least One Read Is Mapped
I might word my initial question like another post, but I really have the opposite meaning, i think:Filtering multiple flags with SAMtoolsI am trying to remove paired-end reads from a .SAM file where...
View ArticleIs There An Elegant Way To Extract Only The Properly-Paired Reads In A...
I know I should be filtering for the following tags: 99,163,83,147 and I know that samtools would work to get all the pairs. For example:samtools view -F 0x99 -b in.bam I was wondering if there was a...
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 ArticleBwa Sampe Segmentation Fault
Hi everyone, I'm running bwa in the sampe mode and, after successfully processing >10M reads it fails with a segmentation fault (as follows) on what appears like a set of poorly-alignable reads. Any...
View ArticlePaired End Vs Single End Detection
Hi, I am wondering if anyone has a subroutine or function (hopefully in Perl or pseudocode) to detect whether a fastq file is a shuffled paired-end or not. I think I've seen a few different syntaxes on...
View ArticleHow To Determine If Paired–End Illumina Rnaseq Reads Are Strand–Specific
I've been provided with more than a billion reads of RNAseq data for a poorly annotated nematode species. They appear to be 2x100 paired-end Illumina reads – I currently know frustratingly little about...
View Article