Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,477,016	0	T	A	21.1%	74.6 / 12.7	39	Q205L (CAG→CTG)	dsdC	dsd operon activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (14/16);  new base A (5/3);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCA  >  NC_000913/2476869‑2477156                                                                                                                                                    |                                                                                                                                             cccATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCGATTCATCCGTACCGGAGTCGTTGCTCCATGCct                                                                                                                                              <  2:423/148‑1 (MQ=255)                         aTCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGt                                                                                                                      >  7:67826/1‑148 (MQ=255)                                          tAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACAGCATAATGTTGCGCCCAACTATGACATTCATCCGTACCGGATTCGTTGCTCCATGCCTGTCTGGAAATGAGCAACGTACAGTGACACAGTTTAATTAcag                                                                                                    >  7:46865/1‑147 (MQ=255)                                             aTCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGACAGACAGGCATGGAGCAAc                                                                                                                             >  8:96853/1‑121 (MQ=39)                                             aTCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGACAGACAGGCATGGAGCAAc                                                                                                                             <  7:96853/121‑1 (MQ=39)                                               cAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTACATGCCTGTCAGTCATGGAGCAAAGTACAGTGAAACCGGTTAAATACCGTGTTg                                                                                               >  4:67219/1‑149 (MQ=255)                                                       tCAAAGCCAATTCCAGAAGATGTTGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAg                                                                                        >  8:35408/1‑148 (MQ=255)                                                       tCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTTCTCCCAGGCGGTGTGTCAAGGAGGAACGGACAGTGACACAGGGTAATTACCgt                                                                                                   >  2:117725/1‑137 (MQ=255)                                                       tCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCgt                                                                                                   <  1:117725/137‑1 (MQ=255)                                                         aaaGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTg                                                                                               <  4:114714/139‑1 (MQ=255)                                                         aaaGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTg                                                                                               >  3:114714/1‑139 (MQ=255)                                                                                gCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCgtg                                                                                                  >  6:86244/1‑113 (MQ=255)                                                                                gCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCgtg                                                                                                  <  5:86244/113‑1 (MQ=255)                                                                                  aaaTTAACCGCATAATGTTGCGCCAAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTATTTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACa                                                            <  8:46865/149‑1 (MQ=255)                                                                                   aaTTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACa                                                            >  7:43343/1‑148 (MQ=255)                                                                                          cGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCa                                                    <  3:67219/149‑1 (MQ=255)                                                                                           gCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAg                                                   <  7:35408/149‑1 (MQ=255)                                                                                             aTAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTc                                                                     <  3:42105/129‑1 (MQ=255)                                                                                             aTAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTc                                                                     >  4:42105/1‑129 (MQ=255)                                                                                                           aaCTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTaaa                                                                                         >  3:18871/1‑95 (MQ=255)                                                                                                           aaCTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTaaa                                                                                         <  4:18871/95‑1 (MQ=255)                                                                                                                    aTTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCAt                                                                                     <  7:74549/90‑1 (MQ=255)                                                                                                                    aTTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCAt                                                                                     >  8:74549/1‑90 (MQ=255)                                                                                                                    aTTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGc                                                                           >  3:153433/1‑100 (MQ=255)                                                                                                                    aTTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGc                                                                           <  4:153433/100‑1 (MQ=255)                                                                                                                    aTTCATCCGTACCGGAGTCGTTGCACCAATCCAGGATGTCATGGATCAAACTAACCTGACACAGGTTAATTACCCTGTTGTTTAAAACAAGACATTGACCGTATTACGAGCTAAATAATTCCACGAATTCTTAATCCAACAGAAAGtg                           >  4:15895/1‑148 (MQ=255)                                                                                                                       cATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtgatg                        <  8:67826/148‑1 (MQ=255)                                                                                                                         tCCGTACCGGAGTCGTTGCTCCATGGCAGTCGGTCATGGAGCAACCTACAGTGCCAAAAGTTAATTACCGTGTTGGGTAAAGCATGTCGTTTAGCGTATTCCGGGCAACATACTTGCAGGATTTCTTAATCCATCAGAAAGTGAtgagt                     >  4:48051/1‑149 (MQ=255)                                                                                                                                  gAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGTTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgag                      <  7:62679/139‑1 (MQ=255)                                                                                                                                  gAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgag                      <  7:133167/139‑1 (MQ=255)                                                                                                                                     tCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGa         <  8:43343/149‑1 (MQ=255)                                                                                                                                       gTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGAc        <  3:15895/148‑1 (MQ=255)                                                                                                                                       gTTGCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtga                          <  1:109699/130‑1 (MQ=39)                                                                                                                                       gTTGCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtga                          >  2:109699/1‑130 (MQ=39)                                                                                                                                         tGCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCa                                                    <  1:173629/102‑1 (MQ=38)                                                                                                                                         tGCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCa                                                    >  2:173629/1‑102 (MQ=38)                                                                                                                                            tCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca  >  5:1265/1‑149 (MQ=255)                                                                                                                                              cATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtga                          >  6:37430/1‑123 (MQ=255)                                                                                                                                              cATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtga                          <  5:37430/123‑1 (MQ=255)                                                                                                                                               aTGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGc                                                                >  6:111277/1‑84 (MQ=255)                                                                                                                                                    |                                                                                                                                             CCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCA  >  NC_000913/2476869‑2477156

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