Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,559,812	G→A	7.5%	G96G (GGC→GGT)	PP_3143 ←	LysR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,559,812	0	G	A	7.5%	154.4 / 2.3	53	G96G (GGC→GGT)	PP_3143	LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (25/24);  new base A (2/2);  total (27/26)
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.47e-01

AAGGTCGAAGTCGCCCGGCGTCAGCGCCTGGGCACTGTCCTCAGCCACCAGGCGCAACTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTC  >  NC_002947/3559665‑3559958                                                                                                                                                         |                                                                                                                                                   aaGGTCGAAGTCGCCCGGCGTCAGCGCCTGGGCACTGTCCTCAGCCACCAGGCGCAACTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGc                                                                                                                                                   >  7:190801/1‑149 (MQ=255)  aGGTCGAAGTCGCCCGGCGTCAGCGCCTGGGCACTGTCCTCAGCCACCAGGCGCAACTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGcc                                                                                                                                                  <  2:269947/149‑1 (MQ=255)       gAAGTCGCCCGGCGTCAGCGCCTGGGCACTGTCCTCAGCCACCAGGCGCAACTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCgg                                                                                                                                              <  5:135871/148‑1 (MQ=255)              cccGGCGTCAGCGCCTGGGCACTGTCCTCAGCCACCAGGCGCAACTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTcccc                                                                                                                                      >  5:4065/1‑149 (MQ=255)                            cTGGGCACTGTCCTCAGCCACCAGGCGCAACTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGctg                                                                                                                          >  3:223295/1‑145 (MQ=255)                                         tCAGCCACCAGGCGCAACTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTc                                                                                                           >  3:102781/1‑149 (MQ=255)                                         tCAGCCACCAGGCGCAACTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTc                                                                                                           >  6:237794/1‑149 (MQ=255)                                           aGCCACCAGGCGCAACTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAGGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGa                                                                                                         <  5:335071/149‑1 (MQ=255)                                           aGCCACCAGGCGCAACTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTc                                                                                                           <  1:264419/147‑1 (MQ=255)                                            gCCACCAGGCGCAACTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGAt                                                                                                        <  1:197804/149‑1 (MQ=255)                                               accaGGCGCAACTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTACCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTc                                                                                                                 <  5:13021/137‑1 (MQ=255)                                               accaGGCGCAACTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTACCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTc                                                                                                                 <  5:17946/137‑1 (MQ=255)                                               accaGGCGCAACTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTACCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTc                                                                                                                 >  6:13021/1‑137 (MQ=255)                                               accaGGCGCAACTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTACCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTc                                                                                                                 >  6:17946/1‑137 (MQ=255)                                                 caGGCGCAACTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGc                                                                                                   >  4:202474/1‑149 (MQ=255)                                                      gcAACTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAg                                                                                              >  4:113999/1‑149 (MQ=255)                                                        aaCTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTa                                                                                               <  2:62388/146‑1 (MQ=255)                                                        aaCTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTa                                                                                               >  1:62388/1‑146 (MQ=255)                                                         aCTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGcc                                                                                           >  7:116257/1‑149 (MQ=255)                                                          cTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGcat                                                                                          >  2:275340/1‑147 (MQ=255)                                                           tCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTc                                                                                         >  6:16086/1‑149 (MQ=255)                                                            cGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCg                                                                                        >  8:38734/1‑149 (MQ=255)                                                             gATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCgg                                                                                       <  6:199045/149‑1 (MQ=255)                                                             gATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCgg                                                                                       <  4:102781/149‑1 (MQ=255)                                                                 tCCGGGTATCGTTGTTGAATACGCGGCAAGCGCGGCGTGAGCCATTTGGCGTGCACTAAGTAGTCGGGTTTCATCGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATAGTGCCGTAGGCCTCg                                                                                        <  1:284139/149‑1 (MQ=255)                                                                   cGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCg                                                                                 >  8:124827/1‑149 (MQ=255)                                                                     ggTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTc                                                                               >  5:135904/1‑149 (MQ=255)                                                                           gttgttGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGc                                                                          >  6:219443/1‑148 (MQ=255)                                                                            ttgttgAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGc                                                                          >  5:354292/1‑147 (MQ=255)                                                                            ttgttgAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCa                                                                        <  8:116257/149‑1 (MQ=255)                                                                             tgttgAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCt                                                                                                                  >  8:174326/1‑106 (MQ=255)                                                                             tgttgAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCt                                                                                                                  <  7:174326/106‑1 (MQ=255)                                                                               ttgAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGc                                                                     <  5:219443/149‑1 (MQ=255)                                                                                    aaCGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGc                                                                >  5:286934/1‑149 (MQ=255)                                                                                     aCGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCACCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGTCGTAGGCCTCGGCCACCGTCAGGGCCAGGCGTTGcc                                                               >  4:55417/1‑149 (MQ=255)                                                                                      cgcgGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGccc                                                              <  3:125326/149‑1 (MQ=255)                                                                                          gCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCTGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCg                                                          <  4:223295/149‑1 (MQ=255)                                                                                            aaGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGgc                                                        <  6:4065/149‑1 (MQ=255)                                                                                                gcggcgTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTc                                                    <  8:190801/149‑1 (MQ=255)                                                                                                   gcgTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGc                                                 <  5:237794/149‑1 (MQ=255)                                                                                                     gTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTc                                               <  6:135904/149‑1 (MQ=255)                                                                                                           cATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCAcgcc                                         <  3:68979/149‑1 (MQ=255)                                                                                                                       aGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCTGATTTCCCCCGCGTCCAGCCgca                                                                                                                        <  7:145678/61‑1 (MQ=255)                                                                                                                       aGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCgca                                                                                                                        >  8:145678/1‑61 (MQ=255)                                                                                                                                acgGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTCGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAg                          >  8:165548/1‑148 (MQ=255)                                                                                                                                   ggagtgAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAAc                      <  8:141467/144‑1 (MQ=255)                                                                                                                                         aaCGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCgg                >  4:74455/1‑149 (MQ=255)                                                                                                                                           cGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTCCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTa              <  7:38734/149‑1 (MQ=255)                                                                                                                                              acacGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGGCCCTCGGCGCTCAGCTCCACGCCGCGCGccc                                 >  3:307013/1‑127 (MQ=255)                                                                                                                                              acacGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGGCCCTCGGCGCTCAGCTCCACGCCGCGCGccc                                 <  4:307013/127‑1 (MQ=255)                                                                                                                                              acacGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGcc           >  3:205113/1‑149 (MQ=255)                                                                                                                                                  gCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAgg       >  3:105985/1‑149 (MQ=255)                                                                                                                                                     aTGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGgcgc    <  3:202474/149‑1 (MQ=255)                                                                                                                                                        tgCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGctc  <  3:113999/148‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                   AAGGTCGAAGTCGCCCGGCGTCAGCGCCTGGGCACTGTCCTCAGCCACCAGGCGCAACTCGATGTCCGGGTAACGTTGTTGAAAACGCGGCAAGCGCGGCGTGAGCCATTTGGC‑‑‑CAGGAA‑‑CGACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTC  >  NC_002947/3559665‑3559958

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