Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	354,282	0	T	G	53.5%	40.3 / 52.3	44	V526G (GTT→GGT)	prpE	propionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (9/11);  new base G (23/0);  total (33/11)
Fisher's exact test for biased strand distribution p-value = 2.92e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-06

GACTGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGC  >  NC_000913/354137‑354430                                                                                                                                                  |                                                                                                                                                     gACTGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTgggg                                                                                                                                                 <  1:232501/151‑1 (MQ=255)  aCTGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTggggg                                                                                                                                                >  1:435488/1‑150 (MQ=255)    tGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTgggggg                                                                                                                                               >  1:869075/1‑148 (MQ=255)    tGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTgggg                                                                                                                                                 >  1:105342/1‑148 (MQ=255)     ggggCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGaa                                                                                                                                             <  2:868276/151‑1 (MQ=255)         cATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAt                                                                                                                                         <  2:701223/151‑1 (MQ=255)         cATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGa                                                                                                                                          >  1:174683/1‑150 (MQ=255)            ccGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATgcg                                                                                                                                      <  1:457677/151‑1 (MQ=255)            ccGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGAAAg                                                                                                                                           >  1:650426/1‑146 (MQ=255)             cGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGAAAGATgcgc                                                                                                                                     >  2:305231/1‑151 (MQ=255)                     tGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGATGCGCCGAAAggg                                                                                                                             >  1:849782/1‑151 (MQ=255)                       aCGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGCGCTGAAAggg                                                                                                                             >  2:4216/1‑149 (MQ=255)                         ggTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAgg                                                                                                                         <  2:844750/151‑1 (MQ=255)                          gTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGGGCTGGAAgggggggg                                                                                                                        >  1:741409/1‑146 (MQ=255)                           ttATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGAAAGATGCGCTGGAAGGGCAgg                                                                                                                         >  1:758678/1‑149 (MQ=255)                                  tttATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGAAAGATGCGGTGAAAGGGTAggtggcgg                                                                                                                   >  2:251285/1‑148 (MQ=255)                                     aTTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGt                                                                                                             <  1:381815/151‑1 (MQ=255)                                        ctcGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGgtt                                                                                                           >  2:187510/1‑150 (MQ=255)                                        ctcGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATTGGCGGAAAGGGCAGGGGGGGGGGGCGtgt                                                                                                           >  2:185563/1‑148 (MQ=255)                                              cgcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGATGGGCTGAAAGGGGAGGTGGCGGGGGGGGTTGTCAtt                                                                                                     >  2:352008/1‑150 (MQ=255)                                                     aCGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAaga                                                                                             <  2:451243/151‑1 (MQ=255)                                                        aTGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAagagag                                                                                          <  1:393445/151‑1 (MQ=255)                                                               tAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGATGCGCTGAAAGGGGAGGTGGCGGGGGGGTTTGTCCTTTCGAAAGAGAGCGAc                                                                                      >  2:73064/1‑148 (MQ=255)                                                                aaCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGGGTTTGTCATTCCGAAAGAGAGCGACAGTc                                                                                  >  2:582067/1‑151 (MQ=255)                                                                  cGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTg                                                                                <  2:92510/151‑1 (MQ=255)                                                                    ttGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCCGGGGGGGGTGGGGGTTGTCATTCCGAAAGAGAGCGACAGGCTgg                                                                               >  1:456807/1‑150 (MQ=255)                                                                          ggCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTTGGGGGGAAGATGCGCCCTAAGGGGAGGGGGGGGGGGCGTTTGTCATTCCCGAAGAGCGCGACaagccggaggggcg                                                                        >  1:739095/1‑138 (MQ=255)                                                                             aTCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGCAGGGGGGGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGGAGAACggg                                                                      >  2:305883/1‑148 (MQ=255)                                                                             aTCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGGAAAATGGGCTGAAAGGGCGGGTGGGGGGGGGGGTTGTTAATCCTGAAAAGAGGGAccggcttggaga                                                                           >  2:859427/1‑134 (MQ=255)                                                                                ggCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGTGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGAtg                                                                   >  2:751815/1‑150 (MQ=255)                                                                                   tGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTgg                                                                >  2:303477/1‑150 (MQ=255)                                                                                          gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGGAAGGGCAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGGCTGGGAGATCGTGATGTGGCGCACTc                                                         >  1:474332/1‑150 (MQ=255)                                                                                           cgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGCAGGTGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGc                                                       >  1:680216/1‑151 (MQ=255)                                                                                           cgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGCAGGTGGCGGGGGCGGTTGTCATTCCGAAAAAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCg                                                        >  2:719843/1‑150 (MQ=255)                                                                                             cgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCg                                                        >  2:725110/1‑148 (MQ=255)                                                                                             cgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCCGAAAGGGCAGGTGGGGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGGAGATCGGGATGTGGCGCACTc                                                         >  1:244888/1‑147 (MQ=255)                                                                                                    gAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGTGGCGGGGGCGTTTGTCATTTCGAAAGAGAGCGGCAGGCTGGAAGATCGTGATGTGGCGCACTCGCAAGaggagg                                              >  2:585963/1‑151 (MQ=255)                                                                                                      agagagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGATGGGGGGAAAGATGCGCTGAAAGGGGAGGGGGGGGGGGGGTTTGTCCTTCCGGAAGAGAGGGAGAGGCTGGGAGATGGTGGTGGGGGGCACTCGCAAGAGAAGGCg                                            >  1:432/1‑151 (MQ=255)                                                                                                               ctcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGgcg                                   >  2:651161/1‑151 (MQ=255)                                                                                                                   aGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtgg                               <  1:113151/151‑1 (MQ=255)                                                                                                                   aGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGAAAAATGCGCTGAAAGGGCAGGTGGCGGGGGCGGTTGTCATTCCGAAAGAGGGGGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAAGCGATTATGGGGCtgg                               >  1:513970/1‑151 (MQ=255)                                                                                                                    gCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGATGCGCTGAAAGGGCAGGGGGCGGTGGCGTTTGTCATTTCGAAAGAGAGCGACAGTCTGGAAGATCGTGGTGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtgg                               >  1:897509/1‑150 (MQ=255)                                                                                                                                  gCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTg                <  2:243424/151‑1 (MQ=255)                                                                                                                                          ggCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACttt        >  1:786157/1‑151 (MQ=255)                                                                                                                                                ggtTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgccgc  >  1:387330/1‑151 (MQ=255)                                                                                                                                                ggtTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgccg   <  1:4216/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GACTGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGC  >  NC_000913/354137‑354430

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: ‑