Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	718,659	0	T	A	44.0%	67.6 / 22.2	51	D601V (GAT→GTT)	speF	ornithine decarboxylase isozyme, inducible
Reads supporting (aligned to +/- strand):  ref base T (5/23);  new base A (22/0);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 9.87e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTG  >  NC_000913/718510‑718804                                                                                                                                                      |                                                                                                                                                    ggCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATc                                                                                                                                                    <  2:42550/151‑1 (MQ=255)   ccTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCAt                                                                                                                                                  <  2:627928/151‑1 (MQ=255)      tCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACa                                                                                                                                                   >  1:120686/1‑147 (MQ=255)          cATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTTCAAATACTGC‑C‑ttt                                                                                                                                          >  2:421400/1‑150 (MQ=255)           aTCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTcc                                                                                                                                        <  1:410163/151‑1 (MQ=255)             ccGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTg                                                                                                                                      <  1:406634/151‑1 (MQ=255)                      gAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGC‑A‑TTTCCTTACACAACTg                                                                                                                             >  2:221303/1‑151 (MQ=255)                        aCCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATCCATACCATGC‑AGTTTTCTGGCACAAcctt                                                                                                                            >  1:952201/1‑148 (MQ=255)                                 aCTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGtat                                                                                                                  <  2:420254/151‑1 (MQ=255)                                    tCACCGCGTAAATATGCGTAGTTGGCGTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTGACACATCTGACGCAGGGtata                                                                                                                 <  1:610245/149‑1 (MQ=255)                                    tCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAAAAAAAATGC‑A‑TTTTCTTACACAAATGACGCCGGGGATaa                                                                                                                >  2:936365/1‑150 (MQ=255)                                              aaTAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCGGACAAAACATGC‑A‑TTTTCTTACACAACCGACGCAGGGGATAACCgggggag                                                                                                       >  1:69513/1‑144 (MQ=255)                                              aaTAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAgcg                                                                                                     <  1:226084/151‑1 (MQ=255)                                               aTAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGGTAAAAAAAATTC‑A‑TTTTCTTACACAAACAACGCCGGGTGTAAACCGCGTAgcgc                                                                                                    >  1:859779/1‑151 (MQ=255)                                                  ggCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCGTACAAAAAATGC‑A‑TTTTCTTACACAAATGACGCAGGGGATAAACCGCGTAg                                                                                                       >  2:196534/1‑145 (MQ=255)                                                   gCGTAGTTGGCTTCTTGCGGATTCATGCTGATGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCtt                                                                                                <  2:607030/151‑1 (MQ=255)                                                   gCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCtt                                                                                                <  1:724361/151‑1 (MQ=255)                                                     gTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCa                                                                                              <  1:822158/151‑1 (MQ=255)                                                         ttGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGAGAAAAAACATGT‑A‑TTTTCTTACACAACTGTCGCAGGGGGTAACCGGGGTAGCGgtcct                                                                                                >  1:460316/1‑143 (MQ=255)                                                                     ggATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATg                                                                               <  1:942298/150‑1 (MQ=255)                                                                            tGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGC‑A‑TTTTCTGACACAAATGACGCAGGGGATAAACGGCGGAGCGCGCTTCATGCCGTTTGTAGATGGAAgg                                                                          >  1:493062/1‑148 (MQ=255)                                                                             gCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGc                                                                      <  1:1006189/151‑1 (MQ=255)                                                                              cTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGc                                                                      <  2:452254/150‑1 (MQ=255)                                                                               tGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAACATGCAA‑TTTCCT‑ACACCACCGACCCAAGGGATAAACAGCGGAGAGCTCTTTCTGCTGTTTGTTGATGGAAGGTAGGAc                                                                    >  1:374512/1‑151 (MQ=255)                                                                                gACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCAACAAAACATGC‑A‑TTTTCTTACAAAACCTAACCAGGGGATAACCGGCGGAGGGCCCTTCCTCCTGTTTTTTGATGGGAgg                                                                          >  1:927580/1‑144 (MQ=255)                                                                                gACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACt                                                                   <  1:879766/151‑1 (MQ=255)                                                                                gACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTATTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCag                                                                    <  1:436873/150‑2 (MQ=255)                                                                                     gTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGC‑A‑TTTTCTGACACAAATGACGCAGGGGATAAAAGGCGTAGCGCTCTTCATGCTGTTTGTTGGTGGGAGGGAGCACtt                                                                  >  2:496294/1‑147 (MQ=255)                                                                                       gggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTAAAAAATATTC‑A‑TTTTCTTACACAAATGACGCAGGGGATAACAGGCGGAGGGCGCTTCCTGCTGCTTGTAGATGGAAGGGAGCACCTCTggt                                                             >  1:182587/1‑148 (MQ=255)                                                                                          aaGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATg                                                         <  2:874541/151‑1 (MQ=255)                                                                                           aGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATgg                                                        <  1:179576/151‑1 (MQ=255)                                                                                            gTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGgc                                                       <  2:407486/151‑1 (MQ=255)                                                                                             tGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAAAATGC‑A‑TTTTCTTACACAACAAACGCAGGGGATAAACCGCGGAGCGCTCCTCATGCTGTTTTTTGATGGAAGGGAGGACTTCTGCTTATTgcg                                                      >  2:983355/1‑151 (MQ=255)                                                                                                tCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAAAAAATATGC‑A‑TTTTCTTACACAAATGACGCAGGGGATAAACCGCGGAAAGCTCTTCCTGCTGTTTGTAGATGGAAGGGAGGACTTCTGCTTATGgcg                                                      >  2:78514/1‑148 (MQ=255)                                                                                                   ttACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGCTACAAATAATGC‑A‑TTTCCTTACAAAACTGACGCAGGGGATAACCGGCGTAGGGCTCTTCTTGCTGGTTTTAGGTGGGAGGGGGCACTTTTTGTAATTGcgcggc                                                  >  2:258386/1‑147 (MQ=255)                                                                                                    tACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAATATTT‑A‑TTTTCTTTACCAACTGACGCAGGGGATAAACCGCGGAGCGCCCTTTCTTCCTTTTTTTGATGGAAGGGAGGACTTCTTCTAATGgggc                                                     >  1:464193/1‑145 (MQ=255)                                                                                                     aCGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCGACAAAAAATGC‑T‑TTTTCTGACACAACTCACGCCGGGGGTAACCCGCGGAGCGCCCTTCATGCTGTTTTTTGATGGGAGGGAGGACCTCTTCTAATGgcg                                                      >  1:238810/1‑143 (MQ=255)                                                                                                      cGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCGTAAAAAACATTC‑A‑TTTTCTTGCACCAACGACGCAGGGGGTAAACGGCGGAGCGCCCTTCTTGCTGTTTTTAGATGGGAGGGAGAACTTTTGCTAATGgcggggcg                                                 >  1:101271/1‑142 (MQ=255)                                                                                                      cGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACt                                             <  2:1005627/151‑1 (MQ=255)                                                                                                         aaCATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCGTACAAAACATGC‑A‑TTTTCTGACACAAATGACGCCGGGGGTAACCGGCGTAGCGCTCTTCCTTCTGTTTTTTGATGGGAGGGAGGACTTCTGCTAATTGCGCGGCGGACTc                                            >  2:437815/1‑149 (MQ=255)                                                                                                          aCATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATTAAAAACATTC‑A‑TTTTCTTGCACAAATTACGCAGGGGATAAACGGCGTCGCGCCCTTTCTTCTGTTTTTTGATGGGAGGGAGGACTTTTGCTAATTTGGCGTCgggc                                              >  1:743256/1‑144 (MQ=255)                                                                                                          aCATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAg                                         <  2:923050/151‑1 (MQ=255)                                                                                                                      gCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGaa                             <  1:811015/151‑1 (MQ=255)                                                                                                                      gCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGaa                             <  1:682165/151‑1 (MQ=255)                                                                                                                       cAGTTGTTTCACGTTGTGGCGGGCATACAAAACAAGC‑A‑TTTTCTGACACCAAATACGCAGGGGATAAACCGCGGAGCGCCCTTCATGCTGTTTGTAGATGGAAGGGAGCAATTCTGCTAATGGCGCGTCGGGGTCAAGCAGgtttttgta                             >  1:944286/1‑146 (MQ=255)                                                                                                                        aGTTGTTTCACGTTGTGGCGGGCATACAAAACATGC‑A‑TTTTCTTACACAAATGACGCAGGGGATAAACCGCGGAGCGCTCTTCATGCTGTTTTTAGATGGAAGGGAGCACTTCTGCTAATGGCGCGTCGGGCTCAAGCAGTTTTTCGAAGc                           >  2:228966/1‑151 (MQ=255)                                                                                                                            gttTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTa                        <  2:884362/150‑1 (MQ=255)                                                                                                                                     tgtgGCGGGCCTACAAAACATGC‑A‑TTTCCTGACCCAACCGACGCCGGGGATAAACCGCGGAGCGCCCTTCATTCTGTTTTTTTATGGAAGGGAGGACTTCTGCTAATGGCGCGTCGGGCTCCAGCAGTTTTTTGAAGCGTACCAACCggg               >  1:158545/1‑150 (MQ=255)                                                                                                                                           gggCATACAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGcaaca          >  2:922218/1‑149 (MQ=255)                                                                                                                                                 aCAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTg  >  2:35181/1‑151 (MQ=255)                                                                                                                                                  cAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTg  >  2:835053/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                    GGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGC‑A‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTG  >  NC_000913/718510‑718804

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