Question: Splitting pair end file
0
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