Skip to content
Snippets Groups Projects
Commit fc3f65ef authored by Blaise Li's avatar Blaise Li
Browse files

Reverting, siRNA can start with any nucleotide.

parent 5bc3c6ce
No related branches found
No related tags found
No related merge requests found
...@@ -355,8 +355,13 @@ def count_sis(annot_infos, write_queue): ...@@ -355,8 +355,13 @@ def count_sis(annot_infos, write_queue):
# Eliminate most common non-siRNA target cases based on annotation name. # Eliminate most common non-siRNA target cases based on annotation name.
# Also eliminate the reads with no signature # Also eliminate the reads with no signature
# (the annot_name should be "no_signature") # (the annot_name should be "no_signature")
if (annot_name in COMMON_NO_TARGET # if (annot_name in COMMON_NO_TARGET
or annot_name[:3] == "no_"): # or annot_name[:3] == "no_"):
# continue
# Actually, we do accept reads that do not start with G
# Provided they are not plain "no_signature"
# (that is: without other annotations)
if annot_name in COMMON_NO_TARGET:
continue continue
# annotations: set of tuples # annotations: set of tuples
# (attrs["gene_biotype"], gtf.strand, gtf.start, gtf.end, attrs["gene_id"]) # (attrs["gene_biotype"], gtf.strand, gtf.start, gtf.end, attrs["gene_id"])
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment