Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	477,435	T→C	T395A (ACT→GCT)	ylaB ←	putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	477,435	0	T	C	100.0%	215.9 / NA	61	T395A (ACT→GCT)	ylaB	putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/33);  total (28/33)

CGACATCCAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATATGCTGC  >  NC_000913/477300‑477553                                                                                                                                        |                                                                                                                       cGACATCCAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAg                                                                                                                     >  1:652538/1‑139 (MQ=255)  gACATCCAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGt                                                                                                                    >  1:13380/1‑139 (MQ=255)        cAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAgc                                                                                                              >  1:73401/1‑139 (MQ=255)         aGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTGTCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAgcg                                                                                                             >  1:994206/1‑139 (MQ=255)          gATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAgcgc                                                                                                            >  1:234614/1‑139 (MQ=255)             ccTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                         <  2:918291/139‑1 (MQ=255)             ccTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                         <  2:652538/139‑1 (MQ=255)             ccTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                         <  2:642782/139‑1 (MQ=255)             ccTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                         <  2:997452/139‑1 (MQ=255)             ccTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                         <  1:1003598/139‑1 (MQ=255)                 tAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAtctgt                                                                                                     >  1:121455/1‑139 (MQ=255)                      aaCTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTggg                                                                                                >  1:891578/1‑139 (MQ=255)                        cTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                         >  1:315076/1‑128 (MQ=255)                        cTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                         <  2:315076/128‑1 (MQ=255)                                tGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGa                                                                                      <  2:73401/139‑1 (MQ=255)                                 gAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGcgcg                                                                                                           <  2:608820/117‑1 (MQ=255)                                 gAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGcgcg                                                                                                           >  1:608820/1‑117 (MQ=255)                                 gAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGTGCGATCTGTCTGGGATTAACCTGAt                                                                                     >  2:520926/1‑139 (MQ=255)                                  aaTACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                         >  2:916194/1‑118 (MQ=255)                                  aaTACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                         <  1:916194/118‑1 (MQ=255)                                  aaTACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATa                                                                                    >  2:490910/1‑139 (MQ=255)                                          gTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGtgattga                                                                            <  2:13380/139‑1 (MQ=255)                                           tACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGAt                                                                                              <  2:224330/120‑1 (MQ=255)                                           tACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGAt                                                                                              >  1:224330/1‑120 (MQ=255)                                           tACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGtgattgat                                                                           >  2:1007407/1‑139 (MQ=255)                                            aCCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATc                                                                          >  2:144320/1‑139 (MQ=255)                                                                     tGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                         <  2:111255/83‑1 (MQ=255)                                                                     tGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                         >  1:111255/1‑83 (MQ=255)                                                                          cGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTcc                                            <  1:490910/139‑1 (MQ=255)                                                                                ccGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTg                                      <  2:234614/139‑1 (MQ=255)                                                                                       cGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                         >  1:374978/1‑65 (MQ=255)                                                                                       cGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                         <  2:374978/65‑1 (MQ=255)                                                                                        gagaAATTGTCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACggg                              <  2:994206/139‑1 (MQ=255)                                                                                        gagaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACggg                              <  2:978705/139‑1 (MQ=255)                                                                                        gagaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACggg                              <  2:697183/139‑1 (MQ=255)                                                                                         agaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                         <  2:814643/63‑1 (MQ=255)                                                                                         agaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                         <  2:184589/63‑1 (MQ=255)                                                                                         agaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                         >  1:814643/1‑63 (MQ=255)                                                                                         agaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                         >  1:184589/1‑63 (MQ=255)                                                                                         agaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACgggg                             >  1:480226/1‑139 (MQ=255)                                                                                                aTCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCaa                      >  1:380321/1‑139 (MQ=255)                                                                                                 tCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAg                     <  2:380321/139‑1 (MQ=255)                                                                                                  cGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGTGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGa                    <  1:520926/139‑1 (MQ=255)                                                                                                  cGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGt                                                                      >  1:1011300/1‑89 (MQ=255)                                                                                                  cGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGt                                                                      <  2:1011300/89‑1 (MQ=255)                                                                                                  cGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGa                    <  2:6211/139‑1 (MQ=255)                                                                                                      gcgcTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACgg                               >  1:118495/1‑124 (MQ=255)                                                                                                      gcgcTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACgg                               <  2:118495/124‑1 (MQ=255)                                                                                                      gcgcTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGa                <  2:194381/139‑1 (MQ=255)                                                                                                      gcgcTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGa                <  2:121455/139‑1 (MQ=255)                                                                                                         cTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                         <  1:1065348/47‑1 (MQ=255)                                                                                                         cTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                         >  2:1065348/1‑47 (MQ=255)                                                                                                         cTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGAt                                                                                              <  1:597967/58‑1 (MQ=255)                                                                                                         cTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGAt                                                                                              >  2:597967/1‑58 (MQ=255)                                                                                                                    aTCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATAtgctgc  <  1:816997/139‑1 (MQ=255)                                                                                                                      cGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGa                                                                                      >  2:753303/1‑53 (MQ=255)                                                                                                                      cGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGa                                                                                      <  1:753303/53‑1 (MQ=255)                                                                                                                       ggCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAAc                                                                                          >  2:708924/1‑48 (MQ=255)                                                                                                                       ggCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAAc                                                                                          <  1:708924/48‑1 (MQ=255)                                                                                                                                   ttCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATAtgct    >  1:221571/1‑122 (MQ=255)                                                                                                                                   ttCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATAtgct    <  2:221571/122‑1 (MQ=255)                                                                                                                                        |                                                                                                                       CGACATCCAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATATGCTGC  >  NC_000913/477300‑477553

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑