Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,776,263	0	T	G	42.6%	81.4 / 13.0	55	V226G (GTG→GGG)	ydiF	putative acetyl‑CoA:acetoacetyl‑CoA transferase: alpha subunit/beta subunit
Reads supporting (aligned to +/- strand):  ref base T (10/21);  new base G (23/0);  total (34/21)
Fisher's exact test for biased strand distribution p-value = 9.33e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGTGCAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCACCCT  >  NC_000913/1776117‑1776410                                                                                                                                                    |                                                                                                                                                     aGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGTGCAg                                                                                                                                                  <  1:992578/151‑1 (MQ=255)           ccAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGGGCAGAAAATGGtt                                                                                                                                         >  1:193660/1‑150 (MQ=255)           ccAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGATACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGTGCAGAAAATGGTTa                                                                                                                                        <  1:450261/151‑1 (MQ=255)            cAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGTGGAGAAAAAgggt                                                                                                                                         >  2:690994/1‑147 (MQ=255)            cAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGTGCAGAAAATGGTTaa                                                                                                                                       <  2:220535/151‑1 (MQ=255)              gATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGGGCAGAAAATgggt                                                                                                                                         >  1:920451/1‑145 (MQ=255)                 atTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGTGCAGAAAATGGTTAAG‑AAAg                                                                                                                                 <  2:966991/151‑1 (MQ=255)                   tGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGGTG‑CAGGTGCAGAAAATGGTTAAG‑aaa                                                                                                                                  >  1:940455/1‑148 (MQ=255)                    gCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGTGCAGAAAATGGTTAAG‑AAAAcc                                                                                                                               >  1:333789/1‑150 (MQ=255)                        tCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGTGCAGAAAATGGTTAAG‑AAAGCCACGc                                                                                                                           <  1:969357/150‑1 (MQ=255)                             cgcgCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGTGCAGAAAATGGTTAAG‑AAAGCCACGCTGCAt                                                                                                                      <  1:526052/150‑1 (MQ=255)                                  tACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGGGCAGAAAATGGTTAAG‑AAAACCACCCTGCATCCTaa                                                                                                                 >  1:597058/1‑150 (MQ=255)                                   accaccTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGGGCAGAAAAAGGTTAAG‑AAAAACACGCTGCATCCTaaa                                                                                                                >  1:359728/1‑150 (MQ=255)                                    ccaccTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGTGCAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATc                                                                                                              <  1:652117/151‑1 (MQ=255)                                     caccTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGATGATG‑CCGGTGCAGAAAAAGGGTAAG‑AAAACCACGCTGCATCCTAAATCt                                                                                                             >  1:695177/1‑151 (MQ=255)                                      accTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGGGCAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATCt                                                                                                             >  1:261310/1‑150 (MQ=255)                                           cGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGTGCAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATCTGTc                                                                                                          <  2:804758/148‑1 (MQ=255)                                                   aaGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGGTGATG‑CAGGGGCAGAAAATGGGTTAG‑AAAAACACGCTGCATCCTAAATCTGTCCCGATTCCggg                                                                                               >  2:971130/1‑151 (MQ=255)                                                    aGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGGGCAGAAAATGGGTAAG‑AAAACCACGCCGCATCCTAAATCTGTCCGTATTCCggg                                                                                               >  2:776632/1‑150 (MQ=255)                                                        tACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGGGCAGAAAATGGGTAAG‑AAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGGTAt                                                                                           >  1:169307/1‑150 (MQ=255)                                                             cACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGTGCAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCtggtg                                                                                     <  1:907587/151‑1 (MQ=255)                                                               cTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGGGCAGAAAATGGTTAAG‑AAAACCACGCTGCATCCTAAATCTTTCCGTATTCCGGGTTATTtggtgg                                                                                    >  2:93705/1‑150 (MQ=255)                                                                ttttGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGGTGGTG‑GAGGGGCAGGAAATGGTTTAG‑AAAGACACCCCGGATCCTAAATCTGTCCGTATTGCggggttttttggggg                                                                                   >  2:40393/1‑138 (MQ=255)                                                                ttttGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGGTGATG‑CAGGGGCAGGAAAAGGTTAAG‑AAAACCACGCTGCATCCTAAAACTGTCCGTATTCCGGGGTATCtggggg                                                                                    >  1:310176/1‑149 (MQ=255)                                                                         tGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGGTGATG‑CAGGGGGAGGAAAAGGGTAAG‑AAAGCCACGCTGCATCCTAAAACTGTCCGTATTTCGGGTTATCTTGTGGGTATTgtggtg                                                                         >  1:908956/1‑151 (MQ=255)                                                                          gAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGTGCAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTgtggtg                                                                         <  1:41150/150‑1 (MQ=255)                                                                              tGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGTGCAGAAAATGGGTAAG‑AAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGGTATCTGGGGGGTATTGTGGTGGTc                                                                      >  1:462384/1‑149 (MQ=255)                                                                               gATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGGATTGAGATGATG‑GAGGGGCAGAAAATGGTTAAG‑GAAACCACGCTGCATCCTCAATCTTTCCGTATTTCGGGGTATTTTGGGGATAttttcgggggcgg                                                                    >  2:422228/1‑139 (MQ=255)                                                                                    aTCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGGTGAAG‑AAAGAGCAGAAAATAGGTAAGAAAAACCACGTTG‑AACCTTAAACTGTGCGTATTCCGGGTTATCtggtgg                                                                                    >  2:94172/1‑129 (MQ=255)                                                                                       tcGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGGTGGTT‑CAGGGGGAGAAAATGGTTAAA‑AAAGACACGCCGGATCCTAAAACTGTCCGGATTCCGGGGTTTCTGGGGGGTATTGTGGGGGGGGGTCCGGATc                                                            >  2:866197/1‑150 (MQ=255)                                                                                           tgCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGGGATGGTG‑GCAGTGCCGGAAAAGGTTAAG‑AAAACCCCGCTGCATCCTAAATCTGTCCGGCTTTACGGGTATCTTGTGGGGATTTTTGTGGTCGATTCGGGTc                                                            >  2:483837/2‑146 (MQ=255)                                                                                            gCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGGGCAGAAAATGGGTAAG‑AAAACCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGGGGGTATTGTGGTGGTCGATCCGGATCAAAccc                                                      >  1:435205/1‑151 (MQ=255)                                                                                             cATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGGGCAGAAAAAGGGTAAG‑AAAACCACGCTGCATCCTAAATCTGTTCGTATTCCGGGGTATTTGGTGGGTATTTTGGGGGTCGATCCGGATCAAAccc                                                      >  1:247713/1‑150 (MQ=255)                                                                                              aTTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGATGATG‑CAGGGGCAGAAAATGGGTAAG‑AAAGCCACGCTGCATCCTAAATCTTTCCGTATTCCGGGTTATCTGGTGGGTATTTTGGTGGTCGATCCGGATCAAAacc                                                      >  2:731392/1‑149 (MQ=255)                                                                                              aTTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGGATTTTGATGGTG‑CAGGGGCAGAAAAAGGGTTAG‑AAAGGCACGCTGGATCCTAAATCTTTCCGTCTTTCGGGGTATCTGGTGGGTATTTTGGTGGTCGGTCCGGGACAAAccc                                                      >  2:336653/1‑149 (MQ=255)                                                                                                 ggTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGATGATG‑CAGGGGCAGAAAAAGGGTAAG‑AAAGCCACGCCGCATCCTAAATCTGTCCGTATTCCGGGTTATTTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAAc                                                   >  1:687023/1‑149 (MQ=255)                                                                                                    tATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGTGCAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTa                                               <  1:445448/150‑1 (MQ=255)                                                                                                    tATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGTGCAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTa                                               <  1:438297/150‑1 (MQ=255)                                                                                                         cccAGGCGGTGCACAATAACGGCGGAATTGTTTTGGTGCCAGCGGCGAAAAAAGGTCATC‑AAA‑ACACGCTGCAACCTTAATCTGTCCGCGGTCCGGGCTATTTTGGGGATATTTTTGTGGGCCATTCTGGTAAAAACCAACTGGATTGCg                                          >  1:832071/1‑150 (MQ=255)                                                                                                           cAGGCGGTGCACAATAACGGCGGTATTTTGATGATG‑CAGGGGCAGAAAATGGGTAAG‑AAAACCACGCTTCATCCTCAATCTGTCCGGCTTCCGGGTTTTCTGGCGGGTATTGTGGTGGTCGATCCGGGTCAAACCCCAATGTATGgggggg                                       >  2:799454/1‑149 (MQ=255)                                                                                                            aGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGTGCAGAAAAAGGTTAAG‑AAAAACACGCTGCATCCTAAAACTGTCCGTATTCCGGGGTATCTGGTGGATATTGTGGGGGTCGGTCCGGATCAAAACCCACTGTATGGCgggg                                       >  1:612107/1‑148 (MQ=255)                                                                                                            aGGCGGTGCAAAATAACGGCGGTATTGTGATGATT‑CAGGGGGAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATCTGTCCGTATTTCCGGGTATCTGGTGGATATTTTGGGGGGCGATCCGGATCAAACCCAACTGTATGGCgggg                                       >  2:89528/1‑148 (MQ=255)                                                                                                             ggCGGTGCACAATAACGGCGGTATTGTGATGATG‑CCGGGGCAGAAAATGGGTAAG‑AAAGCCACACTGCATCCTAAATCTGTCCGCATTCCGGGGTATTTGGTGGGTATTGTGGGGGTCGATCCGGCTTAAACCCAACCGTATGGCGGTg                                       >  1:356476/1‑149 (MQ=255)                                                                                                             ggCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGTGCAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCa                                     <  2:614930/151‑1 (MQ=255)                                                                                                                 gTGCACAATAACGGCGGTATTGTGATGATG‑CAGGTGCAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCgg                                 <  2:887271/151‑1 (MQ=255)                                                                                                                   gCACAATAACGGCGGTATTGTGATGATG‑CAGGTGCAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGtt                               <  1:116952/151‑1 (MQ=255)                                                                                                                       aaTAACGGCGGTATTGTGATGATG‑CAGGTGCAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATCTGTCCGTGTTCCGGGTTGTCTGGTGGATATTGTGGTGGTCGATGCGCATCAAACCCAACTGTATGGCGGTGCACCGGTTa                              <  2:939001/148‑1 (MQ=255)                                                                                                                            cggcggTATTGTGATGATG‑CAGGTGCAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCttt                      <  2:426607/151‑1 (MQ=255)                                                                                                                              gcggTATTGTGATGATG‑CAGGTGCAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTAt                    >  1:655276/1‑151 (MQ=255)                                                                                                                                ggTATTGTGATGATG‑CAGGTGCAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTAttt                  <  2:193643/151‑1 (MQ=255)                                                                                                                                 gTATTGTGATGATG‑CAGGTGCAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTc                 <  2:671576/151‑1 (MQ=255)                                                                                                                                     tgtgATGATG‑CAGGTGCAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTgg              >  1:306018/1‑150 (MQ=255)                                                                                                                                       tgatgaTG‑CAGGTGCAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTg            >  2:665955/1‑150 (MQ=255)                                                                                                                                              g‑CAGGTGCAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCAcc    <  1:825565/151‑1 (MQ=255)                                                                                                                                                 aGGTGCAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCACCCt  <  1:145631/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     AGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATG‑CAGGTGCAGAAAATGGTTAAG‑AAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCACCCT  >  NC_000913/1776117‑1776410

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