Hello,
I have been using dwgsim to stimulate the reads for my masters thesis project.
reads generated by dwgsim are as follows:
@gi|29366675|ref|NC_000866.4|_148598_1_0_1_0_0_7:0:0_0:0:0_0/1
TCAATGTTTAAAATTTTTCTATAAGCCTTTAACTTTATAGAATTATTATTCCAGACTAAATTTTCAGTCTGTTCATCATGTTTATCAATGATATTTAAAATCGAATCAAGCAAGATAAACGTCTCAAACGAAATTATGTTCGATTGCAGAAGTTTAAAAATATAACTTGATTGAACTTTTGAATTATACTCAAAAATTTCTTTAAAAGCAGAAACTTCAACTTTTTTACTAAAATAATAAATGTTGCGAATATCTTCTTCAAACTTAAATTTAATTTGCTCTAAGCGTCCGATATATTCA
+
274243352322242112304101.421322033226312022224227224543221700163322322254222233202222323410323211235002421417.1212/32032144113/34243333142240265131237222310434142423502121331306222223125411125322651421222221301325123622224320232443510122.240104223//211227532272513317.654312223233312046223331322/3520
@gi|29366675|ref|NC_000866.4|_76814_1_1_0_0_0_6:0:0_0:0:0_1/1
TGGACAAGAATTTTTTATTGAAATAAAACCTAAAAAAGAAACACAACCACCAGTTAAACCAGCACATCTAACGACCGCAGCGAAGAAAAGATTTATGAATGAAATTTATACCTGGTCTGTGAACACTGACAAATGGAAAGCAGCACAATCTTTAGCTGAAAAGCGTGGAATAAAATTTAGAAGTCTAACAGAAGATGGATTAGGAGCTCTTGGCTTTAAGGGGGCATAATGGCTATTTTTTAAATAATTAATGAAAGCACTCCCCAATTTCCAAAGGTTAAGCAATCATTAAATGATAAG
+
22522022353212222523315421112234/62214021243134223,43/1223252232230232/3.223202222224233232313252311325252022243342030246412026523234522/22132/22/15341151224513230152142411053242222234-22431351303212121442/22222144322524220/251221301223332102013114322254132.252140334236327221301340224342321213223224
@gi|29366675|ref|NC_000866.4|_107901_1_0_1_0_0_5:0:0_0:0:0_2/1
TGCTCCTAAATTCTTGGAAGATGTGCTTGCAACTGAAATGGCAGATGAAATCAATAAAGACATTCTGCAGTCTATGATTACAGTGTCAAAACGCTATAAAGTTACAGGAATTACTGATAGTGGATTCATCGATTTGAGTTATGCATCTGCTCCTGAAGCTGGTCGTTCATTATACCGAATGGTATGTGAAATGGTTTCGAATATCCAAAAAGAATCAACTTATACAGCAACGTTCTGTGGTGCTTCAGCTCGTGCCGCTGCGATTCTTGCTGCATCAGGCTGGTTAAAACATAAACCAGA
+
000122222226252023235322.522002744/4354861323223512353233012456213104435423/2061742/2303222333236.531/6/22113.2322841222215253302242/04024543251222220232343020061613.32023222362232225221145475/032212202230/44451252162-0423101/3104/64221342214.15423122522032240254130331452422223422/022232323313226243
@gi|29366675|ref|NC_000866.4|_98258_1_1_0_0_0_2:1:0_0:0:0_3/1
TGTCCGGAGATAATAAAGTCATTTTTAATCCTCTTTAATATGCTTTAAAATATTTATACCATTGACATACCATGAGATACTGGAACATACTCAGCAGAATGAACCGAATCCACAAATATAACTGGCGCGTAGTCGTCGCTCATATCCTGAAGCTCTTTTGAAAATACTTCAGATGCTAATCGCATGTCATCTTTATCCGCATAATCAATAAATTTTGACTGCGTTGATAACCATCCAAAAATCACTAAAGACATTCCTAAATCGTCATGATAACCTTCTTCAGCCGCCCAAGACACGCCT
+
2211223513131223225222123224233522355451323223322312222200121/122332322306232222122122420532154611333313662/202536242223136710122332362322223533020225222221/3/42424432025020224312234203210.41425200660222222261/222523444326015421214235422/430203321424132223//52-2222.3321431123251222/25312263422323023
As per the documentation of dwgsim given on sourceforge: in reads name header after 2nd underscore there should be start end 2 (zero-based) but in my case its always 1. I dont know what is that happening. Have anybody faced problem like this? Any help is greatly appreciated!!
I have been using dwgsim to stimulate the reads for my masters thesis project.
reads generated by dwgsim are as follows:
@gi|29366675|ref|NC_000866.4|_148598_1_0_1_0_0_7:0:0_0:0:0_0/1
TCAATGTTTAAAATTTTTCTATAAGCCTTTAACTTTATAGAATTATTATTCCAGACTAAATTTTCAGTCTGTTCATCATGTTTATCAATGATATTTAAAATCGAATCAAGCAAGATAAACGTCTCAAACGAAATTATGTTCGATTGCAGAAGTTTAAAAATATAACTTGATTGAACTTTTGAATTATACTCAAAAATTTCTTTAAAAGCAGAAACTTCAACTTTTTTACTAAAATAATAAATGTTGCGAATATCTTCTTCAAACTTAAATTTAATTTGCTCTAAGCGTCCGATATATTCA
+
274243352322242112304101.421322033226312022224227224543221700163322322254222233202222323410323211235002421417.1212/32032144113/34243333142240265131237222310434142423502121331306222223125411125322651421222221301325123622224320232443510122.240104223//211227532272513317.654312223233312046223331322/3520
@gi|29366675|ref|NC_000866.4|_76814_1_1_0_0_0_6:0:0_0:0:0_1/1
TGGACAAGAATTTTTTATTGAAATAAAACCTAAAAAAGAAACACAACCACCAGTTAAACCAGCACATCTAACGACCGCAGCGAAGAAAAGATTTATGAATGAAATTTATACCTGGTCTGTGAACACTGACAAATGGAAAGCAGCACAATCTTTAGCTGAAAAGCGTGGAATAAAATTTAGAAGTCTAACAGAAGATGGATTAGGAGCTCTTGGCTTTAAGGGGGCATAATGGCTATTTTTTAAATAATTAATGAAAGCACTCCCCAATTTCCAAAGGTTAAGCAATCATTAAATGATAAG
+
22522022353212222523315421112234/62214021243134223,43/1223252232230232/3.223202222224233232313252311325252022243342030246412026523234522/22132/22/15341151224513230152142411053242222234-22431351303212121442/22222144322524220/251221301223332102013114322254132.252140334236327221301340224342321213223224
@gi|29366675|ref|NC_000866.4|_107901_1_0_1_0_0_5:0:0_0:0:0_2/1
TGCTCCTAAATTCTTGGAAGATGTGCTTGCAACTGAAATGGCAGATGAAATCAATAAAGACATTCTGCAGTCTATGATTACAGTGTCAAAACGCTATAAAGTTACAGGAATTACTGATAGTGGATTCATCGATTTGAGTTATGCATCTGCTCCTGAAGCTGGTCGTTCATTATACCGAATGGTATGTGAAATGGTTTCGAATATCCAAAAAGAATCAACTTATACAGCAACGTTCTGTGGTGCTTCAGCTCGTGCCGCTGCGATTCTTGCTGCATCAGGCTGGTTAAAACATAAACCAGA
+
000122222226252023235322.522002744/4354861323223512353233012456213104435423/2061742/2303222333236.531/6/22113.2322841222215253302242/04024543251222220232343020061613.32023222362232225221145475/032212202230/44451252162-0423101/3104/64221342214.15423122522032240254130331452422223422/022232323313226243
@gi|29366675|ref|NC_000866.4|_98258_1_1_0_0_0_2:1:0_0:0:0_3/1
TGTCCGGAGATAATAAAGTCATTTTTAATCCTCTTTAATATGCTTTAAAATATTTATACCATTGACATACCATGAGATACTGGAACATACTCAGCAGAATGAACCGAATCCACAAATATAACTGGCGCGTAGTCGTCGCTCATATCCTGAAGCTCTTTTGAAAATACTTCAGATGCTAATCGCATGTCATCTTTATCCGCATAATCAATAAATTTTGACTGCGTTGATAACCATCCAAAAATCACTAAAGACATTCCTAAATCGTCATGATAACCTTCTTCAGCCGCCCAAGACACGCCT
+
2211223513131223225222123224233522355451323223322312222200121/122332322306232222122122420532154611333313662/202536242223136710122332362322223533020225222221/3/42424432025020224312234203210.41425200660222222261/222523444326015421214235422/430203321424132223//52-2222.3321431123251222/25312263422323023
As per the documentation of dwgsim given on sourceforge: in reads name header after 2nd underscore there should be start end 2 (zero-based) but in my case its always 1. I dont know what is that happening. Have anybody faced problem like this? Any help is greatly appreciated!!
Comment