Best 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 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 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 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 ArticleResampling 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 ArticleOrientation in paired-end sequencing?
I am new to bioinformatics and currently learning how to use Bowtie 2. As written in the manual: A pair that aligns with the expected relative mate orientation and with the expected range of distances...
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 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 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 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 ArticleHow To Rearrange Paired End Bam File?
Hello all,I have a paired end bam file and I want to use bedtools for them. After merging, the paired end read alignments are not lying next to each other. It is making problems in the bedtools...
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 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 ArticleSamtools Count Paired-End Reads
Hi, I used tophat to align paired-end reads from an rna-seq experiment and I obtained an accepted_hits.bam alignment file. Using the accepted_hits.bam I'd like to count the number of properly aligned...
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 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 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 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 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 ArticleNeed help using Shrimp2 on paired end color-space SOLiD data.
Hi, I have SOLiD reads which are paried-end (75bp and 35bp) in .csfasta and .QV.qual format. I would like to use Shrimp2 to align them. So far I have been having trouble using it. I used the following...
View Article