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

ACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925757‑926050                                                                                                                                                |                                                                                                                                                         aCTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAAc                                                                                                                                                   <  1:987632/151‑1 (MQ=255)           aGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCTGGGGGGGGGCTCGAACCCCCGCTa                                                                                                                                           >  2:280652/1‑149 (MQ=38)                   tAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCCAAACGTTGcc                                                                                                                                    >  1:975886/1‑148 (MQ=255)                              aGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTc                                                                                                                     <  2:169753/151‑1 (MQ=255)                                     aTACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGt                                                                                                              <  1:119352/151‑1 (MQ=255)                                      tACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTCGAAACCCCCATACGTTGTCGTATACACA‑CTTTCCAGGCg                                                                                                               >  2:295071/1‑149 (MQ=255)                                       aCTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTg                                                                                                             <  2:971383/150‑1 (MQ=255)                                          ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCt                                                                                                           <  2:151025/149‑1 (MQ=255)                                          ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGTTCGGACCCCCCATAAGTTGCCGTATATACA‑CTTTCCCGGCGTGCTcc                                                                                                         >  2:254586/1‑151 (MQ=21)                                           cACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGATCGAACCCCCCAAAAGCTGCCCGATACTCACCTTTCCCgggggctgc                                                                                                          >  2:389202/1‑143 (MQ=21)                                           cACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTCGAACCCCCCATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTcc                                                                                                         >  1:416142/1‑150 (MQ=255)                                              ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATCGAACCCCCCCTACGTTGCCGTATACACA‑CTTTTCACGCGTGCTCCtt                                                                                                       >  1:98430/1‑149 (MQ=18)                                                    accaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCAAAAGGTGGCGTATACAAA‑CTTTTCCGGCGTGGTCCTTCAGc                                                                                                   >  2:502261/1‑147 (MQ=18)                                                            cACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGacac                                                                                        <  1:638921/150‑1 (MQ=255)                                                                     acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGGTCGGAACCCCCACACGTTTCCGTATACAAA‑CTTTTCACGGGTGCTCCTTCCGCCCCCCGGGCACCCCACC‑caa                                                                              >  2:601769/1‑150 (MQ=14)                                                                       agacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTg                                                                           <  2:423474/151‑1 (MQ=255)                                                                         acagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGtt                                                                         <  1:373042/151‑1 (MQ=255)                                                                             aTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTTTTGc                                                                     <  1:503441/151‑1 (MQ=255)                                                                               aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAAACCCCCCTACGTTTCCCTATACACA‑CTTTCCCGGCGTGCCCCCTCAGCCACCCCGACACCCCACC‑CAAATGtttt                                                                       >  2:982553/1‑147 (MQ=14)                                                                               aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTTTTGCTa                                                                   <  2:162322/151‑1 (MQ=255)                                                                                  aTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCATACGTTGCCGGATAAACC‑CCTTTCCGGCGGGCTCCTTCCGCCACTCGGACACCCCCCC‑AAAATTTTTTTTTAcc                                                                 >  1:250913/1‑150 (MQ=11)                                                                                        aaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTTTTGCTACCAAACCTc                                                          <  1:861727/151‑1 (MQ=255)                                                                                        aaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGATAGGGGGGGGGTCCAAACCCCCCTACGTTGTCGTATACACA‑CTTTCCAGGGGTGGTCCCTCAGCCACACCGACACCCCACC‑AAATTTTTTTTCTACCAAAcc                                                            >  2:390734/1‑149 (MQ=14)                                                                                           aaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAAAGGGGGGGGGTCGAAACCCCCCTACCTTTCCGTTTTAACA‑CTTTTCAGGGGGGCTCCCTCCGCCCCTCCGGAACCCCACCAAAATTTTTTTTTTTCCAAAccc                                                           >  2:38214/1‑147 (MQ=2)                                                                                              aGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCGCTACGTTGCCGTATACACA‑CTTTCCAGGCGGGCTCCTTCAGCCCCTCGGGCACCCCCCC‑AAAATGTTTTTCTACCAAACCTCATggg                                                     >  1:988757/1‑150 (MQ=17)                                                                                                    cACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCCATACGGTGCCGGATACACA‑CCTTTCAGGCGGGGTCCTTCCCCCCCCCCGGCACCCCCCC‑AAAATTTTTTTCTACCAAACCTCCTGGGGGGcca                                               >  2:149370/1‑148 (MQ=9)                                                                                                    cACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCCATACGTTGCCGGATACACA‑CTTTTCAGGCGTGCTCCTTCAGCCACCCGGACACCCCACC‑AAAATGTTTTTCTTCCAAACCTccttggggggaa                                               >  1:530571/1‑139 (MQ=11)                                                                                                      cTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCGATACCTTGCCGTATACACC‑CTTTCCCGGCGGGCTCCTTCAGCCCCTCCGACACCTCCCC‑CAATTGTTTTGCTACCCAAACCCATggggggggacg                                             >  1:807886/1‑145 (MQ=14)                                                                                                        ttCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTTTTGCTACCAAACCTCATGGGTGGCAACggg                                           <  1:858889/150‑1 (MQ=255)                                                                                                                      tcttTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAg                               <  1:199682/148‑1 (MQ=255)                                                                                                                       cttTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCATACGGTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCCGACACCTCACC‑AAAATGTTTTTCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAggggg                           >  1:778267/1‑149 (MQ=18)                                                                                                                       cttTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCGCTACGTTTCCGTATATACA‑CTTTTCCGGGGGGCCCCTTTCTCCCCCCGGGCACCCCCCC‑AAAATTTTTTTTTTCCCAACCCCATGGGGGGGAACCGGGGGCTAATATAtgggg                           >  2:864161/1‑149 (MQ=2)                                                                                                                          tAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGtt                         <  2:230657/150‑1 (MQ=255)                                                                                                                          tAAATATGGCGGCGGGGGGGGGGTTTGAACCCCCCATACGTTTCCGGATAAACA‑CTTTTCAGGCGGGGTCCTTCAGCCCCCCCGACACCTCACC‑AAATTTTTTTGCTACCAAACCTCATGGGTGGGAAAGGGGCGGCACTATAagggggt                         >  2:941599/1‑146 (MQ=255)                                                                                                                           aaaTATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTgg                       <  2:833189/151‑1 (MQ=255)                                                                                                                           aaaTATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTgg                       <  1:68680/151‑1 (MQ=255)                                                                                                                            aaTATGGCGGTGAGGGGGGGGTTCGAACCCCCCATACGTTGCCGTATACACA‑CTTTCCAGGCGTGGTCCTTCAGCCACCCGGACACCCCCCC‑AAAATGTTTTTTTAACAAAACCCAAGGGGGGCAAAGGGGGGCGACCATAggggggtg                        >  2:929359/1‑144 (MQ=9)                                                                                                                            aaTATGGCGGTGAGGGGGGGATTCGAACCCCCCATACGTTGCCGTATACACA‑CTTTTCAGGGGGGCTCCCTCAGCCACCCGGACACCTCACC‑CAAATTTTTTGCTACCAAACCCCATGGGGGGGAAAGGGGGGCCACCATAGGGAggttg                       >  2:376720/1‑146 (MQ=17)                                                                                                                               atGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGt                    <  2:842053/150‑1 (MQ=255)                                                                                                                                   cGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc               >  2:44911/1‑151 (MQ=255)                                                                                                                                   cGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc               <  1:580684/151‑1 (MQ=255)                                                                                                                                      tGAGGGGGCGATTCGAACCCCCGATACGTTGCTGTATACACA‑CTTGCCAGGCGTGCTCCGTCAGCCACTCGGACACCTCACC‑AAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc               <  2:62279/148‑1 (MQ=255)                                                                                                                                       gAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc           >  2:649839/1‑151 (MQ=255)                                                                                                                                         gggggggATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa          >  2:753297/1‑150 (MQ=255)                                                                                                                                          ggggggATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaag        <  1:192928/151‑1 (MQ=255)                                                                                                                                               gATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACCCCTCACC‑AAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAgaagaa   <  2:655608/151‑1 (MQ=255)                                                                                                                                                aTTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  >  1:789848/1‑151 (MQ=255)                                                                                                                                                aTTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  <  2:924693/151‑1 (MQ=255)                                                                                                                                                aTTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  >  2:93051/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                                         ACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925757‑926050

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