Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,912	0	G	C	47.1%	40.5 / 18.5	34	noncoding (60/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base G (7/11);  new base C (16/0);  total (23/11)
Fisher's exact test for biased strand distribution p-value = 1.27e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTAAT  >  NC_000913/925777‑926056                                                                                                                                        |                                                                                                                                                 gAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTCGAAACCCCCCTACGTTGCCGta                                                                                                                                     >  1:397051/1‑149 (MQ=255) gAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGtata                                                                                                                                   <  2:228659/151‑1 (MQ=255)   aaTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATaca                                                                                                                                 <  1:369203/151‑1 (MQ=255)       cacaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAATCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATacacac                                                                                                                              <  2:381968/150‑1 (MQ=255)                 aTACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGt                                                                                                                   <  1:446125/151‑1 (MQ=255)                     tCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAAACCCCCCCACGTTTCCGTATACACACCTTCCAGGCGTg                                                                                                                  >  2:322333/1‑148 (MQ=21)                        aCCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCGATACGTTGCCGGATAAACACTTTCCCGGCGTGCTCCtt                                                                                                            >  2:449684/1‑151 (MQ=255)                        aCCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGTTCGGACACCCCACACGGTTCCGTATAAACACCTTTCCGGGGTGCTCCtt                                                                                                            >  1:292493/1‑151 (MQ=21)                          ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCCCACCTTGCCGTATAAACCCTTTCTCGGCGGGGTcc                                                                                                              >  1:499959/1‑147 (MQ=21)                           cGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCGATACGGTGCCGTATACACACTTTCCAGGCGGGCTccct                                                                                                            >  2:165009/1‑146 (MQ=255)                            gCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGc                                                                                                        <  2:145582/151‑1 (MQ=255)                             cccACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCCATACGTTGCCGTATACACACCTTCCACGCGGGCTCCTTCAGcc                                                                                                       >  2:692459/1‑151 (MQ=255)                               caccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGGTTGAACACCCCCAACCCTGCCCTATACACACTTTTCCCGCGCGCTCCCTTAGCCAc                                                                                                     >  1:486377/1‑151 (MQ=18)                                            gCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCCCCACGTTGCCGGATACAAACTTTTCCGGGGGGCTCCTTCTGCCACTCGGGCAcccccc                                                                                        >  2:429241/1‑147 (MQ=17)                                            gCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGGTTCGAAACCCCCCCAAGTTGCCGCATAACCACCTTCCACGCGGTCTCCCTCAACCACCCCGacac                                                                                             >  1:489519/1‑146 (MQ=17)                                                 acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTGAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCaaa                                                                                    <  1:498588/150‑1 (MQ=255)                                                   agacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAAACCCCCATACGTTTGCGGATAAACACTTTTCAGGCGTGCTCCTTCAGCCCCTCGGACAccccccccaaat                                                                                  >  2:181803/1‑140 (MQ=14)                                                           aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAAGGGGGGGTTCGAAACCCCGAGACGTTGTCGTATAAACACTTTCCAGGCGGGCCCCTTCAGCCACCCGGACAACCCCCCCAATTGtttt                                                                             >  1:90119/1‑147 (MQ=14)                                                                gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAAACCCCCATACGTTGCCGGATACACCCTTTCCCGGCGTGGTCCTTCCGCCACTCGGGCACCTCCCCAAACTGTTTTTCTAccca                                                                     >  1:448728/1‑148 (MQ=14)                                                                         aaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCATACCTTGCCGGATACACACCTTTCAAGCGGGCCCCCTCCGCCACCCCGGCACCCCCCCCAATTTTTTTTTTACCAAAACCCATtgg                                                           >  1:36469/1‑151 (MQ=2)                                                                                 aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCCCTACCTTGTCGGATAAACACTTTTCCGGGGGGGTGCTTCTCCCCCCCCGACACCCCCCCAAAATGTTTTTCTTCCAAAACCCCTgggggg                                                        >  1:626670/1‑146 (MQ=2)                                                                                   tttCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGTTTGGAACCCCCCTACGTTTGCGTATAAACACATTTCCGGGGGGGTCCCTTAGCCCCCCGGACACCTCACCCAAATGTTTTTCTTCCCAAACCCCTGAGTgg                                                        >  1:449823/1‑144 (MQ=2)                                                                                    ttCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCCCTACGTTGCCGTATACACACTTTCCAGGCGGGCTCCTTCAGCCACCCGGACACCTCACCCAAATGTTTTTCTACCAAACCCCCTTGGGGGCAACgggg                                                >  2:384921/1‑151 (MQ=14)                                                                                    ttCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTCGGACCCCCCCTACGTTGCCGTATACACACTTTCCAGGCGGGGTCCTTCCGCCACTCCGACACCTCACCAAAATGTTTTTCTACCAAAACTCATGGGGGGGAACgggg                                                >  2:103534/1‑151 (MQ=14)                                                                                      cGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGgc                                               <  2:315165/150‑1 (MQ=255)                                                                                          cGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGGAACCCCCATACGGTGCCGGATAAACACCTTTCAGGCGGGCTCCCTCAGCCACTCGGACACCCCACCAAAATTTTTTTCTACCAAACCCCCTGGGGGGGAAAgggggg                                              >  2:555101/1‑145 (MQ=9)                                                                                                  tcttTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCGCTACGTTGCCGTATACACACCTTCCCGGCGTGCTCCTTCAGCCACCCGGACACCCCACCAAATTGTTTTTCTACCAAACCTCATGGGGGGCAAAGGGGCGCTACTATAggg                                   >  2:155965/1‑150 (MQ=18)                                                                                                  tcttTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAggg                                   <  2:462730/150‑1 (MQ=255)                                                                                                              gCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaaaa                      <  1:246667/151‑1 (MQ=255)                                                                                                                ggTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACg                    <  2:36465/151‑1 (MQ=255)                                                                                                                  tGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGt                  >  1:141692/1‑151 (MQ=255)                                                                                                                      ggggggATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaag              <  1:169285/151‑1 (MQ=255)                                                                                                                         gggATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaaga             >  2:667062/1‑149 (MQ=255)                                                                                                                                  aCCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTAAt  >  2:463841/1‑151 (MQ=255)                                                                                                                                        |                                                                                                                                                 GAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTAAT  >  NC_000913/925777‑926056

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