Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	354,288	0	T	G	52.3%	53.2 / 71.4	44	V528G (GTG→GGG)	prpE	propionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (3/18);  new base G (22/1);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 2.81e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-08

CATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCG‑CTGAAAGGGCAGGT‑GGCGGTGGCGTTTGT‑CATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCG  >  NC_000913/354145‑354436                                                                                                                                                |                                                                                                                                                        cATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGa                                                                                                                                                   <  2:505543/150‑1 (MQ=255) cATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGa                                                                                                                                                   <  2:353740/150‑1 (MQ=255) aatcggtcgcGCTGACGGTTATCGCTCTATGCTCGGGCGAACTGACGATGTAATTAACGTTGCCGGGCATCGGGTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGa                                                                                                                                                   <  1:61837/141‑1 (MQ=255)     cGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGAAAAATggg                                                                                                                                               >  2:418179/1‑148 (MQ=255)          cGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATgcg‑ccgcaa                                                                                                                                        >  1:465023/1‑146 (MQ=255)                cGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGGAAGATGTG‑GCGGAAGGGccg                                                                                                                                  >  1:153375/1‑149 (MQ=255)                cGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTTGGGGGGAAGAAGCGCCGGAAAGGGcg                                                                                                                                   >  1:668178/1‑150 (MQ=255)                 ggTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGCG‑CTGAAAGGGCCgg                                                                                                                                 >  2:616173/1‑151 (MQ=255)                                              cGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCG‑CTGAAAGGGCAGGT‑GGCGGTGGCGTTTGT‑CATTCCGAAagag                                                                                                  <  2:561414/151‑1 (MQ=255)                                                aTGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTTTGGGGGAAGGTGTG‑CCGTAAAGGGAGGT‑GGCGGGGGGGGTTTT‑CATTCCGGAaaagag                                                                                                >  2:275756/1‑151 (MQ=255)                                                aTGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAAAAGGG‑CTGGAAGGGGCGGG‑GGGGGGGGGGTTTGT‑CATTTCGGAagagag                                                                                                >  1:669244/1‑151 (MQ=255)                                                 tGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCG‑CTGAAAGGGCAGGT‑GGCGGTGGCGTTTGT‑CATTCCGAAAGAGAGc                                                                                               <  2:156545/151‑1 (MQ=255)                                                  gTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGGG‑CTGAAAGGGCCGGG‑GGCGGGGGGGGTTGT‑CCTTTCGAAagagaggg                                                                                              >  2:663899/1‑149 (MQ=255)                                                  gTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGAAGGG‑GTGCAAGGGGAGGG‑GGCGGGGGGGGTTGTGAATTCCGaagaggggg                                                                                               >  2:488975/1‑144 (MQ=255)                                                   tAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTTGGGGGAAAGATGCG‑CTGAAAAGGCAGGT‑GGCGGGGGCGGTTGT‑CATTTCTAAAGGGAGCgg                                                                                             >  2:641876/1‑150 (MQ=255)                                                    aaTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGGGGCGGTGGTTGGGGGGAAAAAAACG‑TCGAGAGGGGCGGG‑GGGGGTGGGGTTTTT‑GATTTCGAAATAGAGCgga                                                                                            >  1:310078/1‑149 (MQ=255)                                                               ccGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCG‑CTGAAAGGGCAGGT‑GGCGGTGGCGTTTGT‑CATTCCGAAAGAGAGCGACAGTCTGGAAg                                                                                  <  1:767301/150‑1 (MQ=255)                                                                cGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCG‑CTGAAAGGGCAGGT‑GGCGGTGGCGTTTGT‑CATTCCGAAAGAGAGCGACAGTCTGGAAGAt                                                                                <  2:465255/151‑1 (MQ=255)                                                                     aTCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCG‑CTGAAAGGGCAGGT‑GGCGGTGGCGTTTGT‑CATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGa                                                                           <  2:807614/151‑1 (MQ=255)                                                                         gCTGGGGACGCGCGAGGTTGAAGAGAGTAGCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCG‑CTGAAAGGGCAGGT‑GGCGGTGGCGTTTGT‑CATTCCGAAAGAGGGCGACAGTCTGGAAGATCGGGAtgt                                                                        <  1:60163/150‑1 (MQ=255)                                                                                aCGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCG‑CTGAAAGGGCAGGT‑GGCGGTGGCGTTTGT‑CATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCAc                                                                 <  2:26043/150‑1 (MQ=255)                                                                                  gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAAAAGCG‑GTGGAAGGGGCGGG‑GGCGGGGGGGGTTTT‑CATTCCGGAAGAGAGCGACAGGCTGGGAGATCGGGATGGGGCGGACTCg                                                              >  1:747677/1‑151 (MQ=255)                                                                                       agaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCG‑CTGAAAGGGCAGGT‑GGCGGTGGCGTTTGT‑CATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAag                                                          <  2:753544/150‑1 (MQ=255)                                                                                       agaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAAAAGGG‑CTGGAAGGGGAGGG‑GGCGGTGGGGTTTGT‑CATTCCGAAAGAGCGGGACAGTCTCGAAGATCGTGGTTGGGCGCCCTCGCCaga                                                         >  1:765128/1‑151 (MQ=255)                                                                                       agaTTGAAGAGAGTATCGCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGTGAACAAGAG‑CTTTAAGAGGAGGTGGGCGGTGG‑GTTTGT‑CCTTCCGAAAGAGAGCCACAGTCTGGCAGATCGTGGTGTTGGGCACttccccga                                                         >  1:68799/1‑144 (MQ=255)                                                                                         aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCG‑CTGAAAGGGCAGGT‑GGCGGTGGCGTTTGT‑CATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAagag                                                        <  2:359154/150‑1 (MQ=255)                                                                                         aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCG‑CTGAAAGGGCAGGT‑GGCGGTGGCGTTTGT‑CATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAagag                                                        <  2:284555/150‑1 (MQ=255)                                                                                         aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCG‑CTGAAAGGGCAGGT‑GGCGGTGGCGTTTGT‑CATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAagag                                                        <  2:104788/150‑1 (MQ=255)                                                                                         aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCG‑CTGAAAGGGCAGGT‑GGCGGTGGCGTTTGT‑CATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAgaga                                                       <  1:407260/151‑1 (MQ=255)                                                                                         aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCG‑CTGAAAGGGCAGGT‑GGCGGTGGCGTTTGT‑CATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAgaga                                                       <  2:572120/151‑1 (MQ=255)                                                                                         aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGATGCG‑CTGAAAGGGCAGGG‑GGGGGGGGGGGTTGT‑CATTCCGAAAGGGAGGGACAGGCTGGAAGATCGGGGTGTGGGGCACCCGCAAaaga                                                       >  2:301363/1‑151 (MQ=255)                                                                                         aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCG‑CTGAAAGGGCAGGG‑GGGGGGGGGGTTTGT‑CATTCCGAAAGAGAGCGACAGGCTGGAAGATCGGGGTGTGGCGCCCCCGCAAaaga                                                       >  2:84327/1‑151 (MQ=255)                                                                                            gAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGAGGTTGGGGGGAAAGATGCG‑CTGAAAGGGCAGGG‑GGGGGGGGGGTTTGT‑CATTCCGAAAGAGAGGGACAGTCTGGAAGATCGTGATGTGGGGCACTCGCAAAAGAAgg                                                    >  2:568476/1‑151 (MQ=255)                                                                                                   gTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTTGGGGGGAAGATGCG‑CTGAAAGGGCAGGG‑GGGGGGGGGGGTTGT‑CATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGGGGCGCACTCGCAAGAGAAGGGGATTAt                                             >  1:50772/1‑151 (MQ=255)                                                                                                    tATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCG‑CTGAAAGGGCAGGT‑GGCGGTGGCGTTTGT‑CATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTAt                                             <  2:334285/150‑1 (MQ=255)                                                                                                     aTCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCG‑CTGAAAGGGCAGGG‑GGCGGGGGGGGTTGT‑CATTCCGAAAGAGAGCGACAGTCTGGAAAATCGTGGGGGGGGGCACCCGCAAGAGAAGGCGGTTAttg                                           >  2:34685/1‑149 (MQ=255)                                                                                                     aTCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGAGGGTGGGGGGAAAGATGAG‑CTGAAAAGGCAGGG‑GGCGGGGGCGTTTGT‑CCTTCCGAAAGAGAGGGACAGTCTGGAAAAACGTGATGTGGCGCACTCGCAAGAGAAAGAGAATATgg                                           >  2:696757/1‑151 (MQ=255)                                                                                                               aTCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCG‑CTGAAAGGGCAGGG‑GGCGGGGGCGTTTGT‑CATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGGTGTGGCGCCCTCGCCAGAGAAAGCGATTATTGCGCtggtggg                                 >  2:499599/1‑150 (MQ=255)                                                                                                                tCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGAAAAATGCG‑CTGGAAGGGGCGGG‑GGGGGGGGGGGTTTT‑GCTTTCGGAAGAGAGGGGCAGTCTGGGAGAACGTGGTGTGGCGCACTCGCCAGCGAAGGCGGGTATTGGGCtgggggga                                >  1:601065/1‑149 (MQ=255)                                                                                                                        ttGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCG‑CTGAAAGGGCAGGT‑GGCGGTGGCGTTTGT‑CATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGAt                        <  1:24691/151‑1 (MQ=255)                                                                                                                              aaGTGGCGGTGGTTGGGGTGAAAGATGCG‑CTGAAAGGGCAGGT‑GGCGGTGGCGTTTGT‑CATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCaa                  <  1:102055/151‑1 (MQ=255)                                                                                                                                             gggTGAAAGATGCG‑CTGAAAGGGCAGGT‑GGCGGTGGCGTTTGT‑CATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCgg    >  2:553061/1‑150 (MQ=255)                                                                                                                                             gggTGAAAGATGCG‑CTGAAAGGGCAGGT‑GGCGGTGGCGTTTGT‑CATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCgg    >  2:464218/1‑150 (MQ=255)                                                                                                                                               gTGAAAGATGCG‑CTGAAAGGGCAGGT‑GGCGGTGGCGTTTGT‑CATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGgcg  >  2:797989/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                        CATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCG‑CTGAAAGGGCAGGT‑GGCGGTGGCGTTTGT‑CATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCG  >  NC_000913/354145‑354436

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