Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,900	0	A	G	41.2%	87.5 / 14.1	52	noncoding (72/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (3/27);  new base G (21/0);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 1.77e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGA  >  NC_000913/925756‑926048                                                                                                                                                  |                                                                                                                                                      cACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGa                                                                                                                                                  <  2:229058/150‑1 (MQ=255)   cTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAAcc                                                                                                                                               >  1:954611/1‑151 (MQ=255)      ccGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTTA‑GGGGGGGGTTCGAAccccc                                                                                                                                            >  1:646548/1‑151 (MQ=255)         tACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGGAccccc                                                                                                                                            >  1:669794/1‑148 (MQ=255)          aCAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGAt                                                                                                                                         <  2:214720/150‑1 (MQ=255)            aGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑c                                                                                                                                      <  2:370802/150‑1 (MQ=255)            aGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑Cg                                                                                                                                     <  1:567192/151‑1 (MQ=255)              gCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑CGtt                                                                                                                                   <  1:271232/151‑1 (MQ=255)              gCAGCTTAGAAATTCACAGGAAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCAAGGGGGGGGTTTCGAAcccccccaggtt                                                                                                                                     >  1:62606/1‑144 (MQ=37)                 gCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGTTCGAACCCCCCACA‑Cggt                                                                                                                                   >  1:610611/1‑146 (MQ=255)                     aGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGtat                                                                                                                            <  1:367759/151‑1 (MQ=255)                       aaaTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGTATac                                                                                                                          <  2:316119/151‑1 (MQ=255)                         aTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGGTTGAACCCCCCCTA‑CGCTGACgcctac                                                                                                                          >  2:824792/1‑144 (MQ=25)                                   aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGTATACACACTTTCCAg                                                                                                               <  2:270814/150‑1 (MQ=255)                                   aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTACGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGTATACACACTTTCCAg                                                                                                               <  2:400921/150‑1 (MQ=255)                                      aTACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAAACCCCCCTA‑CGTTGTCGTATACACACTTTTCAgg                                                                                                              >  1:218066/1‑148 (MQ=21)                                       tACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGGGTTGAACCCCCGACA‑CGCTGCCGTATACACACTTTTCAGGCGTg                                                                                                          >  1:387605/1‑151 (MQ=21)                                            cACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGTATACACACTTTCCAGGCGTGCTCCt                                                                                                     <  1:354746/151‑1 (MQ=255)                                                cGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAAACCCCCATA‑CGTTGGCGGATACACACTTTTCCAGCGTGCTCCtt                                                                                                    >  1:25300/1‑148 (MQ=18)                                                  cccACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGCGCTCTTCGTTAAATATGGCGGTGA‑GGGGGGGATTCGCACCCCCGATA‑CGTTGCCGTATGCACACTTTCCAGGCGTGCTCCtt                                                                                                    <  2:795394/146‑1 (MQ=255)                                                  cccACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGcc                                                                                               <  1:301871/151‑1 (MQ=255)                                                      ccaTGTTCACTGCCGTACAGACAGATAACATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGC‑GGGGGGGGTTTCAACCCCCGATACCTTTACCGT‑TAACCCCTTTTCAGGCGTGcccccttccaccccc                                                                                           >  2:929618/1‑140 (MQ=16)                                                            tCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTTAACCCCCCCAA‑ACTTTCCGTATACAAACCTTTCAGGCGGGCTCCCTCAACCAccc                                                                                           >  1:593259/1‑143 (MQ=17)                                                            tCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGaca                                                                                      <  1:686143/150‑1 (MQ=255)                                                               cTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGTTTGAACCCCCCCTA‑CGTTTCCGTATAAACAACTTCCAGGCGGGCTCCTTCACCCACCCgggcccccc                                                                                  >  2:859776/1‑144 (MQ=17)                                                               cTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACCTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGTATACACGCTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCt                                                                                   <  2:883467/150‑1 (MQ=255)                                                                tGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGTTCGAAACCCCCATA‑CGGTTCCGTATACACACTTTTCCGGCGGGCTCCCTCAGCCCCTccggccccc                                                                                   >  2:187369/1‑141 (MQ=17)                                                                  ccGTACAGACAGATAAAATGCG‑AAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGAAACCCCCCTA‑CGTTGCCGTATACACACATTTCAGGCGGGCTCCTTCCGCCCCCCGGGCAccccccc                                                                               >  1:829165/1‑145 (MQ=17)                                                                        agacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACACCCCATA‑CGGTTCCGTATATAAAATTTTCAGGGGGGGTTCTTTAGACACTCGGACACCCCCCCaa                                                                             >  2:337030/1‑147 (MQ=14)                                                                        agacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTg                                                                         <  1:466565/151‑1 (MQ=255)                                                                            agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGAACCCCCCATA‑CGTTGCCGTATACACACTTTTCAGGCGGGCTCCTTCAGCCACTCGGACACCCCACCAAAATGtttt                                                                     >  1:663497/1‑151 (MQ=14)                                                                              aTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGGGTTGTAACCCCCAGA‑AGTTTCCCTATAAACACTTTTCCGGCGGGGTCCTTCCGCCACTCGGACACCCCCCCCAAAtttttt                                                                     >  1:217753/1‑149 (MQ=9)                                                                                aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGGTTTGAAACCCCCATA‑AGTTTTCGGATAAAACCTTTTCAGGGGGGGTCCTTCTGCCACCCGGAGACCCCACCCAAATGtttt                                                                     >  1:886046/1‑147 (MQ=9)                                                                                           aaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCAAA‑CGTTTTCGTTTAAACACTTTTCAGGCGGGCTCCTTCTGCCCCTCCGGCACCCCACCAAAATGTTTTTCTTCCAAACCTc                                                        >  2:215878/1‑149 (MQ=14)                                                                                                   cGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCa                                              <  1:239102/151‑1 (MQ=255)                                                                                                   cGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCa                                              <  1:591942/151‑1 (MQ=255)                                                                                                     cACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAAc                                            <  1:727524/151‑1 (MQ=255)                                                                                                     cACTTTCGTACGAGCTCTTCTTGAAACATGGCGGTGA‑GGGGGGGATTCGATCCCCCGATA‑CGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAAc                                            <  2:595891/151‑1 (MQ=255)                                                                                                         ttCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACgg                                          <  1:332676/149‑1 (MQ=255)                                                                                                             tACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAACCCCCGATA‑CGATGCCGGATACAAACATTCCAGGCGTGCTCCTTCAACCACTCGGAGACCTCACCAAAATGGTTTGCTAGCAAAACTCATTGGGGGGAAAGGGgcgcg                                    >  1:817093/1‑150 (MQ=14)                                                                                                              aCGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCATA‑CGTTGCCGTATACACACCTTCCAGGCGGGCTCCCTCAGCCACTCCGACACCTCCCCCAAATGTTTTTCTTCCAAACCCCATGGGGGGCAACgggggg                                      >  1:183309/1‑146 (MQ=14)                                                                                                              aCGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGct                                    <  2:856990/150‑1 (MQ=255)                                                                                                               cGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGAACACCCGCCA‑CGTTTCCCTATAAAAACATTCCAGGCGTGCTCCCTTTACCACTCCGGCAAACCCCCCAATTTTTTTTTTAACAAACCTCATGGGGGGGAACcggggg                                      >  2:572876/1‑145 (MQ=2)                                                                                                                          ttAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGa                          <  1:598330/148‑1 (MQ=255)                                                                                                                           tAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGtt                       <  2:183249/150‑1 (MQ=255)                                                                                                                              atatGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAg                   <  1:187429/151‑1 (MQ=255)                                                                                                                                    cGGTGA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaaaa              >  2:21826/1‑150 (MQ=255)                                                                                                                                       tGA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACgg           >  2:131167/1‑150 (MQ=255)                                                                                                                                        gA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGt          <  2:62253/150‑1 (MQ=255)                                                                                                                                            ggggggATTCGAACCCCCGATA‑CGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaag      <  1:383428/151‑1 (MQ=255)                                                                                                                                             gggggATTCGAACCCCCGATA‑CGTTGCCGTATACACACGTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaaga     <  2:148728/151‑1 (MQ=255)                                                                                                                                                 gATTCGAACCCCCGATA‑CGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAagaaga  >  2:831426/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                      CACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATA‑CGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGA  >  NC_000913/925756‑926048

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