AssertionError: not all umis are the same length(!): 4 - 5
0
0
Entering edit mode
3.8 years ago

Hello,

I am trying to run UMI-tools to remove duplicated reads based on UMI. But getting error "AssertionError: not all umis are the same length(!): 4 - 5". Can anyone suggest how to resolve this error?

I am using following command :

umi_tools dedup --stdin=filtered.bam --umi-separator=":" --log=LOGFILE  --output-stats=stats.txt  -S output.bam > OUTFILE

This is the head of my bam file:

A00609:116:H7JCGDSXY:1:2171:13937:29825 16      1       10534   255     3S96M2D21M31S   *       0       0       CGCAGTACCACCGAAATCTGTGCAGAGGAGAACGCAGCTCCGCCCTCGCGGTGCTCTCCGGGTCTGTGCTGAGGAGAACGCAACTCCGCCGTCGCAAAGGCGCCGCGCCGGCGCAGACGCCCCCATGTACTCTGCGTTGATACCACTGCTT FFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFF,FFFFFFFFF,FFFFFFFFFFFFFFFFFF:FFFF:FF:FFFFFFF:FFFFFFF:FFFFF:FFFF  NH:i:1  HI:i:1  AS:i:103        nM:i:3  RE:A:I  xf:i:0  li:i:0  CR:Z:CCGTGAGAGAACGCGT   CY:Z::FF:F:FF:FFF,F:F   CB:Z:CCGTGAGAGAACGCGT-1 UR:Z:TTTTCCGCACTT       UY:ZG:Z:cells:0:1:H7JCGDSXY:1

I am getting following error:

Traceback (most recent call last):
  File "/home/ubuntu/miniconda3/bin/umi_tools", line 11, in <module>
    sys.exit(main())
  File "/home/ubuntu/miniconda3/lib/python3.8/site-packages/umi_tools/umi_tools.py", line 61, in main
    module.main(sys.argv)
  File "/home/ubuntu/miniconda3/lib/python3.8/site-packages/umi_tools/dedup.py", line 329, in main
    reads, umis, umi_counts = processor(
  File "/home/ubuntu/miniconda3/lib/python3.8/site-packages/umi_tools/network.py", line 419, in __call__
    clusters = self.UMIClusterer(counts, threshold)
  File "/home/ubuntu/miniconda3/lib/python3.8/site-packages/umi_tools/network.py", line 367, in __call__
    assert max(len_umis) == min(len_umis), (
AssertionError: not all umis are the same length(!):  4 - 5
umitools RNA-Seq umi • 1.9k views
ADD COMMENT
0
Entering edit mode

Posted at GitHub and answered by the developers:

ADD REPLY

Login before adding your answer.

Traffic: 3112 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6