Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	447,148	Δ1 bp	7.0%	intergenic (+28/‑274)	PP_0368 → / → PP_0369	acyl‑CoA dehydrogenase/diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	447,148	0	G	.	7.0%	181.4 / 6.9	57	intergenic (+28/‑274)	PP_0368/PP_0369	acyl‑CoA dehydrogenase/diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (28/25);  new base . (2/2);  total (30/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 = 8.44e-01

CGCGCTTCTACTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACC  >  NC_002947/447005‑447294                                                                                                                                                |                                                                                                                                                    cgcgcTTCTACTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGccc                                                                                                                                                  <  1:279475/147‑1 (MQ=255)           cTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGc                                                                                                                                     >  7:92506/1‑149 (MQ=255)               cAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCAAAAACCCGCCT‑CCCGGTCAGGGCAGGCGGGt                                                                                                                                 <  3:40/149‑1 (MQ=255)                    cATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGttttt                                                                                                                             >  7:42324/1‑148 (MQ=255)                    cATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGttttt                                                                                                                             >  6:79156/1‑148 (MQ=255)                        cTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCgg                                                                                                                                   >  3:54231/1‑138 (MQ=255)                        cTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCgg                                                                                                                                   <  4:54231/138‑1 (MQ=255)                        cTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACg                                                                                                                         <  8:42324/148‑1 (MQ=255)                         tGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTc                                                                                                                       >  4:159480/1‑149 (MQ=255)                                  cccgtacccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGttttt                                                                                                                             <  2:235846/134‑1 (MQ=255)                                  cccgtacccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGttttt                                                                                                                             >  1:235846/1‑134 (MQ=255)                                   ccgtacccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCa                                                                                                             <  5:79156/149‑1 (MQ=255)                                      tacccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGttttt                                                                                                                             >  5:53853/1‑130 (MQ=255)                                      tacccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGttttt                                                                                                                             <  6:53853/130‑1 (MQ=255)                                                  gACCATGCTGTGGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCg                                                                                                <  5:222754/147‑1 (MQ=255)                                                  gACCATGCTGTAGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCg                                                                                                >  6:222754/1‑147 (MQ=255)                                                   aCCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAg                                                                                             <  5:322088/149‑1 (MQ=255)                                                   aCCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAg                                                                                             <  1:209332/149‑1 (MQ=255)                                                     cATGCTGTCGGGCGCCAACAATCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCCAAACCCGCCTGccc                                                                                                                                                  >  2:88692/1‑95 (MQ=255)                                                     cATGCTGTCGGGCGCCAACAATCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCCAAACCCGCCTGccc                                                                                                                                                  <  1:88692/95‑1 (MQ=255)                                                     cATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGa                                                                                           <  8:92506/149‑1 (MQ=255)                                                        gCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCgg                                                                                                                                   >  7:292808/1‑106 (MQ=255)                                                        gCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCgg                                                                                                                                   <  8:292808/106‑1 (MQ=255)                                                         cTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCgg                                                                                       <  5:33506/149‑1 (MQ=255)                                                            tCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCggg                                                                                      <  2:322497/147‑1 (MQ=255)                                                              gggCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTa                                                                                  >  3:226173/1‑149 (MQ=255)                                                               ggCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGt                                                                                                      >  6:188945/1‑128 (MQ=255)                                                               ggCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGt                                                                                                      <  5:188945/128‑1 (MQ=255)                                                                    caacaaCCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCcgc                                                                             >  1:271837/1‑148 (MQ=255)                                                                                   gATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCctg                                                             >  3:222586/1‑149 (MQ=255)                                                                                      ggACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAaga                                                                        >  4:64006/1‑135 (MQ=255)                                                                                      ggACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAaga                                                                        <  3:64006/135‑1 (MQ=255)                                                                                       gACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTgg                                                         <  4:222586/149‑1 (MQ=255)                                                                                          gAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCa                                                      <  4:226173/149‑1 (MQ=255)                                                                                                   ttCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCAt                                              >  3:156544/1‑148 (MQ=255)                                                                                                               tACTGATACTGGCGCTGTACCCATAACCCGCTTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCAcc                                  <  6:250592/148‑1 (MQ=255)                                                                                                                aCTGATACTGGCGCTGTACCCAAAACCCGCCAGCGCTGACCGAG‑AGGCGGGTTTTATACGTCTTCAGGTTCAGCGTGGACCCTCGGAGGAGCGGGGTTAACCGCGAAGAGTCCGGCCCTCCTGGACAATCTCAATAATCCCTTCACCg                                 >  5:250592/1‑148 (MQ=255)                                                                                                                  tGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAggccg                                                                   >  5:91159/1‑112 (MQ=255)                                                                                                                  tGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAggccg                                                                   <  6:91159/112‑1 (MQ=255)                                                                                                                    aTACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGTCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTga                             >  3:240109/1‑147 (MQ=255)                                                                                                                            gcTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCCTCAGCCGTTCACCGTTGCCATCGCTAc                    >  4:272876/1‑149 (MQ=255)                                                                                                                                  cccAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCAcc                                  <  2:22348/129‑1 (MQ=255)                                                                                                                                  cccAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCAcc                                  >  1:22348/1‑129 (MQ=255)                                                                                                                                  cccAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTAtt              <  4:156544/149‑1 (MQ=255)                                                                                                                                  cccAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTAt               >  1:340234/1‑148 (MQ=255)                                                                                                                                    cAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTcg            >  2:80796/1‑149 (MQ=255)                                                                                                                                     aaaaCCCGCCT‑CCCGGTCAGGGCAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCAACAGCCGTTCACCGTTGCCAGCGCTACCGTATTcgc           >  2:23538/1‑149 (MQ=25)                                                                                                                                      aaaCCCGCCT‑CCCGGTCAGGGCAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCa                           >  6:286985/1‑132 (MQ=25)                                                                                                                                      aaaCCCGCCT‑CCCGGTCAGGGCAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCa                           <  5:286985/132‑1 (MQ=25)                                                                                                                                         cccGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCggg        >  3:96305/1‑148 (MQ=255)                                                                                                                                          ccGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCctgc                                                            <  7:241084/95‑1 (MQ=255)                                                                                                                                          ccGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCctgc                                                            >  8:241084/1‑95 (MQ=255)                                                                                                                                          ccGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGt       >  4:56154/1‑148 (MQ=255)                                                                                                                                          ccGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGt       >  1:268304/1‑148 (MQ=255)                                                                                                                                             ccTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAAc   <  4:96305/149‑1 (MQ=255)                                                                                                                                              cTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAAcc  >  2:84937/1‑149 (MQ=255)                                                                                                                                              cTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAAcc  <  7:32636/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                    CGCGCTTCTACTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGG‑AGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACC  >  NC_002947/447005‑447294

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