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

GCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑C‑GAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAG  >  NC_000913/925760‑926047                                                                                                                                             |                                                                                                                                                      gCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑C‑GAAcccc                                                                                                                                           <  1:162884/151‑1 (MQ=255) gCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑C‑GAAcc                                                                                                                                             <  1:50886/149‑1 (MQ=255)     tACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑C‑GAACCCCCGa                                                                                                                                        <  1:220067/150‑1 (MQ=255)         ggCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑C‑GAACCCCCGATACGt                                                                                                                                   <  2:809330/151‑1 (MQ=255)           cAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑C‑GAACCCCCGATACGTTg                                                                                                                                 <  1:882359/151‑1 (MQ=255)               ttAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑T‑GAACCCCCCATACGTTGcc                                                                                                                               >  1:1017757/1‑149 (MQ=255)                tAGAAATTCACAGGTAACATACTCCACCCGCCCACCATATTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑C‑GAACCCCCGATACGTTGCCGt                                                                                                                             <  2:488159/150‑1 (MQ=255)                     aTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGT‑C‑GAACCCCCCATACGTTGCCGTATaca                                                                                                                        >  2:424262/1‑150 (MQ=255)                        cacaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑C‑GAACCCCCGATACGTTGCCGTATACACACt                                                                                                                    <  1:56101/151‑1 (MQ=255)                                      tCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGAT‑C‑GCACCCCCGATACGTTTTCCTATACACACTTTTCAGGCgag                                                                                                         >  2:320386/1‑146 (MQ=21)                                           ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑C‑GAAACCCCCGTACGGTTCCGTATACACACTTTTCAGGCGGGCTCCtt                                                                                                   >  2:899036/1‑149 (MQ=18)                                             gCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGGT‑C‑GAAACCCCCATACGTTGTCGTATACAAACTTTCCCGGGGTGccccct                                                                                                   >  1:388527/1‑145 (MQ=18)                                             gCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCC‑GAACCCCC‑CTACGTTTCCGTATACACACTTTTCCGGCGGGCCCCCTCAGc                                                                                               >  1:566036/1‑151 (MQ=18)                                               ccaccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑C‑GAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCa                                                                                             <  1:557711/151‑1 (MQ=255)                                                 accaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑C‑GAACCCCCGATACGTTTCCGTATACACACTTTCCCGGCGGGCTCCTTCAGcc                                                                                              >  2:430830/1‑148 (MQ=255)                                                 accaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑C‑GAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACt                                                                                           <  2:535457/151‑1 (MQ=255)                                                  ccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑C‑GGACCCCCGCTACGTTGCCGTATACACACTTTTCCGGCGTGCGCCTTTAGCCAccc                                                                                          >  1:842254/1‑149 (MQ=18)                                                         cACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑C‑GCACCCCCCCTACGTTGCCGTATACACACCTTTCAGGCGGGCTCCCTCCGCCACTCGGACAcc                                                                                   >  1:320360/1‑151 (MQ=17)                                                            tGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑T‑GGAACCCCCAAACGTTTCCGTAAAAACACTTTTCAGGGGTGCTCCTTCCGCCCCCCGGACACCTc                                                                                 >  2:399004/1‑150 (MQ=17)                                                             gCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGTT‑C‑GAACCCCCCCTACGTTGCCGTATACACACCTTTCAGGGGGGCTCCCTCCGCCAccccggacccc                                                                                  >  1:571995/1‑138 (MQ=17)                                                              ccGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑C‑GGACCCCCCATACGTTGCCGTATACACACATTCCCGGCGGGCTCCTTCCGCCACCCCGACACCTCAcc                                                                              >  2:451990/1‑151 (MQ=17)                                                               cGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑T‑GAAACCCCGATAAGTTGCCGTATACACACCTTCCAGGCGTGCTCCTTCAGCCACCCGGACACCTCAcc                                                                              >  2:264773/1‑150 (MQ=17)                                                                  acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑T‑GAAACCCCCATAAGTTGCCGTATACACACTTTTCCGGCGTGCTCCTTCAGCCACCCGGACACCCCACCaaa                                                                           >  2:139809/1‑150 (MQ=14)                                                                   cagacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCGGGGGGGGGGTT‑T‑GAAACCCCGCTAAGCTGCCCGATACACACCTTTCCGGCGCCCTCCTCTCGCCCCCCCGGCACCTCAcc                                                                              >  2:361/1‑146 (MQ=9)                                                                     gacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑C‑GAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGt                                                                       <  2:113689/151‑1 (MQ=255)                                                                      acagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑T‑GAACCCCCCAAACGTTTCCGTATACACACTTTTCCCGCGGGGTCCTTCAGCCCCCCGGGCACCCCCCCAAAATGtt                                                                      >  2:679732/1‑151 (MQ=14)                                                                       cagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑T‑GAACCCCCCATACGTTTCCGTATAAACACTTTCCCGGGGGGGTCCTTCAGCCACCCCGACACCCCCCCAAAATtttt                                                                     >  1:206846/1‑151 (MQ=14)                                                                       cagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑T‑GAACCCCCCATACGTTGCCGTATACACACCTTCCAGGCGGGGTCCCTCCGCCCCCCCGGCACCCCCCCCAAATGttt                                                                     >  2:623608/1‑151 (MQ=14)                                                                       cagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑C‑GAAACCCCCCTACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCCGACACCCCCCCCAAATGttt                                                                     >  1:43493/1‑151 (MQ=14)                                                                            aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑T‑GAAACCCCCCTACGTTTCCGGATAAACACTTTCCCGGCGGGGCCCTTCAGCCACTCCGACAACCCACCAAAATGTTTTTCt                                                                 >  2:924742/1‑150 (MQ=14)                                                                               aTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑C‑CAACCCCCCCTACCTTGCCGCCTACACACCTTCCCCGCGCGCTCCCTCCCCCCCCCCGCCACCCCCCCCAAATGTTTTTCTccc                                                              >  1:516239/1‑150 (MQ=2)                                                                                   gAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑T‑GAAACCCCCATAAGTTGGCGTATACACAATTTTCAGGCGGGGTCCCTCAGCCACTCGGGCACCCCACCCAAATGTTTTGCTACCCAAcc                                                         >  1:342564/1‑151 (MQ=14)                                                                                      aaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑C‑GAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCa                                                      <  2:731467/151‑1 (MQ=255)                                                                                          aaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑C‑GAACCCCCCCTACGTTGCCGTATAAACACTTTCCAGGCGTGGTCCTTCAGCCACCCGGACAACCCACCAAAATGTTTTGCTACCCAACCCCATggg                                                  >  2:54250/1‑151 (MQ=14)                                                                                            gCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑C‑GAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTg                                                <  2:794098/151‑1 (MQ=255)                                                                                             cTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGGTT‑T‑CGACCCCCCATACGTTGCCGGATATACACCTTTCAGGCGTGCCCCCTCAGCCACTCGGACCCCCCCCCAAATTGTTTTTCTACCCCACCACCTgggggg                                               >  1:220383/1‑151 (MQ=2)                                                                                                 cACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑C‑GAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAAc                                           <  2:953990/151‑1 (MQ=255)                                                                                                  aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑C‑GAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACg                                          <  2:923253/151‑1 (MQ=255)                                                                                                      tCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑C‑GAACCCCCGATACGCTGCCGGATACACACCTTCCAGGCGGGCTCCTTCAGCCACTCCGACACCCCACCAAATTTTTTTTCTACCAAAACTCATGGGGGGGAACgggg                                       >  2:150518/1‑150 (MQ=14)                                                                                                       cGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑C‑GAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGgc                                      <  2:401248/150‑1 (MQ=255)                                                                                                        gTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑CGGAACCCCCG‑TACCTTGCCGTATACACACTTTCCAGGGGTGCTCCTTCAGCCACCCGGGCACCCCACCAAAATGTTTTTCTACCAAAACTCCTGGGGGGCAAAGGGgcg                                     >  1:683448/1‑150 (MQ=14)                                                                                                               ctctTCTTTAAATATGGCGGGGAGGGGGGGGTT‑C‑GAACCCCCCATACGTTGCCGTATACACACCTTTCAAGCGCGCTCCCTCCGCCCCTCCGGCACCCCACCCAATTTTTTTTTTTCCCAACCCCATGGGGGGGAAAGGGGGGccacta                               >  1:546341/1‑149 (MQ=2)                                                                                                                      ttAAATATGGCGGTGAGGGGGGGATT‑C‑GAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAg                        <  1:551658/149‑1 (MQ=255)                                                                                                                       tAAATAAGGCGGTGAGGGGGGGATT‑C‑GGACCCCCGATTCGTTGGCGTATACCCCCTTTCCAGGCGTTCCCCCTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGtt                      <  2:718536/150‑1 (MQ=21)                                                                                                                            atGGCGGTGAGGGGGGGATT‑C‑GAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCCAACCTCCTGGGTGGCAACGGGGCGCTACTATAGGGGGTTGGAg                  >  1:54867/1‑149 (MQ=255)                                                                                                                            atGGCGGTGAGGGGGGGATT‑C‑GAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTa                <  2:854477/151‑1 (MQ=255)                                                                                                                            atGGCGGTGAGGGGGGGATT‑C‑GAACCCCCGATACGGTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGGGTTGGAg                  >  1:1027272/1‑149 (MQ=255)                                                                                                                              ggCGGTGAGGGGGGGATT‑C‑GAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaa               <  1:11308/150‑1 (MQ=255)                                                                                                                               gCGGTGAGGGGGGGATT‑C‑GAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaaaa             <  2:393673/151‑1 (MQ=255)                                                                                                                                   tGAGGGGGGGATT‑C‑GAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGt         >  1:257167/1‑151 (MQ=255)                                                                                                                                    gAGGGGGGGATT‑C‑GAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGt         <  2:324673/150‑1 (MQ=255)                                                                                                                                      gggggggATT‑C‑GAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaa      <  1:590624/151‑1 (MQ=255)                                                                                                                                          gggATT‑C‑GAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaag  <  1:559228/151‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                      GCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑C‑GAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAG  >  NC_000913/925760‑926047

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