Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,974,960	G→A	A122V (GCC→GTC)	cheA ←	chemotaxis protein CheA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,974,960	0	G	A	98.2%	180.9 / NA	55	A122V (GCC→GTC)	cheA	chemotaxis protein CheA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (27/27);  minor base C (0/1);  total (27/28)
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

AGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGT  >  NC_000913/1974825‑1975094                                                                                                                                        |                                                                                                                                       aGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTg                                                                                                                                    <  2:679015/140‑1 (MQ=255)  ggCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTgg                                                                                                                                   <  2:153794/140‑1 (MQ=255)       gAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGat                                                                                                                              <  2:748474/140‑1 (MQ=255)             aTAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATc                                                                                                                        <  1:542834/140‑1 (MQ=255)                      cgcgGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGGTTCTAATGGCAGATGATGCAAGACCAGGCACATATAATCGAAGCTggc                                                                                                               >  2:240559/1‑140 (MQ=255)                      cgcgGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTggc                                                                                                               >  1:164098/1‑140 (MQ=255)                      cgcgGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTggc                                                                                                               >  2:111299/1‑140 (MQ=255)                       gcggcgACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTggcg                                                                                                              >  1:142708/1‑140 (MQ=255)                       gcggcgACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTggcg                                                                                                              >  1:537556/1‑140 (MQ=255)                                gactgCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGc                                                                                                     >  2:708016/1‑140 (MQ=255)                                       aCTCTGCTCATCTTGAGGTTCAGTTATGGGAACCACATTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCg                                                                                              <  1:407618/140‑1 (MQ=255)                                        ctctGCTCATCTTGCGGTTCACTTTTGGCAACCCCACTTAATCCGGGCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGCCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGa                                                                                             <  2:785682/140‑1 (MQ=255)                                        ctctGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGa                                                                                             <  1:3485/140‑1 (MQ=255)                                        ctctGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGa                                                                                             <  1:708016/140‑1 (MQ=255)                                        ctctGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGa                                                                                             <  2:629184/140‑1 (MQ=255)                                           cgCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTg                                                                                          >  1:610544/2‑140 (MQ=255)                                            gCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAATACCTGGCAGATAGAATGGAAGCTGGAGGCATCCGGCTCTCGCGACTGt                                                                                         >  2:394809/1‑140 (MQ=255)                                                    tGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCTGCAACCCGCTCTTACGACTGTTTATAAGc                                                                                 >  1:345254/1‑140 (MQ=255)                                                     gCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCg                                                                                >  2:560027/1‑140 (MQ=255)                                                     gCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCg                                                                                >  2:688939/1‑140 (MQ=255)                                                     gCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCg                                                                                <  2:216123/140‑1 (MQ=255)                                                     gCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCg                                                                                >  1:216123/1‑140 (MQ=255)                                                          tCACTTTTGGCAACCTCACTTACTCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAg                                                                           <  1:240559/140‑1 (MQ=255)                                                          tCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATTATCGAAGCGGGCAGCAGCCGGCTCTGGCGACTGTTTATAAGCGTCGAg                                                                           >  2:271495/1‑140 (MQ=255)                                                             cTTTTGGCAACCACATTTAATCAGGTCACTGCGGATGTCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTg                                                                        <  1:271495/140‑1 (MQ=255)                                                             cTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGACTTTCGCTTCTAATGCCAGTTGACGCAAGACCCGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTg                                                                        <  2:6335/140‑1 (MQ=255)                                                                  ggCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCtt                                                                   <  2:574341/140‑1 (MQ=255)                                                                   gCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTg                                                                  <  2:610544/140‑1 (MQ=255)                                                                     aaCCACACTTAATCGTGTCACTACGTATGCCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCa                                                                <  1:404283/140‑1 (MQ=255)                                                                     aaCCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGGGCTGTTCTTGCa                                                                <  2:164098/140‑1 (MQ=255)                                                                      aCCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTTCCACTGTTTAAAAGCGGCGAGCTGCTCTTGCat                                                               >  2:404283/1‑140 (MQ=255)                                                                        cacaCTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCatga                                                             >  1:195093/1‑140 (MQ=255)                                                                        cacaCTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCATGCGACTATTTGTATGCGTCGTGCTGTTCTTGCatga                                                             >  1:410278/1‑140 (MQ=255)                                                                                     gTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCg                                                    <  1:506563/136‑1 (MQ=255)                                                                                     gTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCg                                                    >  2:506563/1‑136 (MQ=255)                                                                                     gTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTc                                                >  1:459328/1‑140 (MQ=255)                                                                                     gTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTc                                                <  2:195093/140‑1 (MQ=255)                                                                                     gTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTc                                                >  2:168939/1‑140 (MQ=255)                                                                                     gTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTc                                                >  2:799714/1‑140 (MQ=255)                                                                                       cACTGCGGATGGCGTTCCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCCACTGTTTAAAAGGGTCCAGCTGTTCTTGCATGATGTCCTTCGTTTCCa                                              >  2:362934/1‑140 (MQ=255)                                                                                         cTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCaaa                                            <  1:799714/140‑1 (MQ=255)                                                                                           gCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCaaaaa                                          >  2:842363/1‑140 (MQ=255)                                                                                             ggATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACa                                        <  1:842363/140‑1 (MQ=255)                                                                                               aTGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGa                                      <  1:563461/140‑1 (MQ=255)                                                                                                      ttCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGCCGAGCTGTACTTGCATCATGTCCTTCGTTTCCAAAAACAGATTGATaa                               >  2:309278/1‑140 (MQ=255)                                                                                                        cGCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATg                             >  1:613070/1‑140 (MQ=255)                                                                                                         gCCTTTCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGt                            <  2:587627/140‑1 (MQ=255)                                                                                                            tttCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCgg                         <  2:537556/140‑1 (MQ=255)                                                                                                            tttCGCTTCTAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCgg                         <  2:142708/140‑1 (MQ=255)                                                                                                                    cTAATGCCAGTTGACTCAAGACCTGCCAGTTATAATCAAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGt                 <  1:309278/140‑1 (MQ=255)                                                                                                                     tAATGCCAGTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGtt                >  1:774598/1‑140 (MQ=255)                                                                                                                             gTTGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCa        <  1:561793/140‑1 (MQ=255)                                                                                                                              ttGACGCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCAc       >  2:522220/1‑140 (MQ=255)                                                                                                                              ttGACGCAAGACCTGGCAGATATAAGCGAGGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCAc       <  1:394809/140‑1 (MQ=255)                                                                                                                                   gCAAGACCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGt  >  1:208082/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                       AGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGT  >  NC_000913/1974825‑1975094

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