This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to take out sequences with barcode?

Hello,

https://www.ncbi.nlm.nih.gov/sra/SRX2791702[accn]

In the above link, the design section has barcode sequences, how do i get all the reads with each particular barcode?

I did try using grep '^<barcode sequence="">' from the fastq file. But as you can see the last column in the link is "no. of sequences" and when i try to count the number by using grep, I'm getting a different number. The number I get is not matching with what they have provided.

Am i using grep incorrectly? what is the position of these barcode sequences?

Thanks!

sequencing fastq illumina barcode
• 3,054 views • link
updated by Ido Tamir Ido Tamir • written by suvratha • suvratha

use GBS tools such as GBSX for extracting reads with defined bar codes. suvratha

• 0 views • link
written by cpad0112 cpad0112

1 answer

You could have been more precise with the difference between your read numbers and the stated numbers. If the stated one is bigger in all samples, then its because often and by default demultiplexing is done with 1 mismatch, which grep can not do.

• 0 views • link
written by Ido Tamir Ido Tamir

grep gives more than the number mentioned. for e.g - one the mentioned numbers there is about 6.1k and grep gives me 13.5k.

• 0 views • link
written by suvratha • suvratha

Log in to answer this question.

More posts like this