Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,704,041	0	G	T	33.3%	54.0 / 10.1	33	H90N (CAT→AAT)	ydgJ	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (14/8);  new base T (0/11);  total (14/19)
Fisher's exact test for biased strand distribution p-value = 4.91e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTAT‑CAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTG  >  NC_000913/1703891‑1704187                                                                                                                                                        |                                                                                                                                                   cTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTAT‑CAACGACCACATg                                                                                                                                                    >  2:50946/1‑151 (MQ=255)    ccGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTAT‑CAACGACCACATGtt                                                                                                                                                  <  1:389852/150‑1 (MQ=255)      gcgAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTAT‑CAACGACCACATGTTTAc                                                                                                                                               >  2:356080/1‑151 (MQ=255)      gcgAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTAT‑CAACGACCACATGTTTAc                                                                                                                                               >  1:79349/1‑151 (MQ=255)                          cAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTAT‑CAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTg                                                                                                                           >  2:58852/1‑151 (MQ=255)                           aaGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTAT‑CAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTgg                                                                                                                          >  2:197958/1‑151 (MQ=255)                           aaGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTAT‑CAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTgg                                                                                                                          >  2:142099/1‑151 (MQ=255)                            aaaaaTCGCCACCCCGGCGACGTTTAGGGAAAACGGACACCACACCCCCCAGGCTTTGTGCCAGCGCTCCCACCGCTCGCGCTTGTGACGGGGTCA‑GGTAAGGGGTTTTTGCACCGCCCAAATGTGTACCCGCCTCAAGCGCCGCTTTGGc                                                                                                                         <  1:339602/149‑1 (MQ=255)                            aGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTAT‑CAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGc                                                                                                                         >  1:26803/1‑151 (MQ=255)                               aaccGCTACCCCAGCGAGGTTTAGGAAATACAGCCACCCCCCCCCCCAGTTTTTTTGCCAGCGCATCCAGCCCTCGCGTTTGTGACAGTTCCAGGGAAAGGGGTTTAT‑CAACGCCCACATGTTTCCCCGCCTCAAGCGCCGCTTTGGCTaa                                                                                                                      <  2:270237/148‑1 (MQ=255)                               aaTCCCTACCCCAGCGGCGGTTAGGGAATACAGACACCCCCCCCCCCAGGTTTTTTGCCGGCCCATCCGGCTCCCGCGTTTGTGACAGTGTCAGGGTAAAGGTTTTAT‑CAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTaa                                                                                                                      <  2:116554/151‑1 (MQ=255)                                  cattaccccGGCGAGGGTTATGAAACACAGACAGCCCCCGCCCCAGGTTTTTCGCCAGCGCACCCAGCCCTCGCGTTTGGGACGGTGCCACGGAAAGGGTTTTAA‑CAACGCCCAATTTTTTACCCGCCTCAAGCGCCGCTTTGGCTAACg                                                                                                                    <  1:179282/144‑1 (MQ=255)                                    caaccccAGCGACGGTTATGGAATACAGACACCCCCCCCCCCAGGTTTTTTGCCAGCGCACCCAGCTCCCGCGTTTGTGACAGTGTCACGGTAAAGGGTTTAC‑AAACGACCAATTTTTTACCCGCCTCAAGCGCCGCTTTGGCTAACggg                                                                                                                  <  2:295065/146‑1 (MQ=255)                                     tATCCCAGCGACGGTTATGGAATACAGACAGCCCACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTAT‑CAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACgg                                                                                                                   >  1:126819/1‑148 (MQ=255)                                              gACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTAT‑CAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATc                                                                                                       >  2:62767/1‑151 (MQ=255)                                                   ttATGAAATACAGACAGCCCCCCCCCCGGGTTTTTTGCCCGCCCATCCAGCTCTCGCGCTTGGGACATTGTCAGGTTAAGGGTTTTAT‑CAACACCCACATTTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTg                                                                                                   <  2:339877/150‑1 (MQ=255)                                                          aTACAGACACCACCCGCCCCAGGCTTTTGGCCAGCACCTCCAGCTCCCGGTTTTGGGACAGTGTCAGGTTAAGGGGTTTAT‑AAGCAACCATATTTTTCACCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAgg                                                                                           <  1:287279/151‑1 (MQ=255)                                                                acacCACCCCCCCCAGGTTTTTTGCCAGCCCACCCAGCTCTCGCGCTTGGGACAGTGTCGCGGTAAGGGTTTTAT‑CAACGACCACATTTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATg                                                                                       <  2:14976/149‑1 (MQ=255)                                                                  agCACACGCCCCAGGCTTTTTGCCAGCACATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTAT‑CAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACa                                                                                    >  2:84359/1‑150 (MQ=255)                                                                           cccGGGTTTTTGCCCACCGCACCCACCTCCCGCGTTTGGGACGGTGTCACGGAAAGGGTTTTAT‑AAACGCCCACTTGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAgg                                                                              <  2:372641/147‑1 (MQ=255)                                                                                 cTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTAT‑CAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTgg                                                                    >  1:45361/1‑151 (MQ=255)                                                                                              cacccAGCTCTCGCGTTTGTGACAGTGTCACGGTAAGGGTTTTAT‑CAACGCCCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACa                                                        <  2:234777/147‑1 (MQ=255)                                                                                                     ctctCGCGTTTGGGACGGTGCCAGGGAAAAGGTTTTAC‑CAACGCCCCCATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCt                                                  <  2:18847/149‑1 (MQ=255)                                                                                                      tctcGCGCTTGTGACAGTGTCACGGTAAAGGGTTTAT‑CAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCg                                               >  1:315596/1‑151 (MQ=255)                                                                                                             cTTGTGACAGTGTCACGGCAAAGGGTTTAT‑CAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATCGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTc                                           >  1:224476/1‑148 (MQ=255)                                                                                                                gtgAAAGTGTCACGGTAAGGGTTTTAT‑CACCGACCACTTTTTTCCCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCt                                            <  2:386899/144‑1 (MQ=255)                                                                                                                   aCAGTGTCACGGTAAAGGGTTTAT‑CAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACaac                                  >  2:8790/1‑151 (MQ=255)                                                                                                                                   gggTTTAT‑CAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGc                     >  1:91773/1‑148 (MQ=255)                                                                                                                                   gggTTTAA‑CAACGCCCACATTTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCa                   <  1:195892/150‑1 (MQ=255)                                                                                                                                      tttAT‑CAAAGACCACTTTTTTGCCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTc                <  1:284853/150‑1 (MQ=255)                                                                                                                                      tttAC‑CAACGCCCACATTTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGt                 <  2:376006/149‑1 (MQ=255)                                                                                                                                       ttAT‑CAACCACCACATTTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCg               <  2:251339/150‑1 (MQ=255)                                                                                                                                                  ccACTTTTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACtt     <  2:228975/150‑1 (MQ=255)                                                                                                                                                     caTGTTTCCCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTg  <  2:51847/150‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                   CTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTAT‑CAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTG  >  NC_000913/1703891‑1704187

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