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

CTGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGG  >  NC_000913/354139‑354417                                                                                                                                                      |                                                                                                                                  cTGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGCTGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTTGGGTg                                                                                                                                  >  2:854166/1‑151 (MQ=255)       cATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAt                                                                                                                            <  2:297712/151‑1 (MQ=255)          ccGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATgcg                                                                                                                         <  1:194009/151‑1 (MQ=255)           cGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATgcgc                                                                                                                        <  2:385878/151‑1 (MQ=255)                   tGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGAAAGATGCGCTGGAAggg                                                                                                                >  2:679114/1‑151 (MQ=255)                     aCGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAAATGCGCTGGAAggg                                                                                                                >  2:678731/1‑149 (MQ=255)                      cGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAg                                                                                                             <  1:474222/151‑1 (MQ=255)                       ggTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAgg                                                                                                            <  2:293577/151‑1 (MQ=255)                        gTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggt                                                                                                           <  1:237589/151‑1 (MQ=255)                         ttATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGGTGAAAGGGCAgg                                                                                                            >  2:297971/1‑149 (MQ=255)                          tATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGAGCGGGAAGGGGAggggg                                                                                                         >  2:479983/1‑151 (MQ=255)                           aTCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtggc                                                                                                        <  2:272301/151‑1 (MQ=255)                            tCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAAGCGCTGAAAGGGCAggggg                                                                                                         >  1:102569/1‑149 (MQ=255)                                 ttATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtggcggtgg                                                                                                   <  2:455852/150‑1 (MQ=255)                                   aTTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGt                                                                                                <  1:416793/151‑1 (MQ=255)                                     tctcGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAAGGCAGGGGGGGGGGGCGttt                                                                                              >  2:319721/1‑151 (MQ=255)                                        cGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGttt                                                                                              <  2:457406/148‑1 (MQ=255)                                             gcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGGAAGGGCAGGTGGGGGGGGCGGTTGTCATTcc                                                                                      >  2:401954/1‑151 (MQ=255)                                             gcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGGAGGGGGCGGGGGCGGTTGTCATTcc                                                                                      >  2:937985/1‑151 (MQ=255)                                                  gACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAag                                                                                 <  2:656867/151‑1 (MQ=255)                                                       tGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAAGCGCTGAAAGGGCAGGGGGGGGGGGGGTGTGTTATTCCGGAaga                                                                                >  1:440625/1‑147 (MQ=255)                                                                cGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTg                                                                   <  1:514162/151‑1 (MQ=255)                                                                 gTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACATTCTgg                                                                  <  1:781586/151‑1 (MQ=255)                                                                   tGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTgg                                                                  >  2:601891/1‑149 (MQ=255)                                                                    gCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGaa                                                                <  1:378575/150‑1 (MQ=255)                                                                      cGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTgg                                                                  >  2:109729/1‑146 (MQ=255)                                                                        ggCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCg                                                           <  1:725143/151‑1 (MQ=255)                                                                        ggCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTgg                                                                  >  2:314816/1‑144 (MQ=255)                                                                                 tGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGGTGCGCTTAAAGGGCAGGTGGCGGGGGCGGTTGTCATTCCGGAAGAGAGCGACAGTCTGGGAGATCGGGGTGTgg                                                   >  1:306700/1‑150 (MQ=255)                                                                                     gACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGTTTGTCATTCCGCAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCAc                                              >  1:908398/1‑151 (MQ=255)                                                                                        gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCg                                           <  1:679037/151‑1 (MQ=255)                                                                                            gagaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAag                                       <  1:323714/151‑1 (MQ=255)                                                                                             agaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAag                                       <  1:542272/150‑1 (MQ=255)                                                                                              gaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAag                                       <  2:683573/149‑1 (MQ=255)                                                                                               aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAag                                       >  2:121780/1‑148 (MQ=255)                                                                                               aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAgaga                                    <  1:950061/151‑1 (MQ=255)                                                                                               aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTTGGGTGGAAGATGCGCTGAAAGGGCAGGTGGCGGGGGCGTTTGTCATTTCGAAAGAGAGCGACAGTCTTGAAGAGCGTGATGTGGCGCACTCGCaa                                        >  2:471152/1‑147 (MQ=255)                                                                                                             ctcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGTGGGGGGGGCGCTTGTCATTCCGAAAGAGAGAGACAGTCTGGAAGATCGTGGTGTGGGGCACTCGCAAAAGAAGGCGAGTAtag                        >  1:809120/1‑147 (MQ=255)                                                                                                              tcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGGGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGGGATGTGGCGCACTCGCAAGAGAAAGCGATTATGgcg                      >  1:68692/1‑150 (MQ=255)                                                                                                                cAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAAATGCGCTGAAAGGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGgcg                      >  2:926458/1‑148 (MQ=255)                                                                                                                     aTCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGGAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGGTTATGGCGCtggggg               >  1:363739/1‑150 (MQ=255)                                                                                                                        cGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGGGATTATGGCGCtggggg               >  1:679141/1‑147 (MQ=255)                                                                                                                           gCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCa        <  1:479932/151‑1 (MQ=255)                                                                                                                                gCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTg   <  1:758359/151‑1 (MQ=255)                                                                                                                                 ccGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGCAGGGGGGGGGGGCGTTTGCCATTTCCAAAGAGAGCGACAGTCTGGAAGATCGTGGTTTGGCGCCCCCGCCAGAGAAGGGGGATTTTGCGCTGGTGGACAGCCCGAATgg  >  1:595667/1‑151 (MQ=255)                                                                                                                                   gAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGCAGGTGGGGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGGGATGTGGGGCACTCGCCAGAGAAGGCGATTATTGCGCTGGTGGGCCGCCAGaattg  >  1:937726/1‑145 (MQ=255)                                                                                                                                                      |                                                                                                                                  CTGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGG  >  NC_000913/354139‑354417

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