Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	1,843,613	T→C	N408S (AAC→AGC)	gacS ←	sensor protein GacS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,843,613	0	T	C	98.0%	151.4 / NA	49	N408S (AAC→AGC)	gacS	sensor protein GacS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (24/24);  minor base A (0/1);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCAGGTTG  >  NC_002947/1843467‑1843759                                                                                                                                                   |                                                                                                                                                   gtgtGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTgg                                                                                                                                                 >  2:64205/1‑150 (MQ=255)   gtgCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTg                                                                                                                                               <  2:49555/150‑1 (MQ=255)    tgCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGa                                                                                                                                              <  1:530947/150‑1 (MQ=255)        cacaTCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGAt                                                                                                                                          <  2:411951/150‑1 (MQ=255)        cacaTCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGAt                                                                                                                                          >  2:262200/1‑150 (MQ=255)         acaTCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCCAGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTGGCTGACCAGTCCGGTGAGGATc                                                                                                                                         >  2:253899/1‑150 (MQ=255)                      aCAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCa                                                                                                                            >  1:148108/1‑150 (MQ=255)                        aGGCTGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAAc                                                                                                                          <  1:364229/150‑1 (MQ=255)                        aGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAAc                                                                                                                          >  2:122591/1‑150 (MQ=255)                             gATGCCGGTGTCCTGCACGCTGATGGCCAAGAGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGGACTTGATGGCGTTGATGACGAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTc                                                                                                                     >  1:173761/1‑150 (MQ=255)                                     tgtCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCa                                                                                                             <  1:64205/150‑1 (MQ=255)                                          tGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAg                                                                                                                                 >  1:282340/1‑125 (MQ=255)                                          tGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAg                                                                                                                                 <  2:282340/125‑1 (MQ=255)                                          tGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGaa                                                                                                        >  1:175655/1‑150 (MQ=255)                                          tGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGaa                                                                                                        <  1:253899/150‑1 (MQ=255)                                          tGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGaa                                                                                                        >  2:622753/1‑150 (MQ=255)                                             aCGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCa                                                                                                     >  1:46626/1‑150 (MQ=255)                                                 tGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTAGGTCGCGGTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGgtg                                                                                                 >  1:350543/1‑150 (MQ=255)                                                   gtcttcAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGgtgtg                                                                                               >  1:489876/6‑150 (MQ=255)                                                            tGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAg                                                                                      >  1:530170/1‑150 (MQ=255)                                                                cgcAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCa                                                                                  <  2:148108/150‑1 (MQ=255)                                                                cgcAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCa                                                                                  <  1:451236/150‑1 (MQ=255)                                                                          gTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCaa                                                                        <  1:262200/150‑1 (MQ=255)                                                                            gCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAAtt                                                                      >  1:383883/1‑150 (MQ=255)                                                                              tCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTcc                                                                    >  2:350201/1‑150 (MQ=255)                                                                                        gCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTgg                                                          >  1:115638/1‑150 (MQ=255)                                                                                           tGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCGGGCTCGGGAATTCCAGCTGCTTGGCGt                                                       <  2:384949/150‑1 (MQ=255)                                                                                             gcgcgcACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTgg                                                     <  2:46626/150‑1 (MQ=255)                                                                                                  cACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGggcggc                                                <  1:327449/150‑1 (MQ=255)                                                                                                    cgacgaTGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCg                                              <  1:84212/150‑1 (MQ=255)                                                                                                     gacgaTGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCgg                                             >  1:486291/1‑150 (MQ=255)                                                                                                     gacgaTGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCgg                                             >  2:172376/1‑150 (MQ=255)                                                                                                          tGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGa                                        <  1:449287/150‑1 (MQ=255)                                                                                                           ggTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAg                                       <  1:210148/150‑1 (MQ=255)                                                                                                              gCCTTCGCGGGTGAACTTGATGGCGTTGATGACCAGGCTGGCGAGGACTTGGCGCGGGCGTAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGAt                                    <  2:173761/150‑1 (MQ=255)                                                                                                                catCGCGGGAGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGACGGCCGGAGCGAGGACgg                                  >  2:40593/3‑150 (MQ=255)                                                                                                                cTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATgg                                  >  2:13573/1‑150 (MQ=255)                                                                                                                     cgGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAgg                             <  2:115638/150‑1 (MQ=255)                                                                                                                        gTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAgc                                                                 >  1:365231/1‑111 (MQ=255)                                                                                                                        gTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAgc                                                                 <  2:365231/111‑1 (MQ=255)                                                                                                                           aaCTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTcc                       <  1:122591/150‑1 (MQ=255)                                                                                                                            aCTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCt                      <  2:410538/150‑1 (MQ=255)                                                                                                                            aCTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCt                      >  2:446076/1‑150 (MQ=255)                                                                                                                             cTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTg                     >  2:614794/1‑150 (MQ=255)                                                                                                                              ttGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTgg                    <  1:446076/150‑1 (MQ=255)                                                                                                                               tGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGa                   <  1:172376/150‑1 (MQ=255)                                                                                                                                       ttGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTcc           <  1:614794/150‑1 (MQ=255)                                                                                                                                               cAGGATGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCTGGCTCAGCAATTCCAGCTGCTTGGCGTGTGCGGCTGGAGCGAGGATTGTCAGGGTGTCCTGGATCAGGTCCCGCAGttt   <  1:40593/150‑1 (MQ=255)                                                                                                                                                aGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGCTCAGGTCCCGCAGGTTg  >  2:186403/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCAGGTTG  >  NC_002947/1843467‑1843759

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑