Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,900	0	A	G	45.2%	57.8 / 15.8	42	noncoding (72/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (4/19);  new base G (19/0);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 2.08e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCA‑CTCGGA‑CACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925760‑926050                                                                                                                                             |                                                                                                                                                         gCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAAcccc                                                                                                                                                <  2:70687/151‑1 (MQ=255)    gTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGa                                                                                                                                             <  1:937210/151‑1 (MQ=255)         ggCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGGGGGAAcccccccaa                                                                                                                                           >  2:222651/1‑144 (MQ=37)            aGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGGTCGAAcccccccaacaatg                                                                                                                                      >  2:634191/1‑141 (MQ=37)                  gAAATTCACAGGTAACATACTCCACCCGCTCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGtata                                                                                                                               <  2:144462/151‑1 (MQ=255)                  gAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGGTTGGAcccccccaaaggtgccggat                                                                                                                                >  2:1099955/1‑135 (MQ=37)                   aaaTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCCCTACGTTGCCGtat                                                                                                                                >  1:830468/1‑149 (MQ=255)                           aGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTc                                                                                                                      <  2:802387/151‑1 (MQ=255)                            ggTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTcc                                                                                                                     <  2:986753/151‑1 (MQ=255)                               aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAgg                                                                                                                  <  1:1005471/151‑1 (MQ=255)                                   tACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAAACCCCCCCACGTTGCCGTATACAAAATTTTCAGGCg                                                                                                                >  2:1044524/1‑149 (MQ=21)                                      tCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGGGTTGGACCCCCCACACGTTGGCGTATACACACCTTTCCGGCggg                                                                                                              >  2:1008255/1‑146 (MQ=21)                                        cACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAAACCCCCCGACGTTTCCCGATACACACTTTTCCCGCGTGCTcc                                                                                                          >  1:246750/1‑150 (MQ=21)                                             gCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGGACCCCCCCTACGTTGCCGTATAAACACTTTCCCAGGGGGCTCCTTc                                                                                                       >  2:836229/1‑148 (MQ=18)                                               ccaccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGAGTCGGAACCCCCATACGTTGCCGTATAAAAACTTTCCCCGGGTGCTCCTTCCGcc                                                                                                   >  1:879062/1‑150 (MQ=18)                                                  ccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTGAACCCCCCATACGGTGCCGTATTAAAACCTTCCCGGCGGGGTCCTTCcacca‑ccc                                                                                              >  2:764136/1‑149 (MQ=18)                                                                gTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTGCACCCCCCCCACGTTGCCGGTTACACACCTTCCCCGCGGGCTCCTTCCGCCC‑CTCGGG‑CACCCCAcc                                                                                 >  2:701357/1‑149 (MQ=17)                                                                  acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGTTTGGACCCCCCATACCTTGGCGGATACACACTTTTCACGGGTGCTCCTTCCGCCC‑CCCGGA‑CAcccccccca                                                                               >  2:856315/1‑142 (MQ=14)                                                                           tAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTTAACCCCCCACACGTTTGCCTATACACACCTTTCAAGCGGGCTCCCTCCGCCC‑CCCGGACCCCCCCCCCAAATtttttttt                                                                     >  2:441230/1‑149 (MQ=9)                                                                           tAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTGGACCCCCCCAACGGTGCCGTATACACACATTTCAGGCGTGGTCCTTTACCCA‑CCCGGA‑CACCCCACCACAATGTTTTgtt                                                                    >  2:468528/1‑149 (MQ=14)                                                                            aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCA‑CTCGGA‑CACCTCACCAAATTGTTTTGCTa                                                                   <  2:252009/151‑1 (MQ=255)                                                                             aaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGTTCGAACCCCCCCTACGTTGCCGTATAAACACCTTTCCGGCGGGCGCCTTCCGCCC‑CTCGGG‑CACCCCCCCCAAATGTTTTGCt                                                                    >  2:74673/1‑149 (MQ=14)                                                                                      aaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATAGGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCA‑CTCGGA‑CACCTCACCAAATTGTTTTGCTACCAAACCTCa                                                         <  1:630359/151‑1 (MQ=255)                                                                                        aaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTCGAACCCCCCCCACGTTGCCGTATACACACTTTTCCGGCGTGCGCCTTCCCCCACCCCGGC‑CCCCCCCCCAAATT‑TTTTTTTTCCCAACCTCAt                                                        >  2:694475/1‑150 (MQ=2)                                                                                         aaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCA‑CTCGGA‑CACCTCACCAAATTGTTTTGCTACCAAACCTCATgg                                                      <  2:54243/151‑1 (MQ=255)                                                                                              tCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCCTACGTTTCCGGATACACACTTTTCCGGCGTGCTCCCTCAACCA‑CTCGGG‑CACCCCACCAAAATGTTTTGCTACCAAACCTCCTgggggg                                                  >  1:1005334/1‑150 (MQ=14)                                                                                               cGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCA‑CTCGGA‑CACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGc                                                 <  1:1171532/150‑1 (MQ=255)                                                                                                 cACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCA‑CTCGGA‑CACCTCACCAAATTGTTTTGCTACCAACCCTCATGGGTGGCAAc                                              <  2:339345/151‑1 (MQ=255)                                                                                                  aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCATACGTTTCCGTATAAACACTTTCCCGGCGGGCTCCCTCCACCC‑CTCCGG‑CACCCCCCCCAAATGTTTTTTTTCCaaaaccccttggggg                                                  >  2:1086309/1‑134 (MQ=2)                                                                                                  aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTTGGACCCCCCCTAAGTTGCCGTATACACACTTTTCCGGCGTGGCCCTTCAGCCC‑CTCCGA‑CACCCCACCAAAATGTTTTTCTACCAAACCTCATGGGTGGCaa                                               >  1:513802/1‑149 (MQ=17)                                                                                                       cGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCA‑CTCGGA‑CACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGgc                                         <  1:24451/150‑1 (MQ=255)                                                                                                            gAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGGACCCCCCATACGTTGCCGTATACAAACTTTTCCGGGGGGCCCCCTCAGCCC‑CTCGGA‑CACCCCACCAAAATTTTTTTCTACCAAAACTCCTGGGGGGGAACGGGGGGctac                                    >  1:788509/1‑150 (MQ=9)                                                                                                            gAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCA‑CTCGGA‑CACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGctact                                   <  2:240940/151‑1 (MQ=255)                                                                                                                 cttcttTAAATATGGCGGTGAGGGGGGGGGTCGCACCCCCCATACCTTGCCCTATACACACTTTTCAGGCGTGCCCCCTCTGCCA‑CTCCGA‑CACCTCCCCCAAATTTTTTTTTTCCCAACCTCATTGGGGGCAACgggg                                          >  2:259904/1‑139 (MQ=9)                                                                                                                        aaaTATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCA‑CTCGGA‑CACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTg                        <  2:359663/150‑1 (MQ=255)                                                                                                                        aaaTATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCA‑CTCGGA‑CACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTg                        <  1:856321/150‑1 (MQ=255)                                                                                                                        aaaTATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCA‑CTCGGA‑CACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTg                        <  1:348734/150‑1 (MQ=255)                                                                                                                              ggCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCA‑CTCGGA‑CACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaaa                 <  1:37767/151‑1 (MQ=255)                                                                                                                              ggCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCA‑CTCGGA‑CACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaa                  >  1:210976/1‑150 (MQ=255)                                                                                                                                 ggTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCA‑CTCGGA‑CACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACg              <  1:830698/151‑1 (MQ=255)                                                                                                                                       ggggggATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCA‑CTCGGA‑CACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa          >  2:611517/1‑149 (MQ=255)                                                                                                                                          gggATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCA‑CTCGGA‑CACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaag     <  2:408252/151‑1 (MQ=255)                                                                                                                                             aTTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCA‑CTCGGA‑CACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  <  1:33002/151‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                         GCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCA‑CTCGGA‑CACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925760‑926050

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