How to build a reference genome given an excel file of sequences
Hello,
I have sequencing data from a library prep that I've removed the adapter sequences from, leaving a set of ~20 nt barcodes generated from the reads (about 2e6). I have a file listing all of the possible barcodes (~5000, they are also ~20 nt) and would like to know the number of times each barcode appears among the reads.
I know this is probably a more powerful tool than I need, but I was thinking of using Bowtie 2 to bin the reads, using a reference genome where each barcode is its own chromosome. However, I'm not sure how to convert the table of all the barcodes into the correct format to feed into Bowtie 2.
Any help would be much appreciated. Thanks!
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Unexpected read structure in Ribo-seq libraries prepared with the QIAseq miRNA kit: missing 5′ adap…
written by GiovannaPrezia •
-
Demultiplexing asymmetrically barcoded reads from multiplex amplicon nanopore sequencing
written by michael.flower.14 •
-
Split large fastq based on value in description line
written by mmarcell •
-
Picard demultiplexing with random barcodes
written by RNG_Daemon •
-
Can I use FLEXBAR to detect and remove barcodes out of my short reads without demultiplexing ?
written by giulia.trauzzi •
-
Trimming Illumina Seq with TRIMMOMATIC
written by giulia.trauzzi •
-
Unrecognized Sequence in NGS reads
written by jer364 •
-
Count pooled sequencing barcode in fastq
written by vinaykusuma •
-
Untrimmed reads in small-RNA seq datasets
written by russell.stewart.j •
-
Processing Barcoded Hiseq Data
written by bbio
If your reads and barcodes are both 20 bp and you are looking for exact matches then simply using UNIX
grepusing full word match mode (or ideallyseqkit grep(LINK) may be the easiest option.