Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,786,567	0	T	G	51.7%	31.2 / 35.8	29	G436G (GGT→GGG)	gpmM	phosphoglycero mutase III, cofactor‑independent
Reads supporting (aligned to +/- strand):  ref base T (6/8);  new base G (15/0);  total (21/8)
Fisher's exact test for biased strand distribution p-value = 7.00e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-07

TCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACT  >  NC_000913/3786419‑3786716                                                                                                                                                     |                                                                                                                                                      tCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCgtgggtgg                                                                                                                                                     <  1:218925/151‑1 (MQ=255)   aaGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGAc                                                                                                                                                   <  1:133907/151‑1 (MQ=255)          gCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAActgc                                                                                                                                             <  1:144011/150‑1 (MQ=255)          gCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAActgc                                                                                                                                             <  1:153582/150‑1 (MQ=255)              aTACGACACCAGCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCAGTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATc                                                                                                                                        <  1:450190/151‑1 (MQ=255)                        atcatcTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGGAGCGCTGGGTCACTGCGTGGGAGAAGGCGCGAAAACGGTTGAGTCCGTGGGGGGGCAACTGCTGATCACCGCCGAc                                                                                                                               >  1:405757/1‑150 (MQ=255)                        atcatcTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGAcc                                                                                                                              <  1:357066/151‑1 (MQ=255)                           atcTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGGTGAGTCCGTGGGGGGACAACTGCTGATCACCGCTGACCACg                                                                                                                           >  1:441372/1‑151 (MQ=255)                           atcTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGGTGAAGCGCTGGATCACTGCGTGGGAGAAGTCGCGAAAGCGGGTGAGTCCGTGGGGGGACAACTGCTGATCACCGCTGACCACg                                                                                                                           >  1:99118/1‑151 (MQ=255)                           atcTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGGTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGGGGACAACTGCTGATCACCGCTGACCACg                                                                                                                           >  1:46350/1‑151 (MQ=255)                               gTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGGTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAACGGGTGAGGCCGTGGGGGGACAACTGCTGATCACCGCTGAACACGGTaa                                                                                                                       >  1:14216/1‑151 (MQ=255)                                     aTCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGGTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGGTGAGTCCGTGGGGGGACAACTGCTGATCACCGCTGACCACGGTAACGCTgg                                                                                                                 >  1:194870/1‑150 (MQ=255)                                                   aTGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGGTGAAGCGCTGGGTCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGGGGGCAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGAt                                                                                                    >  2:23794/1‑150 (MQ=255)                                                      gTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGGAGAAGTCGCGAAAACGGGTGAGTCCGTGGGGGGACAACTGCTGATCACCGCTGACCCCGGGAACGCTGAGCAGGTGCGCGGTCCgg                                                                                                >  1:352153/1‑151 (MQ=255)                                                      gTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGGGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCgg                                                                                                >  1:128129/1‑151 (MQ=255)                                                      gTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAACGGTTGAAGCGCTGGAACACTGCGTGGGAGAAGGCGCGCAAGCGGGTGAGTCCGGGGGGGGGCAACTGCTGATCACCGCCGAACCCGGGAACGCTGAGCAGATGCGCGATCCgg                                                                                                >  1:50291/1‑151 (MQ=255)                                                        aGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGGTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGGGGACAACTGCTGATCACCGCTGACCACGGGAACGCTGAGCAGATGCGCGATCCGGCa                                                                                              >  2:324979/1‑151 (MQ=255)                                                          gTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGGAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGGGGACAACTGCTGATCACCGCTGACCACGGGAACGCTGAGCCGATGCGCGATCCGGCaaa                                                                                            >  2:262365/1‑150 (MQ=255)                                                          gTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCaa                                                                                             <  2:170413/150‑1 (MQ=255)                                                              cacCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGGCGCGAAAGCGGGTGAGTCCGTGGGGGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACgggg                                                                                        >  1:333207/1‑150 (MQ=255)                                                              cacCGGGGTGATGGAAGCGGCGGTTAAAGCGGGTGAAGCGCTGGGACACTGCGTGGAAGAAGTCGCGAAAACGGGTGAGTCCGGGGGGGGGCAACTGCTGGTCACCGCCGACCCCGGGAACGCTGAGGAGATGCGCGATCCGGGAACgggg                                                                                        >  1:320934/1‑150 (MQ=255)                                                                       gATGGAAGCGGCGGTTAAAGCGGGTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGGTGAGTCCGTGGGGGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAggcaca                                                                                >  1:313343/1‑150 (MQ=255)                                                                               cggcggTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAggcacacacggcac                                                                        <  1:305740/150‑1 (MQ=255)                                                                                ggcggTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAggcacacacggcaca                                                                       >  2:170504/1‑150 (MQ=255)                                                                                                   gcgcTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCg                                                    >  2:293136/1‑150 (MQ=255)                                                                                                         gATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGAttt                                             >  1:299134/1‑151 (MQ=255)                                                                                                             aCTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGt                                         <  2:405768/151‑1 (MQ=255)                                                                                                                     gaagaaGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGGGATa                                 >  1:75310/1‑151 (MQ=255)                                                                                                                                          gAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAg            >  1:308951/1‑151 (MQ=255)                                                                                                                                            gTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAgg           >  1:193050/1‑150 (MQ=255)                                                                                                                                                gtgggtggACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCa      >  1:297658/1‑151 (MQ=255)                                                                                                                                                    gtggACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAAc   <  2:78120/150‑1 (MQ=255)                                                                                                                                                     tggACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACt  <  2:121178/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      TCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACT  >  NC_000913/3786419‑3786716

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