Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	718,669	0	C	T	50.0%	20.2 / 15.2	28	E598K (GAA→AAA)	speF	ornithine decarboxylase isozyme, inducible
Reads supporting (aligned to +/- strand):  ref base C (3/11);  new base T (14/0);  total (17/11)
Fisher's exact test for biased strand distribution p-value = 3.39e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCC  >  NC_000913/718526‑718814                                                                                                                                                |                                                                                                                                                  cAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAACAAGCATTTTCTCacac                                                                                                                                            >  1:199784/1‑151 (MQ=255)      gAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTg                                                                                                                                       <  2:718027/151‑1 (MQ=255)       aaCCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCCTTTTCTTTCACaa                                                                                                                                          >  2:297017/1‑147 (MQ=255)         ccAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTAAAAAATATTTATTTTCTTACACaa                                                                                                                                          >  1:71367/1‑145 (MQ=255)          cAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAAATGCATTTTCTTACACAACTGACg                                                                                                                                    >  1:278524/1‑150 (MQ=255)                      aCCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATGCAAAACATTCATTTTTTGCCACACAAGACGGAgggggtaa                                                                                                                          >  1:143294/1‑143 (MQ=255)                             aaaTAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAg                                                                                                                 <  1:43159/150‑1 (MQ=255)                                     gTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTTGCAAAACCGACGCAGGGGAAAAACGGCGGAGCGCGCtt                                                                                                          >  1:52426/1‑149 (MQ=255)                                       aGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCAt                                                                                                       <  1:479731/150‑1 (MQ=255)                                             cttcttGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACAAGCATTTTCTGACACAAATGACGCAGGGGATAAACCGCGGAGCGCTCTTCATGCTGttt                                                                                                >  1:555101/1‑151 (MQ=255)                                                cttGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTg                                                                                               <  2:406579/149‑1 (MQ=255)                                                 ttGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGGTAAAAAAAATGTATTTTCTGACACAAAAGAAGCCGGGGAGAAACCGCGGCGCGCGCTTCTTGCTGtttt                                                                                               >  2:228632/1‑147 (MQ=255)                                                    cGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATg                                                                                         <  1:26530/151‑1 (MQ=255)                                                    cGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAAATTCATTTTCTTACCCAACTGACGCCGGGGATAAACGGCGGAGGGctccttcttcttgtttt                                                                                              >  2:710039/1‑132 (MQ=255)                                                      gATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGa                                                                                       <  1:743957/151‑1 (MQ=255)                                                          cATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAACATGCATTTTCTTACACAACTTACGCCGGGGGTAAACCGCGTAAAGCTCTTCATGCTGTTTTTTGTTGGAAg                                                                                     >  1:696303/1‑149 (MQ=255)                                                                   gcgTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCt                                                                          <  2:120701/151‑1 (MQ=255)                                                                     gTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCGACAAAACAATCCTTTTCTGTCACCACTGACGCCGGGGATAAACCGCGGAGGGCCCCTTTTTCTGTTTTTAGATGggagggagacctttt                                                                          >  2:522662/1‑142 (MQ=255)                                                                       gggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCt                                                                       <  1:50830/150‑1 (MQ=255)                                                                       gggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAAATGCATTTTCTGACACAACAGACGCAGGGGATAACCGGCGGAGCGCTCTTCATGCTGTTTTTTGATGGAAGGGAGCACTtcttct                                                                       >  2:105033/1‑150 (MQ=255)                                                                               cTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTTACACCAATGACGCAGGGGATAAACCGCGGAGCGCTCTTCCTTCTGTTTTTAGATGGGAGGGAGCACTTTTGCTTATGGcgcg                                                              >  1:402886/1‑151 (MQ=255)                                                                                      cGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACt                                                       <  2:246191/151‑1 (MQ=255)                                                                                                ctTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAAAAAAAATGCTTTTTCTGGCAAAACATACGCAGGGGATAAACGGCGAAGCGCTCTTCATGCTGTTTGTTGAGGGGAGGGGGCAAATCTTCTAATAGCGCGTCGGGCTCAAGCAGttt                                             >  2:138540/1‑151 (MQ=255)                                                                                                            gttTCACGTTGTGGCGGGGATAAAAATAATGCCTTTCTTTGCACAAATCAAGCAGGGGGTAAACGGCGGAGCGCTCCTCATGTTGTTTGTAGATGGAAGGTAGCACTTCTTCTAATGGGGCGTGGGGGCCAAGCCGgttttttatgcgt                                   >  2:311807/1‑141 (MQ=255)                                                                                                                            ggCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGt                 >  1:427677/1‑151 (MQ=255)                                                                                                                                    aaTCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGa           >  1:257311/1‑149 (MQ=255)                                                                                                                                       cATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGtt       <  2:745575/150‑1 (MQ=255)                                                                                                                                           cATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGcc  <  1:574418/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                  CAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCC  >  NC_000913/718526‑718814

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