how to exclude sequence from sam file?
1
0
Entering edit mode
9.6 years ago
Medhat 9.7k

The issue that I have Illumina mate pairs sequence, I filtered them using skewer it results in two fastq files forward and reverse.

There are some sequences that was removed from one of the files "did not pass the filtering process"

when I use bowtie2 it gives this error:

Error, fewer reads in file specified with -1 than in file specified with -2
terminate called after throwing an instance of 'int'
bowtie2-align died with signal 6 (ABRT)

Is there any script or software to exclude those sequences from the two files or I need to write my own?

sequencing mapping illumina sam mate-pairs • 2.4k views
ADD COMMENT
2
Entering edit mode
9.6 years ago

There's a convenient tool in BBtools for this. Otherwise, see How To Sort Two Mate Pair (Fastq) Files So That The Order Of The Identifiers Is The Same?.

ADD COMMENT
1
Entering edit mode

I tried BBtools worked fine and have a very good and helpful options thanks a lot

ADD REPLY

Login before adding your answer.

Traffic: 1771 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6