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

CGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTG  >  NC_000913/718509‑718804                                                                                                                                                       |                                                                                                                                                   cGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAt                                                                                                                                                    <  1:1010236/151‑1 (MQ=255)          gCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCA‑ttt                                                                                                                                          >  1:64551/1‑151 (MQ=255)            aTCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGCTACAAAACATGCC‑ttt                                                                                                                                          >  2:547886/1‑149 (MQ=255)            aTCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCA‑tttt                                                                                                                                         >  1:12841/1‑149 (MQ=255)              ccGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCA‑TTTCCt                                                                                                                                       <  1:344484/150‑1 (MQ=255)                ggCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAACATACA‑TTTcc                                                                                                                                        >  1:406190/1‑147 (MQ=255)                   aGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCA‑TTTCCTGAcaca                                                                                                                                 <  2:243196/151‑1 (MQ=255)                      cGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCGAACAAAAAATGCA‑TTTTCTGTCACaa                                                                                                                                >  2:409551/1‑149 (MQ=255)                               tCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAgg                                                                                                                      <  1:418156/150‑1 (MQ=255)                                     tCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAAAATGCA‑TTTTCTCACACAAATGACGCCGGGGATaa                                                                                                                >  1:65020/1‑150 (MQ=255)                                     tCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGtata                                                                                                                 <  1:98876/149‑1 (MQ=255)                                      cACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAAc                                                                                                               <  2:119921/150‑1 (MQ=255)                                       aCCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCA‑TTTTCTTACACAACTGACGCAGGGGATaa                                                                                                                >  1:509557/1‑148 (MQ=255)                                        ccGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCgg                                                                                                            <  1:627511/151‑1 (MQ=255)                                          gcgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCg                                                                                                          <  1:971959/151‑1 (MQ=255)                                                   ggCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGctct                                                                                                 <  2:938955/151‑1 (MQ=255)                                                   ggCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAAAAAACATGCA‑TTTTCTTACACAACTGACGCAGGGTATAAACGGCGTAgcgccc                                                                                                  >  1:16079/1‑148 (MQ=255)                                                      gTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAAAATGGA‑TTTTCTTTAACAAATGAGGCAGGGGATAACCGGCGGAGCGctc                                                                                                  >  1:67687/1‑147 (MQ=255)                                                           tGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTAAAAAAAATGCA‑TTTCCTTGCACAACCGACGCAGGGGGTAAACGGCGTCGCGCTCTTCATGc                                                                                           >  2:308844/1‑149 (MQ=255)                                                            ggCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGt                                                                                        <  2:941968/151‑1 (MQ=255)                                                              cttcttGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAAAAAATATTCA‑TTTTCTGTCACAAACGACGCCGGGGATAAACGGCGGAGCGCCCTTCATGCTGttt                                                                                      >  2:1011587/1‑151 (MQ=255)                                                                   tGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTa                                                                                   <  2:64549/149‑1 (MQ=255)                                                                      ggATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATAAAAAAAATGCA‑TTTTCTGGCACCAAAGACGCAGGGGATAAACGCCGGAGCGCTCTTCCTTCTGttttttgaagg                                                                              >  2:657229/1‑143 (MQ=255)                                                                      ggATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATgg                                                                              <  2:382053/151‑1 (MQ=255)                                                                        aTTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGaa                                                                            <  1:954231/151‑1 (MQ=255)                                                                            aTGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTa                                                                        <  2:113293/151‑1 (MQ=255)                                                                                tGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGCTACAAATCATGCA‑TTTTCTCACACAAATGACGCAGGGGGGAAACGGCCTAGGGCTCTTCATTCTGTTTTTAGATGGAAGGGAGCa                                                                     >  2:170765/1‑150 (MQ=255)                                                                                tGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGAAAAAAAAAATGCA‑TTTTCTTACACAAATTACGCAGGGGGTAACAAGGGGAGCGGTCTTTCTGGTGTTTGTTGATGGAggg                                                                          >  2:549310/1‑145 (MQ=255)                                                                                       tGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAATCATGCA‑ATTTCTGACACAAAAGACGCAGGGGAAAAACGGCGTGGGGCCCTTCCTGCTGTTTTTAGATGGAAGGGAGCACTTCTg                                                               >  1:179417/1‑149 (MQ=255)                                                                                       tGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAAAATTCA‑TTTCTTGTCACAACTGACGCAGGGGGTAACCGGCGTAGGGGTCTTCTTGCTGTTTTTAGATGGGAGGGAGCACCTc                                                                 >  2:352131/1‑147 (MQ=255)                                                                                         ggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATAAAAAACATTCA‑TTTTCTTACACAACCGAAGCAGGGGATAAACCGCGTAGCGGTCTTCATGCTGTTTGTTGAAGGGAGGGAGGACtt                                                                  >  1:252655/1‑144 (MQ=255)                                                                                           aaGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATg                                                         <  2:196601/151‑1 (MQ=255)                                                                                            aGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGTA‑TTTTCTGACACAACAGACGCAGGGGATAAACAGCGTAGCGCGCTTCCTGCTGTTTTTTGATTGGAGGGAGCACtt                                                                  >  2:803200/1‑141 (MQ=255)                                                                                             gTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCA‑TTTTCTGACACAACACACGCACGGGAAAAAAGGCGTAGCGCTCTTCATGCTGTTTTTTGATGGAAGGGAGCACTTCTGCt                                                             >  1:306206/1‑145 (MQ=255)                                                                                                       cGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCGGAAAAAAAATTCC‑ATTTCTGGCAAAAACAAAGCAGGGGATAAACCGCCCAGCGCCCCCTATTCTGTTTGTAGGTGGAAGGGAGGACTTTTGccgatagggcgccgg                                                >  2:83353/1‑134 (MQ=255)                                                                                                       cGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACt                                             <  1:128232/151‑1 (MQ=255)                                                                                                       cGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATTCA‑TTTTCTTTCACAACCGACGCAGGGGGTAAACCGCGGAGGGCTCTTCATGCTGTTTTTTGATTGGAGGGAGCACTTCTTCTAATTGGGCGTCgggc                                              >  2:315982/1‑148 (MQ=255)                                                                                                         gAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAATATGCA‑TTTTTTTTCACAACTTACGCCGGGGGTAACCGGCGCAGGGCTCTTCATTCTGTTTTTTGATGGAAGGGAGCACTTCTGCTTATGGCGCGGCGGGCt                                             >  2:996787/1‑149 (MQ=255)                                                                                                          aaCATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCaa                                          <  1:921707/151‑1 (MQ=255)                                                                                                                tctTTTTGCAGTTGTCTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTTGCCCCAATTG‑CGCAGGGGATAAACGGCGTAGCGCTCCTCATTCTGTTTGTTGATGGAAGGGAACACTTCTGCTAATTGCGCGGCGGGCTCAAACAGt                                     >  1:92721/1‑150 (MQ=255)                                                                                                                     ttGCAGTTGTTTCACGTTGTGGCGGGGCTACAAAACATTCC‑TTTTCTGACACAAACAACGCCGGGGATAAACAGCGGAGCGCCCTTTCTGCTGTTTTTTGATGGAAGGGAGCACTTCTTCTTATTGCGCGGCGGGCCCAAACAGgtttt                                 >  1:901357/1‑149 (MQ=255)                                                                                                                      tGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAAATGCA‑TTTCCTTAAACAAATGACGCAGGGGATAAACCGCGGAGCGCTCTTCCTGCTGTTTGTAGATGGAAGGGAGCACTTCTGCTAATGGCGCGGCGGGCTCAAGCAGgtttt                                 >  1:386324/1‑148 (MQ=255)                                                                                                                         aGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCA‑TTTCCTTACACAACTTACGCAGGGGATAAACGGCGGAGCGCTCTTCATGCTGTTTGTAGATGGGAGGGAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGGTTTTTGaa                             >  1:1002698/1‑149 (MQ=255)                                                                                                                             gttTCACGTTGTGGCGGGCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTa                        <  2:262346/150‑1 (MQ=255)                                                                                                                             gttTCACGTTGTGGCGGGCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTa                        <  2:818372/150‑1 (MQ=255)                                                                                                                                 cACGTTGTGGCGGGCATACAAAAAATGCA‑TTTTCTGACACAACTGACGCCGGGGATAAACCGCGGAGCGCCCTTCCTTCTTTTTGTTGATGGAAGGGAGAACTTCTGCTAATGGCGCGTCGGGCTCAAACAGATTTTTTAAGCGTAc                       >  2:848860/1‑147 (MQ=255)                                                                                                                                       gtgGCGGGCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGc              <  2:518082/150‑1 (MQ=255)                                                                                                                                         ggCGGGCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGcaac           <  1:285167/151‑1 (MQ=255)                                                                                                                                           cGGGCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGcaacaa         <  1:111810/151‑1 (MQ=255)                                                                                                                                            gggCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTTCAGCAGTTTTTCGAAGCGTACTAGCAGGGcaaca          >  1:206385/1‑149 (MQ=255)                                                                                                                                            gggCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAg        <  2:10807/151‑1 (MQ=255)                                                                                                                                                 tACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCt   >  2:203348/1‑151 (MQ=255)                                                                                                                                                  aCAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTg  <  2:892808/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                   CGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCA‑TTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTG  >  NC_000913/718509‑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: ‑