Resampling Fastq Sequences Without Replacement
Hello, I want to extract a random sample (without replacement) of 7.5 million fastq sequences from illumina sequencing data that contains about 30 million sequences each in of the reads. I want to...
View ArticleDealing With Read Counts Under Pe And Se Scenarios
Hi, I am unsure how to deal with this case to go about analysing RNA-seq data. Suppose that you have a control and treatment setup with 4 biological replicates each. However, two in control and two in...
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 ArticleClik here to view.
Trinity Question
Hi， I 'd like to use Trinity to analyze my strand-specific pair-end sequencing data (dUTP) and I am very confused about how to choose the argument. I got two reads files, _1.fa and _2.fa. For...
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 ArticleMerging Illumina Paired End Reads
Dear All,I have fastq a dataset containing forward and reverse sequences obtained through paired end module of Illumina platform. I am trying to merge these paired end reads. I have a query which I...
View ArticlePair End Sequencing Problem
Dear all,I have a quick question on pair end sequencing. I used to work with Illumina without the pair end reads and I have dificulties to understand how the pair end reads work.In the "old" system you...
View ArticleIdentifying mutations from Paired-End Sequencing data
Hello! I'm trying to get mutations from paired-end sequenced reads aligned with BWA Â using SamTools. Coverage is about 16,000. Generally it works fine, but in one fragment (TGGGC) i see that in reads...
View ArticleWarnings In Bowtie Mapping
Hello, I am trying to use bowtie on small synthetic data for short read mapping. My command is ./bowtie -p 8 -t -S hg19 -1 synthetic_sample1.fq -2 synthetic_sample2.fq > bowtie.sam. The alignment...
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 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 ArticleX And Y Chromsome Crossover Position Alignments
How do I find the genotypes on the X chromosome which match the Y SNPs listed in raw data from 23andMe, Ancestry or FTDNA
View ArticleTool Recommendation Wanted For Cleaning Fasta/Fastq Files To Remove Unpaired...
Hi Everyone, I've been digging around the web trying to find a tool that would allow me to clean-up my paired-end Illumina data before mapping. My pipeline thus far has been to:1) FASTQC - my R1 file...
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 ArticleBasic Paired-End Sam Questions
I can't seem to find answers to 2 very basic questions.For a paired-end sam, is there a separate sam line for mate1's and mate2's? If so, how do find the mate of a given sam line?Thanks
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 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 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 ArticleAmos With Paired End Data
Hi all, how do I take advantage of paired end data with AMOS? I have a 2x100bp fastq file that I would like to use AMOScmp-shortreads with. I see that when I google for it, there might be a .mates file...
View ArticleEstimating Insert Size From Paired End Data.
Hi,I have paired end data from illumina. To estimate the insert size in silico ( from scratch ), I have aligned the reads as single end reads to the genome ( mouse ). Now I have the two alignment files...
View Article