Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	354,282	0	T	G	60.0%	35.0 / 56.3	45	V526G (GTT→GGT)	prpE	propionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (2/16);  new base G (27/0);  total (29/16)
Fisher's exact test for biased strand distribution p-value = 2.37e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-07

GGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAT‑GCGCTG‑AAAGGGCAGGTGGCGGTGGCGTTTGTCA‑TTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCG  >  NC_000913/354142‑354429                                                                                                                                             |                                                                                                                                                       gggCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTgggggggaa                                                                                                                                              >  1:74148/1‑146 (MQ=255)    cATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGa                                                                                                                                            <  2:557811/150‑1 (MQ=255)    cATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGa                                                                                                                                            <  1:827118/150‑1 (MQ=255)            aCGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAAAA‑gcgccgga                                                                                                                                  >  1:61586/1‑147 (MQ=255)                   cGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGCGATTGAAGAGAGAATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGGCGTGG‑AAAGGGcg                                                                                                                           >  2:682185/1‑150 (MQ=255)                   cGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAT‑GCGCTG‑GAAGGGGAg                                                                                                                          >  2:692176/1‑151 (MQ=255)                           aCTTTATTCTCGGGCCCACTGACGATGTAATTAACGTTGCCGGGCCTCTGCTGGGGACTCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGAGGGTGGGGGGAAAAAT‑ACCCCG‑GACGGGGAggtggcggg                                                                                                                  >  2:92842/1‑150 (MQ=255)                            cTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTTGGGGGGAAGAT‑GGGCTG‑GAAGGGGGgggggcgggg                                                                                                                 >  1:686944/1‑149 (MQ=255)                            cTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAA‑GGGCCG‑AAAGGGGAgggggggggg                                                                                                                 >  1:666009/1‑149 (MQ=255)                                             cTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGTGAAAAAT‑GCGCTG‑AAAGGGCAGGGGGGGGGGGGGGTTGTCA‑GTCCgga                                                                                               >  2:797757/1‑149 (MQ=255)                                                aCGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAT‑GCGCTG‑AAAGGGCAGGTGGCGGTGGCGTTTGTCA‑TTCCGAAag                                                                                             <  1:1057484/150‑1 (MQ=255)                                                  gATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGAAAGAA‑GGGCTG‑GAAGGGCCGGGGGGGGGGGGGGTTGTTA‑TTCCGAAagagg                                                                                          >  2:1031106/1‑150 (MQ=255)                                                   aTGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAAAA‑GCGCTG‑GAAGGGGGGGGGGGGGGGGGGGTTGTGG‑TTTCGAAagagag                                                                                         >  2:499196/1‑151 (MQ=255)                                                     gaaaTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGAT‑GCGCTG‑AAAGGGGCGGGGGGGGGGGCGTTTTGCT‑TTCCGAAagagaggg                                                                                       >  2:796511/3‑149 (MQ=255)                                                     gTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTTGGGTGGAAAAT‑GAGCTGTAAAGGCCGGGGGGCGGGGGGGTTTGTcttctgtaaaaaggag                                                                                         >  1:1011938/1‑136 (MQ=255)                                                     gTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGGGGGGGGGAAAAA‑GCGCCG‑GAAGGGCGGGGGGGGGGGGGGGTTGGCATTTCCGAAagggggg                                                                                        >  1:658126/1‑146 (MQ=255)                                                       aaTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAAAT‑GCGCTG‑GAAGGGGCGGGGGGGGGGGGGGTTGTCC‑TTTCGGAagagagggg                                                                                      >  2:61737/1‑147 (MQ=255)                                                               ttGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAAAT‑GGGCTG‑AAAGGGGAGGGGGGGGGGGGGGTTGTCA‑TTCCGCAAGAGAGGGACAGGCTGGa                                                                             >  2:927394/1‑151 (MQ=255)                                                                  ccGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAT‑GCGCTG‑AAAGGGCAGGTGGCGGTGGCGTTTGTCA‑TTCCGAAAGAGAGCGACAGTCTGGAAg                                                                           <  2:706447/150‑1 (MQ=255)                                                                  ccGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAT‑GCGCTG‑AAAGGGCAGGTGGCGGTGGCGTTTGTCA‑TTCCGAAAGAGAGCGACAGTCTGGAAg                                                                           <  2:55069/150‑1 (MQ=255)                                                                   cGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAT‑GCGCTG‑AAAGGGCAGGTGGCGGTGGCGTTTGTCA‑TTCCGAAAGAGAGCGACAGTCTGGAAGAt                                                                         <  2:572829/151‑1 (MQ=255)                                                                      gCATCGGCTGGGGCCGCGGGAGATTGGCGAGAGTATATGCAGCAATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAT‑GCGCTG‑AAAGGGCAGGTGGCGGTGGCGTTTGTCA‑TTCCGAAAGAGAGCGACAGTCTGGAAGATCg                                                                       <  1:92833/150‑1 (MQ=255)                                                                      gCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAT‑GCGCTG‑AAAGGGCAGGTGGCGGTGGCGTTTGGCA‑TTCCGAAAGAGAGCGACAGTCTGGAAGATCg                                                                       <  2:910506/150‑1 (MQ=255)                                                                           ggCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAT‑GCGCCG‑AAAGGGGAGGTGGGGGGGGCGGTTGTCA‑TTCCGAAAGAGAGCGACAGTCTGGGAGATCGGGgtgt                                                                 >  1:452495/1‑151 (MQ=255)                                                                                 ggACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGTGAAAAAT‑GCGCTG‑AAAGGGCGGGGGGGGGGGGCGGTTGTGA‑TTTCGAAAGAGAGGGACAGTCTGGGAGAGCGGGGTGGGgcgc                                                            >  2:554566/1‑150 (MQ=255)                                                                                   aCGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAT‑GCGCTG‑AAAGGGCAGGTGGCGGTGGCGTTTGTCA‑TTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCAc                                                          <  1:239920/150‑1 (MQ=255)                                                                                          agaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAT‑GCGCTG‑AAAGGGCAGGTGGCGGTGGCGTTTGTCA‑TTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCGCTCGCAag                                                   <  2:795548/150‑1 (MQ=255)                                                                                          agaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGAAAGAA‑GCGCTG‑AAAAGGCAGGTGGGGGGGGGGTTTGTCA‑TTCCGAAAAAGGGAGACAGTCCGGAAGAACGTGAGGGGGCGCGCTCGCCagg                                                  >  1:1029029/1‑150 (MQ=255)                                                                                          agaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGCGGGTGGGGGGAAAAAT‑GCGCTG‑AAAGGGCCGGGGGCGGGGGGGGTTGTCA‑TTCCGCAAGAGAGCGACAGGCTGGGAGATCGTGGTGTGGCGGACTCGCCaga                                                  >  1:611443/1‑151 (MQ=255)                                                                                            aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAT‑GCGCTG‑AAAGGGCAGGTGGCGGTGGCGTTTGTCA‑TTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAgaga                                                <  2:156969/151‑1 (MQ=255)                                                                                            aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAT‑GCGCTG‑AAAGGGCAGGTGGCGGTGGCGTTTGTCA‑TTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAgaga                                                <  1:691920/151‑1 (MQ=255)                                                                                            aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGAAAGAT‑TCGCTG‑AAAGGGGAGGGGGGGGGGGGGGTTGTGC‑TTTCGGAAGAGAGGGACAGGCTGGGAGAACGGGGTGTGGGGCACTCGCCAgaga                                                >  2:623517/1‑151 (MQ=255)                                                                                            aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGAT‑GGGCTG‑GAAGGGGAGGGGGGGGGGGCGTTTGTCA‑TTTCGGAAGAGAGGGCCCGGCTGGGAGGTCGGGATTTGGCGCACTCCCCAgaga                                                >  2:1018347/1‑151 (MQ=255)                                                                                                 agagagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGAT‑GCGCTG‑GAAGGGCAGGGGGCGGTGGCGTTTGTCA‑TTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCg                                           >  1:860020/1‑151 (MQ=255)                                                                                                    gcgTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGGAAGAC‑GCGCCG‑GAAGGGGAGGGGGGGGCGGCGGTTGGCC‑TTCCGAAAGAGAGTGACAGTCTGGGAGGTCGGGCTGTGGCGCACCCGCaacgggaggggggt                                        >  1:42702/3‑139 (MQ=255)                                                                                                    gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGGAAAAA‑GCGCTG‑AAAGGGGAGGGGGGGGGGGCGGTTGTCA‑TTCCGAAAGAGAGCGGCAGTCTGGAAGATCGGGATGTTGCGCACTCGCAAGAGAAGGCGAtt                                        >  1:984494/1‑151 (MQ=255)                                                                                                    gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGAT‑GCGCTG‑AAAGGGGAGGGGGGGGGGGGGGTTGTCA‑TTCCGAAAAAGAGGGGAAGTCTGGGAGAGCGGGAGGTGGGGCCCTCGCAAAAGAAgggggtt                                        >  2:399844/1‑146 (MQ=255)                                                                                                        aTCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAAAT‑GCGCTG‑AAAGGGCAGGGGGCGGGGGGGTTTGTCC‑TTCCGAAAGAGAGGGACAGGCTGGGAAATCGGGAGGTGGGGCCCTCGCAAGAGAAGGGGGTTATgg                                    >  2:1006993/1‑151 (MQ=255)                                                                                                        aTCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTTGGGGGAAAGAT‑GCGCCG‑AAAGGGCAGGTGGCGGGGGGGGTTGTGA‑TTCCGAAAGAGAGCGACAGTCTGGAAGATCGGGATGTGGCGCACTCGCAAGAGAAGGCGATTATgg                                    >  1:67553/1‑151 (MQ=255)                                                                                                         tctcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGAA‑GCGCTG‑AAAGGGGAGGGGGCGGTGGCGTTTGTGA‑TTCCGAAAGAGAGCGACAGTCTGGGAGATCGTGTTGTGGCGCACCACCCAGAGAAGGCGGTTATgg                                    >  1:388949/1‑150 (MQ=255)                                                                                                           tcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAT‑GCGCTG‑AAAGGGCAGGTGGCGGTGGCGTTTGTCA‑TTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGgcgc                                 <  2:986356/151‑1 (MQ=255)                                                                                                                ccATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAT‑GCGCGG‑AAAGGGCAGGTGGCGGTGGCGTTTGTCA‑TTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtggtg                            <  2:227328/151‑1 (MQ=255)                                                                                                                  aTCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAT‑GCGCTG‑AAAAGGCGGGGGGGGGGGGGGGTTGTCG‑TTTCGGAAGAGGGCGGCAGTCTGGGAGATCGTGGTGTTGCGCCCCCGCCAGAGAAAGCGGATATGGCGCcggtggg                          >  1:505384/1‑150 (MQ=255)                                                                                                                     cGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAT‑GCGCTG‑AAAGGGCAGGTGGCGGTGGCGTTTGTCA‑TTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAg                       <  1:208336/151‑1 (MQ=255)                                                                                                                       ggCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAT‑GCGCTG‑AAAGGGCAGGTGGCGGTGGCGTTTGTCA‑TTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGc                      <  2:679897/150‑1 (MQ=255)                                                                                                                                         gtggtTGGGGTGAAAGAT‑GCGCTG‑ACAGGGCAGGTGGCGGTGGCGTTTGTCA‑TTCCGAAAGAGAGCGACAGTCTGGAAGATAGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgc    >  2:836915/1‑150 (MQ=255)                                                                                                                                          tggtTGGGGTGAAAGAT‑GCGCTG‑AAAGGGCAGGTGGCGGTGGCGTTTGTCA‑TTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgccg  >  1:919093/1‑151 (MQ=255)                                                                                                                                           ggtTGGGGTGAAAGAT‑GCGCTG‑AAAGGGCAGGTGGCGGTGGCGTTTGTCA‑TTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgccg  >  2:136306/1‑150 (MQ=255)                                                                                                                                             |                                                                                                                                                       GGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAT‑GCGCTG‑AAAGGGCAGGTGGCGGTGGCGTTTGTCA‑TTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCG  >  NC_000913/354142‑354429

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