Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,755,533	A→C	6.9%	D325A (GAC→GCC)	PP_3319 →	GGDEF domain‑containing sensory box protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,755,533	0	A	C	6.9%	172.7 / 4.8	58	D325A (GAC→GCC)	PP_3319	GGDEF domain‑containing sensory box protein
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base C (2/2);  total (29/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 = 6.06e-01

TGGTCAATGGCGAGATCGTCGGGGTATTCGGCATCGCGCGGGATATCAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCT  >  NC_002947/3755398‑3755679                                                                                                                                        |                                                                                                                                                     tGGTCAATGGCGAGATCGTCGGGGTATTCGGCATCGCGCGGGATATCAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCg‑g                                                                                                                                                 >  3:184493/1‑140 (MQ=255) tGGTCAATGGCGAGATCGTCGGGGTATTCGGCATCGCGCGGGATATCAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCg‑g                                                                                                                                                 <  4:184493/140‑1 (MQ=255)   gTCAATGGCGAGATCGTCGGGGTATTCGGCATCGCGCGGGATATGAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGc                                                                                                                                     <  6:179143/149‑1 (MQ=255)     cAATGGCGAGATCGTCGGGGTATTCGGCATCGCGCGGGATATCAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCgt                                                                                                                                   >  8:31424/1‑149 (MQ=255)      aaTGGCGAGATCGTCGGGGTATTCGGCATCGCGCGGGATATCAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCgt                                                                                                                                   <  7:31424/148‑1 (MQ=255)                tcgtcgGGGTATTCGGCATCGCGCGGGATATCAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCagcagc                                                                                                                        <  1:313088/149‑1 (MQ=255)                   tcgGGGTATTCGGCATCGCGCGGGATATCAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGt                                                                                                                     >  5:305107/1‑149 (MQ=255)                      gggTATTCGGCATCGCGCGGGATATCAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAgg                                                                                                                                      <  6:5699/129‑1 (MQ=255)                      gggTATTCGGCATCGCGCGGGATATCAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAgg                                                                                                                                      >  5:5699/1‑129 (MQ=255)                              ggCATCGCGCGGGATATCAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGt                                                                                                          >  3:128530/1‑149 (MQ=255)                              ggCATCGCGCGGGATATCAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCCGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGt                                                                                                          >  7:137379/1‑149 (MQ=255)                                cATCGCGCGGGATATCAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGtt                                                                                                         >  7:134994/1‑148 (MQ=255)                                 aTCGCGCGGGATATCAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTc                                                                                                        >  2:261634/1‑148 (MQ=255)                                     cgcgGGATATCAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTaaa                                                                                                            <  8:286387/140‑1 (MQ=255)                                     cgcgGGATATCAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTaaa                                                                                                            >  7:286387/1‑140 (MQ=255)                                        gggATATCAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCg                                                                                                <  7:17226/149‑1 (MQ=255)                                           atatCAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGt                                                                                                               >  1:106500/1‑131 (MQ=255)                                           atatCAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGt                                                                                                               <  2:106500/131‑1 (MQ=255)                                              tCAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATg                                                                                          >  5:178851/1‑149 (MQ=255)                                               cAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGAGGTCATg                                                                                          >  4:289998/1‑148 (MQ=255)                                                 gCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCg‑g                                                                                                                                                 <  2:181378/92‑1 (MQ=255)                                                 gCAAACAGAAGCATTACGAGCGCCAGTTCAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCg‑g                                                                                                                                                 >  1:181378/1‑92 (MQ=255)                                                      cAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGa                                                                                  >  1:266579/1‑149 (MQ=255)                                                       aGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGAc                                                                                 <  6:178851/149‑1 (MQ=255)                                                       aGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGAc                                                                                 >  4:1230/1‑149 (MQ=255)                                                       aGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGAc                                                                                 >  1:302900/1‑149 (MQ=255)                                                        gAATCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGAcc                                                                                <  4:256160/149‑1 (MQ=255)                                                         aaGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCt                                                                               <  4:128530/149‑1 (MQ=255)                                                         aaGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCCGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCt                                                                               <  8:137379/149‑1 (MQ=255)                                                             aTTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCtggat                                                                           <  3:75498/149‑1 (MQ=255)                                                                    gcgcAAGTTGAGTATCAATGCAAGCATTGACGGGCTTAGGGGTTTGATTAACGGCCTCTCTCTGGAAGACCGCTATCTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTggatggat                                                                       <  1:261634/148‑1 (MQ=255)                                                                     cgcCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCaa                                                                   <  3:1230/149‑1 (MQ=255)                                                                       ccAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAAt                                                                 <  6:305107/149‑1 (MQ=255)                                                                         aGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATc                                                                >  3:20811/1‑148 (MQ=255)                                                                             gAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGt                                                                                         <  4:52700/119‑1 (MQ=255)                                                                             gAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGt                                                                                         >  3:52700/1‑119 (MQ=255)                                                                                    aaTGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTc                                                    >  7:159798/1‑149 (MQ=255)                                                                                     aTGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATcc                                                               >  5:187618/1‑137 (MQ=255)                                                                                     aTGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATcc                                                               <  6:187618/137‑1 (MQ=255)                                                                                         cagccaTGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCg                                               >  7:283461/1‑149 (MQ=255)                                                                                         cagccaTGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCg                                               >  1:141124/1‑149 (MQ=255)                                                                                          gttgaTGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCgg                                              <  3:231779/145‑1 (MQ=255)                                                                                              aTGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCac                                          <  8:283461/149‑1 (MQ=255)                                                                                                 aCGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCactac                                       >  5:101811/1‑149 (MQ=255)                                                                                                   gggCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACtt                                     <  2:141124/149‑1 (MQ=255)                                                                                                         aCGGGTTTGCTTAACCGCCTGTCTCTGGAAGCCCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAAt                                                                 >  7:166131/1‑115 (MQ=255)                                                                                                         aCGGGTTTGCTTAACCGCCTGTCTCTGGAAGCCCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAAt                                                                 >  7:170650/1‑115 (MQ=255)                                                                                                         aCGGGTTTGCTTAACCGCCTGTCTCTGGAAGCCCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAAt                                                                 <  8:170650/115‑1 (MQ=255)                                                                                                         aCGGGTTTGCTTAACCGCCTGTCTCTGGAAGCCCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAAt                                                                 <  8:166131/115‑1 (MQ=255)                                                                                                         aCGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGc                                 <  6:101811/147‑1 (MQ=255)                                                                                                          cGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGAt                              <  8:134994/149‑1 (MQ=255)                                                                                                              attGCTTACCCGCCGGTCTGTGGAAGACCG‑GGT‑TGTCCTGGCTTGTCGCAGTAGCGGTCAAGTTAAAGTTCATGTGGCGGTCTTGTACATCGACCTGGTTGGATACAAATCCTTCAACGACTCCATCGGGCACTACTTCGGCGATCaa                           <  6:212311/147‑1 (MQ=255)                                                                                                                 gCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTg                        <  2:266579/148‑1 (MQ=255)                                                                                                                 gCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGACCAGGCGTGACGCAGCAGCCGTCAACGTAAAGTTCAACTGGCCGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAgag                        >  5:212311/1‑146 (MQ=255)                                                                                                                   ttAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGCGCTg                     >  1:205621/1‑149 (MQ=255)                                                                                                                                    gAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCa                                           <  5:259183/110‑1 (MQ=255)                                                                                                                                    gAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCa                                           >  6:259183/1‑110 (MQ=255)                                                                                                                                      aGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCt  >  5:198599/1‑149 (MQ=255)                                                                                                                                        |                                                                                                                                                     TGGTCAATGGCGAGATCGTCGGGGTATTCGGCATCGCGCGGGATATCAGCAAACAGAAGCATTACGAGCGCCAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCG‑GCT‑TGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCT  >  NC_002947/3755398‑3755679

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