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

CATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATA  >  NC_000913/718519‑718817                                                                                                                                                        |                                                                                                                                                      cATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTc                                                                                                                                                       <  2:889361/151‑1 (MQ=255)       gCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTATAAATA‑ATGCATTTTCt                                                                                                                                                     >  2:857202/1‑147 (MQ=255)            cGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACt                                                                                                                                            <  1:374561/151‑1 (MQ=255)             gAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAC‑AAACATTTTCTGTCACAACTg                                                                                                                                           >  2:683236/1‑151 (MQ=255)              aaCCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGa                                                                                                                                          <  2:848459/151‑1 (MQ=255)               aCCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATA‑ATCGCTTTTCTTACACaa                                                                                                                                              >  2:592526/1‑146 (MQ=255)                     tCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑gg                                                                                                                                   <  2:634866/150‑1 (MQ=255)                           tCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATaa                                                                                                                             <  1:960905/150‑1 (MQ=255)                                gcgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATAAAAAAC‑ATGAATTTTCTGTCCCAAAAGACGCC‑GGGGATAAACgggg                                                                                                                       >  1:533521/1‑149 (MQ=255)                                  gTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAGCATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTa                                                                                                                     <  1:190627/151‑1 (MQ=255)                                  gTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTa                                                                                                                     <  2:146569/151‑1 (MQ=255)                                  gTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTa                                                                                                                     <  2:795239/151‑1 (MQ=255)                                     aaTAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAC‑ATTCATTTTCTTACACAAATGACGCA‑gggggtaaaccgcg                                                                                                                       >  2:634349/1‑135 (MQ=255)                                        aGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAAC‑ATTCATTTTCTGACACAACTGACGCA‑AGGGATAAACCGCGCAgcgc                                                                                                                 >  2:252362/1‑149 (MQ=255)                                        aGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGctc                                                                                                               <  2:194725/151‑1 (MQ=255)                                        aGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGct                                                                                                                <  2:832461/150‑1 (MQ=255)                                              aGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGCTCTTCAt                                                                                                          <  1:328057/150‑1 (MQ=255)                                                  ggCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGCTCTTCATGCTg                                                                                                      <  2:379701/150‑1 (MQ=255)                                                   gCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAAC‑ATTCAATTTCTTGCACCAAATTCGCA‑CGGGATAACCCGCGTAGCGCCCTTCCTGCtttt                                                                                                    >  2:260953/1‑151 (MQ=255)                                                   gCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGCTCTTCATGCTGtt                                                                                                    <  2:165397/151‑1 (MQ=255)                                                        ttGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAC‑ATGCATTTTCTTACACAACTGACGCC‑GGGGGTAAACGGCGTAGGGCCCTTCCTTCTGttttt                                                                                                 >  2:158001/1‑147 (MQ=255)                                                          gCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATAAAAAACAAAGCTTTTCCTG‑CCCAACACACGCA‑GGGGGTAACCCGCGGAGCGCTCTTTCTGCTGttt                                                                                                   >  1:177111/1‑145 (MQ=255)                                                              aTTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAAAAATA‑ATGCATTTTCTTACACAACTAACGCA‑GGGGGTAACCCGCGTAGGGCTCTTCCTTCTGttttttggtgg                                                                                           >  1:564737/1‑141 (MQ=255)                                                                  aTGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATAAAAAAA‑ATCCCTTTTCTTGCACCACACACGCA‑GGGGGTAAACAGCGGAGCGCCCTTCATTCTGTTTTTTGATGGGAgg                                                                                       >  1:588968/1‑149 (MQ=255)                                                                    gCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAg                                                                                    <  1:764055/150‑1 (MQ=255)                                                                         cgcgTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGACAGACAAATC‑ATGCCTTTTATTACACAACAGACGCT‑GGGGTGAAAAGGCGGAGCGCTCCTTTTCCTGTTTTTTGATGGAAGGGAGCAc                                                                                 >  1:128218/1‑148 (MQ=255)                                                                          gcgTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCt                                                                             <  1:558529/151‑1 (MQ=255)                                                                           cgTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCGTACAAAAA‑ATTCATTTTCCGACACAACTGAAGCA‑GGGGATAAACCGCCGAGCGCCCTTCCTGCTGTTTTTAGAAGGAAGGAGGCAATTCt                                                                             >  2:936760/1‑150 (MQ=255)                                                                           cgTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTTCTGACACAACTGACGCA‑GGGGATAACCGGCGGAGCGCTCTTCATGCTGTTTTTTGATGGGAGGGAGCACTTCt                                                                             >  2:83963/1‑150 (MQ=255)                                                                              gggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTa                                                                         <  2:564748/151‑1 (MQ=255)                                                                               ggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATTCAAAAA‑AATCATTTTCTTAAACAAATCACGCA‑GGGGATAAACGGGGGAGCGCTCTTTATTCTGTTTGTAGGTGGAAGGAAGCACATTTGCTa                                                                         >  1:751983/1‑150 (MQ=255)                                                                               ggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTaa                                                                        <  1:719301/151‑1 (MQ=255)                                                                                 aaGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATg                                                                      <  2:178129/151‑1 (MQ=255)                                                                                   gTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGgc                                                                    <  2:269058/151‑1 (MQ=255)                                                                                    tGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAA‑ATTCCTTTTCTGCCACAAACAACGCC‑GGGGGTAACCCGCGGCGGGCCCTTCCTTCTGTTTTTTGATGGAAGGGAGCACTTCTGCTAATgg                                                                     >  1:904901/1‑149 (MQ=255)                                                                                       tCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAACAATC‑ATGCAATTTCTGTCACCACTTACGCA‑GGGGATAAACCGCGTATCGGTCTTCATGATTTTTTTGGATGGAAGGAAGGCCTTTTGCTAAATGcgcg                                                                 >  2:733345/1‑150 (MQ=255)                                                                                        ccTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGcgcg                                                                 <  1:447992/149‑1 (MQ=255)                                                                                             cGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAAC‑ATGCATTTCCTGACACAACTGACGCA‑GGGGATAACCGGCGGAGCGCTCTTCATGCTGTTTTTAGATGGAAGGGAGCACTTCTGCTAATGGCGcggcgg                                                             >  2:416915/1‑148 (MQ=255)                                                                                             cGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAT‑ATTCATTTTCTGACACCAATGACGCA‑GGGGGGTAACGGCGGAGGGGTCTTCCTTCTGGTTGTTGATGGGAGGGAGGACTTTTTCTAATGgggccgcgg                                                             >  1:11249/1‑143 (MQ=255)                                                                                                   aTCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAC‑AAGCATTTTCTTACACAAAAGAAGCA‑GGGGATAAACAGCGTAGCGCCCCTCCTGCTGTTTGTAGATGGGAGGGAGTACTTTTGCTAATGGCGCGGCGGGCTcca                                                       >  2:247454/1‑146 (MQ=255)                                                                                                     ctctTTTTGCAGTTGTTTCACGTTGTGGCGGGCACACAAAAC‑ATACATTTTTTTACACCAATCAAACAGGGGAATAACC‑GCCTAGAGCTCTTCATTCTGTTTGTGGATGGAAGGTAGAACTTCTGCTCAATGAGCGGCGGGCTCCAGCa                                                    >  2:555426/1‑149 (MQ=255)                                                                                                         ttttGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGttt                                                <  2:840467/149‑1 (MQ=255)                                                                                                              cAGTTGTTTCACGTTGTGGCGGGCATACAAAAA‑AAGCATTTTCTTACACAACAGAAGCA‑GGGGATAAACGGCGTAGCGCTCTTCATGCTGTTTTTTGATGGAAGGGAGGAATTCTGCTAATGGCGCGTCGGGGTCAAACAGgttttt                                             >  2:116170/1‑146 (MQ=255)                                                                                                                   gttTCACGTTGTGGCGGGCATACAAAAC‑AATCCTTTTCTTACACAACTCACGCA‑GGGGGTAACCCGCGGAGCGCTCTTCCTTCTGTTTTTTTATGGAAGGGAGCACCTCTTCTAATGGCGCGTCGGGGTCCAGCAGgttttttaag                                         >  1:377833/1‑141 (MQ=255)                                                                                                                    tttCACGTTGTGGCGGGCATACAAAAC‑ATTCATTTTCTTACACAAATGACGCA‑GGGGATAAACCGCGTAGCGCTCTTCATGCTGTTTTTAGATGGAAGGGAGCACTTCTGCTAATTGCGCGTCGGGCTCCAGCAGGTTTTTGAAGCGTAcc                                   >  2:334144/1‑151 (MQ=255)                                                                                                                       cACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACcagc                                <  2:81758/151‑1 (MQ=255)                                                                                                                                     cATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTg                  >  2:743300/1‑151 (MQ=255)                                                                                                                                        aCAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTg               >  2:962768/1‑151 (MQ=255)                                                                                                                                         cAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTg               >  2:9492/1‑150 (MQ=255)                                                                                                                                                ttGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTg        >  1:798359/2‑151 (MQ=255)                                                                                                                                                  gCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGTACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTgg       <  1:247451/150‑1 (MQ=255)                                                                                                                                                  gCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTgg       <  1:857185/150‑1 (MQ=255)                                                                                                                                                       tCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCata  <  1:698585/150‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                      CATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCA‑GGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATA  >  NC_000913/718519‑718817

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