Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	4,916,414	C→T	6.7%	G54S (GGC→AGC)	cheW ←	Chemotaxis protein CheW

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	4,916,414	0	C	T	6.7%	196.2 / 3.4	60	G54S (GGC→AGC)	cheW	Chemotaxis protein CheW
Reads supporting (aligned to +/- strand):  ref base C (28/28);  new base T (2/2);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01

CGACCACTTGCTTGTCGGCTTCGATGATGACGATACGGGTGTTGTCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACTCGTTGTCCAGACGGAAGGTAACCCACTGCAGGATCGGATCTTCGGA  >  ACB143/4916288‑4916549                                                                                                                               |                                                                                                                                        cGACCACTTGCTTGTCGGCTTCGATGATGACGATACGGGTGTTGTCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACg                                                                                                                                  >  1:192312/1‑134 (MQ=255) cGACCACTTGCTTGTCGGCTTCGATGATGACGATACGGGTGTTGTCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACg                                                                                                                                  >  5:8748/1‑134 (MQ=255) cGACCACTTGCTTGTCGGCTTCGATGATGACGATACGGGTGTTGTCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACg                                                                                                                                  >  4:47595/1‑134 (MQ=255)      aCTTGCTTGTCGGCTTCGATGATGACGATACGGGTGTTGTCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCg                                                                                                                           >  4:91108/1‑136 (MQ=255)      aCTTGCTTGTCGGCTTCGATGATGACGATACGGGTGTTGTCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCg                                                                                                                           >  4:38894/1‑136 (MQ=255)      aCTTGCTTGTCGGCTTCGATGATGACGATACGGGTGTTGTCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCg                                                                                                                           >  4:34601/1‑136 (MQ=255)      aCTTGCTTGTCGGCTTCGATGATGACGATACGGGTGTTGTCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCg                                                                                                                           <  4:314254/136‑1 (MQ=255)          gcttgTCGGCTTCGATGATGACGATACGGGTGTTGTCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGgc                                                                                                                         <  1:146430/134‑1 (MQ=255)           cttgTCGGCTTCGATGATGACGATACGGGTGTTGTCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGgcgc                                                                                                                       <  6:274459/135‑1 (MQ=255)           cttgTCGGCTTCGATGATGACGATACGGGTGTTGTCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGgcgc                                                                                                                       <  1:300154/135‑1 (MQ=255)                cGGCTTCGATGATGACGATACGGGTGTTGTCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCgg                                                                                                                 >  8:224602/1‑136 (MQ=255)                 ggCTTCGATGATGACGATACGGGTGTTGTCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGt                                                                                                                <  3:193635/136‑1 (MQ=255)                    ttCGATGATGACGATACGGGTGTTGTCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTAcc                                                                                                             >  8:23431/1‑136 (MQ=255)                       gatgatGACGATACGGGTGTTGTCTGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCgg                                                                                                           <  5:14907/135‑1 (MQ=255)                       gatgatGACGATACGGGTGTTGTCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCgg                                                                                                           <  7:152639/135‑1 (MQ=255)                           atgaCGATACGGGTGTTGTCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCaa                                                                                                      <  5:34984/136‑1 (MQ=255)                                gATACGGGTGTTGTCGGTCACTTCGGTCGGCATCAGTCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCTGCGCACCCGGTACCGGGTCAATCTCg                                                                                                 <  1:13227/136‑1 (MQ=255)                                gATACGGGTGTTGTCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCg                                                                                                 <  6:165178/136‑1 (MQ=255)                                    cGGGTGGTGTCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGgtgt                                                                                             <  7:216838/136‑1 (MQ=255)                                        tgttgtCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAgcg                                                                                         <  2:100307/136‑1 (MQ=255)                                             tCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGa                                                                                    >  7:73996/1‑136 (MQ=255)                                             tCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGa                                                                                    >  7:225468/1‑136 (MQ=255)                                             tCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGCTGTAGCGCAGGa                                                                                    >  4:65789/1‑136 (MQ=255)                                                    cgtcggtcggCATCAGGCCAAAACGCTGACGCGTGGCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTcctg                                                                             <  7:172133/134‑1 (MQ=255)                                                    cTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTcctg                                                                             <  7:191785/136‑1 (MQ=255)                                                    cTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTcctg                                                                             <  4:121693/136‑1 (MQ=255)                                                          tcggCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTcctgcacctg                                                                       <  3:204862/136‑1 (MQ=255)                                                          tcggCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTcctgcacctg                                                                       <  1:90770/136‑1 (MQ=255)                                                            ggCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTcctgcacctgca                                                                     >  2:264011/1‑136 (MQ=255)                                                            ggCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTcctgcacctgca                                                                     >  8:193312/1‑136 (MQ=255)                                                              cATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATc                                                                   <  7:101017/136‑1 (MQ=255)                                                                  aGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCTCGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCTGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACAt                                                               <  5:90417/136‑1 (MQ=255)                                                                  aGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACAt                                                               >  6:176305/1‑136 (MQ=255)                                                                  aGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACAt                                                               <  3:10291/136‑1 (MQ=255)                                                                                aCGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGAc                                                 <  6:233693/136‑1 (MQ=255)                                                                                aCGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGAc                                                 <  5:42235/136‑1 (MQ=255)                                                                                aCGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGAc                                                 <  6:127936/136‑1 (MQ=255)                                                                                 cgcgTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACt                                                >  6:258188/1‑136 (MQ=255)                                                                                 cgcgTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACt                                                >  2:126004/1‑136 (MQ=255)                                                                                 cgcgTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACt                                                >  1:115628/1‑136 (MQ=255)                                                                                 cgcgTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACt                                                >  2:93447/1‑136 (MQ=255)                                                                                 cgcgTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACt                                                >  2:105720/1‑136 (MQ=255)                                                                                     tgtCGATCACCGTCACCACGTTGGCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACTCGtt                                            <  3:114249/136‑1 (MQ=255)                                                                                        cGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCTCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCg                                                                                                 >  1:247030/1‑80 (MQ=255)                                                                                        cGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCTCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCg                                                                                                 >  1:251522/1‑80 (MQ=255)                                                                                        cGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCTCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCg                                                                                                 <  2:247030/80‑1 (MQ=255)                                                                                        cGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCTCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCg                                                                                                 <  2:251522/80‑1 (MQ=255)                                                                                        cGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACTCGTTGTc                                         >  7:290257/1‑136 (MQ=255)                                                                                             aCCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACTCGTTGTCCAGAc                                    >  6:202945/1‑136 (MQ=255)                                                                                             aCCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACTCGTTGTCCAGAc                                    >  7:84744/1‑136 (MQ=255)                                                                                                            ccGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACTCGTTGTCCAGCCGGAAGGTAACCCACt                     >  1:298920/1‑136 (MQ=255)                                                                                                            ccGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACTCGTTGTCCAGACGGAAGGTAACCCAc                      <  5:308635/135‑1 (MQ=255)                                                                                                              gcgcAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACTCGTTGTCCAGACGGAAGGTAACCCACTg                    >  3:90327/1‑135 (MQ=255)                                                                                                                  aaGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACTCGTTGTCCAGACGGAAGGTAACCCACTGCAgga               >  2:240120/1‑136 (MQ=255)                                                                                                                  aaGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACTCGTTGTCCAGACGGAAGGTAACCCACTGCAgga               >  7:296639/1‑136 (MQ=255)                                                                                                                  aaGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACTCGTTGTCCAGACGGAAGGTAACCCACTGCAgga               >  5:264030/1‑136 (MQ=255)                                                                                                                  aaGTTGATGACGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACTCGTTGTCCAGACGGAAGGTACCCCACTGCAgga               >  3:228244/1‑136 (MQ=255)                                                                                                                   aGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACTCGTTGTCCAGACGGAAGGTAACCCACTGCAggat              <  7:123582/136‑1 (MQ=255)                                                                                                                     ttGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACTCGTTGTCCAGACGGAAGGTAACCCACTGCAggatcg            <  3:6192/136‑1 (MQ=255)                                                                                                                           atgCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACTCGTTGTCCAGACGGAAGGTAACCCACTGCAGGATCGGATCt       <  1:310193/135‑1 (MQ=255)                                                                                                                               ccAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACTCGTTGTCCAGACGGAAGGTAACCCACTGCAGGATCGGATCTTCGGa  <  1:142421/136‑1 (MQ=255)                                                                                                                               |                                                                                                                                        CGACCACTTGCTTGTCGGCTTCGATGATGACGATACGGGTGTTGTCGGTCACTTCGGTCGGCATCAGGCCAAAACGCTGACGCGTGTCGATCACCGTCACCACGTTGCCGCGCAAGTTGATGATGCCCAGCACGTAGCTCGGCGCACCCGGTACCGGGGCAATCTCGGTGTAGCGCAGGACTTCCTGCACCTGCATCACATTGATGCCGTAGGACTCGTTGTCCAGACGGAAGGTAACCCACTGCAGGATCGGATCTTCGGA  >  ACB143/4916288‑4916549

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