Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	718,669	0	C	T	43.9%	52.4 / 15.8	41	E598K (GAA→AAA)	speF	ornithine decarboxylase isozyme, inducible
Reads supporting (aligned to +/- strand):  ref base C (4/19);  new base T (18/0);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 5.49e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGG  >  NC_000913/718523‑718812                                                                                                                                                   |                                                                                                                                                cGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGa                                                                                                                                             <  2:247194/151‑1 (MQ=255)   gCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGaca                                                                                                                                           <  1:154746/151‑1 (MQ=255)     aGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGAcaca                                                                                                                                         <  2:430593/151‑1 (MQ=255)     aGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAAAAAACATGCATTTTCTGaca                                                                                                                                           >  2:272946/1‑149 (MQ=255)      gACGAACCAGTTCCACTTCACCGCGTAAATGGTCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTGCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAGCATGCATTTCCTGacac                                                                                                                                          <  1:587199/149‑1 (MQ=255)       aCGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAAc                                                                                                                                       <  2:556145/151‑1 (MQ=255)             cAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTTTTGACACAAACGACGc                                                                                                                                 >  2:510681/1‑151 (MQ=255)               gTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAATCATTCATTTTCTGACACAACTGTCGc                                                                                                                                 >  2:645405/1‑149 (MQ=255)                  ccACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAACATTCATTTTCTTACACAAATGACGCAgggg                                                                                                                            >  2:642772/1‑150 (MQ=255)                  ccACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTTCTGACACAAATGACGCAgggg                                                                                                                            >  2:518496/1‑150 (MQ=255)                   cACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAACATGCATTTTCTTACACCACTGACGCAggg                                                                                                                             >  1:200291/1‑149 (MQ=255)                    aCTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGta                                                                                                                           <  1:365876/150‑1 (MQ=255)                        cACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAATATGCATTTTCTGTCACAAATGACGCAGGGGATAAcc                                                                                                                      >  1:594632/1‑151 (MQ=255)                         aCCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAAATGCATTTTCTGACACAACTGACGCAGGGGATaa                                                                                                                        >  2:622947/1‑148 (MQ=255)                            gcgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCg                                                                                                                  <  1:359052/151‑1 (MQ=255)                                 aaTAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAACATGCATTTTCTTTCACAACTGACGCAGGGGATAAACCGCGGAg                                                                                                               >  1:604997/1‑149 (MQ=255)                                  aTAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATTCATTTTCTGACACAAATGACGCAGGGGATAACCAGCGGAgcg                                                                                                             >  1:498364/1‑150 (MQ=255)                                    aGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGct                                                                                                           <  1:195737/150‑1 (MQ=255)                                                      gCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGAt                                                                                        <  1:276964/151‑1 (MQ=255)                                                       cGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATg                                                                                       <  1:284722/151‑1 (MQ=255)                                                              aTGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTTACACAAATGACGCAGGGGATAAACCGCGGAGCGCTCTTCATGCTGTTTGTAGATGGAAgg                                                                                  >  2:538259/1‑149 (MQ=255)                                                                   gACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTTACACAAATGACGCAGGGGATAAACCGCGGAGCGCTCTTTATGCTGTTTTTAGATGGGAgg                                                                                  >  1:228588/1‑144 (MQ=255)                                                                      gcgTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCt                                                                        <  2:523623/151‑1 (MQ=255)                                                                      gcgTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACAAGCATTTTCTGACACAACTGAAGCAGGGTATAAACGGGGTAGCGCTCCTCATGCTGTTTGTTGTTGGAAGGGAGGAcct                                                                          >  1:19177/1‑147 (MQ=255)                                                                        gTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGc                                                                      <  1:430826/151‑1 (MQ=255)                                                                         tGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTTCTGACACAAATGACGCAGGGGATAACCGGCGTAGCGCTCTTCCTGCTGTTTTTAGATGGGAGGGAGCACTTCt                                                                        >  2:574568/1‑148 (MQ=255)                                                                          gggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTa                                                                    <  1:98136/151‑1 (MQ=255)                                                                          gggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAAATGCATTTTCTGACACAAATTACGCCGGGGATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGGAGCACTtcttct                                                                     >  2:491112/1‑150 (MQ=255)                                                                                tGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTCCTGACACAACTGACGCAGGGGATAAACCGCGGAGCGCTCTTCATGCTGTTTGTAGATGGAAGGGAGCACTTCTGCTAATGgcg                                                              >  2:537704/1‑151 (MQ=255)                                                                                    ccTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCAAGCATTTTCTGACACAAATGACGCAGGGGATAAACGGCGGAGCGCTCTTCATGCTGTTTGTTGATGGAAGGGAGGACTTCTGCTAATgg                                                                >  2:536388/1‑145 (MQ=255)                                                                                       tACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCg                                                         <  2:476209/149‑1 (MQ=255)                                                                                          ggAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTc                                                    <  1:225319/151‑1 (MQ=255)                                                                                            aaCATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAACATATCTTTTCTTTCCCAACTCACGCAGGGGATAAACGGCGGAGCGCCCTTCATTCTGTTTTTAGATGGGAGGGAGCACTTCTGCttaagg                                                                >  1:452984/1‑132 (MQ=255)                                                                                              cATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGc                                                <  2:342282/151‑1 (MQ=255)                                                                                               aTCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAATATGCATTTTCTGACACAACAGACGCAGGGGATAAACGGCGTAGCGCTCTTTCTGCTGTTTTTAGTAGGAAGGGAGCACATTTGCTAATGGGGCGGCGGGCTCaa                                                  >  1:367733/1‑148 (MQ=255)                                                                                                       ttGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTc                                        <  2:19183/150‑1 (MQ=255)                                                                                                            gttgttTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCg                                  <  2:246755/151‑1 (MQ=255)                                                                                                            gttgttTCACGTTGTGGCGGGCATACAAAACATGCATTTCCTGACACAACTGACGCAGGGGATAACCCGCGGAGCGCTCTTCCTTCTGTTTGTAGATGGGAGGGAGCACTTCTGCTAATGGCGCGTCGGGCTCAAGCAGTTTTTTGAAGCg                                  >  2:16907/1‑151 (MQ=255)                                                                                                              tgttTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGt                                 <  2:228606/150‑1 (MQ=255)                                                                                                                                      aaaTCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGa         >  1:57375/1‑150 (MQ=255)                                                                                                                                            tGCATTTCGTGACACAAATGACGCAGGGTATAATCGGCGTAGCGCTATTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACTAGCAGGGCAACAAGTTGCTGAAGTTTgg  >  1:188154/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                CGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGG  >  NC_000913/718523‑718812

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