Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	535,920	0	C	T	71.4%	129.1 / ‑4.3	52	R71C (CGC→TGC)	hyi	hydroxypyruvate isomerase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (4/31);  minor base G (14/0);  total (20/32)
Fisher's exact test for biased strand distribution p-value = 4.53e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCA  >  NC_000913/535775‑536065                                                                                                                                                  |                                                                                                                                                  gaAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGt                                                                                                                                              <  2:627983/151‑1 (MQ=255)    aGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTAtt                                                                                                                                           <  1:503482/151‑1 (MQ=255)    aGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAGTGGGCGGGGGGGGGGGGCgg                                                                                                                                               >  1:361535/1‑147 (MQ=255)            aGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTAt                                                                                                                                   <  2:259934/151‑1 (MQ=255)                  gTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTgg                                                                                                                             <  1:57900/151‑1 (MQ=255)                       cgcgGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGCGGCGGGGGGGGGCGGGGTTTGCTGTATTCCttggcggg                                                                                                                        >  1:54137/1‑144 (MQ=255)                        gcggcgTTGAATTTATGTTTCCTTATGACTACGACAGTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGt                                                                                                                         <  2:619086/149‑1 (MQ=255)                         cggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTgaa                                                                                                                      <  1:394181/151‑1 (MQ=255)                            cgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGGGGGGGGGGGGGGATTTCCTGTATTTCTGGcc                                                                                                                           >  1:634645/1‑143 (MQ=255)                                 aaTTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGtt                                                                                                              <  1:112994/151‑1 (MQ=255)                                 aaTTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGCGGGGGGGGGGGGGGGGATTTCCTTTATTCCTGGCCGGGTAGAagag                                                                                                                >  1:709828/1‑149 (MQ=255)                                        gTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGAt                                                                                                       <  2:241153/151‑1 (MQ=255)                                           tCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATgg                                                                                                     <  2:174675/150‑1 (MQ=255)                                             cTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCGGGGGGGCGGGGGAGGGGGGGGTTGGCTGTATTTCTGGCCGTGAAGAAGAGATTTGGGGTgggg                                                                                                   >  2:87458/1‑148 (MQ=255)                                              ttATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGCGGGGGGGGGCGGGATTGCCTGTATTCCTGGCCGTGAAGAAGAGGTTTGGGGTggg                                                                                                    >  2:314964/1‑147 (MQ=255)                                               tATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGt                                                                                                  <  2:255466/149‑1 (MQ=255)                                                aTGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTag                                                                                                <  1:51889/150‑1 (MQ=255)                                                      aCGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTagcagcag                                                                                          <  1:671147/150‑1 (MQ=255)                                                       cGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCg                                                                                        <  2:95598/151‑1 (MQ=255)                                                         aCATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGGGGGGGGGGCGGGATTGCCTGTATTTCTTGCCCGGGAGAAGAGGTTTGGGGTGGCGGAGCAGCAGCGAt                                                                                      >  2:644936/1‑151 (MQ=255)                                                                agaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTAt                                                                               <  1:411511/151‑1 (MQ=255)                                                                 gaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATg                                                                              <  1:850146/151‑1 (MQ=255)                                                                    ttAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGTGCGGTATTGTCTGTATTTCTGGCCGGGGAGAAGAGGTTCGGGGTGGCGGAGCAGCAGCGATTCGTTTTGccc                                                                           >  2:683037/1‑151 (MQ=255)                                                                     tAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGccc                                                                           <  1:614834/150‑1 (MQ=255)                                                                         aCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGCGGGGGGGTGCGGTATTGCCTGTATTCCTGGCCGGGAAGAAGAGTTTCGGGATGGCGTAGCAGCCGCGAATCGTTATGCCCGGgcg                                                                      >  2:126082/1‑151 (MQ=255)                                                                           aTGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCt                                                                    <  1:839123/151‑1 (MQ=255)                                                                             gtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTgg                                                                  <  2:127425/151‑1 (MQ=255)                                                                              tgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGCTGGGGGGGGGGGGGGGGGGGGATTTGCTGTTTTTCTGGGCCGGGAGAAGAGGTTCGGGGGGGGGGAGCAGGAGCGGTTCGTTTTGCCCGTGCGGTggg                                                                 >  1:695003/1‑151 (MQ=255)                                                                              tgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGATTGGGGGCGGGGGGGGGGGGGATTGACTTTATTTCTTGCGCGGGAGAAGAGGTTTGGGGTGGCGGAGGAGCAACGATTTGTTATTGCCGt                                                                         >  1:402172/1‑143 (MQ=255)                                                                                   gCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGtaata                                                            <  2:574257/151‑1 (MQ=255)                                                                                      aGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGGGGGGGGGGCGGGATTGGCTGTATTCCTGGCCGTGAAGAAGAGTTTTGGGGTGGGGTAGCAGCAGCGATTCGTTATTTCCGTGCGCTGGGGAATaaa                                                          >  2:414195/1‑150 (MQ=255)                                                                                       gTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATaaaaa                                                        <  2:683916/151‑1 (MQ=255)                                                                                          ataaACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAAt                                                     <  2:695046/151‑1 (MQ=255)                                                                                            aaaCTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTa                                                   <  1:723508/151‑1 (MQ=255)                                                                                             aaCTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGGGGGGGAGTGCGGGATTGCCTGTATTCCTTGCCGTGGAGAAGAGGTTCGGGGTGGCGGAGCAGCAGCGATTCGTTTTTCCCGTGCGCTGGGGAATAAAAAAAATa                                                   >  1:522929/1‑150 (MQ=255)                                                                                             aaCTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTaa                                                  <  1:700768/151‑1 (MQ=255)                                                                                             aaCTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGGGGGGGAGAGGGGGGTTGCCTGTTTTCCTGGCCGGGGAGAAGAGGTTTGGGGTGGGGTAGGCGCAGCGGTTCGTTATGCCCGGGCGCTGGGGaaaaaaaaaaat                                                    >  1:772673/1‑147 (MQ=255)                                                                                               cTCGAACATACGCTGCACAATTTACCGGCGGGTGGCTGGGCGGCGGGGGGGGGCGGGATTGCCTGTATTCCTGGCCGTGAAGAAGAGGTTCGGGGTGGCGGAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGGAATAAAAAAATTaa                                                  >  2:207310/1‑149 (MQ=255)                                                                                               cTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGCGGCGGGGGGGTGGGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGGTTTGGGGTGGGGTAGCAGCAGCGAATCGTTATGGCCGTGCGCTGGGGAAAAAAAAAATTAACt                                                >  2:493108/1‑151 (MQ=255)                                                                                                   aaCATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCt                                            <  2:634685/151‑1 (MQ=255)                                                                                                    aCATACGCTGCACAATTTACCGGCGGGGGGCTTGGGGGCGGGGGGGAGGGGGATTGACTTTTTTTCTGGCCGTGAAGAAGAGGTTTGGGGTGGCGGAGCAGCAGCGATTCGTTATTGCCGTGCGCTGGGTAATAAAAAAAATAAATGTCt                                            >  1:320657/1‑150 (MQ=255)                                                                                                    aCATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGGGCGGTATTTCCTGTATTTCCGGCCGTGAAGAAGAGGTTTGGGGTGGCGGAGCAGCAGCGGTTCGTTTTGCCCGGGCGCTGGGGAATAAAAAAATTAACCTTCTg                                           >  2:565358/1‑151 (MQ=255)                                                                                                        aCGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTc                                        <  2:233324/150‑1 (MQ=255)                                                                                                             gCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTaa                                   <  1:252496/150‑1 (MQ=255)                                                                                                               acaATTTACCGGCGGGTGACTGGGCGGGGGGGGGGGGGGGGATTGCCTTTATTTCTGGGCGGGGAGGAGAGGTTCGGGGGGGGGGGGGAGCAGCGGTTCGTTTTTCCCGGGCGCTGGGGAATAAAAAAATTAAATGTCTGGGCGGTaaaa                                 >  1:618457/1‑150 (MQ=255)                                                                                                                 aaTTTACCGGCGGGGGACTGGGGGGCGGGGGGGGGGGGGATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTTGGGGTGGCGGAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAAATAACTGTCTGGTCGGTaaaa                                 >  2:467685/1‑148 (MQ=255)                                                                                                                      aCCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCt                         <  1:81424/151‑1 (MQ=255)                                                                                                                           cggGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGttt                    <  1:506435/151‑1 (MQ=255)                                                                                                                            gggTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGttt                    <  2:136605/150‑1 (MQ=255)                                                                                                                             ggTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTca                  <  1:632057/151‑1 (MQ=255)                                                                                                                                           gggggAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGAtt    <  1:123404/151‑1 (MQ=255)                                                                                                                                             gggAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCa  <  2:709871/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  GAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCA  >  NC_000913/535775‑536065

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