Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,233,578 A→G 7.8% V496A (GTC→GCC) ‡ PP_3711 ← diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,233,5780AG 7.8% 140.3 / 2.3 51V496A (GTC→GCC) ‡PP_3711diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (24/23);  new base G (2/2);  total (26/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 = 9.20e-01

ATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTGGTAGTGCAACTGCAACTGCCCCAGGTGCAGCGCTGTGCGCAGGTCGGTTTCCACCAGGCGGC  >  NC_002947/4233431‑4233722
                                                                                                                                                   |                                                                                                                                                
atcaaGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGAc                                                                                                                                                 >  2:146812/1‑149 (MQ=255)
 tcaaGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGAc                                                                                                                                                 <  3:18010/148‑1 (MQ=255)
 tcaaGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGAc                                                                                                                                                 >  8:15355/1‑148 (MQ=255)
 tcaaGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTAGTGGGTTTTGAGACTGAAGGGGTATGTGTCTAGATAATGAAGTCCATCGGGATGag                                                                                                                                                 >  4:18010/1‑147 (MQ=255)
      agCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGcc                                                                                                                                           <  2:196248/149‑1 (MQ=255)
          tCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCAcgg                                                                                                                                       <  2:158389/149‑1 (MQ=255)
                ccccaaaggCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCAcgggtcggg                                                                                                                                 <  4:179993/142‑1 (MQ=255)
                ccTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCAcgggtcggg                                                                                                                                 <  2:102972/149‑1 (MQ=255)
                  tttGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGgtg                                                                                                                               >  1:203516/1‑149 (MQ=255)
                          tttGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGGTCATGccgatgg                                                                                                                                      >  3:157520/1‑129 (MQ=255)
                          tttGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGGTCATGccgatgg                                                                                                                                      <  4:157520/134‑6 (MQ=255)
                                cTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCgca                                                                                                                    >  4:193644/1‑146 (MQ=255)
                                cTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCgca                                                                                                                    <  3:193644/146‑1 (MQ=255)
                                              ggagcaggagcaggagcaggagcaggagcaggagcagTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCg                                                                                                   <  2:203516/149‑1 (MQ=255)
                                              ggagcaggagcaggagcaggagcaggagcaggagcagTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCg                                                                                                   >  5:169062/1‑149 (MQ=255)
                                                 gcaggagcaggagcaggagcaggagcaggagcagTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGt                                                                                                 >  8:240130/1‑148 (MQ=255)
                                                        caggagcaggagcaggagcaggagcagTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACgctgc                                                                                         <  8:68812/149‑1 (MQ=255)
                                                         aggagcaggagcaggagcaggagcagTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACgctgct                                                                                        >  1:158552/1‑149 (MQ=255)
                                                           gagcaggagcaggagcaggagcagTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAgcgc                                                                                                                     <  6:228146/118‑1 (MQ=255)
                                                           gagcaggagcaggagcaggagcagTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCAATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAgcgc                                                                                                                     >  5:228146/1‑118 (MQ=255)
                                                              caggagcaggagcaggagcagTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCt                                                                                   <  6:169062/149‑1 (MQ=255)
                                                               aggagcaggagcaggagcagTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCtg                                                                                  <  1:146812/149‑1 (MQ=255)
                                                                 gagcaggagcaggagcagTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCtgat                                                                                >  6:173611/1‑149 (MQ=255)
                                                                 gagcaggagcaggagcagTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCtgat                                                                                <  7:240130/149‑1 (MQ=255)
                                                                     aggagcaggagcagTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATc                                                                            <  2:158552/149‑1 (MQ=255)
                                                                     aggagcaggagcagTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATc                                                                            >  6:38842/1‑149 (MQ=255)
                                                                             gagcagTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATc                                                                            <  5:154442/141‑1 (MQ=255)
                                                                             gagcagTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATc                                                                            >  6:154442/1‑141 (MQ=255)
                                                                                 agTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTg                                                                <  5:38842/149‑1 (MQ=255)
                                                                                   tCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTgg                                                               >  1:307578/1‑148 (MQ=255)
                                                                                     aaaaCCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACgctgc                                                                                         <  7:66645/120‑1 (MQ=255)
                                                                                     aaaaCCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACgctgc                                                                                         >  8:66645/1‑120 (MQ=255)
                                                                                     aaaaCCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTGGTa                                                             <  1:190451/148‑1 (MQ=255)
                                                                                           ttGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTaa                                                                                               >  2:223230/1‑108 (MQ=255)
                                                                                           ttGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTaa                                                                                               <  1:223230/108‑1 (MQ=255)
                                                                                                  ttttGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTGGTAGTgcaactgcaact                                               <  5:173611/149‑1 (MQ=255)
                                                                                                  ttttGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTGGTAGTgcaactgcaact                                               >  7:252533/1‑149 (MQ=255)
                                                                                                       acacTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTGGTAGTGCAACTGCAACTGccc                                           >  5:41259/1‑148 (MQ=255)
                                                                                                               ggggTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTGGTAGTGCAACTGCACCTGGGGCAGTTGCAg                                  <  7:15355/149‑1 (MQ=255)
                                                                                                                           tAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTGGTAGTGCAACTGCAACTGCCCCAGGTGCAGCGCTGTGCGCAg                      >  2:88668/1‑149 (MQ=255)
                                                                                                                           tAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTGGTAGTGCAACTGCAACTGCCCCAGGTGCAGCGCTGTGCGCAg                      >  5:135120/1‑149 (MQ=255)
                                                                                                                           tAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTGGTAGTGCAACTGCAACTGCCCCAGGTGCAGCGCTGTGCGCAg                      >  2:92786/1‑149 (MQ=255)
                                                                                                                             gAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTGGTAGTGCAACTGCAACTGCCCCAGGTGCAGCGCTGTGCGCa                       <  6:41259/146‑1 (MQ=255)
                                                                                                                               aaaTGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTGGTAGTGCAACTGCAACTGCCCCAGGTGCAGCGCTGTGCGCAGGTCg                  <  6:189199/149‑1 (MQ=255)
                                                                                                                                aaTGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTGGTAGTGCAACTGCAACTGCCCCAGGTGCAGCGCTGTGCGCAGGTCgg                 >  8:254746/1‑149 (MQ=255)
                                                                                                                                   gAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTGGTAGtgca                                                        >  4:175257/1‑107 (MQ=255)
                                                                                                                                   gAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTGGTAGtgca                                                        <  3:175257/107‑1 (MQ=255)
                                                                                                                                          cgtggCATGGTCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTGGTAGTGCAACTGCAACTGCCCCAGGTGCAGCGCTGTGCGCAGGTCGGTTTCCACCAg       >  8:127633/4‑149 (MQ=255)
                                                                                                                                           gtggCATGGTCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTGGTAGTGCAACTGCAACTGCCCCAGGTGCAGCGCTGTGCGCAGGTCGGTTTCCACCAgg      <  7:127633/147‑1 (MQ=255)
                                                                                                                                              gCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTGGTAGTGCAACTGCAACTGCCCCAGGTGCAGCGCTGTGCGCAGGTCGGTTTCCACCAggcgg   >  3:30809/1‑149 (MQ=255)
                                                                                                                                              gCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTGGTAGTGCAACTGCAACTGCCCCAGGTGCAGCGCTGTGCGCAGGTCGGTTTCCACCAggcgg   >  4:316086/1‑149 (MQ=255)
                                                                                                                                               cATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTGGTAGTGCAACTGCAACTGCCCCAGGTGCAGCGCTGTGCGCAGGTCGGTTTCCACCAggcggc  >  8:156648/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCCCAGCGCAGCAGCGCCTCATAACCGGTAACGCTGCTGGTCTGATGATCGACTACCGGCTGGTAGTGCAACTGCAACTGCCCCAGGTGCAGCGCTGTGCGCAGGTCGGTTTCCACCAGGCGGC  >  NC_002947/4233431‑4233722

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