Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	354,288	0	T	G	45.3%	72.2 / 18.9	53	V528G (GTG→GGG)	prpE	propionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (7/22);  new base G (24/0);  total (31/22)
Fisher's exact test for biased strand distribution p-value = 3.50e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGG  >  NC_000913/354138‑354434                                                                                                                                                        |                                                                                                                                                    aCTGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGTTGGGGt                                                                                                                                                     <  2:423126/151‑1 (MQ=255)    ggggCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGTTGGGGTGaa                                                                                                                                                  <  2:437137/151‑1 (MQ=255)    ggggCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGTTGGGGTGa                                                                                                                                                   <  1:401465/150‑1 (MQ=255)       gCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGTTGGGGTGAAAGa                                                                                                                                               <  1:62253/151‑1 (MQ=255)          tCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTTGGGGGGAAAgt                                                                                                                                               >  2:470208/1‑148 (MQ=255)               gACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGGTGGGGGGAAAGATGCGCTg                                                                                                                                        >  2:90559/1‑150 (MQ=255)                 cGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGTTGGGGTGAAAGATGCGCTGaaa                                                                                                                                     <  2:153019/151‑1 (MQ=255)                       cGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGGGGTTGGGGTGAAAAATGCGCTGAAAggg                                                                                                                                  >  2:613902/1‑148 (MQ=255)                         gTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggt                                                                                                                             <  1:747369/151‑1 (MQ=255)                          ttATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGGGGGTGGGGGGGAAGATGCGCTGAAAGGGCAgggg                                                                                                                            >  2:386128/1‑149 (MQ=255)                              cACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtggcg                                                                                                                         <  1:312970/150‑1 (MQ=255)                                    aTTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGGTGGGGGGGAAGATGCGCCGGAAGGGGCggtggggggggggg                                                                                                                  >  1:290096/1‑148 (MQ=255)                                          gggCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGGTGCCGAAGTGGC‑GGGGGGTGGGGGGGAAGATGGGCTGGAAAGGCAggggggggggggggtt                                                                                                                >  2:437691/1‑142 (MQ=255)                                            gcgcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCAtt                                                                                                         <  2:634195/151‑1 (MQ=255)                                                aCTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGa                                                                                                     <  1:674416/151‑1 (MQ=255)                                                    aCGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGGTTGGGGGAAAAAAGGGCTGGAAGGGGAGGGGGGGGGGGGGTTTTGTC‑TTCCGaaa                                                                                                   >  1:38326/1‑149 (MQ=255)                                                      gATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGGGGGTGGGGGGGAAGATGCGCTGGAAGGGCAGGGGGGGGGGGCG‑GTTGTCATTTCGGAagag                                                                                                >  2:11317/1‑150 (MQ=255)                                                       aTGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGGTGGGGGGAAAGATGCGCTGGAAGGGGAGGGGGGGGGGGGG‑GTTGTCATTCCGaaa                                                                                                   >  1:231867/1‑146 (MQ=255)                                                       aTGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGGGGGTGGGGGGGAAGATGCGCTGAAAGGGCCGGGGGGGGGGGGG‑TTTGTCATTCCGGAaaag                                                                                                >  2:464040/1‑149 (MQ=255)                                                         gTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGGGGGTGGGGGGGAAGAAGCGCTGGAAGGGCAGGGGGGGGGGGCG‑GTTGTCATTCCGCAAGAGAGCg                                                                                            >  1:324217/1‑151 (MQ=255)                                                               aaCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGTGGCGGGGGCG‑TTTGTCATTCCGAAAGAGAGCg                                                                                            >  2:143301/1‑145 (MQ=255)                                                                   ttGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGCAGTGGC‑GGTTGGTGGGGGGAAAAAAACGCGGAAACGGCAGGGGGGGGGGGCG‑GTTGGCTTTCCGTAAGAGAGCGTCAGTCTg                                                                                    >  1:532666/1‑149 (MQ=255)                                                                    tGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGGTTGGGGGGAAGATGCGCTGAAAAGGCAGGTGGGGGGGGCG‑GTTGTCATTCCGAAAGAGAGCGACAGTCTgg                                                                                   >  2:181178/1‑149 (MQ=255)                                                                           cATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGGTGGGGGGGAAGATGCGCCGGAAGGGGAGGTGGGGGGGGCG‑GTTGTTGTTCCGGAAGAGAGCGACAGGCCGGgaga                                                                               >  1:591183/1‑146 (MQ=255)                                                                            aTCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGGGGGTGGGGGGGAAGAAGGGCTGGAAGGGGAGGGGGGGGGGGGG‑GTTTTGATTTCGGAAGAGAGCGGCaagttggga                                                                                 >  1:57298/1‑135 (MQ=255)                                                                            aTCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGGGGGTGGGGGGGAAGAAGCGCCGGAAGGGGAGGGGGGGGGGGGG‑GTTTTTATTTCGGAAGAGAGGGGCAGGCTGGGAGAACg                                                                            >  1:803852/1‑148 (MQ=255)                                                                                     ggACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGGGGTTGGGGGGAAAGATACACGGGAAGGGCAGGGGGGGGAGGGG‑GTTGTTTTTTCGAAAGAGAAGGACAGACTTGAAGATCGTgg                                                                         >  2:194394/1‑141 (MQ=255)                                                                                         gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCg                                                            <  1:111739/151‑1 (MQ=255)                                                                                         gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTAGGTGGGGGGGAAGAAGAGCTGCAAGGGCAGGAGGCGGGGGCG‑GTTTTCATTCCGAAAGAGAGGGGCAGTCTGGAAGATccggacgttggcaac                                                               >  1:699417/1‑134 (MQ=255)                                                                                         gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGGGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGCGGCGGGGGGG‑GTTGTCATTCCGAAAGAGAGCGACAGACTGGAAGATCGGGATGTGGCGCACTc                                                             >  1:615281/1‑150 (MQ=255)                                                                                           gcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGGTGGGGGGAAAGATGCGCTGGAAGGGGAGGTGGGGGGGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCCCTc                                                             >  1:219721/1‑148 (MQ=255)                                                                                             gagaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAag                                                        <  1:505114/151‑1 (MQ=255)                                                                                             gagaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAag                                                        <  1:366590/151‑1 (MQ=255)                                                                                                          gTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTAt                                           <  1:482343/151‑1 (MQ=255)                                                                                                             tctcCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGGGGGGGGGAAAAATGCGCTGAAAGGGCAGGGGGGGGGGGCG‑TTTGTCATTCCGAAAGAGAGGGACAGTCTGGGAGATCGTGGTGGGGCGCACTCGCAAGAGAAAGGGGTTATgg                                         >  1:473829/1‑150 (MQ=255)                                                                                                              ctcCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAg                                                                                <  1:193101/110‑1 (MQ=255)                                                                                                              ctcCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGGTGGGGTGGAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTgg                                                                                   >  2:193122/1‑107 (MQ=255)                                                                                                               tcCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGGTGGGGGGAAAAATGCGCTGAAAGGGCAGGTGGCGGGGGGG‑TTTGTCCTTCCGAAAGAGAGCGACAGGCTGGAAGATCGGGAGGTGGCGCACTCGCAAGAGAAGGCGGTTATGgcg                                       >  1:348227/1‑150 (MQ=255)                                                                                                               tcCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGGGGGTGGGGGGAAAGAAGCGCTGAAAGGGCAGGGGGCGGTGGCG‑TTTGTAATTCCGAAAAAGGGCGACAGGCCGGAAGATCGTTATGTTGCGCACTTGCAACAGAAgggggctatagcg                                       >  2:147535/1‑139 (MQ=255)                                                                                                                ccAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGgcg                                       <  2:351983/149‑1 (MQ=255)                                                                                                                ccAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGgc                                        <  1:6962/148‑1 (MQ=255)                                                                                                                ccAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCt                                     <  1:613868/151‑1 (MQ=255)                                                                                                                     cATCCGGGCGTTGCCGAAGTGGC‑GGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtggtgg                                <  1:783591/151‑1 (MQ=255)                                                                                                                      aTCCGGGCGTTGCCGAAGTGGC‑GGTGGGTGGGGGGGAAAATGAGCCGAAAGGGCCGGGGGGGGGGGCG‑GTTGTCATTTCGAAAAAAGGCGACAGTCTGGGAGATCGTGGGGGGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtggggg                                >  1:363603/1‑150 (MQ=255)                                                                                                                                 gCCGAAGTGGC‑GGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTg                    <  2:535394/151‑1 (MQ=255)                                                                                                                                     aaGTGGC‑GGGGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGAGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGAAAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAAAGGAGGCGATTATGGGGCTGGAGGACAGCCAGAtt                     >  1:170665/1‑146 (MQ=255)                                                                                                                                        tGGC‑GGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACtt             >  1:790748/1‑151 (MQ=255)                                                                                                                                         ggC‑GGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACttt            <  1:650366/151‑1 (MQ=255)                                                                                                                                          gC‑GGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTg           <  1:24454/151‑1 (MQ=255)                                                                                                                                           c‑GGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTgg          >  1:704827/1‑151 (MQ=255)                                                                                                                                             ggtggtTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgc         >  1:258631/1‑151 (MQ=255)                                                                                                                                              gtggtTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgcc        <  2:126769/151‑1 (MQ=255)                                                                                                                                                 gtTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgccgcc     <  1:245503/151‑1 (MQ=255)                                                                                                                                                    ggggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCgg  >  2:701040/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                    ACTGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGC‑GGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCG‑TTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGG  >  NC_000913/354138‑354434

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