Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,297,192 G→C 7.7% R57R (CGC→CGG)  TpbB ← diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,297,1920GC 7.7% 160.4 / 4.3 52R57R (CGC→CGG) TpbBdiguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (22/26);  new base C (2/2);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

GCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCG  >  NC_002947/5297045‑5297335
                                                                                                                                                    |                                                                                                                                               
gcgTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCAC‑TTCTTCGGTACGGGCTATCATCGCTAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCAGAGCCGCTTCTACGGTGTAGTTGATCGAg                                                                                                                                                 <  1:94151/148‑1 (MQ=255)
  gTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCgg                                                                                                                                              >  6:79199/1‑149 (MQ=255)
   ttGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATTGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCggg                                                                                                                                             <  3:315299/149‑1 (MQ=255)
     gCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACATTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGc                                                                                                                                            >  1:79710/1‑149 (MQ=255)
               cAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGCCCGCTCGATCAACTACa                                                                                                                                 <  3:256191/149‑1 (MQ=255)
                aGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAg                                                                                                                                <  5:143355/149‑1 (MQ=255)
                aGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAg                                                                                                                                <  6:22747/149‑1 (MQ=255)
                    gCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATc                                                                                                                                        >  6:130643/1‑137 (MQ=255)
                    gCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATc                                                                                                                                        <  5:130643/137‑1 (MQ=255)
                     ctctCCCTCGTTGTTGAAGACCTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCTTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTAt                                                                                                                           <  2:277621/149‑1 (MQ=255)
                     cTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATc                                                                                                                                        <  7:228229/136‑1 (MQ=255)
                     cTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGAGTAGTTGATCGAGCGGGCGATc                                                                                                                                        >  8:228229/1‑136 (MQ=255)
                       cGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATgg                                                                                                                         >  7:237400/1‑149 (MQ=255)
                        gCCCTCGTTGTTGAAGACCTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATgg                                                                                                                         <  5:79199/148‑1 (MQ=255)
                         cccTCGTTGTTGAAGACCTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGtt                                                                                                                       <  2:103444/149‑1 (MQ=255)
                           cTCGTTGTTGAAGACCTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTgg                                                                                                                     <  8:237400/149‑1 (MQ=255)
                              gttgttGAAGACCTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCa                                                                                                                   >  4:58765/1‑148 (MQ=255)
                                    gAAGACCTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGacac                                                                                                            >  1:288211/1‑149 (MQ=255)
                                    gAAGACCTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGaca                                                                                                             <  2:140775/148‑1 (MQ=255)
                                      aGACCTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAg                                                                                                                                >  3:140107/1‑127 (MQ=255)
                                      aGACCTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAg                                                                                                                                <  4:140107/127‑1 (MQ=255)
                                          cTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATTGTTGGCATAGACAAGCAggg                                                                                                      >  5:195959/1‑149 (MQ=255)
                                          cTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAggg                                                                                                      >  5:168080/1‑149 (MQ=255)
                                          cTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAggg                                                                                                      >  4:46032/1‑149 (MQ=255)
                                                     ccAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGccc                                                                                           <  3:46032/149‑1 (MQ=255)
                                                      cAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCca                                                                                          <  5:85227/149‑1 (MQ=255)
                                                          ttcttcGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGAt                                                                                                                                         <  2:39566/99‑1 (MQ=255)
                                                          ttcttcGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGAt                                                                                                                                         >  1:39566/1‑99 (MQ=255)
                                                               cGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGt                                                                                                                        >  7:45969/1‑111 (MQ=255)
                                                               cGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGt                                                                                                                        <  8:45969/111‑1 (MQ=255)
                                                                 gTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAg                                                                                <  2:79710/149‑1 (MQ=255)
                                                                 gTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAg                                                                                >  2:29430/1‑149 (MQ=255)
                                                                       gCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTg                                                                                                                      >  5:317546/1‑105 (MQ=255)
                                                                       gCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTg                                                                                                                      <  6:317546/105‑1 (MQ=255)
                                                                               cGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCa                                                                  <  5:280109/149‑1 (MQ=255)
                                                                                  cAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGc                                                               <  3:190869/149‑1 (MQ=255)
                                                                                    gCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAg                                                                                                                >  4:180790/1‑98 (MQ=255)
                                                                                    gCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAg                                                                                                                <  3:180790/98‑1 (MQ=255)
                                                                                           gTTGGAGGCGGCGCTGTCATCGAACACCACAGCCTCTTCTACGTTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGgcaa                                                      <  6:195959/149‑1 (MQ=255)
                                                                                               gcggcggcgCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGgcaagc                                                    >  2:61213/1‑147 (MQ=255)
                                                                                                                gAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTgg                                  <  1:209368/148‑1 (MQ=255)
                                                                                                                     ccacAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACgg                            <  1:61213/149‑1 (MQ=255)
                                                                                                                      cacaGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGc                           >  6:329729/1‑149 (MQ=255)
                                                                                                                       acaGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGcc                          >  6:280815/1‑149 (MQ=255)
                                                                                                                          gCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGTGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAg                       <  2:288211/149‑1 (MQ=255)
                                                                                                                                  tACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGgcaa                                                      >  1:222168/1‑110 (MQ=255)
                                                                                                                                  tACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGgcaa                                                      <  2:222168/110‑1 (MQ=255)
                                                                                                                                      gtgtAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAg           <  6:168080/149‑1 (MQ=255)
                                                                                                                                            ttGATCGCCCGCTCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAgg                                                                >  1:102862/1‑90 (MQ=255)
                                                                                                                                            ttGATCGCCCGCTCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAgg                                                                <  2:102862/90‑1 (MQ=255)
                                                                                                                                            ttGATCGCCCGCTCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCgg     >  7:120407/1‑149 (MQ=255)
                                                                                                                                            ttGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCgg     >  3:74808/1‑149 (MQ=255)
                                                                                                                                               aTCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGgcg  >  3:982/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCAC‑TTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCG  >  NC_002947/5297045‑5297335

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