Query Regarding Generating Reads With Genome Coverage Values Of 1X, 2X Etc
I have a genome length of length approx 3 million basepairs. I need to generate 100 bp random reads so that the coverage is 1X.
Ideally how many reads should I output in the forward strand and how many in the reverse orientation. Should these numbers be equal.
1 answer
1x coverage of a 3,000,000 bp genome is 3,000,000bp. Since you want 100bp reads, you want to generate 3000000/100 = 30000 reads. You should simulate the true sequencing process by randomly choosing the strand for each read. You will end up with approximately equal reads.
I would suggest not writing your own simulator however and use a tool like this: https://github.com/jstjohn/SimSeq
Log in to answer this question.
More posts like this
-
Cut&Run Matrix Generation
written by Mohammed •
-
Reads mapped on wrong strand
written by qstefano •
-
Primmer trimming 4C seq
written by amie.e.siemonek •
-
Number of reads per lane
written by kamel •
-
what read should i use for classifying?
written by life99945 •
-
Paired end sequencing with different read length: better to trim everything to a short length or us…
written by Lalla •
-
How does IGB provide strand specific output
written by nagamalaimohan •
-
Presence of both forward and reverse motifs in transcription factor binding sites
written by jjl5sw •
-
RNA seq FASTX quality trimming
written by Rahul •
-
Determining Orientation Of Fastq Files
written by dan79 •
