Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,764,771	0	T	A	45.5%	37.1 / 17.0	33	intergenic (‑385/+163)	sufA/ydiH	Fe‑S cluster assembly protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (16/2);  new base A (0/15);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 1.47e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGG‑CACTCTCAGTGGCTTAGCTCTTGAAGGGGC‑GATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGC  >  NC_000913/1764628‑1764921                                                                                                                                                   |                                                                                                                                                       aTGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGG‑CACTCTCAGTGGCTTAGCTCTTGAAGGGGC‑GATAAGaat                                                                                                                                                  >  2:159218/1‑150 (MQ=255)      tATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGG‑CACTCTCAGTGGCTTAGCTCTTGAAGGGGC‑GATAAGAATAAtctc                                                                                                                                            >  1:672877/1‑151 (MQ=255)           aaGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGG‑CACTCTCAGTGGCTTAGCTCTTGAAGGGGC‑GATAAGAATAATCTCATaaa                                                                                                                                       <  2:17151/151‑1 (MQ=255)             gagaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGG‑CACTCTCAGTGGCTTAGCTCTTGAAGGGGC‑GATAAGAATAATCTCATAAAGc                                                                                                                                     >  1:567510/1‑151 (MQ=255)                aTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGG‑CACTCTCAGTGGCTTAGCTCTTGAAGGGGC‑GATAAGAAGAATCTCATAAAGCTaa                                                                                                                                  <  1:276986/151‑1 (MQ=255)                 tAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGG‑CACTCTCAGTGGCTTAGCTCTTGAAGGGGC‑GATAAGAATAATCTCATAAAGCTAAc                                                                                                                                 >  1:198133/1‑151 (MQ=255)                 tAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGG‑CACTCTCAGTGGCTTAGCTCTTGAAGGGGC‑GATAAGAATAATCTCATAAAGCTAAc                                                                                                                                 >  2:217630/1‑151 (MQ=255)                  aGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGG‑CACTCTCAGTGGCTTAGCTCTTGAAGGGGC‑GATAAGAATAATCTCATAAAGCTAAcc                                                                                                                                >  2:145437/1‑151 (MQ=255)                      aaaaaaGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGG‑CACTCTCAGTGGCTTAGCTCTTGAAGGGGC‑GATAAGAATAATCTCATAAAGCTAACccgcc                                                                                                                            >  2:78243/1‑151 (MQ=255)                              gctgCGTTCAGAAATAAGAAAACGC‑TTAAGTCTCTGCGACACAGGCTTAAGGGATTCTACCCCATCCGGCGCTTATCTCCGG‑CACTCTCAGTGGCTTAGCTCTTGACGCGGC‑TATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAac                                                                                                                    >  1:176847/4‑151 (MQ=255)                                    ttCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGG‑CACTCTCAGTGGCTTAGCTCTTGAAGGGGC‑GATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAAc                                                                                                              >  1:279940/1‑151 (MQ=255)                                       aGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGG‑CACTCTCAGTGGCTTAGCTCTTGAAGGGGC‑GATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGc                                                                                                           >  2:351374/1‑151 (MQ=255)                                         aaaaaaaaaaCCCC‑TAACGTCTGGGCAACACTGGTTTAAGTTTTCCTCCCCCACCGGCCGCTTAACCCCCC‑CCCCCTCATTGTCTTCCCTCTTGAGGGCCC‑CAAAATAATAACATCATAAATCAAACCCGCCGTTTTAACACAAACTGCGa                                                                                                         <  1:18032/151‑1 (MQ=255)                                              aaaaaaCCCGTTAAGTCTGTGCAACCCAGGTTTATGGTTTTCTCCCCCCCCCGGCGTTTATCCCCGG‑CTCTCTCATGGTTTTAGCTTTTGATGGAGC‑AAAAATAATAATCTAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTa                                                                                                      <  1:538564/148‑1 (MQ=255)                                                        ttAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGG‑CACTCTCAGTGGCTTAGCTCTTGAAGGGGC‑GATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttttt                                                                                           >  2:189312/1‑151 (MQ=255)                                                          aaGTCTGTGCGACACAGGCTTAAGGTTTTCTCCCCCACCCGGCGCTTACCCCCGG‑CCCTCTCAGTGGTTTAGCTCTGGAAGGGGC‑GAAAAAAATAATCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttttt                                                                                           <  2:14956/149‑1 (MQ=255)                                                          aaGTCTGGGCGACACAGGTTAAGGGGTTTCTCCCCCACCCGGCGTTTATCCCCGG‑CACTCTCAGTGGTTTACCTTTTGAGGGGGC‑AAAAAGAAAAATCTAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTg                                                                                          <  2:150902/150‑1 (MQ=255)                                                                   ctacacAGTTTTAGGGTTTCCTCCCCCATCCGGC‑TTTATCCCCGGCCCCTCACAG‑GGCTTTACCTTTGAAAGGGCAAAAAAGAATATTATAAAAATCTTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAata                                                                                   <  1:425859/146‑1 (MQ=255)                                                                   cGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGG‑CACTCTCAGTGGCTTAGCTCTTGAAGGGGC‑GATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATc                                                                                 >  2:94428/1‑150 (MQ=255)                                                                       acgggCTTAGGGTTTTCTACCCCATCCGGCGTTTATCTCCGG‑CACTCTCAGTGTCTTAGCTCTTGAAGGGGC‑GAAAAGAATAATCTAATAAATCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCg                                                                            <  1:534094/148‑1 (MQ=255)                                                                          gttttAAGGGTTTCCCCCCCACCGGGCGTTTACCTCCGG‑CTCTCTCGGGGGTTAAGCTCTTAAGGGGGC‑AAAAAAAATATTCCAAAAAAGCAAGCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTa                                                                         <  2:672546/148‑1 (MQ=255)                                                                          ggCTTAAGGGTTTCTACCCCATCCGGCGCTTATCCCCGG‑CACTCTCAGTGGCTTAGCTCTTGAAGGGGC‑GATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGt                                                                          >  2:410756/1‑150 (MQ=255)                                                                             tgaaGGTTTTCACCCCCACCCGGCGCTTATCCCCGG‑CACTCTAAGTGGATGTGTTCTGGATGGGCA‑AAAAAGAATAATCTAAAAAATCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCg                                                                            <  1:412516/143‑1 (MQ=255)                                                                                       ctccccccaccGGCGCTTATCTCCGC‑CCCCCTAGGTGGCTTTCTTCTTGAAGGGGG‑AAAAAGAAAAATTTAATAAGGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTa                                                              <  1:145664/140‑1 (MQ=255)                                                                                          ccccATCCGGCGCTTATCCCCGG‑CCCTCTCAGTGGTTTAGCTTTTGAAGGGGC‑AAAAAAAATAATCTAATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTa                                                              <  2:418596/146‑1 (MQ=255)                                                                                            caaccgggCCCTTATCCCCGA‑CCCTCGCAGTGGTTTATCTCTTGAAGGGAC‑AAAAAAAATAATATAATAAAGATAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTaagaa                                                       <  1:66484/145‑1 (MQ=255)                                                                                                 cGGCGTTTATCCCCGG‑CCCTCTCAGTGGCTTAGCTCTTGAAGGGGC‑GAAAAAAATAATCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTc                                                   <  2:793/150‑1 (MQ=255)                                                                                                  ggCGCTTATCTCCGG‑CACTCTCAGTGGCTTAGCTCTTGAAGGGGC‑GATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCAt                                                 >  2:74132/1‑151 (MQ=255)                                                                                                        tATCTCCGG‑CACTCTCAGTGGCTTAGCTCTTGAAGGGGC‑GATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAAcc                                           >  2:185541/1‑151 (MQ=255)                                                                                                             ccag‑CACCCTAAGTTGGTTCGCTGTTGAGGGGGC‑AAAAAAAATATTTTAAAAAATAGAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAATAAGGTCGGCGAACCTAGAc                                      <  1:157432/148‑1 (MQ=255)                                                                                                                                    tctTAAGGGGGC‑AAAAAGAATAATCAAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTc                     <  2:552897/146‑1 (MQ=255)                                                                                                                                         aaGGGGA‑GAAAAAAAAATTCTCAAAAAGTTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGAttt                <  1:41520/146‑1 (MQ=255)                                                                                                                                            gggC‑GATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACt        >  1:184168/1‑151 (MQ=255)                                                                                                                                                   tAAGTATAATCTCATAAAGCTAACACGACGTTTTAACACAAACTGCGACTAGTATTATTTTTGAACAATCTCAGGCGGTAGATAAGGAGTATGAAGAAGGTCTTCGAACCTGGAAGGAGGTTAATCCAGGTAGATTTGACGAGCTTGctg   >  1:35677/1‑148 (MQ=255)                                                                                                                                                   tAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGTACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGc  >  2:236607/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                       ATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGG‑CACTCTCAGTGGCTTAGCTCTTGAAGGGGC‑GATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGC  >  NC_000913/1764628‑1764921

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