RNASeqReadSimulator is a set of scripts generating simulated RNA-Seq reads. It provides users a simple tool to generate RNA-Seq reads for research purposes, and a framework to allow experienced users to expand functions. RNASeqReadSimulator offers the following features:
- It allows users to randomly assign expression levels of transcripts and generate simulated single-end or paired-end RNA-Seq reads.
- It is able to generate RNA-Seq reads that have a specified positional bias profile.
- It is able to simulate random read errors from sequencing platforms.
- The simulator consists of a few simple Python scripts. All scripts are command line driven, allowing users to invoke and design more functions.
Hey I tried using the tools But I faced some errors.
ReplyDeleteI would be very thankful if you good help me sort out these errors at the earliest.
I hereby provide you fresh outputs.
TAIR10_GFF3_genes.sorted.BED
Chr1 0 30427671 chromosome . .
Chr1 3630 3759 five_prime_UTR . +
Chr1 3630 3913 exon . +
Chr1 3630 5899 gene . +
Chr1 3630 5899 mRNA . +
Chr1 3759 3913 CDS . +
Chr1 3759 5630 protein . +
Chr1 3995 4276 CDS . +
Chr1 3995 4276 exon . +
Chr1 4485 4605 CDS . +
>genexplvprofile.py -e -4,4 out.bed
Mean and variance for lognormal distribution: -4.0,4.0
#ID Length Dir Exons Position GroupID NIsoformInGroup Explv
> genexplvprofile.py --statonly TAIR10_GFF3_genes.sorted.BED
#ID Length Dir Exons Position GroupID NIsoformInGroup
> gensimreads.py --stranded -n 100 -l 88 -o out.bed -f 12 --stranded TAIR10_GFF3_genes.bed
Read count: 100
Read length: 88
Output bed file: out.bed
Field: 12
Assigning weights...
Sampling...
Error: incorrect number of fields (should be 12)
Error: incorrect number of fields (should be 12)
Error: incorrect number of fields (should be 12)
Error: incorrect number of fields (should be 12)
Error: incorrect number of fields (should be 12)
Error: incorrect number of fields (should be 12)
.
.
.
Error: incorrect number of fields (should be 12)
Total 590264 lines...
Total 0 reads...
What exactly is the error related to field I have changed the value of fields n number of times form 7 to 24 including 12 but everytime I get the same error
And as for the generation of expression profile
all I get it single line..nothing else.
How exactly does this work.Could you clarify.
Regards,
Ganga Jeena
Can this tools also generate fastq reads.
DeleteIf Bed files do not have quality information fastq can not be generated from bed file by bedtools.
Regards,
Ganga Jeena
I'm 15 years old. I was born with HIV my mother passed away because of the HIV infection And I regret why i never met Dr Itua he could have cured my mum for me because as a single mother it was very hard for my mother I came across Dr itua healing words online about how he cure different disease in different races diseases like HIV/Aids Herpes Copd Diabetes Hepatitis even Cancer I was so excited but frighten at same time because I haven't come across such thing article online then I contacted Dr Itua on Mail drituaherbalcenter@gmail.com I also chat with him on what's app +2348149277967 he tells me how it works then I tell him I want to proceed I paid him so swiftly Colorado post office I receive my herbal medicine within 4/5 working days he gave me guild lines to follow and here am I living healthy again can imagine how god use men to manifest his works am I writing in all articles online to spread the god work of Dr Itua Herbal Medicine Yes can ask me anything about this on my twitter @ericamilli or text 7205992850.He's a Great Man.
ReplyDelete