Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,145,313 C→A 100% S683Y (TCC→TAC)  yegE → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,145,3130CA100.0% 207.2 / NA 58S683Y (TCC→TAC) yegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (26/32);  total (26/32)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAA  >  NC_000913/2145175‑2145448
                                                                                                                                          |                                                                                                                                       
ggAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCa                                                                                                                               <  1:14423/149‑1 (MQ=255)
 gAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCAc                                                                                                                              <  4:633619/149‑1 (MQ=255)
  aGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACt                                                                                                                             <  4:203188/149‑1 (MQ=255)
      gcggcAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACg                                                                                                                         >  3:476870/1‑149 (MQ=255)
        ggcAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCa                                                                                                                       >  1:309565/1‑149 (MQ=255)
         gcAGTTACGACGTGCATTACAGTATTACGACGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCAt                                                                                                                      >  2:15611/1‑149 (MQ=255)
              tACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGc                                                                                                                 <  3:65111/149‑1 (MQ=255)
               acgacgTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGcc                                                                                                                <  3:675413/149‑1 (MQ=255)
                         ttACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCa                                                                                                      >  1:439109/1‑149 (MQ=255)
                             aGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGttt                                                                                                  >  2:206917/1‑149 (MQ=255)
                              gTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGtttt                                                                                                 <  1:694247/149‑1 (MQ=255)
                               tATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTg                                                                                                >  2:224886/1‑149 (MQ=255)
                                aTTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTga                                                                                               <  4:58895/149‑1 (MQ=255)
                                aTTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTga                                                                                               <  3:354320/149‑1 (MQ=255)
                                    cGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGaaa                                                                                           <  2:440515/149‑1 (MQ=255)
                                     gCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAaac                                                                                          >  4:247303/1‑149 (MQ=255)
                                        gTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAaacaac                                                                                       >  3:597124/1‑149 (MQ=255)
                                           aaGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGc                                                                                    <  1:285725/149‑1 (MQ=255)
                                              tACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTa                                                                                 >  4:672587/1‑149 (MQ=255)
                                               aCTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTAt                                                                                <  3:15170/149‑1 (MQ=255)
                                                 tctGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATcc                                                                              >  3:57171/1‑149 (MQ=255)
                                                 tctGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATcc                                                                              >  2:517506/1‑149 (MQ=255)
                                                     gACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCctgct                                                                          <  2:628951/149‑1 (MQ=255)
                                                             caATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACgg                                                                  >  1:27172/1‑149 (MQ=255)
                                                              aaTATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGt                                                                 >  2:373961/1‑149 (MQ=255)
                                                                    ggttcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAg                                                           <  1:657565/149‑1 (MQ=255)
                                                                    ggttcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGAACTGACGCAGCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAg                                                           <  2:409667/149‑1 (MQ=255)
                                                                      ttcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTa                                                         >  3:292940/1‑149 (MQ=255)
                                                                      ttcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTa                                                         <  2:398816/149‑1 (MQ=255)
                                                                      ttcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTa                                                         >  4:414702/1‑149 (MQ=255)
                                                                      ttcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTa                                                         >  1:625737/1‑149 (MQ=255)
                                                                      ttcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCCCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTa                                                         <  2:280244/149‑1 (MQ=255)
                                                                                 tGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGa                                              >  3:374814/1‑149 (MQ=255)
                                                                                    ttCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACAt                                           <  1:224711/149‑1 (MQ=255)
                                                                                            gTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGgt                                   >  3:87653/1‑149 (MQ=255)
                                                                                               aCCGAATCACGCAAAATGCTGCGCAAGCTGATCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGtt                                >  1:612515/1‑149 (MQ=255)
                                                                                                 cGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTa                              <  3:32452/149‑1 (MQ=255)
                                                                                                      cACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCGACCAAGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGAt                         <  4:445069/149‑1 (MQ=255)
                                                                                                       aCGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATc                        >  3:398424/1‑149 (MQ=255)
                                                                                                       aCGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATc                        <  4:292940/149‑1 (MQ=255)
                                                                                                       aCGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATc                        <  2:27172/149‑1 (MQ=255)
                                                                                                        cGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCt                       >  2:98083/1‑149 (MQ=255)
                                                                                                            aaaTGCTGCGCCAGCTGAGCTACAGCGCCTACCCTGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCtt                      >  3:141900/1‑146 (MQ=255)
                                                                                                            aaaTGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCtt                      <  4:141900/146‑1 (MQ=255)
                                                                                                                 ctgcGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCtt              <  2:124374/149‑1 (MQ=255)
                                                                                                                  tgcGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCttt             <  1:247716/149‑1 (MQ=255)
                                                                                                                  tgcGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCttt             <  3:10825/149‑1 (MQ=255)
                                                                                                                   gcgcCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGAt                         >  3:226830/1‑136 (MQ=255)
                                                                                                                   gcgcCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGAt                         <  4:226830/136‑1 (MQ=255)
                                                                                                                     gcCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTaaa          >  1:529755/1‑149 (MQ=255)
                                                                                                                     gcCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGATTTaaa          <  4:289662/149‑1 (MQ=255)
                                                                                                                       cAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGc        <  2:625737/149‑1 (MQ=255)
                                                                                                                        aGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCg       <  4:312796/149‑1 (MQ=255)
                                                                                                                          cTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGt     <  1:224886/149‑1 (MQ=255)
                                                                                                                            gAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCGGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGa   >  1:646030/1‑149 (MQ=255)
                                                                                                                             aGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGaa  <  4:57171/149‑1 (MQ=255)
                                                                                                                                tACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGt                                                                 <  1:308975/83‑1 (MQ=255)
                                                                                                                                tACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGt                                                                 >  2:308975/1‑83 (MQ=255)
                                                                                                                                          |                                                                                                                                       
GGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAA  >  NC_000913/2145175‑2145448

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