Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,535,785	0	A	G	45.2%	52.4 / 14.1	42	F25S (TTC→TCC)	yddE	PhzC‑PhzF family protein
Reads supporting (aligned to +/- strand):  ref base A (2/21);  new base G (19/0);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 9.40e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGG  >  NC_000913/1535643‑1535931                                                                                                                                               |                                                                                                                                                   gCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGaaaaaa                                                                                                                                              >  1:431704/1‑149 (MQ=255)    aaTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCg                                                                                                                                         <  1:206965/151‑1 (MQ=255)         gCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGaaaaaa                                                                                                                                              >  2:254062/1‑141 (MQ=255)          cACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAAcccg                                                                                                                                        >  2:545240/1‑144 (MQ=255)              tCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAACAAAACCCGccggaat                                                                                                                                 >  2:435721/1‑143 (MQ=255)               cAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGcc                                                                                                                              <  2:193698/151‑1 (MQ=255)                 aCCGGTGGCGTCAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAAGAACAACGCCaaaaatttcgt                                                                                                                            >  2:195640/1‑144 (MQ=255)                         cGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAACAAAACCCGCAGAATTGCCGCGAAACggg                                                                                                                    >  1:38611/1‑150 (MQ=255)                          gTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAAACCGCAGAATTTGCGCGCAAAgggt                                                                                                                   >  2:359312/1‑149 (MQ=255)                                tAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGCAGAAATGCCGCGAAACGGGTGTgggg                                                                                                             >  1:505074/1‑151 (MQ=255)                                 aGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAgg                                                                                                             <  2:591494/150‑1 (MQ=255)                                 aGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGt                                                                                                            <  1:517875/151‑1 (MQ=255)                                  gCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTa                                                                                                           <  1:560799/151‑1 (MQ=255)                                        gcgcACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGc                                                                                                      <  2:76214/150‑1 (MQ=255)                                          gcACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATc                                                                                                   <  2:307840/151‑1 (MQ=255)                                           cacaTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCg                                                                                                  <  1:185275/151‑1 (MQ=255)                                             caTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGAc                                                                                                <  2:581326/151‑1 (MQ=255)                                                cGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGCAAGATTGGCGCGCAAAGGGTGTGGGGTgaaaaaaaccaccgg                                                                                             >  2:7910/1‑140 (MQ=255)                                                    gtcgtcGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTa                                                                                          <  1:584122/150‑1 (MQ=255)                                                      cgtcgCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGCAGAAATTCCGCGAAAAGGGTGTGGGGGAAAAgaaacaacgggg                                                                                            >  1:615025/1‑135 (MQ=255)                                                         cgCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTg                                                                                    <  1:131415/151‑1 (MQ=255)                                                          gCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGc                                                                                   <  2:159052/151‑1 (MQ=255)                                                             gtgCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGt                                                                                <  1:342166/151‑1 (MQ=255)                                                               gcagcaGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGtt                                                                               <  1:385287/150‑1 (MQ=255)                                                                    agAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAAACGGGAGAATTGCCCCGGAACGGGTGTGGGGGGAAAAAATCCACGGGGGACACCCGCGGgtt                                                                              >  2:514968/1‑146 (MQ=255)                                                                     gAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTaa                                                                        <  2:3605/151‑1 (MQ=255)                                                                       aaaGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGCAGAATTGCCGCGGAACGGGTGTGGGGGGAAAACAACGACGCGGGACACCTGCGGGTTTAttta                                                                        >  2:122653/1‑147 (MQ=255)                                                                          gCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTtc                                                                   <  2:630386/151‑1 (MQ=255)                                                                            aGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAACAAAAACCGCAGAATTGCCCCCGAACGGGTGTGAGGTGAAAGAAACAACGCGGTACACCTCCGGGTTCCTTAACGTtct                                                                  >  2:434950/1‑150 (MQ=255)                                                                                 ccGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGCAGAAATGCCGCGCAACGGGTGTGGGGGAAAAACAACGACGTGGGACACCCCCGGGTTTTTTTACGtttt                                                                  >  2:498203/1‑143 (MQ=255)                                                                                     aTGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAACAAAAACCGCCGAAATGCCGCCGAACGGGTGTGAGGGAAAAACATCGACGTGGGAAACCCCCGGGTTTATTAAAGGTCTCCtgt                                                             >  2:334196/1‑146 (MQ=255)                                                                                           tAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAAACCGCAGAATTGGCGCGAAACGGGTGTGAGGGAAAAACAACGACGGGGGACACCCGCGGGTTTATTAAAGGTCTCCTGTGACTGGGGaa                                                   >  2:325776/1‑150 (MQ=255)                                                                                                cgcgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAAACCGCCGAAATGCCGCGGAAAGGGTGTGGGGGAAAAAGAACGACGGGGGAGAACTGCGGGTTCATTAAAGTTTTCCTGGGACTgggggagta                                                >  2:306004/1‑141 (MQ=255)                                                                                                  cgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATa                                           <  1:557342/151‑1 (MQ=255)                                                                                                   gggCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGGAAAAAACCCGGCCGAAATTCCGCGAAAAAGGTGTGGGGGAAAAGCATACACGTGGTACACCTGCGGTTTCCTTAAAGATTTCCTTGGTCTGGAGAACTATCa                                             >  2:560360/1‑148 (MQ=255)                                                                                                          aaGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAACACCCGCAGAAATGCCGCGAAACGGGTGTGGGGGGAAAACAACGACGTGGGACACCCGCGGGTTTATTAAAGTTCTCCCGTGACTGGGGAAATATCATAAcc                                        >  1:220588/1‑146 (MQ=255)                                                                                                          aaGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCaa                                   <  1:545251/151‑1 (MQ=255)                                                                                                               gCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGc                              <  2:628873/151‑1 (MQ=255)                                                                                                                       gcTTCGCTAAGATTATCGGCGGGGAAAACAACACCCCGAGAAAATCCGCGAAACGGGGTTGGGGGAAAAAGAACCACGCGGGACAACCCCGGGTTTCTTAAACATCTCCTCTGACTGGGGAACTATCATAACCTCCACGAGGACggg                          >  2:192441/1‑146 (MQ=255)                                                                                                                        cTTCGCTAAGATTATCGGCGGGGGGAAACACACCCGCCGAATTGTCCCGGAAAGGGTGTGGGGGTAAAACAACGACCCGGGACACCTGCCGGTTTATTTAAGTTCTCCTGTGGCCGGGGAACCATCCTATCCTTCCAGTTGcc                             >  2:378860/1‑143 (MQ=255)                                                                                                                                    tATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATcc          <  1:345650/150‑1 (MQ=255)                                                                                                                                           ggggAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAggg  <  2:106807/151‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                   GCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGG  >  NC_000913/1535643‑1535931

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