Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,402,000 A→G I18T (ATC→ACC)  dgcM ← diguanylate cyclase DgcM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,402,0000AG88.9% 63.9 / 2.8 27I18T (ATC→ACC) dgcMdiguanylate cyclase DgcM
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (10/14);  total (12/15)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.96e-01

GTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGC  >  NZ_CP009273/1401917‑1402085
                                                                                   |                                                                                     
gTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGa                                                                                 >  2:768684/1‑90 (MQ=255)
 tAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGAc                                                                                >  1:564532/1‑90 (MQ=255)
      gACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGAc                                                                           <  1:790782/90‑1 (MQ=255)
      gACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGAc                                                                           >  2:635160/1‑90 (MQ=255)
      gACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGAc                                                                           <  1:564677/90‑1 (MQ=255)
              gTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGt                                                                   <  2:710259/90‑1 (MQ=255)
               tGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                                                  <  1:59116/90‑1 (MQ=255)
                    tCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTcc                                                             <  1:488690/90‑1 (MQ=255)
                      ttGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAg                                                           <  2:727239/90‑1 (MQ=255)
                      ttGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAg                                                           <  1:578640/90‑1 (MQ=255)
                      ttGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAg                                                           <  1:768684/90‑1 (MQ=255)
                             aaCAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTg                                                    >  2:116179/1‑90 (MQ=255)
                                        cgcgCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGt                                         <  2:564532/90‑1 (MQ=255)
                                        cgcgCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGt                                         <  1:269728/90‑1 (MQ=255)
                                         gcgcTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTa                                        >  1:502269/1‑90 (MQ=255)
                                         gcgcTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTa                                        >  2:48508/1‑90 (MQ=255)
                                         gcgcTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTa                                        <  2:473/90‑1 (MQ=255)
                                                      aaaTTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATcc                           >  2:261438/1‑90 (MQ=255)
                                                      aaaTTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATcc                           >  2:438412/1‑90 (MQ=255)
                                                           aaCTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                                                  >  2:260322/1‑46 (MQ=255)
                                                           aaCTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                                                  <  1:260322/46‑1 (MQ=255)
                                                                ttCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGc                                           >  2:323352/1‑64 (MQ=255)
                                                                ttCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGc                                           <  1:323352/64‑1 (MQ=255)
                                                                 tCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTa                >  1:349662/1‑90 (MQ=255)
                                                                         gggCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAAtttt        <  1:116179/90‑1 (MQ=255)
                                                                           gCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTttg      >  1:262321/1‑90 (MQ=255)
                                                                               aaCGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGc  <  2:502269/90‑1 (MQ=255)
                                                                                   |                                                                                     
GTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGC  >  NZ_CP009273/1401917‑1402085

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