Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,145,777	G→A	G838S (GGC→AGC)	dgcE →	putative diguanylate cyclase DgcE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,145,777	0	G	A	91.9%	89.6 / ‑1.5	38	G838S (GGC→AGC)	dgcE	putative diguanylate cyclase DgcE
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base A (17/17);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

TCACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAATGGTGGCCGGGGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATC  >  NC_000913/2145638‑2145914                                                                                                                                             |                                                                                                                                           tCACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAATGGTggccggagcc                                                                                                                                         >  1:633856/1‑143 (MQ=255) tCACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCAGTTATGCCTCC‑AAAAATGGTggccggggcc                                                                                                                                         >  2:283270/1‑143 (MQ=255)    cTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAATGGTGgccggagccgg‑g                                                                                                                                     >  1:518837/1‑143 (MQ=255)      ttATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAATGGTGGCCGGAGCCGG‑GTg                                                                                                                                   >  1:713415/1‑143 (MQ=255)        atatGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGACTCC‑AAAAATGGTGGCCGGAGCCGG‑GTGAc                                                                                                                                 <  1:505236/143‑1 (MQ=255)           tGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAATGGTGGCCGGGGCCGG‑GTGACGGt                                                                                                                              <  1:406100/143‑1 (MQ=255)               aaGGCCGTGTGCATCGGGTAGGTGGCAGTGACGGGATAACCGGGATGGAGGAAGACAATCAGCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAATGGTGGCCGGAGCCGG‑GTGACGGTTTAc                                                                                                                          <  1:163570/143‑1 (MQ=255)               aaGGCCGTCTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAAATGGTGGCCGGAGCCGG‑GTGACGGTTTa                                                                                                                           <  2:549773/143‑1 (MQ=255)                        tgCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAATGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCgcag                                                                                                                 <  1:170124/143‑1 (MQ=255)                        tgCAGCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAATGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCgcag                                                                                                                 <  1:564555/143‑1 (MQ=255)                          cATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAATGGTGGCCGGAGCCGGCGCGCCGG‑TTCCGAACCgcagca                                                                                                               >  2:728583/1‑143 (MQ=255)                                gaaagggCCGGTGCCGGGCTCACCTTGATTGATGACAACAAACATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAACGGTGGCCGGAGCCGG‑GTGACGGTTTACTAACCGCAGCAAGCTga                                                                                                         <  1:664216/137‑2 (MQ=255)                                    gtgcgggTGCCGGGATGACCTTGATGGATGACGACAAGCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAACTGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGca                                                                                                     <  2:754581/137‑1 (MQ=255)                                         aGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAATGGTGGCCGGCGCCGG‑GTGACGGTTTACGAACCGCAGCCAGCTGCCGCACATAg                                                                                                >  1:429812/1‑143 (MQ=255)                                             ccGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGGTATCGCCTGTTATGCCTCC‑AAAAATGGTGGCCGAAGCCGG‑GTGACGGTTTTCCAACCGCAGCCAGCTGCCGCACGTAGCGaa                                                                                            >  2:591636/1‑143 (MQ=255)                                                       cTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAATGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGa                                                                                  >  1:234120/1‑143 (MQ=255)                                                       cTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCAGTTATGCCTCC‑AAAAATGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGGAGATGGCGAAGGggcggggg                                                                                  >  1:8364/1‑140 (MQ=255)                                                          gattgatGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAATGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGGGATGt                                                                               >  2:774631/1‑143 (MQ=255)                                                              gatgaCAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAATGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCt                                                                           <  1:267005/143‑1 (MQ=255)                                                               atgaCAACAGTCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAATGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCtt                                                                          <  1:553784/143‑1 (MQ=255)                                                                tgaCAACAATCATCAGGCGGAAGAAGTGATGTAGCAGGCTGATATCGCCTGCTATGCCTCC‑AAAAATGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCGCACCAAGCTGCCGCACACAGCGAACGGGCGGCGCTGTCGCCtg                                                                         >  2:417202/1‑143 (MQ=255)                                                                  acaacaATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAATGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTtgat                                                                       <  1:393969/143‑1 (MQ=255)                                                                      caaTCATCGGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAACGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAAc                                                                   <  2:45420/143‑1 (MQ=255)                                                                               gccgggAGAAGTGATGTCGCGGGCTGATATCGCCTGTTATGCCTCC‑AAAGAGGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTggagga                                                          <  2:325482/141‑1 (MQ=255)                                                                                    cAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAATGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGgatgagt                                                     >  1:632293/1‑141 (MQ=255)                                                                                        aGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAATGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAag                                                 <  1:411085/143‑1 (MQ=255)                                                                                           gATGTCGCCGGCTGATATCGCCTGTTATGCCGCC‑AAAAATGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAgaga                                              <  2:467296/143‑1 (MQ=255)                                                                                           gATGTCGCAGGCTGATAGCGCCGGTTAGGCCTCC‑AAAAATGGTGGCCGGAGCCGG‑GTGACGGGTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAgaga                                              <  1:76935/143‑1 (MQ=255)                                                                                            aTGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAATGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGaa                                             >  1:385015/1‑143 (MQ=255)                                                                                            aTGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAATGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCACAGCAAGCTGCCGCACATAGCGAACGGGAGGCGATGACGCCGGAGGAGCAGTGGCGGAAGAAGAAagggga                                             >  2:142089/1‑139 (MQ=255)                                                                                               tCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAACTGGTGGCCGGGGCCGG‑GTGACGGTTTATGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCa                                          <  1:12872/143‑1 (MQ=255)                                                                                                       tGATATCGCCTGTTATGCCTCC‑AAAAATGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGACTCAGTtgatga                                  >  1:386333/1‑143 (MQ=255)                                                                                                              gCCTGTTATGCCTCC‑AAAAATGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGGACAGTGGCGGATGATTAAAGAGAGTCGGGTGGTGATGCTCGc                           >  1:692485/1‑143 (MQ=255)                                                                                                               ccTGTTATGCCTCC‑AAAAATGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGGTGATTAAAGAGAAGCAGTTGCTGATGCTCGcc                          >  1:227621/1‑143 (MQ=255)                                                                                                                 gggcaTGCCGCC‑AACAAAGGTGGCCGGAGCCGG‑CGTACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTTGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCa                        <  2:129231/139‑1 (MQ=255)                                                                                                                     aTGCCTCC‑AAAAATGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGACTAAAGAGAATCAGTTGATGAAGCTCGCCCACGgt                    >  2:535634/1‑143 (MQ=255)                                                                                                                      tGCCTAC‑AAAAATGGTGGTCGGAGCAGG‑GCGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACCGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGgtg                   <  1:136535/143‑1 (MQ=255)                                                                                                                         cTCC‑AAAAATGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCg                >  1:176604/1‑143 (MQ=255)                                                                                                                              aaaaaTGGTGGCCGGGGCCGG‑GCGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTc            <  1:283270/143‑1 (MQ=255)                                                                                                                              aaaaaTGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTc            >  1:480186/1‑143 (MQ=255)                                                                                                                              aaaaaTGGTGGCCGGAGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTc            >  1:339705/1‑143 (MQ=255)                                                                                                                                       ggccggggccgg‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCCTGATGAACAGTGGCGGATGATTAAAGAGAATCAGGTGATGATGCTCGCCCACGGTGCGCCCTCGCCACGGag   >  2:54311/1‑142 (MQ=255)                                                                                                                                        tccggagccgg‑cGGACCGTCTACGGCCCGCCGCAGGCTGCAGCACATAGCGAACGGGCGGCGATGTCGCTTGAGGAACAGTGGCGGATGATTAAAGAGAATCAGTGCATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATc  <  2:116950/142‑1 (MQ=255)                                                                                                                                             |                                                                                                                                           TCACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCC‑AAAAATGGTGGCCGGGGCCGG‑GTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATC  >  NC_000913/2145638‑2145914

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑