Gsnap: 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 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 ArticleMacs Raises Error: No Such File Or Directory
I am trying to run macs14 with sam files from paired-end data + control. Macs14 returns "No such file":sb7904313:line2 $ macs14 -t /Volumes/Data/G6L2_G6L3/s5/clean_paired_sample.sam -c...
View ArticlePaire-End Alignments And Scripture
Hello, I am slowly getting my hands on Scripture and have stumbled on yet another unclear (at least to me )detail. I have paired-end data and the two pairs are in a sigle TopHat-derived *.bam file. I...
View ArticleJoining Paired-End Illumina Raw Reads
I recently have amplicons sequenced (Illumina PE250) to investigate microbial community. I want to know what quality detection, and what kind of trims of raw reads by what programs should be performed...
View ArticleBWA MEM mate pair rescue
Greetings,Can someone spell out what this option means?  I have several guesses, but would rather ask than guess. -P     In the paired-end mode, perform SW to rescue missing hits only but do not try...
View ArticleForum: Mapping Of Ngs Short Reads
This is a simple explanation of how the mapping of short reads works !http://www.youtube.com/watch?v=1ZyoI-4ObSA&feature=related see the first 16 min ! It helped me a lot to understand the basic...
View ArticlePaired-End Bam Files
Hi,Having two BAM files from NGS data, how can one check if they are the BAM files (left and right) from a paired end mapping of the same sample? Thanks for the help.
View ArticleIn Paired-End Data, If One Read Of The Pair Is Unmapped, Is That Pair...
Hello everyone, I have a simple question about generic paired-end Illumina data. If one read of a pair is unmapped, does this automatically mean that the pair is improper and is (un)flagged in a SAM...
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 ArticleDoes Every Read In Paired-End Sam File Have The 0X0001 Flag?
Hello everyone, I have a simple, generic, question about SAM format flags in paired-end Illumina data. Does every read in a SAM file from a paired-end sequencing run automatically have the 0x0001 flag?...
View ArticleBest Cnv Software?
Hello,I know there are many reviews out there, but I can't seem to find exactly what I like.I found this software called mrCaNaVaR, and it's great since it has it's own aligner which is supposed to be...
View ArticleWhat Does Requirebothendsmapped From Rsubread Package Means?
Hi,I am using the featureCounts from the Rsubread package. And I am trying to understand what does the requireBothEndsMapped option do. The manual says"logical indicating if both ends from the same...
View ArticleCounting Reads On Paired-End Strand-Specific Rnaseq Data
Hello everybody, I've a strand-specific paired-end library on which I'd like to perform some standard DGE analysis. However I'm quite unclear about how to go about counting reads. Usually when I have a...
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 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 ArticlePicard Matequery Slows Process To A Crawl
I'm looking to iterate through an indexed BAM file using picard and perform various tests on both a read and it's mate. For some I would need the full SAMRecord for the mate so I can't just use the...
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 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 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 Article