Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,900	0	A	G	52.3%	39.4 / 63.1	44	noncoding (72/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (4/17);  new base G (23/0);  total (27/17)
Fisher's exact test for biased strand distribution p-value = 8.72e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-07

GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACAC‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925751‑926050                                                                                                                                                      |                                                                                                                                                         gAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGa                                                                                                                                                          <  2:669578/150‑1 (MQ=255)  aGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGtt                                                                                                                                                        >  2:21191/1‑151 (MQ=255)   gTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGTTc                                                                                                                                                       >  2:638821/1‑151 (MQ=255)   gTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGAtt                                                                                                                                                        <  2:98764/150‑1 (MQ=255)    ttCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGGtt                                                                                                                                                        >  2:707889/1‑149 (MQ=255)             gTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTTAGGGGGGGGGTCCCAccccc                                                                                                                                               >  1:258846/1‑149 (MQ=255)                   gCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAAGGGGGGGTTTGAAccccccc                                                                                                                                             >  1:528661/1‑143 (MQ=255)                         tAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGTGGGGGGGGTTCGAAACCCCGATACggt                                                                                                                                       >  1:452175/1‑143 (MQ=255)                         tAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCg                                                                                                                                   <  1:436363/149‑1 (MQ=255)                                  acaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCCATACGTTGTCGTATaaacac                                                                                                                          >  2:578937/1‑149 (MQ=255)                                    aGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGTTTGGAACCCCCATACGTTTCCCGTTACACAC‑att                                                                                                                      >  1:658503/1‑150 (MQ=25)                                          cATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCCCAAAGGTGCCGTATAAACAC‑TTTCCAGGCg                                                                                                               >  2:419644/1‑151 (MQ=21)                                                 cACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGGGTCGAAACCCCCACACGTTTCCGTATAAACAC‑TTTTCAGGCGTGGTcc                                                                                                         >  1:558321/1‑150 (MQ=21)                                                   cccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATCTGGCGGCCGGGGGGGGGGCAGCAACCCCGCCACGTTGTCGTATAAACAC‑CTTCCAGGCGGGccccct                                                                                                       >  1:564684/1‑148 (MQ=21)                                                   cccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACAC‑TTTCCAGGCGTGCTCCtt                                                                                                       <  1:120716/150‑1 (MQ=255)                                                    ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTTAACCCCCCCTACGTTTCCGTATACACAC‑TTTTCAGGCGTGCTcc                                                                                                         >  1:407611/1‑147 (MQ=21)                                                    ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAAACCCCGAAACGTTTCCGGATACAAAC‑TTTTCAGGCGGGGTCCtt                                                                                                       >  2:99403/1‑149 (MQ=18)                                                             aTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGATTCGAACCCCCCCGACATTTGCCTATACACAC‑ATTCCCGGCGGGCCCCCTCAcccccc                                                                                               >  1:84983/1‑145 (MQ=18)                                                                ttCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACAC‑TTTCCAGGCGTTCTCCTTCAGCCACTCGGac                                                                                          <  1:351050/150‑1 (MQ=255)                                                                 tCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACAC‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGaca                                                                                         <  2:133119/150‑1 (MQ=255)                                                                    cTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAAACCCCCATAAGGTGGCGGATAAAAAC‑TTTTCAGGCGGGGCCCTTTAGCCCCCCGGgcaac                                                                                       >  2:633771/1‑147 (MQ=17)                                                                          tACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCCATACGTTTCCGTTTACACAA‑CTTCCCGGCGTGCTCCCTCAGCCACCCCGACAccccccc                                                                                  >  2:496308/1‑143 (MQ=17)                                                                              gacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGGACCCCCCCTACGTTGCCGGCTAAACACTTTTCCAGG‑GGGCTCCCTCAGCCCCTCCGACACCC‑CACCaaa                                                                              >  2:386387/1‑147 (MQ=14)                                                                              gacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTTCAAACCCCGATACCTTTCCCTATACACAC‑CTTCCCGGCGGGCTCCCTCAGCCACCCCGACCCCTCCCCCAAAttt                                                                           >  1:27454/1‑150 (MQ=14)                                                                                cagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGCACCCCCCCTCCTTTTCCGTACAAACAA‑TTTTCCGGGGGGCTCCTTCACCCACCCGGACCCCT‑CACCAAAtttttt                                                                        >  1:208494/1‑151 (MQ=14)                                                                                 agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTGAACCACCCATACGTTTGCGTATACACAC‑ATTCCAGGGGTGCGCCCTCTGCCCCTCCGGCACCT‑ccacacaatttttt                                                                       >  1:432151/1‑138 (MQ=9)                                                                                  gaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACAC‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CACCAAATTGTTTTg                                                                      <  1:592908/151‑1 (MQ=255)                                                                                          gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACAC‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CACCAAATTGTTTTGCTACCaaa                                                              <  1:637249/151‑1 (MQ=255)                                                                                             aaaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCCCTCCGTTTGCGTTTAAACAC‑TTTTCCGGCGGGGTCCTTCCGCCACCCGGACCCCC‑CACCAAAATTTTTTTCTACCcaa                                                              >  1:574683/1‑148 (MQ=2)                                                                                                aaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACAC‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CACCAAATTGTTTTGCTACCAAACCTCAt                                                        <  1:412788/151‑1 (MQ=255)                                                                                                aaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGTTCGAAACCCCCATACGTTTCCGTATAAACAA‑CTTCCAGGCGTGGTCCTTCAGCCCCCCCGACACCC‑CCCCAAAATGTTTTTCTACCAAAcc                                                            >  1:368298/1‑147 (MQ=14)                                                                                                        cGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACAC‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCa                                                <  2:572397/151‑1 (MQ=255)                                                                                                             tttCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATCCCACCCCCGCGTCGATTCCGGATACACCC‑TTTTCAGGCGAGCGCCTTCACCCACCCGGACACCA‑CCCCAAAATTGTTTGCTAACAAACCCCCCCGGTGGCAACggg                                           >  2:641854/1‑151 (MQ=2)                                                                                                               tCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACAC‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACgggg                                          <  1:228246/150‑1 (MQ=255)                                                                                                                      aGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTCGAAACCCCGCTACTTTTCTGTATAAACAC‑ATTCCCGGCGGGCTCCTTTAACCACCCCGACACCT‑CCACCAACTATTTTGTTACCAAAACTCCTTGGGGGCAACgggggggcac                                    >  1:319580/1‑143 (MQ=9)                                                                                                                       gCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTCGAAACCCCCCTACGTTGCCGTATAAACAC‑CTTTCAGGGGGGGTCCTTCAGCCACTCGGGCACCC‑CCCCAAATTGTTTTTTTACCCAACCTCATGGGGGGCAACGGGGCGCCACtat                                 >  1:279835/1‑151 (MQ=17)                                                                                                                            tcttTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACAC‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAg                               <  2:284843/148‑1 (MQ=255)                                                                                                                               ttAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACAC‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGt                          <  1:340772/150‑1 (MQ=255)                                                                                                                                  aaTATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACAC‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGa                      <  2:461401/151‑1 (MQ=255)                                                                                                                                   atatGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACAC‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAg                     <  1:21190/151‑1 (MQ=255)                                                                                                                                      tGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACAC‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CACCAAATTGTTTTGCTACCAAACCTCCTGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTa                   >  2:218376/1‑150 (MQ=255)                                                                                                                                        gCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACAC‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaaaa                <  2:407843/151‑1 (MQ=255)                                                                                                                                           gTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACAC‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CACCAAATTGTTTTGCTACTAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACgg             <  1:643821/151‑1 (MQ=255)                                                                                                                                               gggggggATTCGAACCCCCGATACGTTGCCGTATACACAC‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa          >  2:535562/1‑150 (MQ=255)                                                                                                                                                 gggggATTCGAACCCCCGATACGTTGCCGTATACACAC‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaaga       <  2:311814/151‑1 (MQ=255)                                                                                                                                                      aTTCGAACCCCCGATACGTTGCCGTATACACAC‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  <  2:454608/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                         GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACAC‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925751‑926050

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