Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • northbio
    Junior Member
    • Sep 2010
    • 4

    One error in SSPACE

    I've encounter a strange problem when using SSPACE v3.0. I found that reads dealed by SSPACE were as follows in one fasta file in standard_output/reads :

    >read24410302
    GCCTTATCGCCCAACCTCAGTGATAAACGTTACTAATGCTCCTGTGGCTGAATGGGAGCAAATCCCTGCAGCCAGGTTACAAAATCTTGTTGAAAGCCTTCCAAGAAGAGCGGACGCTGTTATAGCTGTGTATTAATGCC
    >read24410302
    GAGTATGTCCTCAAACCAGAGCTTTCAAACTGTGAGGCGTGATGAGGTAAAGATGTTGTGTGAGTTGAAAGGGCCAGTTATTGTAGTATGGCCCCACTGATTGACACAGTCAGCAGTAGTAGCAGCAGTAGTGGCACACA
    >read24410303
    CAGGACAGAATCTGAACAGAAATGAGTCCAGAAGGAGATGGCTGTCAAAGTGTTAGTATTTTCTGACAAGATTAGTGCTGATATGGTTAACTGATTAAATCAATACAATGAAAAATAATCAATTCCAAATATGATGGTCA
    >read24410303
    TGCAGACTTGGAAAGCTCATCCCAGGGAAGGAGTATGTTTTCCGTATCCGTGCTGAGAACCGCTTTGGCGTGTCAGACCCCATTCAGTCTGACAGAATGGTTGCCAAGTTCCCCTTTGGTAAGTTTAATTGAATTATGTC
    >read24410304
    >ST-E00291:98:HVL7GCCXX:4:1101:5731:38878/1
    >read24410304
    GGGCTTTAAATCAAAAATAACCTGTATCCATTTCTCTGGGTCACGTACTGTAAATAATGCATTTTTTGATTGCATTCCTGCAACCTGGAAAATGTTATTTGTAAAAGGAAATTCTTGAACAAAATACACTACCTCAGATT
    >read24410305
    >ST-E00291:98:HVL7GCCXX:4:1101:5629:38878/1
    >read24410305
    AAAGCTGCAGGCGACTGTGCGTCTGCAACATGCAGCACATACACATGCATCAAAGTGTACACGCTAATACACCT
    >read24410306
    >ST-E00291:98:HVL7GCCXX:4:1101:28260:11488/1
    >read24410306
    TTGCCTTTTCTTCCAGGCTTCTAAATTTAAAAAACACAGGACAAAATGTGCTAAATAGGGAGCACACGGTACCAAGAGAGGACCCTGCCTCAGAAGTTAATCTGCCGAGGGCTTTAGAACAGAGCAACACACTGCAGAGC
    >read24410307
    >ST-E00291:98:HVL7GCCXX:4:1101:6380:38878/1
    >read24410307
    ACTGAGGACAGTATCTGGTATCAAGTTAATCAGTAAGACCATTTTACTCGCTGCACTCTCATTAG

    As you can see, suddenly the fasta file gave wrong format of fasta file. At this time, I checked the reads file I gave in libraries.txt. I found the origin of this error:

    >ST-E00291:98:HVL7GCCXX:4:1101:5609:38878/1
    GCCTTATCGCCCAACCTCAGTGATAAACGTTACTAATGCTCCTGTGGCTGAATGGGAGCAAATCCCTGCAGCCAGGTTACAAAATCTTGTTGAAAGCCTTCCAAGAAGAGCGGACGCTGTTATAGCTGTGTATTAATGCC
    >ST-E00291:98:HVL7GCCXX:4:1101:27245:11488/1
    CAGGACAGAATCTGAACAGAAATGAGTCCAGAAGGAGATGGCTGTCAAAGTGTTAGTATTTTCTGACAAGATTAGTGCTGATATGGTTAACTGATTAAATCAATACAATGAAAAATAATCAATTCCAAATATGATGGTCA
    >ST-E00291:98:HVL7GCCXX:2:1101AGTCCCTCTGACCCTCATGTGTGGGTTTTCTAACTGTACTGTGAATAAG
    >ST-E00291:98:HVL7GCCXX:4:1101:5731:38878/1
    ATGAAACACCCAGAGTACAATACTGGTTAATACGATACTGAATACTGAGTATCACAACACAAACTGAACTTCAAACACAGAATTTGAAGCCCACCTGTGACCTGTTACCTGAATCACACAGATACATAGTACCATCAAGA
    >ST-E00291:98:HVL7GCCXX:4:1101:5629:38878/1
    GTACAGACACAGCGTGGGACAAAGCAACGTTAAACCGAGTTCAATTAAAGAAAATGTCAGTCAGTTATGGGCCCAAAGCTTCAAATGTTGCATTTATCTTAAAACGATGAACGACTGGAGATCTTAAGATGAACTTTCGC
    >ST-E00291:98:HVL7GCCXX:4:1101:28260:11488/1
    TATTATTATTACTATTATTTTATTAATATTATTACCATTACCATTACTATTATTTTAAATCTGTATGTGAAATTTGCATTGTGCTACTGTATGCCCCCCCCTCAACCCAACTGGTCGAGGCAGATGGCTGCCCACATTGA

    The red sequence is the point.......
    But the real real point I want to say is this fasta is produced by bowtie2, which I used --un-conc to output these unmapped reads.
    So, Be careful to use it ......

Latest Articles

Collapse

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, Today, 10:09 AM
0 responses
8 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, Yesterday, 08:59 AM
0 responses
14 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 06-02-2026, 12:03 PM
0 responses
22 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 06-02-2026, 11:40 AM
0 responses
19 views
0 reactions
Last Post SEQadmin2  
Working...