This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to splice the genome into many 700bp fragments ramdonly

To test the performance of my model ,I have splice the genome into many 700bp fragments randomly,which

could be used to construct the negtive and positive examples.I have used the R splice the genome

into many 700bp fragments orderly rather than randomly!!!

r genome sequence next-gen
• 2,498 views • link
updated by Brian Bushnell Brian Bushnell • written by jinxyhucy521 • jinxyhucy521

2 answers

Check out bedtools random command http://bedtools.readthedocs.io/en/latest/content/tools/random.html

• 1 views • link
written by microfuge microfuge

Thanks for your reply very much!!This is a very useful tool!!!

• 1 views • link
written by jinxyhucy521 • jinxyhucy521

Using the BBMap package:

randomreads.sh in=genome.fasta out=fragments.fasta length=700 reads=1000000 adderrors=f

That will give randomly-distributed sequences... alternatively:

shred.sh in=genome.fasta out=fragments.fasta length=700

That will produce evenly-spaced sequences; you can adjust the spacing with the "overlap" flag.

• 1 views • link
written by Brian Bushnell Brian Bushnell

2_chr1_1_1481410_1482109_1473410_NC_000853.1 Thermotoga maritima MSB8 chromosome, complete genome Dear Brian Bushnell,I have used your tools,and the result is the upper,I set the length=700,I have read your readme file,but I can't find the meaning of chr1_1_1481410_1482109_1473410...Could you tell me this meaning?

• 1 views • link
written by jinxyhucy521 • jinxyhucy521

"1473410" is the start position, and the part after (NC_000853.1 Thermotoga maritima MSB8 chromosome, complete genome) is the name of the sequence it came from. The parts before that are internal BBMap coordinates and not very useful. So, that sequence was 700 bp from "NC_000853.1 Thermotoga maritima MSB8 chromosome, complete genome" starting at position 1473410 (zero-based). Meaning, the first base in the sequence would be numbered 0.

• 1 views • link
written by Brian Bushnell Brian Bushnell

I think if the 1473410 is the starting,which one is the ending?And I have another question:does the

result has the fragments from complementary strand?

• 1 views • link
written by jinxyhucy521 • jinxyhucy521

Log in to answer this question.

More posts like this