Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,900	0	A	G	55.6%	31.0 / 46.0	36	noncoding (72/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (4/12);  new base G (20/0);  total (24/12)
Fisher's exact test for biased strand distribution p-value = 1.45e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-06

GCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGA  >  NC_000913/925760‑926045                                                                                                                                              |                                                                                                                                                  gCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGATTCGAAcc                                                                                                                                           <  2:313848/149‑1 (MQ=255)       cAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGGTTCGAAccccc                                                                                                                                        >  2:169283/1‑146 (MQ=255)           cAGCGTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTTAGGTGGGGGATGCAAAcccccacacctt                                                                                                                                 >  1:228655/1‑143 (MQ=37)                  gAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGGGTCGAAACCCCCCTACGTTGCCGta                                                                                                                          >  1:397051/1‑149 (MQ=255)                  gAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGATTCGAACCCCCGATACGTTGCCGtata                                                                                                                        <  2:228659/151‑1 (MQ=255)                    aaTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGATTCGAACCCCCGATACGTTGCCGTATaca                                                                                                                      <  1:369203/151‑1 (MQ=255)                        cacaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAATCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGATTCGAACCCCCGATACGTTGCCGTATacacac                                                                                                                   <  2:381968/150‑1 (MQ=255)                                  aTACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGt                                                                                                        <  1:446125/151‑1 (MQ=255)                                      tCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGGTTCGAAACCCCCCCACGTTTCCGTATACACACCTTCCAGGCGTg                                                                                                       >  2:322333/1‑148 (MQ=21)                                         aCCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGGTTCGAACCCCCGATACGTTGCCGGATAAACACTTTCCCGGCGTGCTCCtt                                                                                                 >  2:449684/1‑151 (MQ=255)                                         aCCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGG‑GGGGGGTTCGGACACCCCACACGGTTCCGTATAAACACCTTTCCGGGGTGCTCCtt                                                                                                 >  1:292493/1‑151 (MQ=21)                                           ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGGTTTGAACCCCCCCCACCTTGCCGTATAAACCCTTTCTCGGCGGGGTcc                                                                                                   >  1:499959/1‑147 (MQ=21)                                            cGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGGTTTGAACCCCCGATACGGTGCCGTATACACACTTTCCAGGCGGGCTccct                                                                                                 >  2:165009/1‑146 (MQ=255)                                             gCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGc                                                                                             <  2:145582/151‑1 (MQ=255)                                              cccACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGGTTCGAACCCCCCATACGTTGCCGTATACACACCTTCCACGCGGGCTCCTTCAGcc                                                                                            >  2:692459/1‑151 (MQ=255)                                                caccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGG‑GGGGGGGTTGAACACCCCCAACCCTGCCCTATACACACTTTTCCCGCGCGCTCCCTTAGCCAc                                                                                          >  1:486377/1‑151 (MQ=18)                                                             gCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGGTTCGAACCCCCCCCACGTTGCCGGATACAAACTTTTCCGGGGGGCTCCTTCTGCCACTCGGGCAcccccc                                                                             >  2:429241/1‑147 (MQ=17)                                                             gCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGG‑GGGGGGTTCGAAACCCCCCCAAGTTGCCGCATAACCACCTTCCACGCGGTCTCCCTCAACCACCCCGacac                                                                                  >  1:489519/1‑146 (MQ=17)                                                                  acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTGAAATATGGCGGTGAGG‑GGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCaaa                                                                         <  1:498588/150‑1 (MQ=255)                                                                    agacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGGTTCGAAACCCCCATACGTTTGCGGATAAACACTTTTCAGGCGTGCTCCTTCAGCCCCTCGGACAccccccccaaat                                                                       >  2:181803/1‑140 (MQ=14)                                                                            aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAAG‑GGGGGGTTCGAAACCCCGAGACGTTGTCGTATAAACACTTTCCAGGCGGGCCCCTTCAGCCACCCGGACAACCCCCCCAATTGtttt                                                                  >  1:90119/1‑147 (MQ=14)                                                                                 gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGATTCGAAACCCCCATACGTTGCCGGATACACCCTTTCCCGGCGTGGTCCTTCCGCCACTCGGGCACCTCCCCAAACTGTTTTTCTAccca                                                          >  1:448728/1‑148 (MQ=14)                                                                                          aaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGGTTTGAACCCCCCATACCTTGCCGGATACACACCTTTCAAGCGGGCCCCCTCCGCCACCCCGGCACCCCCCCCAATTTTTTTTTTACCAAAACCCATtgg                                                >  1:36469/1‑151 (MQ=2)                                                                                                  aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGGTTTGAAACCCCCCTACCTTGTCGGATAAACACTTTTCCGGGGGGGTGCTTCTCCCCCCCCGACACCCCCCCAAAATGTTTTTCTTCCAAAACCCCTgggggg                                             >  1:626670/1‑146 (MQ=2)                                                                                                    tttCGTACGAGCTCTTCTTTAAATATGGCGGAGAGG‑GGGGGGTTTGGAACCCCCCTACGTTTGCGTATAAACACATTTCCGGGGGGGTCCCTTAGCCCCCCGGACACCTCACCCAAATGTTTTTCTTCCCAAACCCCTGAGTgg                                             >  1:449823/1‑144 (MQ=2)                                                                                                     ttCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGGTTCGAACCCCCCCTACGTTGCCGTATACACACTTTCCAGGCGGGCTCCTTCAGCCACCCGGACACCTCACCCAAATGTTTTTCTACCAAACCCCCTTGGGGGCAACgggg                                     >  2:384921/1‑151 (MQ=14)                                                                                                     ttCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGGGTCGGACCCCCCCTACGTTGCCGTATACACACTTTCCAGGCGGGGTCCTTCCGCCACTCCGACACCTCACCAAAATGTTTTTCTACCAAAACTCATGGGGGGGAACgggg                                     >  2:103534/1‑151 (MQ=14)                                                                                                       cGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGgc                                    <  2:315165/150‑1 (MQ=255)                                                                                                           cGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGGTTCGGAACCCCCATACGGTGCCGGATAAACACCTTTCAGGCGGGCTCCCTCAGCCACTCGGACACCCCACCAAAATTTTTTTCTACCAAACCCCCTGGGGGGGAAAgggggg                                   >  2:555101/1‑145 (MQ=9)                                                                                                                   tcttTAAATATGGCGGTGAGG‑GGGGGGTTCGAACCCCCGCTACGTTGCCGTATACACACCTTCCCGGCGTGCTCCTTCAGCCACCCGGACACCCCACCAAATTGTTTTTCTACCAAACCTCATGGGGGGCAAAGGGGCGCTACTATAggg                        >  2:155965/1‑150 (MQ=18)                                                                                                                   tcttTAAATATGGCGGTGAGG‑GGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAggg                        <  2:462730/150‑1 (MQ=255)                                                                                                                               gCGGTGAGG‑GGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaaaa           <  1:246667/151‑1 (MQ=255)                                                                                                                                 ggTGAGG‑GGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACg         <  2:36465/151‑1 (MQ=255)                                                                                                                                   tGAGG‑GGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGt       >  1:141692/1‑151 (MQ=255)                                                                                                                                       g‑gggggATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaag   <  1:169285/151‑1 (MQ=255)                                                                                                                                           gggATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaaga  >  2:667062/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                  GCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGG‑GGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGA  >  NC_000913/925760‑926045

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