Pile up reads without reference?
Hi everyone,
I have multiple groups of reads (count < 1000) which derived from a number of novel junctions, with each reads contain minimum ~10bp overhang on the other side (i.e., exon). Because the number of random links between exon-exon can be quite large (hundreds), it is difficult to construct artificial reference for every case.
I would like to ask, is the any simple way to pile up reads (fasta or fastq) and visualize single consensus sequence?. I thought of assembling them to make contigs, then map them using contigs as reference, and visualize with IGV, but there should be a simpler way to just pile them up and visualize it?
Thank you in advance for your suggestion!
1 answer
Log in to answer this question.
More posts like this
-
De novo assembly for phylogenetic analysis
written by isag •
-
Mapping a small number of contigs to a reference subsequence
written by alexandru.bologa.marian •
-
bowtie2: Print only Alignment Summary?
written by jer364 •
-
Phylogenetic tree for very large numbers of bacterial assemblies
written by yairgatt •
-
Mapping softwares setting minimum matching read length
written by sumaluzn •
-
Identify reads that span intron-exon junctions in RNA-Seq
written by oriolebaltimore
-
14 contigs recieved from SMRTportal using HGAP.3
written by komwit.s •
-
Assemble fully overlapping reads into contigs that may have several clones
written by Iayork •
-
Scaffolding with SSPACE returned scaffold the same as initial contigs set
written by pbigbig
-
Selecting and de-selecting exon-exon junctions based on depth or coverage from RNA-seq mapping
written by mjoyraj
How was the alignment done (has there been any alignment done?)?
Thank you, I used kmer-based searching for target sequence.
Did you extract reads containing certain k-mers from a FASTQ file? If you don't have any shared coordinate space among the different reads then you either need to go via assembly or alignment to line and pile up similar reads
Yes, I did extracted reads using certain k-mer patterns.