Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,027,785 A→G 7.7% L78P (CTG→CCG)  yedQ ← putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,027,7850AG 7.7% 150.6 / 2.4 52L78P (CTG→CCG) yedQputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (25/23);  new base G (2/2);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAACTGGCAGAGTGCGTCGTCGGTTGAGTTCGATTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATA  >  NC_000913/2027638‑2027928
                                                                                                                                                   |                                                                                                                                               
ttAACTGGGAGTGTGCGTCGTGGGTTGAGTTCGATTTTGCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGTGGTGATAT‑GAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAATTAGATCAGt                                                                                                                                               <  2:134181/149‑1 (MQ=255)
ttAACTGGCAGAGTGCGTCGTCGGTTGAGTTCGATTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAAT‑CAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGt                                                                                                                                               <  2:137921/149‑1 (MQ=255)
       gcatattGCGTCGTCGGTTGAGTTCGATTTTCCAGGCGTGCTCTTCGCGATGCTGTTCTAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCTCGCATGCCATTACGTAAAAAGATCAGTTTGTc                                                                                                                                          <  5:105250/142‑1 (MQ=255)
               cgtcgtCGGTTGAGTTCGATTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTtata                                                                                                                                 <  6:115208/149‑1 (MQ=255)
                             ttCGATTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTAt                                                                                                                   >  2:4034/1‑149 (MQ=255)
                             ttCGATTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTAt                                                                                                                   >  6:78098/1‑149 (MQ=255)
                                             tGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAAc                                                                                                   <  2:184722/149‑1 (MQ=255)
                                                  gtcgCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGtttt                                                                                              >  7:208447/1‑149 (MQ=255)
                                                      cgATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTc                                                                                             >  6:98883/1‑146 (MQ=255)
                                                          gCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAgatga                                                                                      >  4:169004/1‑149 (MQ=255)
                                                              ttCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTAc                                                                                  >  5:174119/1‑149 (MQ=255)
                                                                   aCTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCGGTTTGTCGACGTTATACTGc                                                                                                                             >  1:205086/1‑101 (MQ=255)
                                                                   aCTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCGGTTTGTCGACGTTATACTGc                                                                                                                             <  2:205086/101‑1 (MQ=255)
                                                                    cTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACgctggc                                                                            <  7:64604/149‑1 (MQ=255)
                                                                              gCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGata                                                                    >  6:39706/1‑147 (MQ=255)
                                                                              gCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGAtata                                                                  <  5:78098/149‑1 (MQ=255)
                                                                               cATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCTAGATGATTACGCTGGCTGGAGATATAg                                                                 <  6:175220/149‑1 (MQ=255)
                                                                                   aCGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGACATTACGTAAAAAGATAAGTTGGTCGCCGTTATACTGCAAATGCTTATCGCGCGCGTTGGCACCGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCAt                                                             >  7:173535/1‑149 (MQ=255)
                                                                                     gCAGAGAGGTGAAATCCAGTGGCGCTAGGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCt                                                           >  4:220088/1‑149 (MQ=255)
                                                                                     gCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCt                                                           >  1:221839/1‑149 (MQ=255)
                                                                                        gagagGTGAATTCCATTAGCGCTACGAGAGCTTCGCTCATTCCATTACGTAATAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCTCGTTGGCAACGTTTTCCAGATGATTACGCTTTCTTGAGATATAGGCATCTTcc                                                        <  6:179136/149‑1 (MQ=255)
                                                                                              tGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGcac                                                   >  6:162085/1‑148 (MQ=255)
                                                                                                aaaTCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGcaccac                                                <  6:202733/149‑1 (MQ=255)
                                                                                                   tCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGcaccaccac                                             >  4:227346/1‑149 (MQ=255)
                                                                                                   tCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGcaccaccac                                             <  8:208447/149‑1 (MQ=255)
                                                                                                     cAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACtt                                           <  5:162085/149‑1 (MQ=255)
                                                                                                     cAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACtt                                           >  3:14536/1‑149 (MQ=255)
                                                                                                     cAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACt                                            >  2:71810/1‑148 (MQ=255)
                                                                                                       gTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCa                                         <  2:221839/149‑1 (MQ=255)
                                                                                                           cgcTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGcc                                     <  5:39706/149‑1 (MQ=255)
                                                                                                                cGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGt                                 >  6:236414/1‑148 (MQ=255)
                                                                                                                       ttCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAgg                         >  4:161259/1‑149 (MQ=255)
                                                                                                                          gcgcATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTg                       >  2:50568/1‑148 (MQ=255)
                                                                                                                            gcATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGaa                    <  3:161259/149‑1 (MQ=255)
                                                                                                                              aTGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGaaaa                  <  7:251559/149‑1 (MQ=255)
                                                                                                                                 ccATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGaaaaaa                <  5:236414/148‑1 (MQ=255)
                                                                                                                                 ccATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGaaaaaa                <  3:227346/148‑1 (MQ=255)
                                                                                                                                  cATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTcc                                                        >  7:45737/1‑107 (MQ=255)
                                                                                                                                  cATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTcc                                                        <  8:45737/107‑1 (MQ=255)
                                                                                                                                   aTTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCa             <  5:98883/149‑1 (MQ=255)
                                                                                                                                    ttACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGcta            <  6:174119/149‑3 (MQ=255)
                                                                                                                                         tAAAAAGATCGGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCa       >  8:227520/1‑149 (MQ=255)
                                                                                                                                          aaaaaGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAgg                                  >  6:209898/1‑121 (MQ=255)
                                                                                                                                          aaaaaGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAgg                                  <  5:209898/121‑1 (MQ=255)
                                                                                                                                          aaaaaGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATAacac      >  7:153624/1‑149 (MQ=255)
                                                                                                                                          aaaaaGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCAc      >  7:8181/1‑149 (MQ=255)
                                                                                                                                          aaaaaGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCAc      >  1:106961/1‑149 (MQ=255)
                                                                                                                                           aaaaGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCAc      <  8:8181/148‑1 (MQ=255)
                                                                                                                                             aaGATCGGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGAt   <  7:227520/149‑1 (MQ=255)
                                                                                                                                             aaGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGa    >  3:209488/1‑148 (MQ=255)
                                                                                                                                                aTCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATa  >  3:119461/1‑147 (MQ=255)
                                                                                                                                                aTCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGGGGAAAAAAGCAATACCACGATa  >  8:205045/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TTAACTGGCAGAGTGCGTCGTCGGTTGAGTTCGATTTTCCAGGCGTGCTCGTCGCGATGCTGTTCGAACTCGGTTACAGCATCACGCAGAGAGGTGAAATCCAGTGGCGCTACGAGAGCTTCGCGCATGCCATTACGTAAAAAGATCAGTTTGTCGACGTTATACTGCAAATGCTTATCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATA  >  NC_000913/2027638‑2027928

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