Hi Gopi,
If I understand you correctly, you are trying to extract the LCBs from your
alignments. I used the stripsubsetsLCB command to do this with the genomes I am
working with. The only issue I had was that if you ran your alignment in the
Windows GUI and not using Linux then the path names won’t recognize the input
files. Below is the link to the linux commands that I used. Hope this is
helpful.
http://sourceforge.net/p/mauve/mailman/mauve-users/thread/51f26b88.2000...@helsinki.fi/
-Susan
From: Gopi Nath [mailto:itsmegopina...@gmail.com]
Sent: Wednesday, February 03, 2016 1:58 PM
To: mauve-users@lists.sourceforge.net
Subject: [Mauve-users] How to extract the alignment sequence
Hi,
I have aligned 5 bacterial genomes using the progressiveMauve. I would like to
extract the part of the alignment having sequence in all the genome (eg: second
and fourth line positions) using the alignment(.xmfa) file and
backbone(.backbone) file.
seq0_leftend seq0_rightend seq1_leftend seq1_rightend seq2_leftend
seq2_rightend seq3_leftend seq3_rightend seq4_leftend seq4_rightend
0 0 23127 23463 20107 20443 32207 32543 29528 29864
-28145 -28315 23464 23633 20444 20613 32544 32713 29865 30034
0 0 23634 23814 20614 20794 32714 32894 30035 30215
-27824 -28003 23815 23997 20795 20977 32895 33077 30216 30398
Is there any utility to achieve this task ? How to extract only the lines
having sequence from all the genomes(eg: 2nd, 4th) ?
Thanks.
------------------------------------------------------------------------------
Site24x7 APM Insight: Get Deep Visibility into Application Performance
APM + Mobile APM + RUM: Monitor 3 App instances at just $35/Month
Monitor end-to-end web transactions and take corrective actions now
Troubleshoot faster and improve end-user experience. Signup Now!
http://pubads.g.doubleclick.net/gampad/clk?id=272487151&iu=/4140
_______________________________________________
Mauve-users mailing list
Mauve-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/mauve-users
