Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,963,802 T→A 12.8% L5Q (CTA→CAA) ‡ yegE → phosphodiesterase domain‑containing diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,963,8020TA12.8% 122.5 / 8.7 47L5Q (CTA→CAA) ‡yegEphosphodiesterase domain‑containing diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (20/21);  new base A (3/3);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.08e-01

GCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTCATCGCCGTCGTGGTCCTTCA  >  NC_002947/1963656‑1963934
                                                                                                                                                  |                                                                                                                                    
gCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTAc                                                                                                                                    >  8:219598/1‑149 (MQ=255)
  ccGCGAACACGGGCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACg                                                                                                                                   <  7:219598/148‑1 (MQ=255)
  ccGCGAACACGGGCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGc                                                                                                                                  >  3:251899/1‑149 (MQ=255)
  ccGCGAACACGGGCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGc                                                                                                                                  >  6:51554/1‑149 (MQ=255)
  ccGCGAACACGGGCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGc                                                                                                                                  <  3:103733/149‑1 (MQ=255)
   cgcgAACACGGGCGTAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCt                                                                                                                                 <  7:144391/149‑1 (MQ=255)
        acacGGGCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCt                                                                                                                            >  6:162028/1‑149 (MQ=255)
            gggCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTc                                                                                                                        >  5:41139/1‑149 (MQ=255)
                        tggcttCCACCGCGGTGTCTTCTTCGCGGGTTAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATTTCCTACTTTGTGCTTGTCTGGGTCGATTTa                                                                                                            <  7:41309/144‑1 (MQ=38)
                          ccATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACTCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTAc                                                                                                           <  6:210571/148‑1 (MQ=255)
                                     ggTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGct                                                                                               <  5:51554/149‑1 (MQ=255)
                                           ttcttcGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGTAGGCCATTTGCCATCGCtgtggagg                                                                                                                          >  2:258874/1‑112 (MQ=37)
                                           ttcttcGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGTAGGCCATTTGCCATCGCtgtggagg                                                                                                                          <  1:258874/116‑5 (MQ=37)
                                             cttcGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGa                                                                                         <  2:129853/147‑1 (MQ=255)
                                              ttcGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAg                                                                                      >  5:114514/1‑149 (MQ=255)
                                                cgcgGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACa                                                                                       <  4:251899/146‑1 (MQ=255)
                                                     gTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTTTTAATTTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTcc                                                                               <  1:225079/149‑1 (MQ=255)
                                                           ccGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGaa                                                                          <  7:45981/148‑1 (MQ=255)
                                                           ccGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGCCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATg                                                                            >  8:45981/1‑146 (MQ=255)
                                                            cGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAgg                                                                        <  5:162028/149‑1 (MQ=255)
                                                               tCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGa                                                                      >  8:12564/1‑148 (MQ=255)
                                                               tCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGa                                                                      >  8:37541/1‑148 (MQ=255)
                                                                ccTACAATAGCCAGCGCATACTCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTAAGCCTTATATGCTATGTGACAGGGACTCCATGAAg                                                                         >  6:148078/1‑144 (MQ=255)
                                                                ccTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAg                                                                         <  5:148078/144‑1 (MQ=255)
                                                                 cTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACAt                                                                   <  6:41139/149‑1 (MQ=255)
                                                                  tACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGTAGGCCATTTGCCATCGCtgtggag                                                                                                                           >  1:132813/1‑89 (MQ=25)
                                                                  tACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGTAGGCCATTTGCCATCGCtgtggag                                                                                                                           <  2:132813/92‑4 (MQ=25)
                                                                                 aTACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAgcgc                                                   >  6:141367/1‑149 (MQ=255)
                                                                                         tACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGaa                                                                          >  1:27458/1‑118 (MQ=255)
                                                                                         tACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGaa                                                                          <  2:27458/118‑1 (MQ=255)
                                                                                             aaTTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCg                                       <  5:141367/149‑1 (MQ=255)
                                                                                                tCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCATACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCa                                    >  7:114491/1‑149 (MQ=255)
                                                                                                   ccGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCa                                 >  3:146438/1‑149 (MQ=255)
                                                                                                      aTTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCt                              >  4:18199/1‑149 (MQ=255)
                                                                                                                 tAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTCATc                   <  7:12564/149‑1 (MQ=255)
                                                                                                                 tAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTCATc                   >  3:191857/1‑149 (MQ=255)
                                                                                                                         cccTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTCATCGCcgtcgt           >  8:235269/1‑149 (MQ=255)
                                                                                                                          ccTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTCATCGCCGTCGTg          >  2:87310/1‑149 (MQ=255)
                                                                                                                          ccTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTCATCGCCGTCGTg          <  3:18199/149‑1 (MQ=255)
                                                                                                                           cTCCAAATCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGTTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTCATCGCCGTCGTgg         <  8:114491/149‑1 (MQ=255)
                                                                                                                             gcacaGCGTAGGCCATTTGCCATCGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGc                           <  8:8609/128‑1 (MQ=255)
                                                                                                                             gcacaGCGTAGGCCATTTGCCATCGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGc                           >  7:8609/2‑129 (MQ=255)
                                                                                                                                  gCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTCATCGCCGTCGTGGTCCTTCa  <  7:235269/149‑1 (MQ=255)
                                                                                                                                          aaaTGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCtt                        >  6:139115/1‑119 (MQ=255)
                                                                                                                                          aaaTGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCtt                        <  5:139115/119‑1 (MQ=255)
                                                                                                                                               gCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTCATCGCCGTCGTGGTCCTTc   <  7:214697/135‑1 (MQ=255)
                                                                                                                                               gCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTAGCTCGGCATCACCTGGACTTTCATCGCCGTCGTGGTCCtta   >  8:214697/1‑134 (MQ=255)
                                                                                                                                                  |                                                                                                                                    
GCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTCATCGCCGTCGTGGTCCTTCA  >  NC_002947/1963656‑1963934

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: ‑