Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,097,581	0	A	G	46.9%	69.4 / 19.9	50	noncoding (72/88 nt)	serX	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (1/25);  new base G (23/0);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 4.11e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGAAAA  >  NC_000913/1097436‑1097730                                                                                                                                                  |                                                                                                                                                       cACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑Cg                                                                                                                                                   <  2:472162/150‑1 (MQ=31)  aCTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑CGa                                                                                                                                                  <  2:533160/150‑1 (MQ=31)     gCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑CGAAc                                                                                                                                                <  2:223660/149‑1 (MQ=33)      cGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGTT‑TGAGcccc                                                                                                                                             >  1:77168/1‑151 (MQ=255)              gTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑CGAACCCCCGATACGt                                                                                                                                     <  2:356797/151‑1 (MQ=40)                   tCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑TGCACCCCCCATACggttcc                                                                                                                                 >  1:140768/1‑145 (MQ=11)                    ccAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑TGATcccctcccccgcctcgc                                                                                                                                >  1:642612/1‑137 (MQ=11)                                  ccAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTcc                                                                                                                  <  1:175686/150‑1 (MQ=255)                                   cAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAg                                                                                                                <  2:666679/151‑1 (MQ=255)                                     gTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGT‑CGGACCCCCCATACGCTTCCGtataaccacctt                                                                                                                    >  1:678205/1‑137 (MQ=11)                                       tGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGGTT‑CGAACCCCCCCCACAGTTTCGTATTCACAATTTCCACGc                                                                                                              >  1:329658/1‑149 (MQ=11)                                          cGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGc                                                                                                          <  2:118121/150‑1 (MQ=255)                                             aCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGT‑TGAACCCCCCCTACCGTGCCGTATACACACCTTTCAGGCGGGCTcc                                                                                                       >  1:355767/1‑150 (MQ=14)                                                     cGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCa                                                                                               <  1:461146/150‑1 (MQ=255)                                                      ggggCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACt                                                                                             <  2:77188/151‑1 (MQ=255)                                                            tctcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGa                                                                                         <  2:305910/149‑1 (MQ=255)                                                               cATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑TGAACCCCCCATAAATTTCCGTTTACACAATTTTCAGGGGGGGTCCTTCAGCCACTCCGGCAcc                                                                                     >  1:283674/1‑150 (MQ=17)                                                                      ccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑TGAACCCCCCGTACGTTGCCGTATACACACTTTTCCGGCGTGCTCCCTCCGCCCCCCCGGCACCCCAcc                                                                                >  2:69527/1‑148 (MQ=17)                                                                        ggTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑TGAACCCCCCATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACCCGGACACCTCACCAAAtt                                                                           >  2:219557/1‑151 (MQ=255)                                                                        ggTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGT‑TGTAACCCCGATACGTTGCCGTATAAACACCTTCCAGGCGGGCTTCTTCAGCCACCCGGACACCTcccc                                                                                >  2:279669/1‑146 (MQ=14)                                                                        ggTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAAtt                                                                           <  1:121238/151‑1 (MQ=255)                                                                         gtgtgtGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑CGAACCCCCCATACGTTTACCTATACACACTTTTCAGGGGGGCTCCTTCCGCCACCCGGACACCCCCCCaaact                                                                           >  2:301201/1‑148 (MQ=14)                                                                             gtgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCg                                                                       <  1:491241/150‑1 (MQ=255)                                                                              tgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGGT‑TGAACCCCCCGTACGTTTCCCTATACAAACTTTCCAGGCGGGCTCCTTCAGCCACTCCGACACCCCACCaaaatattgtt                                                                     >  1:121503/1‑143 (MQ=14)                                                                               gCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑TGAAACCCCCATACGTTTCCGTATACACACTTTTCAGGCGTGCTCCTTTCGCCCCTCGGACACCCCCCCAAAtttttggt                                                                     >  2:254440/1‑144 (MQ=14)                                                                                       cGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTg                                                            <  1:578754/151‑1 (MQ=255)                                                                                             aaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑CGAACCCCCCATACGTTGCCGTATACACACCTTTCCGGCGGGCTCCTTCAGCCACTCGGGCACCCCACCCAATTTTCGTTCCTGTCtt                                                             >  2:150809/1‑144 (MQ=14)                                                                                               aaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑TGAAACCCCCGTACGTTTTCGTATAAACACTTTTCCGGCGGGGTCCCTCCGCCACCCCGGCACCCCCCCCAAATTTTGTTTCTGTCTTTCTgg                                                        >  1:37433/1‑147 (MQ=2)                                                                                                  cccGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑CGAAACCCCCCTACGTTGCCGTTTACACACTTTTCAGGCGTGCTCCCTCAGCCACCCGGACACCTCACCCAATTTTTGTTCCTTTCTTGCTGGGACggg                                                  >  2:356075/1‑150 (MQ=14)                                                                                                  cccGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGCGGGGGGGGTT‑TGGAACCCCCCTACTCTGCCGTATAAAAAATTTTCCGGCGGGGGCCCTCTCCCACTCGGGGAACCCACCCAAATATCGTTTCTGTTCTGCTGGaa                                                      >  1:349942/1‑146 (MQ=2)                                                                                                       aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGTT‑TTGACCCCCCCTACCTTTCCGTTTTAACACTTTCCCGGCGTGCTCCTTTCCCCCCCCCGACACCCCACCAAAATGTCGTTTCTGTCTTTCTGGGACCGGCGCTa                                             >  1:537558/1‑150 (MQ=9)                                                                                                       aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGG‑TT‑TGAACCCCCCACACGTTGCCCGATAAACACTTTCCAGGGGTGCTCCTTCAGCCACCCGGACACCCCACCAAAATGTCGTTTCTGTTTTTCTTGGACGGGGGCCaa                                            >  2:254105/1‑150 (MQ=11)                                                                                                         tttCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTT‑CGCACACCCGCTACGTTGCCGTATACACCATTTTCCGGGGCGGTCCTTTCACCACCCGGACACCCCACCAAAATTTTGTTCCTGTCTTTCTgggaac                                                    >  2:96156/1‑137 (MQ=2)                                                                                                          ttCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAAtt                                          <  1:143612/150‑1 (MQ=255)                                                                                                           tCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAAttt                                         <  1:665211/150‑1 (MQ=255)                                                                                                            cGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTa                                        <  2:152602/150‑1 (MQ=255)                                                                                                             gTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGTT‑TGAAACCCCCATACCTTGTCGTATACACACTTTCCCGGGGGGCTCCTTCCGCCCCCCGGGCCCCCCCCCAAAATTTCTTTTCTGTCTTTCTTGGACGGGCGCCAAttt                                         >  1:660130/1‑148 (MQ=9)                                                                                                              tACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGT‑CGAAACCCCCATACGTTTCCGTATAAACACTTTTCAGGGGTGCTCCTTCACCCACTCCGGCACCCCACCAAAATTTTGTTCCTGTCTTTCTGGGACGGGGGCCAATTTTggg                                     >  1:601016/1‑151 (MQ=9)                                                                                                               aCGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGTTCCGAACCCCCGA‑ACGTTGCCCTATACACACTTTTCAGGGGTGCTCCTTCAGGCACCCCGGCACCTCACCCAAATGTCGTTCCTGTCTTTTTGGAACAGGCGCCAATTTAggg                                     >  2:214206/1‑150 (MQ=17)                                                                                                               aCGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAggg                                     <  2:338712/150‑1 (MQ=255)                                                                                                                   gcgcttcttTAAATATGGCGGTGAGGGGGGGATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGaaa                                  <  2:424487/146‑1 (MQ=255)                                                                                                                      cttcttTAAATATGGCGGTGAGGGGGGGATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCa                               <  1:531576/149‑1 (MQ=255)                                                                                                                        tcttTAAATATGGCGGTGAGGGGGGGATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCAt                              <  1:175680/148‑1 (MQ=255)                                                                                                                          tttAAATATGGCGGTGAGGGGGGGGGT‑TGTACCCCCCATACATTTTCGTATTCACAATTTTCAGGGGGGGTCCTTCCGCCACCCCGACACCCCCCCAAAATTTTGTTCCCGTCTTGCTGGAACGGGCGCCAATTTATGGaaaaaattaac                          >  1:166592/1‑142 (MQ=12)                                                                                                                            tAAATATGGCGGTGAGGGGGGGGTT‑TGAACCCCCCCTACGTTTCCGTGTACACACTTTCCCGGCGTGCTCCTTCAGCCCCTCCGACACCCCACCAAATTTTCGTTCCTTTCTTTCTTGAACGGGCGCTCATTTAGGGAAATCATGAcc                          >  2:124301/1‑148 (MQ=18)                                                                                                                             aaaTATGGCGGTGAGGGGGGGATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGa                       <  2:69943/150‑1 (MQ=255)                                                                                                                                 atGGCGGTGAGGGGGGGATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTca                  <  2:224341/151‑1 (MQ=255)                                                                                                                                     cGGTGAGGGGGGGATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTcaacaa              <  2:601057/151‑1 (MQ=255)                                                                                                                                             gggggATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTg      <  1:530295/151‑1 (MQ=255)                                                                                                                                                 gATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGaaaa  <  1:275094/151‑1 (MQ=255)                                                                                                                                                  aTT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGaaaa  >  2:116277/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                       CACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATT‑CGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGAAAA  >  NC_000913/1097436‑1097730

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