Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,535,785	0	A	G	46.2%	40.7 / 15.9	39	F25S (TTC→TCC)	yddE	PhzC‑PhzF family protein
Reads supporting (aligned to +/- strand):  ref base A (4/17);  new base G (18/0);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 1.81e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAAC‑AC‑CGGCAGAATTGCCG‑CGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAA  >  NC_000913/1535635‑1535935                                                                                                                                                       |                                                                                                                                                          gcgTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGa                                                                                                                                                           <  1:435321/151‑1 (MQ=255)         gCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAacaac‑a                                                                                                                                                  <  1:48528/151‑1 (MQ=255)            aaTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAAAAAC‑Ac‑c                                                                                                                                               >  1:391091/1‑150 (MQ=255)              tCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGaaaaaaaaac                                                                                                                                                 >  2:170304/1‑144 (MQ=255)                 gCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGGAaaaac‑aac                                                                                                                                                >  2:469374/1‑143 (MQ=255)                     ttCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGTGAGAACAAC‑Ac‑cccccg                                                                                                                                          >  1:779671/1‑141 (MQ=255)                      tCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGaaaaaaaccccccgcga                                                                                                                                          >  2:82892/1‑136 (MQ=255)                         aCCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAA‑AA‑CGGCAGAATTGcc                                                                                                                                   >  1:625539/1‑149 (MQ=255)                         aCCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAAC‑AC‑CGGCAGAATTGCcg‑c                                                                                                                                <  1:506090/151‑1 (MQ=255)                          ccGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAAC‑AC‑CGGCAGAATTGCcg‑cg                                                                                                                               <  1:370447/151‑1 (MQ=255)                                   tAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAA‑AA‑CGGCAGaaat                                                                                                                                      >  1:951067/1‑134 (MQ=255)                                    aaaGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAC‑AC‑CCGCCGGATTTGCG‑CGGAAAgggt                                                                                                                       >  2:204178/1‑147 (MQ=255)                                           cGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAAC‑AC‑CGGCAGAATTGCCG‑CGAAACGGTTGTGAGGTaa                                                                                                              <  1:369298/151‑1 (MQ=255)                                             gATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAAC‑AC‑CGGCAGAATTGCCG‑CGAAACGGTTGTGAGGTaaaa                                                                                                            <  1:149806/151‑1 (MQ=255)                                                  gcACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGGAAAAAA‑AC‑CCGCCGACATGCCGCCGAAAC‑GTTGTGGGGTgaa                                                                                                             >  2:587705/1‑145 (MQ=255)                                                  gcACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAAC‑AC‑CGGCAGAATTGCCG‑CGAAACGGTTGTGAGGTAAAAGCATc                                                                                                       <  2:451467/151‑1 (MQ=255)                                                   cacaTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAAC‑AC‑CGGCAGAATTGCCG‑CGAAACGGTTGTGAGGTAAAAGCATCg                                                                                                      <  2:567000/151‑1 (MQ=255)                                                      aTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAC‑ACACGGC‑GGATTGACG‑CGAAACGGGTGGGAGGGGAAAACAtgca                                                                                                     >  1:815505/1‑146 (MQ=255)                                                            gtcgtcGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAC‑AA‑CAGACGGATTTGAG‑CGGAAAGGGTGTGAGGGaaaaaaaactacgtg                                                                                                 >  2:355330/1‑136 (MQ=255)                                                                 cgCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAAC‑AC‑CGGCAGAATTGCCG‑CGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTg                                                                                        <  1:47787/151‑1 (MQ=255)                                                                       gcagcaGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAAC‑AC‑CGGCAGAATTGCCG‑CGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGtt                                                                                   <  2:388776/150‑1 (MQ=255)                                                                         agcagAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAAC‑AC‑CGGCAGAATTGCCG‑CGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCa                                                                                <  2:841353/151‑1 (MQ=255)                                                                                    aGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAACAACAAC‑CGCCGGAATGGCC‑‑CGGAACGGGTGGGGGGGAAAAACAACGAAGTGGGACACCTGCGGGTTTATTAAAATtct                                                                      >  1:760772/1‑150 (MQ=255)                                                                                           gAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAAC‑AC‑CGGCAGAATTGCCG‑CGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGAc                                                              <  1:89553/151‑1 (MQ=255)                                                                                            aaTGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAA‑AC‑CGGCAGAAATTCCG‑CGGAAAGGGTGTGGGGGAAAAACAACGACGTGGGACACCTGCGGGTTTATTAAAATTCTCCtgt                                                                 >  1:433447/1‑147 (MQ=255)                                                                                             aTGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAA‑AC‑CCGCAGAATTGCCG‑CGAAACGGGTGTGGGGGAAAAAAATCGACGTGGGACACCTGCGGGTTTATTAAAGTTCTCCTGTGACTg                                                            >  1:264397/1‑151 (MQ=255)                                                                                                gCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAA‑AC‑CCCCAGAATTGCCG‑CGCCACCGGGGTGGGGGGAAAGCATAGACGTGGTACACCTGCGGCGTTATTAAAATTCTCCTCTGACTgg                                                           >  1:204740/1‑149 (MQ=255)                                                                                                gCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAAC‑AC‑CGGCAGAATTGCCG‑CGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGgag                                                         <  2:632819/151‑1 (MQ=255)                                                                                                     aCTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAAC‑AC‑CGGCAGAATGGCCG‑CGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGATAACt                                                     <  1:730127/150‑1 (MQ=255)                                                                                                        cgcgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAACAAC‑AC‑CCGCCGAATTGCCG‑CGAAACGGGTGTGAGGGAAAAACAACGACGTGGGACACCCGCGGGTTTATTAAAGTTCTCCTGTGACTGGGGaa                                                       >  1:128582/1‑145 (MQ=255)                                                                                                              cGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAA‑AC‑CCGCAGAAATTTCC‑CGGAAAGGTTGGTGGGGAAAAACAAAGACGGGGGACACCCGCGGGTTTATTAAAATTCTCCTGTGGCTGGGGGAATATCATAAcc                                            >  2:316937/1‑150 (MQ=255)                                                                                                               gATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAAC‑AC‑CGGCAGAATTGCCG‑CGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTg                                          <  2:747958/151‑1 (MQ=255)                                                                                                                 tAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAA‑AA‑CGGCAGAAATGCCG‑CGGAACGGGTGTGAGGGAAAAACATCGACGGGGGACACCTGCGGGTTTATTTAAGTTCTCCTGTGACTGGAGAAATATCATAAcc                                            >  1:303970/1‑147 (MQ=255)                                                                                                                         aTTTGCGCTTCGCTAAGATTATCGGCGGGGGAAACAAC‑AC‑CGGCAGGATTGTCG‑CGAAAAGGGTGTGGGGGAAAAACATCGACGGGGGACACCTGCGGGTTCATTAAAGGTCTCCTGTGACTGGAGAAATATCATAGCCCGCAAGGGGCCg                                >  2:117077/1‑151 (MQ=255)                                                                                                                           ttGCGCTTCGCTAAGATTATCGGCGGGGAAAACAAC‑AC‑CGGCAGAATTGCCG‑CGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCgg                               <  1:82898/150‑1 (MQ=255)                                                                                                                           ttGCGCTTCGCTAAGATTATCGGCGGGGAAAACAAA‑AC‑CCGCCGAAATGCCG‑CGGAAAAGGTGTGAGGGAAAAAAAACGACGGGGGACACCCGCGGGTTCATTAAAGATCTCCTGTGGCTGGGGGACAATAATAACCTGCCAGTGGGCgg                               >  2:111883/1‑150 (MQ=255)                                                                                                                                 ttCGCTAAGATTATCGGCGGGGGAAAAAAA‑AC‑CGGCAGAAATGCCG‑CGAAACGGGTGGGGGGGAAAAACATCCACGCGGGACACCCCCGGGTTCATTAAAAATCTCCTTTGGCTGGGGAAATATCATAACCCGCAAGGGGCCGGggag                           >  1:794992/1‑148 (MQ=255)                                                                                                                                                 gcggGGAAAACAAC‑AC‑CGGCAGAATTGCCG‑CGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGcc          >  2:794542/1‑149 (MQ=255)                                                                                                                                                 gcggGGAAAACAAC‑AC‑CGGCAGAATTGCCG‑CGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAg        <  2:669651/151‑1 (MQ=255)                                                                                                                                                     ggAAAACAAC‑AC‑CGGCAGAATTGCCG‑CGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACGCCTGCGGTTTCATTAACGTTCGCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTg    <  2:557639/151‑1 (MQ=255)                                                                                                                                                       aaaaCAAC‑AC‑CGGCAGAATTGCCG‑CGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGaa  <  1:53285/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                          GCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAAC‑AC‑CGGCAGAATTGCCG‑CGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAA  >  NC_000913/1535635‑1535935

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