Question: Splitting pair end file
14 months ago by
jramo033 • 20
jramo033 • 20 wrote:
This is the type of file for pair end DNA sequence I have:
@MONK:275:C1Y6RACXX:1:1101:1169:1943/1 NAATAAATACTCTAAGACAAATCTATCACTTTAAAAGAGAGAAATCATATTCTAGGTTTTCATTAATGTCAGACTTAGCAGTGAAGATTCTCAGGTTCCT + #4=DFFFFHHHHHJJJJJJIIIJJJJJIJJJJIJJJJIJGJGIJJJJIIJJJJIJJHIJJJJJJJIIJIGIIIJJHIIFIJHHHHHHFFFFFFFEEEEED @MONK:275:C1Y6RACXX:1:1101:1169:1943/2 CAAATAAATACTCTAAGACAAAACTTGAATATAGAAAGTTGTAGAGGGTTTTTAGAAATTGATTTATAATACTTACAAATGATGAGTCTGAAAGAAGCCAAG + CB@@DDFFFHHHHHJIIIJJJJJJJJJFHJGIHHHGIHGHIHIHHIHJ?GHDGHIIIJIJJIIJIJJIIJJIJIJIJIIHFHHGEFCEFFFFEEDEEDDCDA
I need help to split the file in right and reverse reading. The number of reading is not the same (100 vs 102) so I was thinking in using Manipulate FASTQ. Can you help me how to write the code so that at the end I can have two separate files?
ADD COMMENT • link •modified 14 months ago by Jennifer Hillman Jackson ♦ 25k • written 14 months ago by jramo033 • 20