Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,139,649 A→G E968G (GAA→GGA)  dgcE → diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,139,6490AG95.5% 147.7 / ‑1.1 44E968G (GAA→GGA) dgcEdiguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (14/28);  total (15/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.85e-01

ACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGAAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCCGCGATC  >  NZ_CP009273/2139514‑2139780
                                                                                                                                       |                                                                                                                                   
aCCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGgaaa                                                                                                                                  <  2:749152/139‑1 (MQ=255)
   gccgGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCTACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAg                                                                                                                               <  1:764617/139‑1 (MQ=255)
   gccgGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAg                                                                                                                               <  2:847022/139‑1 (MQ=255)
   gccgGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAg                                                                                                                               <  1:125791/139‑1 (MQ=255)
       ggTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCct                                                                                                                           >  1:1000239/1‑139 (MQ=255)
           ttCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTac                                                                                                                       <  2:41415/139‑1 (MQ=255)
               aCGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTaccacc                                                                                                                   <  2:74448/139‑1 (MQ=255)
                  aaTTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACg                                                                                                                <  2:558006/139‑1 (MQ=255)
                  aaTTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACg                                                                                                                <  2:318939/139‑1 (MQ=255)
                   aTTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGgaaa                                                                                                                                  <  2:514958/120‑1 (MQ=255)
                   aTTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGgaaa                                                                                                                                  >  1:514958/1‑120 (MQ=255)
                                  cGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGat                                                                                                 >  2:278154/1‑138 (MQ=255)
                                                ccAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAg                                                                                   <  1:901408/138‑1 (MQ=255)
                                                ccAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAg                                                                                   <  2:1236130/138‑1 (MQ=255)
                                                ccAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAg                                                                                   <  2:1012507/138‑1 (MQ=255)
                                                ccAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAg                                                                                   <  1:156795/138‑1 (MQ=255)
                                                ccAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAg                                                                                   <  1:442876/138‑1 (MQ=255)
                                                          cATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGAt                                                                        <  1:622648/139‑1 (MQ=255)
                                                          cATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGAt                                                                        <  2:401279/139‑1 (MQ=255)
                                                          cATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGAt                                                                        <  1:278154/139‑1 (MQ=255)
                                                          cATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGAt                                                                        <  1:930474/139‑1 (MQ=255)
                                                              aGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACg                                                                    >  2:1089243/1‑139 (MQ=255)
                                                                  tctcCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTaccacc                                                                                                                   <  2:312738/88‑1 (MQ=255)
                                                                  tctcCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTaccacc                                                                                                                   >  1:312738/1‑88 (MQ=255)
                                                                   ctccctccCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTaccac                                                                                                                    <  1:663257/86‑1 (MQ=255)
                                                                   ctccctccCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTaccac                                                                                                                    >  2:663257/1‑86 (MQ=255)
                                                                      cctcccCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAg                                                             <  1:1081647/138‑1 (MQ=255)
                                                                                      gCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGc                                            >  1:237403/1‑139 (MQ=255)
                                                                                              gAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCggg                                     <  2:237403/138‑1 (MQ=255)
                                                                                              gAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCggg                                     <  1:494608/138‑1 (MQ=255)
                                                                                                   gcgcCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTc                                <  1:1089243/138‑1 (MQ=255)
                                                                                                     gcCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAg                                                                                   >  2:227409/1‑85 (MQ=255)
                                                                                                     gcCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAg                                                                                   <  1:227409/85‑1 (MQ=255)
                                                                                                     gcCACGCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGa                             >  2:822579/1‑139 (MQ=255)
                                                                                                          gCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTg                        <  2:744707/139‑1 (MQ=255)
                                                                                                          gCTGGTGAATGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTg                        <  1:822579/139‑1 (MQ=255)
                                                                                                                 aaTGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGttat                                                                                                           <  2:771817/49‑1 (MQ=255)
                                                                                                                 aaTGATCTGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGttat                                                                                                           >  1:771817/1‑49 (MQ=255)
                                                                                                                        tGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTggg          >  1:302633/1‑139 (MQ=255)
                                                                                                                        tGCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTggg          >  2:1019954/1‑139 (MQ=255)
                                                                                                                        tGCTTGAGCAGCTGGAAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTgg                                         <  2:474523/108‑1 (MQ=255)
                                                                                                                        tGCTTGAGCAGCTGGAAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTgg                                         >  1:474523/1‑108 (MQ=255)
                                                                                                                         gCTTGAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGc         >  1:505037/1‑139 (MQ=255)
                                                                                                                             gAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCggg                                     <  2:513062/107‑1 (MQ=255)
                                                                                                                             gAGCAGCTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCggg                                     >  1:513062/1‑107 (MQ=255)
                                                                                                                              agcagcTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCcgcg     <  1:714595/138‑1 (MQ=255)
                                                                                                                                cagcTGGGAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCCGCGATc  >  1:115089/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                   
ACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTTGCCGGTTTGAGTAGCGCCACGCTGGTGAATGATCTGCTTGAGCAGCTGGAAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCCGCGATC  >  NZ_CP009273/2139514‑2139780

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: