Entering edit mode
8 months 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
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.