Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,622,250	0	A	G	42.0%	67.3 / 17.4	51	E161G (GAG→GGG)	yjjJ	putative protein kinase
Reads supporting (aligned to +/- strand):  ref base A (6/23);  new base G (21/0);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 5.11e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCC  >  NC_000913/4622103‑4622388                                                                                                                                                    |                                                                                                                                            cGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCggggaggg                                                                                                                                          <  1:143149/151‑1 (MQ=255)  gCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCggggggggg                                                                                                                                         >  1:393619/1‑151 (MQ=255)      ggCTTTTTAGGGCGCGCGTGTGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTGTGCGCCCTGACCGTATTTAACGGTGAGTATACTGGAGGTTGGTTGGTCGGGGAGGGGAAtt                                                                                                                                     <  1:146836/151‑1 (MQ=255)      ggCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAAtt                                                                                                                                     <  1:903212/151‑1 (MQ=255)         tttttAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGGAtt                                                                                                                                     >  2:110285/1‑148 (MQ=255)           tttAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCa                                                                                                                                 <  1:856599/150‑1 (MQ=255)              aGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGa                                                                                                                             <  2:383793/151‑1 (MQ=255)               gggCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGAt                                                                                                                            <  2:833420/151‑1 (MQ=255)                    cgcgTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGAtt                                                                                                                       <  2:596426/151‑1 (MQ=255)                       gTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGGTGGGGGGGGGGATTATTAGCGATGGGtt                                                                                                                       >  2:635062/1‑148 (MQ=255)                         ggggCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTg                                                                                                                   <  1:232701/150‑1 (MQ=255)                           ggCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACTCCCTGCCCGTATTTAACGGTGAATATACTGGCGGTTTGGTGGTGGGGGTGGGGAGATATTATCGCTGGaaaaccgcc                                                                                                                 >  2:290351/1‑142 (MQ=255)                             cAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGGTGGGGGGGGGGAATATCAGCGGTGGATTAATTcac                                                                                                                >  2:279615/1‑149 (MQ=255)                               ggAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCacaaca                                                                                                            <  2:228357/151‑1 (MQ=255)                                 aaGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGGATTATCAGCGATGGATTACATCACAACAcc                                                                                                          >  1:451430/1‑151 (MQ=255)                                   gTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCt                                                                                                        <  1:279612/151‑1 (MQ=255)                                     tAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTTGGCGGGGAGGGGGGTTATCAGCGGTGGAGTACTGGACAACAccc                                                                                                         >  1:765191/1‑148 (MQ=255)                                      aGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑g                                                                                                    <  1:10423/151‑1 (MQ=255)                                         cgcgCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGaa                                                                                                 <  2:350721/151‑1 (MQ=255)                                           cgcAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTGGGGGGGGGGGATTATTAGCGATGGGTTATTGCCAAAAACCCCGC‑GGaaa                                                                                                >  1:813345/1‑150 (MQ=255)                                                cTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGATTAATGCACAACACCCCGC‑GGAAATTtctc                                                                                          >  2:217452/1‑151 (MQ=255)                                                  gAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCctc                                                                                          <  2:497942/149‑1 (MQ=255)                                                    aTCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGGTGGATTACTGCACAACACCCCGC‑GGGAATTtctc                                                                                          >  2:578997/1‑147 (MQ=255)                                                     tCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAAATCCTCTTGAt                                                                                     >  2:940483/1‑151 (MQ=255)                                                       tgactgaTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGGTTATCACCGATGGATTACTGGACAACACCCCGC‑GGAAATTCccctcggt                                                                                     >  1:197326/1‑145 (MQ=255)                                                        gactgaTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGAt                                                                                     <  2:807545/148‑1 (MQ=255)                                                             atgatATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACt                                                                             <  2:123754/151‑1 (MQ=255)                                                               gatATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTTGGCGGGGGGGGGAATTATTAGCGATGGATTACCGGACAAAACCCTGC‑GGAAATTCCTCCGGATCaaaaaat                                                                             >  1:520859/1‑147 (MQ=255)                                                                  ataCGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCAccc                                                                        <  2:28617/151‑1 (MQ=255)                                                                     cGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCAtt                                                                     <  1:544895/151‑1 (MQ=255)                                                                          ctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATTAGCGATGGATTACTGCACAACACCCTGC‑GGGAATTCCTCTGGATCAAAAAATCACCCATTACGa                                                                 >  1:617126/1‑150 (MQ=255)                                                                               aGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATTAGCGGTGGGTTAATGGACAAAACCCTGCGGGAAATTCC‑CTGGGACAAAAACTCACCCAATACGaa                                                                >  1:397186/1‑146 (MQ=255)                                                                                    gaagaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTTGGTGGTCGGGGGGGGGGATTATCAGCGATGGGTTACTGCACAAAACCCTGC‑GGGAAATTCTCTGGGTCAAAAAATCACCCCTTACGAACAACTgg                                                         >  1:758413/1‑148 (MQ=255)                                                                                          gtgCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGATTATCAGCGATGGATTTACGCACAACACCCCTC‑GGGAATTTCTCTTGGTTAAAAAACCACCCATTACGAAAAGCTGGcca                                                      >  2:856092/1‑143 (MQ=255)                                                                                              tctACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGGTTAATGCACAAAACCCCGC‑GGGAATTCCTCTGGATCAAAAAATCACCCATTACGAAAAACTGGCAAatgatg                                                >  2:910870/1‑147 (MQ=255)                                                                                               ctACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTTGGGGGGGGGGATTATAAGCGATGGGTTAATGCACAACACCCTGC‑GGGAATTCCTCTGGGTCAAAAACACACCCATTACGAACAGCTGGCAAGTGATGCACTg                                           >  2:66789/1‑151 (MQ=255)                                                                                                  cGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGGATTATCAGCGATGGATTACTGCACAACACCCCGC‑GGAAATTCCTCTGGATCAAAAAATCACCCATTACGAACAGCTGGCAAGTGATGCACTgg                                          >  2:526638/1‑149 (MQ=255)                                                                                                   gCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGATTATCAGCGATGGATTAATGCACAAAACCCCGC‑GGGAAATTCTCTGGATCAAAAAATCACCCATTACGAAAAACCGGGAAGGGATGCAcctg                                          >  1:204780/1‑145 (MQ=255)                                                                                                     ccTGACCGTATTTAACGGTGAATATACTGGCGGTTTGGTGGTCCGGGGGGGGAATTATCAGCGATGGATTTCTGCACAACACCCTGC‑GGAAAATCCTCTGGGTCAAAAACTCACCCAATACCAACAGCTGGCCAGGGATGCACTGGCAgg                                      >  1:488598/1‑150 (MQ=255)                                                                                                         aCCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGaaa                                 <  2:553858/151‑1 (MQ=255)                                                                                                           cGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAACACCCCGC‑GGAAATTCCTCTGGATCAAAAAATCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGaaaat                               >  2:36855/1‑149 (MQ=255)                                                                                                            gTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCCGC‑GGGAAATCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGGGAAAtt                               >  2:192417/1‑150 (MQ=255)                                                                                                              aTTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATtgt                             <  1:742456/150‑1 (MQ=255)                                                                                                              aTTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATtgt                             <  2:945137/150‑1 (MQ=255)                                                                                                              aTTTAACGGTGAATATACTGGCGGCTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGGTTACTGCACAAAACCCTGC‑GGGAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGGGAAATtgt                             >  2:250894/1‑150 (MQ=255)                                                                                                                  aaCGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGt                         <  2:489248/150‑1 (MQ=255)                                                                                                                  aaCGGTGAATATACTGGCGGTTGGGTTGGCGGGGGGGGGGAATATCAGCGATTGATTACTGCCCAACACCCCGC‑GGGAAATCCTCTGGGGCAAAAAACCACCCATTAAGAAAAACCGGCAAGAGGTTCACCGGGAGGAGAAAATGTggg                          >  1:355457/1‑149 (MQ=255)                                                                                                                   aCGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGttc                       <  1:567090/151‑1 (MQ=255)                                                                                                                          atatACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTgcg                <  1:400687/151‑1 (MQ=255)                                                                                                                                 ggCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGgcggcg         >  1:536596/1‑151 (MQ=255)                                                                                                                                    ggttggttGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGagc      <  1:631846/151‑1 (MQ=255)                                                                                                                                     gttggttgGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGagc      >  1:875323/1‑150 (MQ=255)                                                                                                                                        ggttggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGcc  >  1:186691/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                            CGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGC‑GGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCC  >  NC_000913/4622103‑4622388

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