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

GGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCGC  >  NC_000913/354141‑354437                                                                                                                                                    |                                                                                                                                                      ggggCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGaa                                                                                                                                                    <  2:899589/151‑1 (MQ=255)        ccGTGACGCTGTCGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATgc                                                                                                                                              <  2:436393/150‑1 (MQ=255)         cGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATgcgc                                                                                                                                            <  2:934532/151‑1 (MQ=255)            gACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTg                                                                                                                                          <  1:345714/150‑1 (MQ=255)                  gACGGTTATCACTTTATTCTCGGGCGCGCTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAgg                                                                                                                                     <  1:751652/149‑1 (MQ=255)                  gACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGAGTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAgg                                                                                                                                     <  2:1070035/149‑1 (MQ=255)                    cGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGGGCTGAAAggg                                                                                                                                    >  1:815923/1‑148 (MQ=255)                      gTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAgg                                                                                                                                <  1:299117/150‑1 (MQ=255)                        tATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtg                                                                                                                              <  2:1073433/150‑1 (MQ=255)                         aTCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtggc                                                                                                                            <  2:116545/151‑1 (MQ=255)                              tttATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGAAGCGCTGAAAGGGCAggggggggggg                                                                                                                       >  2:673291/1‑151 (MQ=255)                                tATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtggcggtggc                                                                                                                      <  2:728827/150‑1 (MQ=255)                                  ttCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGGAggggggggggggggt                                                                                                                   >  1:435684/1‑147 (MQ=255)                                  ttCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAAGCGCTGGAAGGGGAggggggggggggg                                                                                                                     >  2:890904/1‑147 (MQ=255)                                   tctcGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCGGGAAGGGGAggggggggggggggtt                                                                                                                  >  1:877338/1‑146 (MQ=255)                                    ctcGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGAGGCGCTGGAAGGGGAGGGGGCGGGGGCGgtt                                                                                                                  >  2:310149/1‑150 (MQ=255)                                       gggCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGGAGGGGGCGGGGGCGTTTGt                                                                                                                >  1:627837/1‑149 (MQ=255)                                               tGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGTGGTGGGGGGAAAGAGGCGCTGAAAGGGCAGGGGGGGGGGGGGTTTGTCTTTTCGaa                                                                                                       >  1:404149/1‑150 (MQ=255)                                                gACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAag                                                                                                     <  1:128880/151‑1 (MQ=255)                                                  cGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGGAAAGGCGCTGGAAGGGGAGGGGGGGGGGGGGGTTGTCCTTTCGaaa                                                                                                      >  2:758911/1‑148 (MQ=255)                                                    aTGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGCGGGTGGGGGGGAAGAAGGGCTGAAAGGGGAGGGGGCGGGGGGGTTTGTTAtttgg                                                                                                         >  1:528150/1‑140 (MQ=255)                                                      gTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGc                                                                                                <  2:703368/150‑1 (MQ=255)                                                      gTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCg                                                                                               <  1:96160/151‑1 (MQ=255)                                                       tAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGGAGCGCGGGAAAGGCAGGGGGGGGGGGGGGTTttccttcccgaa                                                                                                      >  2:616558/1‑139 (MQ=255)                                                       tAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGAGGGTGGGGTGGAAGATGCGCTGGAAGGGCAGGCGGGGGGGGGGGTTGTCCTTCCGgaa                                                                                                      >  2:1092773/1‑143 (MQ=255)                                                         aTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACa                                                                                            <  2:541573/151‑1 (MQ=255)                                                            aaCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGt                                                                                          <  2:1036845/150‑1 (MQ=255)                                                            aaCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGATGCGCTGAAAGGGGAGGGGGGGGGGGCGGTTGTTCTTTCCAAAGAGAGCggc                                                                                             >  1:895635/1‑145 (MQ=255)                                                                 tGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGAAGGGCTGGAAGGGGAGGGGGGGGGGGGGGTTGTCATTCCGGAAGAGAGGGACAGTCTggg                                                                                     >  2:1020668/1‑149 (MQ=255)                                                                   ccGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGATGCGCTGAAAGGGGAGGTGGCGGGGGCGTTTGTTATTCCGAAAGAGAGCGACAGTCTgg                                                                                      >  1:457855/1‑147 (MQ=255)                                                                   ccGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGGAAGGGCAGGGGGGGGGGGGGGTTGTTATTCCCAAAGAGAGCGACCGGCTgg                                                                                      >  2:491287/1‑147 (MQ=255)                                                                   ccGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGCAGGGGGGGGGGGGGGTTGTCATTTCGAAAGAGGGCGACAGTCTgg                                                                                      >  1:13617/1‑147 (MQ=255)                                                                    cGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGAt                                                                                 <  1:89464/151‑1 (MQ=255)                                                                     gggCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATc                                                                                <  2:486569/151‑1 (MQ=255)                                                                        cATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTg                                                                             <  2:26925/151‑1 (MQ=255)                                                                         aTCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGGGGGGGGGAAGATGGGCTGAAAAGGGCGGGGGGGGGGGGGGTTGTCATTCCCCAAAAGAGAGACAGGCTGGGAGAttg                                                                               >  1:84489/1‑146 (MQ=255)                                                                         aTCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAAATGCGGTGAAAGGGGAGGGGGGGGGGGGGTTTTTCATTTCGAAAGAGAGCGACAGTCTGGGAGAACGTg                                                                             >  1:947423/1‑150 (MQ=255)                                                                             gCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTg                                                                             <  2:734326/146‑1 (MQ=255)                                                                                ggggACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGgc                                                                      <  1:372541/150‑1 (MQ=255)                                                                                      gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCg                                                               <  2:568698/151‑1 (MQ=255)                                                                                        gcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGTGGGGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGGGCACTCGc                                                              >  1:386444/1‑150 (MQ=255)                                                                                          gagaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAAATGCGCTGAAAGGGCAGGGGGCGGGGGCGGTTGTCATTTCGCAAAAGAGCGACAGTCTGGGAGAACGTGGTGTGGGGCACTCGc                                                              >  1:1007593/1‑148 (MQ=255)                                                                                            gaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAagag                                                         <  1:106222/151‑1 (MQ=255)                                                                                             aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAgaga                                                        <  1:190150/151‑1 (MQ=255)                                                                                                    agagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGGGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGGTGTGGCGCACTCGCAAGAGAAgg                                                     >  1:1004493/1‑147 (MQ=255)                                                                                                      agTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTa                                               <  1:28002/151‑1 (MQ=255)                                                                                                       gTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTAt                                              <  2:199866/151‑1 (MQ=255)                                                                                                         aTCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGAAAGAGGCCCTGAAAAGGCAAGGGGCGGGGGGGTTTGGTCTTCCGGAAAAGAGCGGAACTCTGGGAGAACGGGGGGTGGCGGCCTCGCAAGAGAAAGCGAATTTgg                                            >  2:569190/1‑151 (MQ=255)                                                                                                                          cGTTGCCGAAGTGGCGGTGGGTTGGGGGAAAAATGCGCTGAAAGGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGAACGGGATGTGGCGCACTCGCCAGAGAAAGCGATTATGGCGCTGGTGGACAGcc                             >  2:1063534/1‑149 (MQ=255)                                                                                                                             tGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGAt                         <  2:237085/150‑1 (MQ=255)                                                                                                                                            ggtTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgccg          <  2:635693/150‑1 (MQ=255)                                                                                                                                               tGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGccc       <  2:521028/150‑1 (MQ=255)                                                                                                                                                    tGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGgcgc  >  2:935323/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCGC  >  NC_000913/354141‑354437

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