Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,233,856	G→C	7.3%	A403A (GCC→GCG)	PP_3711 ←	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,233,856	0	G	C	7.3%	169.5 / 3.4	55	A403A (GCC→GCG)	PP_3711	diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (24/27);  new base C (2/2);  total (26/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 = 5.26e-01

TTTCCACCAGGCGGCGCTGGCGGGCGGCAACGTCCAGTTCGACATCGAAGCATTCATAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTGCGCGGTATTGCACGCGGCATCCCGCCCCTCCAGGCCCGTCAGCACCACGGCAAACTCGTC  >  NC_002947/4233708‑4233996                                                                                                                                                     |                                                                                                                                              tttCCACCAGGCGGCGCTGGCGGGCGGCAACGTCCAGTTCGACATCGAAGCATTCATAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGGCAACAGGCgcgatg                                                                                                                                               >  8:178166/1‑149 (MQ=255)     caccaGGCGGCGCTGGCGGGCGGCAACGTCCAGTTCGACATCGAAGCATTCATAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCgcgatggcga                                                                                                                                           <  1:241103/149‑1 (MQ=255)     caccaGGCGGCGCTGGCGGGCGACAACGTCCAGTTCGACATCGAAGCATTCATAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCgcgatggc                                                                                                                                             >  7:62706/1‑147 (MQ=255)            cggcgCTGGCGGGCGGCAACGTCCAGTTCGACATCGAAGCATTCATAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGcgatgc‑cgat                                                                                                                                   <  8:62706/149‑1 (MQ=255)                gcTGGCGGGCGGCAACGTCCAGTTCGACATCGAAGCATTCATAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTc                                                                                                                               <  2:375867/149‑1 (MQ=255)                 cTGGCGGGCGGCAACGTCCAGTTCGACATCGAAGCATTCATAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCa                                                                                                                              <  7:233797/149‑1 (MQ=255)                          ggcAACGTCCAGTTCGACATCGAAGCATTCATAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCAc                                                                                                                             <  3:339051/141‑1 (MQ=255)                          ggcAACGTCCAGTTCGACATCGAAGCATTCATAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACcg                                                                                                                     >  1:386958/1‑149 (MQ=255)                          ggcAACGTCCAGTTCGACATCGAAGCATTCATAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTAAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCAc                                                                                                                             >  4:339051/1‑141 (MQ=255)                              aCGTCCAGTTCGACATCGAAGCATTCATAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGaa                                                                                                                 <  6:300934/149‑1 (MQ=255)                              aCGTCCAGTTCGACATCGAAGCATTCATAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGa                                                                                                                  >  5:300934/1‑148 (MQ=255)                               cGTCCAGTTCGACATCGAAGCATTCATAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAAc                                                                                                                <  6:14860/149‑1 (MQ=255)                                         gACATCGAAGCATTCATAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCtt                                                                                                      <  2:386958/149‑1 (MQ=255)                                          aCATCGAAGCATTCATAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATg                                                                                                           >  6:27242/1‑143 (MQ=255)                                          aCATCGAAGCATTCATAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATg                                                                                                           <  5:27242/143‑1 (MQ=255)                                                  gCATTCATAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGTGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCg                                                                                               <  1:308006/147‑1 (MQ=255)                                                  gCATTCATAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCg                                                                                               >  2:308006/1‑147 (MQ=255)                                                        aTAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCtt                                                                                                      >  7:112731/1‑134 (MQ=255)                                                         tAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGa                                                                                                             <  1:393199/126‑1 (MQ=255)                                                         tAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGa                                                                                                             >  2:393199/1‑126 (MQ=255)                                                         tAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAAcggcgg                                                                                      >  4:133586/1‑149 (MQ=255)                                                           tcGATTACGTCCGTTAAGCTTGGCAGTGTACAGCGCCATGTCATCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCTATCCGGAACGGCggg                                                                                     <  4:244492/147‑1 (MQ=255)                                                                   ggCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGcc                                                                            <  7:364034/149‑1 (MQ=255)                                                                    gCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCt                                                                           <  5:155713/149‑1 (MQ=255)                                                                        ttACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGAt                                                                       <  1:238677/149‑1 (MQ=255)                                                                          aCGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACAGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGAt                                                                       >  5:211591/1‑147 (MQ=255)                                                                           cGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCa                                                                     <  3:133586/148‑1 (MQ=255)                                                                               tGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGgcgc                                                                >  7:115561/1‑149 (MQ=255)                                                                               tGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATCGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGgcgc                                                                >  2:297009/1‑149 (MQ=255)                                                                               tGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATCGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGACTTGATCAGgcgc                                                                >  2:302512/1‑149 (MQ=255)                                                                                 gCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGc                                                                                    >  2:319438/1‑127 (MQ=255)                                                                                 gCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGc                                                                                    <  1:319438/127‑1 (MQ=255)                                                                                   cTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTgcg                                                            >  8:269026/1‑149 (MQ=255)                                                                                    tGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTgcgc                                                           >  3:192725/1‑149 (MQ=255)                                                                                                 tGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTGCGCGGTATTGCAcgcg                                              <  4:192725/149‑1 (MQ=255)                                                                                                        gTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACcg                                                                                                                     <  4:366908/71‑1 (MQ=255)                                                                                                        gTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACcg                                                                                                                     >  3:366908/1‑71 (MQ=255)                                                                                                         tAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTGCGCGGTATTGCACGCGGCATCCCg                                      <  8:115561/149‑1 (MQ=255)                                                                                                            cccAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGCGCGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTGCGCGGTATTGCACGCGGCATCCCGcc                                    >  1:96510/1‑149 (MQ=255)                                                                                                             ccAGCAGTTGCTCTGCCTGGTCATGGTCTTCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTGCGCGGTATTGCACGCGGCATCCCGccc                                   >  1:228308/1‑148 (MQ=255)                                                                                                             ccAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTGCGCGGTATTGCACGCGGCATCCCGccc                                   >  6:66508/1‑148 (MQ=255)                                                                                                             ccAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTGCGCGGTATTGCACGCGGCATCCCGccc                                   <  2:228308/148‑1 (MQ=255)                                                                                                              cagcagTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTGCGCGGTATTGCACGCGGCATCCCGCCCCt                                 >  1:194074/1‑149 (MQ=255)                                                                                                                 cagTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATCGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTGCGCGGTATTGCACGCGGCATCCCGCCCCTCCa                              <  1:297009/149‑1 (MQ=255)                                                                                                                 cagTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATCGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGACTTGATCAGGCGCTGCGCGGTATTGCACGCGGCATCCCGCCCCTCCa                              <  1:302512/149‑1 (MQ=255)                                                                                                                    ttGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTGCGCGGTATTGCACGCGGCATCCCGCCCCTCCAGGc                           >  6:105205/1‑149 (MQ=255)                                                                                                                       ctctGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTGCGCGGTATTGCACGCGGCATCCCGCCCCTCCAGGCCCg                        <  7:269026/149‑1 (MQ=255)                                                                                                                        tctGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGa                                                                                 <  8:93006/91‑1 (MQ=255)                                                                                                                        tctGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGa                                                                                 >  7:93006/1‑91 (MQ=255)                                                                                                                              tggtcatggtcaTCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTGCGCGGTATTGCACGCGGCATCCCGCCCCTc                                >  1:409264/1‑134 (MQ=255)                                                                                                                              tggtcatggtcaTCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTGCGCGGTATTGCACGCGGCATCCCGCCCCTc                                <  2:409264/134‑1 (MQ=255)                                                                                                                                gtcatggtcatCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTGCGCGGTATTGCACGCGGCATCCCGCCCCTCCAGGCCCGTc                      >  7:377327/1‑142 (MQ=255)                                                                                                                                gtcatggtcatCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTGCGCGGTATTGCACGCGGCATCCCGCCCCTCCAGGCCCGTc                      <  8:377327/142‑1 (MQ=255)                                                                                                                                   atggtcatCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTGCGCGGTATTGCACGCGGCATCCCGCCCCTCCAGGCCCGTCAGCACCACg             <  2:194074/148‑1 (MQ=255)                                                                                                                                           caGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTGCGCGGTATTGCACGCGGCATCCCGCCCCTCCAGGCCCGTCAGCACCACGGCACACtcg    >  4:99134/1‑149 (MQ=255)                                                                                                                                             ggCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTGCGCGGTATTGCACGCGGCATCCCGCCCCTCCAGGCCCGTCAGCACCACGGCAAACtcgtc  <  5:66508/149‑1 (MQ=255)                                                                                                                                              gcgcGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTGCGCGGTATTGCACGCGGCATCCCGc                                     <  7:390705/113‑1 (MQ=255)                                                                                                                                              gcgcGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTGCGCGGTATTGCACGCGGCATCCCGc                                     >  8:390705/1‑113 (MQ=255)                                                                                                                                                     |                                                                                                                                              TTTCCACCAGGCGGCGCTGGCGGGCGGCAACGTCCAGTTCGACATCGAAGCATTCATAACGATTACGGCCGTTACGCTTGGCCTGGTACAGCGCCATGTCAGCGTAGCCCAGCAGTTGCTCTGCCTGGTCATGGTCATCAGGCGCGATGGCGATGC‑CGATGCTCACCCCCACCGCGAACTGATGCCCTTCGATCCGGAACGGCGGGCTGATGGCCTTGATCAGGCGCTGCGCGGTATTGCACGCGGCATCCCGCCCCTCCAGGCCCGTCAGCACCACGGCAAACTCGTC  >  NC_002947/4233708‑4233996

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