This will only serve a small niche of people doing RNA-seq, but figured I'd post it here anyway.
I wrote a script that essentially implements the counting described here:
If you have an n-mer of random nucleotides at the end of your RNA-seq adapter, you can use that n-mer to count the number of molecules of RNA. Basically, if you have.
This script takes a FASTQ file and collapses it based on the barcode. For example, if you have the reads
GAGAGAGAGGGGGGGGGG
GAGAGAGAGGGGGGGGGG
CACACACAGGGGGGGGGG
and the barcode is the first 8 nucleotides, it gets collapsed down to
GGGGGGGGGG
GGGGGGGGGG
since it is likely that the first two barcodes arose from PCRing the same molecule.
Requires: Python 2.7, numpy, and biopython.
Best,
David
-
The recent pandemic caused worldwide health, economic, and social disruptions with its reverberations still felt today. A key takeaway from this event is the need for accurate and accessible tools for detecting and tracking infectious diseases. Timely identification is essential for early intervention, managing outbreaks, and preventing their spread. This article reviews several valuable tools employed in the detection and surveillance of infectious diseases.
...11-27-2023, 01:15 PM -
Microbiome research has led to the discovery of important connections to human and environmental health. Sequencing has become a core investigational tool in microbiome research, a subject that we covered during a recent webinar. Our expert speakers shared a number of advancements including improved experimental workflows, research involving transmission dynamics, and invaluable analysis resources. This article recaps their informative presentations, offering insights...11-09-2023, 07:02 AM
| Topics | Statistics | Last Post | ||
|---|---|---|---|---|
|
0 responses
11 views
0 likes
|
by seqadmin
Today, 08:26 AM
|
|||
|
0 responses
10 views
0 likes
|
by seqadmin
Today, 08:12 AM
|
|||
|
0 responses
19 views
0 likes
|
by seqadmin
11-27-2023, 08:12 AM
|
|||
|
1 response
60 views
0 likes
|
by VilliamPast
11-25-2023, 02:47 AM
|