Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	814,298	G→A	8.1%	P29S (CCC→TCC)	PP_0700 ←	transmembrane sensor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	814,298	0	G	A	8.1%	134.5 / 2.6	49	P29S (CCC→TCC)	PP_0700	transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base G (26/19);  new base A (2/2);  total (28/21)
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.42e-01

TTGCCAGAACGTCTCCAGCTCGGCGTAGGCGT‑GGGCGTTCTCTTGTACCTGGCACCAGGCTGCAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCATTGGCAGAGGATGCTAATTATATTGAGAACCATTTTCAAGTACGTGGCGACAGGGGGATAAAAAA  >  NC_002947/814155‑814446                                                                                                                                                 |                                                                                                                                                     ttGCCAGAACGTCTCCAGCTCGGCGTAGGCGT‑GGGCGTTCTCTTGTACCTGGCACCAGGCTGCAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCa                                                                                                                                                 <  6:248971/149‑1 (MQ=255)    ccAGAACGTCTCCAGCTCGGCGTAGGCGTGGGGCGTTCTCTTGTACCTGGCACCAGGCTGCAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAgc                                                                                                                                               >  1:236432/1‑149 (MQ=255)      aGAACGTCTCCAGCTCGGCGTAGGCGT‑GGGCGTTCTCTTGTACCTGGCACCAGGCTGCAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAgctgc                                                                                                                                            >  2:62243/1‑149 (MQ=255)      aGAACGTCTCCAGCTCGGCGTAGGCGT‑GGGCGTTCTCTTGTACCTGGCACCAGGCTGCAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGACAGCGCTTCCTGTTGTGGATCGGGCTGCAgctgc                                                                                                                                            >  7:258660/1‑149 (MQ=255)              tcCAGCTCGGCGTAGGCGT‑GGGCGTTCTCTTGTACCTGGCACCAGGCTGCAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGa                                                                                                                                    >  7:220189/1‑149 (MQ=255)                cAGCTCGGCGTAGGCGT‑GGGCGTTCTCTTGTACCTGGCACCAGGCTGCAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACa                                                                                                                                  >  3:292515/1‑149 (MQ=255)                          tAGGCGT‑GGGCGTTCTCTTGTACCTGGCACCAGGCTGCAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAgc                                                                                                                                               <  1:36589/126‑1 (MQ=255)                          tAGGCGT‑GGGCGTTCTCTTGTACCTGGCACCAGGCTGCAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAgc                                                                                                                                               >  2:36589/1‑126 (MQ=255)                                   ggCGTTCTCTTGTACCTGGCACCAGGCTGCAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCgg                                                                                                                >  2:48321/1‑149 (MQ=255)                                   ggCGTTCTCTTGTACCTGGCACCAGGCTGCAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCg                                                                                                                 <  1:62243/148‑1 (MQ=255)                                            ttGTACCTGGCACCAGGCTGCAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCAc                                                                                                       >  8:221481/1‑149 (MQ=255)                                                 ccTGGCACCAGGCTGCAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATa                                                                                                   <  6:294104/148‑1 (MQ=255)                                                 ccTGGCACCAGGCTGCAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAg                                                                                                  <  5:106578/149‑1 (MQ=255)                                                  cTGGCACCAGGCTGCAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCCCGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAgg                                                                                                 >  7:21986/1‑149 (MQ=255)                                                  cTGGCACCAGGCTGCAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAgt                                                                                                 >  7:307482/1‑148 (MQ=255)                                                     gCACCAGGCTGCAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAgggg                                                                                               <  1:48321/148‑1 (MQ=255)                                                       aCCAGGCTGCAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCg                                                                                            >  1:209866/1‑149 (MQ=255)                                                        ccAGGCTGCAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGt                                                                                           >  8:77689/1‑149 (MQ=255)                                                                cAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCt                                                                                    >  2:181164/1‑148 (MQ=255)                                                                cAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCt                                                                                    >  8:72634/1‑148 (MQ=255)                                                                 aaaTGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGTGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTa                                                                                   >  4:352653/1‑148 (MQ=255)                                                                 aaaTGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTaa                                                                                  >  8:57183/1‑149 (MQ=255)                                                                    tGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGa                                                                               <  6:129121/149‑1 (MQ=255)                                                                      ccTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGAt                                                                              <  7:77689/148‑1 (MQ=255)                                                                          gcgcTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACa                                                                          <  4:322604/148‑1 (MQ=255)                                                                               gggTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGTGGCGGTCTGCTCACGATAGTGGCCGTTTGGCCTAACGATGACAGGAGt                                                                     >  4:31255/1‑148 (MQ=255)                                                                                  tttCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCAtt                                                                 >  8:58497/1‑149 (MQ=255)                                                                                     cgtcgCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCATTgg                                                               <  7:72634/148‑1 (MQ=255)                                                                                       tcgCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGa                                                                               >  2:339063/1‑130 (MQ=255)                                                                                       tcgCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGa                                                                               <  1:339063/130‑1 (MQ=255)                                                                                       tcgCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCATTgg                                                               >  3:331562/1‑146 (MQ=255)                                                                                            aCCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGTTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCATTGGCAGAGGAt                                                       <  4:292515/149‑1 (MQ=255)                                                                                              cGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTaa                                                                                  >  4:57442/1‑120 (MQ=255)                                                                                                       gcagcaACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCATTGGCAGAGGATGCTAATTatat                                            <  7:58497/149‑1 (MQ=255)                                                                                                          gcaACGAAATCCACGCCAGCGCTGCCTGTGGTGGAGCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCATTGGCAGAGGATGCTAATTATATTga                                         <  5:169363/149‑1 (MQ=255)                                                                                                          gcaACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCATTGGCAGAGGATGCTAATTATATTga                                         <  2:209866/149‑1 (MQ=255)                                                                                                           caACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCATTGGGAGAGGATGCTAATTATATTgag                                        <  7:57183/149‑1 (MQ=255)                                                                                                             aCGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCATTGGCAGAGGATGCTAATTATATTGAGaa                                      >  5:297621/1‑149 (MQ=255)                                                                                                                     cACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCATTGGCAGAGGATGCTAATTATATTGAGAACCAtttt                               >  4:49719/1‑148 (MQ=255)                                                                                                                           aGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCATTGGCAGAGGATGCTAATTATATTGAGAACCATTTTCAAGTAc                        >  2:16918/1‑149 (MQ=255)                                                                                                                              gcTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCATTGGCAGAGGATGCTAATTATATTGAGAACCATTTTCAAGTACGTg                     >  1:141757/1‑149 (MQ=255)                                                                                                                               cTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGTGGCCGTTTGGCCTAACGATGACAGGAGTCATTGGCAGAGGATGCTAATTATATTGAGAACCATTTTCAAGTACGTgg                    <  1:181164/149‑1 (MQ=255)                                                                                                                                    tgttgtGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCATTGGCAGAGGATGCTAAGTATATTGAGAACCATTTTCAAGTACGTGGCGa                 >  1:243536/1‑147 (MQ=255)                                                                                                                                        gtgGATCGGACTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCATTGGCAGAGGATGCTAATTATATTGAGAACCATTTTCAAGTACGTGGCGACAgggg           >  5:104165/1‑149 (MQ=255)                                                                                                                                        gtgGATCGGACTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCATTGGCAGAGGATGCTAATTATATTGAGAACCATTTTCAAGTACGTGGCGACAgggg           >  5:109423/1‑149 (MQ=255)                                                                                                                                         tgGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCATTGGCAGAGGATGCTAATTATATTGAGAACCATTTTCAAGTACGTGGCGACAggggg          <  2:126182/149‑1 (MQ=255)                                                                                                                                         tgGATCGGACTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCATTGGCAGAGGATGCTAATTATATTGAGAACCATTTTCAAGTACGTGGCGACAgggg           <  6:104165/148‑1 (MQ=255)                                                                                                                                         tgGATCGGACTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCATTGGCAGAGGATGCTAATTATATTGAGAACCATTTTCAAGTACGTGGCGACAgggg           <  6:109423/148‑1 (MQ=255)                                                                                                                                            aTCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCATTGGCAGAGGATGCTAATTATATTGAGAACCATTTTCAAGTACGTGGCGACAGGGGGATa       <  2:141757/149‑1 (MQ=255)                                                                                                                                                 gCTGCATCGGCTGTGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCATTGGCAGAGGATGCTAATTATATTTAGAACCATTTTCAAGTACGTGGCGACAGGGGGATaaaaaa  <  8:258660/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     TTGCCAGAACGTCTCCAGCTCGGCGTAGGCGT‑GGGCGTTCTCTTGTACCTGGCACCAGGCTGCAAATGCCTGGCGCTGGGTTTCGTCGCAACCGGGCTGGCGCAGCAACGAAAACCACGCCAGCGCTTCCTGTTGTGGATCGGGCTGCAGCTGCTTGGCGGACATCGGGCTCATGGGGGCGGTCTGCTCACGATAGGGGCCGTTTGGCCTAACGATGACAGGAGTCATTGGCAGAGGATGCTAATTATATTGAGAACCATTTTCAAGTACGTGGCGACAGGGGGATAAAAAA  >  NC_002947/814155‑814446

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑