Extract Umi sequence to readname
0
0
Entering edit mode
5 weeks ago
Smilesky • 0
@BH04640:4:22CNJJLT3:5:1101:27120:1048 1:N:0:CAAGGTGA+AGTCGCGA
GNACTGTTATCTATTGAAGGACATGTGCATTACCTTATACAGGAAGCTACTGATGAAAACTTACTATGCCAGATGTATCTTGGTTGGACTCCATATATGTGAAATGAAATTATGTAAAAGAATATGTTAAAGTAGAGATCGGAAGAGCAC
+
I#IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII          

My sequence looks like this. I am trying to extract the umi sequence using FGBIO - FastqToBam tool. can anyone tell me How to Extract the UMI in ubam. Because when i run the FastqToBam the umis are not extracted

BH04640:4:22CNJJLT3:5:1101:27009:1048:CAAGGTGA+AGTCGCGA 4   *   0   0   *   *   0   0   CNAGTAGCTAGAGTGGTGGCACATGCCTGTAGTCCCAGCTACTCGGAAGGCTGAGGTGGGAGAATCACTTGAACTCAGGAGGCGGAGGTTGCAGTGAGCCGAGATTGCGCCACTGCACTCCAGCCTGGGAGACAGAGCAAGACTCTGCCT  I#IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII  RG:Z:A.

I also want to know the coressponding read structure for the sample given

Sequence FgBio UMI • 230 views
ADD COMMENT
0
Entering edit mode

Can you post your entire FastqToBam command?

Are you telling the program the structure of your UMI's as noted here: http://fulcrumgenomics.github.io/fgbio/tools/latest/FastqToBam.html with --read-structures option.

ADD REPLY

Login before adding your answer.

Traffic: 1484 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