Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,649,236 A→T 8.0% C653S (TGC→AGC)  PP_4959 ← diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,649,2360AT 8.0% 143.6 / 2.8 50C653S (TGC→AGC) PP_4959diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (21/25);  new base T (2/2);  total (23/27)
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.86e-01

CGGGAAGCCGCTAGCGTCAGTCGCGGTCATCGGATAAGCCACCAGGAAACCCTGAACGTATTCGCAACCATTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGCGCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCC  >  NC_002947/5649093‑5649384
                                                                                                                                               |                                                                                                                                                    
cGGGAAGCCGCTAGCGTCAGTCGCGGTCATCGGATAAGCCACCAGGAAACCCTGAACGTATTCGCAACCATTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGc                                                                                                                                                 <  5:131338/149‑1 (MQ=255)
 gggAAGCCGCTAGCGTCAGTCGCGGTCATCGGATAAGCCACCAGTAAACCCTGAACGTATTCGCAACCATTTGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCTTGATACGGTAGTCGGCGCACAGGcc                                                                                                                                                <  7:202350/149‑1 (MQ=255)
     aGCCGCTAGCGTCAGTCGCGGTCATCGGATAAGCCACCAGGAAACCCTGAACGTATTCGCAACCATTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATc                                                                                                                                            <  3:116227/149‑1 (MQ=255)
                          tCATCGGATAAGCCACCAGGAAACCCTGAACGTATTCGCAACCATTGGCCTTTAGCCAGGCGGCCTGCGCCGGTGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACgcagca                                                                                                                       <  8:38914/149‑1 (MQ=255)
                            aTCGGATAAGCCACCAGGAAACCCTGAACGTATTCGCAACCATTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGACGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATc                                                                                                                     >  6:39312/1‑149 (MQ=255)
                                gATAAGCCACCAGGAAACCCTGAACGTATTCGCAACCATTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTAcc                                                                                                                 <  2:77556/149‑1 (MQ=255)
                                      ccaccaGGAAACCCTGAACGTATTCGCAACCATTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGaa                                                                                                           >  6:147012/1‑149 (MQ=255)
                                       caccaGGAAACCCTGAACGTATTCGCAACCATTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAAt                                                                                                          >  2:68479/1‑149 (MQ=255)
                                                       aCGTATTCGCAACCATTGGCCTTTAGCCAGGCGGCCTGCTCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGccac                                                                                          <  2:300132/149‑1 (MQ=255)
                                                         gTATTCGCAACCATTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGACCCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGccac                                                                                          >  5:317051/1‑147 (MQ=255)
                                                         gTATTCGCAACCATTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGACGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGccacca                                                                                        <  5:39312/149‑1 (MQ=255)
                                                           aTTCGCAACCATTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGccac                                                                                          >  7:46279/1‑145 (MQ=255)
                                                            ttCGCAACCATTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTc                                                                                                              <  4:124301/124‑1 (MQ=255)
                                                            ttCGCAACCATTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTc                                                                                                              >  3:124301/1‑124 (MQ=255)
                                                            ttCGCAACCATTCGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGa                                                                                      <  8:109873/148‑1 (MQ=255)
                                                                 aaCCATTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACgg                                                                                >  2:320174/1‑149 (MQ=255)
                                                                    cATTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTc                                                                              >  2:130349/1‑148 (MQ=255)
                                                                    cATTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTc                                                                              <  1:130349/148‑1 (MQ=255)
                                                                              tAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAgg                                                                    <  5:147012/148‑1 (MQ=255)
                                                                              tAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAg                                                                     <  6:195688/147‑1 (MQ=255)
                                                                               aGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGCTGACGGTCAAGCTTGAGGgt                                                                  >  1:52336/1‑149 (MQ=255)
                                                                                 ccAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGc                                                                          <  1:875/139‑1 (MQ=255)
                                                                                 ccAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGc                                                                          >  2:875/1‑139 (MQ=255)
                                                                                 ccAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAg                                                                     >  7:136067/1‑144 (MQ=255)
                                                                                 ccAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAg                                                                     <  8:136067/144‑1 (MQ=255)
                                                                                 ccAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGgtgt                                                                >  3:162504/1‑149 (MQ=255)
                                                                                   aGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTc                                                               >  4:296111/1‑148 (MQ=255)
                                                                                         ccTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGc                                                        <  1:68479/149‑1 (MQ=255)
                                                                                                      ctcCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCggg                                                                                                      >  2:15430/1‑90 (MQ=255)
                                                                                                      ctcCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCggg                                                                                                      <  1:15430/90‑1 (MQ=255)
                                                                                                       tcCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAg                                                      <  2:35461/137‑1 (MQ=255)
                                                                                                       tcCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAg                                                      >  1:35461/1‑137 (MQ=255)
                                                                                                         cacaCCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAgg                                                                    >  6:240997/1‑121 (MQ=255)
                                                                                                         cacaCCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAgg                                                                    <  5:240997/121‑1 (MQ=255)
                                                                                                         cacaCCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGCGCAGCAt                                          >  3:207410/1‑147 (MQ=255)
                                                                                                           cacCCTCGGCGATCACCGTGATACGGTAGTCGGCGCTCAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGCGCAGCATgcg                                       >  3:241086/1‑148 (MQ=255)
                                                                                                           cacCCTCGGCGATCACCGTGATACGGTAGTCGGCGCTCAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGCGCAGCATgcg                                       >  3:236846/1‑148 (MQ=255)
                                                                                                           cacCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGCGCAGCATgcgc                                      <  3:296111/149‑1 (MQ=255)
                                                                                                            acCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGCGCAGCATgcgc                                      >  1:82561/1‑148 (MQ=255)
                                                                                                              ccTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGCGCAGCATGCGCAAc                                   >  7:203751/1‑149 (MQ=255)
                                                                                                              ccTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGCGCAGCATGCGCAAc                                   <  1:320174/149‑1 (MQ=255)
                                                                                                                  ggCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGCGCAGCATGCGCAACGaa                                <  8:46279/148‑1 (MQ=255)
                                                                                                                   gCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGCGCAGCATGCGCAACGAACa                              <  2:52336/149‑1 (MQ=255)
                                                                                                                    cGATCACCGTGATACGGTAGTCGGCGCTCAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGCGCAGCATGCGCAAcctgcct                             <  4:241086/149‑7 (MQ=255)
                                                                                                                    cGATCACCGTGATACGGTAGTCGGCGCTCAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGCGCAGCATGCGCAAcctgcct                             <  4:236846/149‑7 (MQ=255)
                                                                                                                          ccGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGCGCAGCATgc                                        <  7:279008/132‑1 (MQ=255)
                                                                                                                          ccGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGCGCAGCATgc                                        >  8:279008/1‑132 (MQ=255)
                                                                                                                                        tCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGCGCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATc         <  2:82561/149‑1 (MQ=255)
                                                                                                                                           gcgcACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGc                                                        <  2:210111/99‑1 (MQ=255)
                                                                                                                                           gcgcACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGc                                                        >  1:210111/1‑99 (MQ=255)
                                                                                                                                               acaGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGCGCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGcc  >  8:28307/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CGGGAAGCCGCTAGCGTCAGTCGCGGTCATCGGATAAGCCACCAGGAAACCCTGAACGTATTCGCAACCATTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGCGCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCC  >  NC_002947/5649093‑5649384

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