SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   Demultiplexing (http://seqanswers.com/forums/showthread.php?t=62812)

pingu 09-19-2015 01:45 AM

Demultiplexing
 
Hello,
I have to demultiplex my fastq file, I have seen that about 200 reads are not assigned (total is 17000), there is a tool that can check barcode on 5' and 3'?
Thank you very much

GenoMax 09-19-2015 08:30 AM

Are you referring to an internal barcode that is present in the read itself? Do you know what those barcodes are?

pingu 09-19-2015 09:25 AM

Thank you for your help, I know the barcodes and they are at the of the reads (last 10bp)
For example, this read is not assigned:
@JBMBIZC01CKLPI
CACTGTAGAAGACTCGGCAGCATCTCCATGAGCCAGTATTGAAAATGTTGAAGATCAAAAAA
CACTAAGTTTTTCCAAAGTTAATATCCAATGTAAAAGATAGCAAATGCATACCCACAAACTGTA
AATGAAGATATTTGCGTTGAGGAACTTGTGACTAGCTCTTCACCCTGCAAAAATAAAAATGCA
GCCATTAAATTGTCCATATCTAATAGTAATAATTTTGAGGTAGGGCCACCTGCATTTAGGATAG
CCAGTGGTAAAATCGTTTGTGTTTCACATGAAACATTAAAAAAGTGAAAGACATATTTACAGA
CAGTTTCAGTAAAGTAATTAAGGAAAACAACGAGAATAAATCAAAAATTTGCCAAACGAAAA
TTATGGCAGGTTGTTTGGAGAACAGTGACGATCGCCTACAGTGCT

The barcodes of the reads is: AGCACTGTAG and it is a complementary reverse: CTACAGTGCT,
I have MID at the start of the reads (5') e at the end (3'), the barcode splitter tool checks MID at 5' and I have allowed 2 mismatches.

sklages 09-21-2015 02:44 AM

Looks like 454 data. If you have the original SFF file you could use Roche's sfffile tool to demultiplex your data. It recognizes 5' and 3' barcodes/MIDs.


All times are GMT -8. The time now is 07:52 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.