Difference Between "Mate Pair" And "Pair-End"
Just as the title , I can't tell the difference between those two conception. :) waiting for your help.
View ArticlePaired-End Protocol For Micrornaseq
In another post, a guy wanted to know how to analyze paired-end data and use them to predict microRNAs.I never heard about a paired-end protocol for miRNAseq and would be interested in some more...
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 ArticleSamtools Mpileup And Overlapping Paired-End Reads
I am using samtools mpileup to generate a pileup for paired-end RNA sequencing data. I am curious about how samtools handles pair mates whose read mappings overlap. With regard to the simple example...
View ArticlePaired-End Reads Alignment For Variant Calling ?
I'm trying to do variant calling (SNPs, Indels) from exome-sequencing data, and the sequencing was done with paired end reads. I would like to use BWA for mapping/alignment, followed by PiCard and GATK...
View ArticleFastq Sort By Sequence
Dear all,I have some questions about manipulating fastq files. Its the first time I do this, so I want advice, how to do it. Its a special case, because I only want to keep reads when the following...
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 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 ArticleTake A Subset Of A Fastq Paired-End Sample
Hi,I have two paired-end fastq compressed files coming from HiSeq RNA-SEq experiment, ie., pair.1.fastq.gz and pair.2.fastq.gz.The files are very large so I wanted to just take a few million/thousand...
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 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 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 ArticleWgsim Mutations In Output After Setting Everything To 0
I was just wondering, is there any useful information on wgsim? Tutorial? Anything? I have been stuck with it for the last 2 weeks. I'm really not sure how to use it. I need it for a project of mine....
View ArticleMapping Trimmed Paired End Reads With Mapsplice
Hi! I'm trying to map illumina paired end using MapSplice. I used script named sickle to trim the bad quality tails and the output looks like this:Paired end 1@TUPAC_0006:1:1:3062:1473#0/1...
View ArticleHow To Assemble Genome Generated With Bac Clones?
I have 2 fastq files from illimina with reads length 250b. Sequences from one file obtained by sequencing from "right" and from "left" in another. This is paired end sequencing. As it is whole genome...
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 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 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 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 ArticleDownload Large Paired-End Rna-Seq And Microarray Data Of The Same Sample (>...
Hi,I was trying to find if there is any large size paired-end RNA-Seq and microarray data of the same sample, but I wasn't able to find as much data as I wanted.Could somebody please point out where I...
View Article