SEQanswers

Go Back   SEQanswers > Search Forums


Showing results 1 to 5 of 5
Search took 0.00 seconds.
Search: Posts Made By: GenoMax
Forum: General 09-17-2019, 03:03 PM
Replies: 14
Views: 1,303
Posted By GenoMax
Then we have been on wrong track all along. You...

Then we have been on wrong track all along. You should have said you need a consensus fasta file from the BAM alignment.

What you need is: https://www.biostars.org/p/367960/

I had posted this...
Forum: General 09-17-2019, 12:00 PM
Replies: 14
Views: 1,303
Posted By GenoMax
Output fasta reads were interleaved (put in one...

Output fasta reads were interleaved (put in one file) since you had paired-end reads to begin with. To separate them into two files you should use the command below.


for i in *.fa; do...
Forum: General 09-12-2019, 08:23 AM
Replies: 14
Views: 1,303
Posted By GenoMax
Good catch SNPsaurus. I have updated my original...

Good catch SNPsaurus. I have updated my original command. There was an additional ".bam" in there.

@Manuelly: Use for i in *.bam; do reformat.sh in=$i out=${i}.fa; done

If you are trying to...
Forum: General 08-28-2019, 08:30 AM
Replies: 14
Views: 1,303
Posted By GenoMax
Yes you can convert multiple files at once. If...

Yes you can convert multiple files at once. If your data came from paired-end reads then be aware that the "out=" files will contain data in interleaved format. You will need to separate R1/R2 reads....
Forum: General 08-27-2019, 03:52 PM
Replies: 14
Views: 1,303
Posted By GenoMax
Use reformat.sh from BBMap suite. Take a look at...

Use reformat.sh from BBMap suite. Take a look at BAM processing options in in-line help to decide if you want to keep primary reads etc.

for i in *.bam; do reformat.sh in=${i} out=${i}.fa...
Showing results 1 to 5 of 5

 


All times are GMT -8. The time now is 05:44 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO