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_000913477,4350TC100.0% 137.7 / NA 40T395A (ACT→GCT) ylaBputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/20);  total (20/20)

ATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCG  >  NC_000913/477310‑477542
                                                                                                                             |                                                                                                            
aTCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGcgcg                                                                                                 >  1:651570/1‑139 (MQ=255)
  ccTGTAAATAACTTAAACTTGAATACCCCTTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                               <  2:497351/139‑1 (MQ=255)
     gTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAtctg                                                                                            >  2:129797/1‑139 (MQ=255)
     gTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAtctg                                                                                            >  2:607389/1‑139 (MQ=255)
      tAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAtctgt                                                                                           >  2:653094/1‑139 (MQ=255)
      tAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAtctgt                                                                                           >  2:315163/1‑139 (MQ=255)
      tAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAtctgt                                                                                           >  2:533368/1‑139 (MQ=255)
      tAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAtctgt                                                                                           >  1:698279/1‑139 (MQ=255)
      tAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAtctgt                                                                                           >  1:358897/1‑139 (MQ=255)
       aaaTAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAtctgtc                                                                                          <  2:651570/139‑1 (MQ=255)
       aaaTAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAtctgtc                                                                                          <  1:533368/139‑1 (MQ=255)
                   cTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCt                                                                             <  2:362796/139‑1 (MQ=255)
                       aaTACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATa                                                                         <  2:453139/139‑1 (MQ=255)
                        aTACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGCTCAAGCGCGATCTGTCTGGGATTAACCTGATAg                                                                        >  1:185378/1‑139 (MQ=255)
                            cccGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGtgat                                                                    >  1:432324/1‑139 (MQ=255)
                                tACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTGAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGa                                                                           <  1:515689/128‑1 (MQ=255)
                                tACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGa                                                                           >  2:515689/1‑128 (MQ=255)
                                tACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGtgattgat                                                                >  2:410687/1‑139 (MQ=255)
                                 aCCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                               <  1:588577/108‑1 (MQ=255)
                                 aCCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                               >  2:588577/1‑108 (MQ=255)
                                     aacgcatcaAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTGCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGt                                                           >  2:496142/5‑139 (MQ=255)
                                           tcaAGATAAGTTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCATCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGa                                                                           <  2:164473/117‑1 (MQ=255)
                                           tcaAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGa                                                                           >  1:164473/1‑117 (MQ=255)
                                            caAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                               <  1:444933/97‑1 (MQ=255)
                                            caAGATAAATTTCATGGCCCGCCTCCCGGAAGCGAGAAATTATCGAGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                               >  2:444933/1‑97 (MQ=255)
                                                   aaTTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGc                                             <  2:432324/139‑1 (MQ=255)
                                                          tGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGAttt                                      >  1:270676/1‑139 (MQ=255)
                                                                    cccGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGt                            <  1:374162/139‑1 (MQ=255)
                                                                           gCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGCTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACg                     <  2:185378/139‑1 (MQ=255)
                                                                             gagaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTTCGAGGTGAGCACggg                   <  1:653094/139‑1 (MQ=255)
                                                                             gagaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACggg                   <  1:129797/139‑1 (MQ=255)
                                                                                        tggtctcTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGGTCAAGCGGGAGCTTCAGAGGTGTAACCTGGTAGTGATTGATTATGTCACGCAGCAACTGCAGGATTTTTACCGAGGTGAGCACGGGCAATTCGAGAt        <  2:28189/133‑1 (MQ=255)
                                                                                           gcgcTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGtt                                                                                                         >  2:103681/1‑40 (MQ=255)
                                                                                           gcgcTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGtt                                                                                                         <  1:103681/40‑1 (MQ=255)
                                                                                           gcgcTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTgg                                                                                      >  1:656346/1‑58 (MQ=255)
                                                                                           gcgcTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTgg                                                                                      <  2:656346/58‑1 (MQ=255)
                                                                                            cgcTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGAt    <  2:259736/139‑1 (MQ=255)
                                                                                              cTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCg  <  1:438042/139‑1 (MQ=255)
                                                                                                           cGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGAtt                                                                                  >  1:685234/1‑46 (MQ=255)
                                                                                                           cGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGAtt                                                                                  <  2:685234/46‑1 (MQ=255)
                                                                                                                             |                                                                                                            
ATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCG  >  NC_000913/477310‑477542

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