Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,900	0	A	G	43.1%	81.3 / 12.9	51	noncoding (72/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (5/24);  new base G (22/0);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 8.69e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGC‑AACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAA  >  NC_000913/925761‑926049                                                                                                                                               |                                                                                                                                                         ccGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGATTCGAAccccc                                                                                                                                               <  1:832613/151‑1 (MQ=255)  cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGG‑G‑AGGGGGGGGGTTCAAAccccc                                                                                                                                               >  2:301765/1‑151 (MQ=255)     aCAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGGTTCGAAccccc                                                                                                                                               >  2:26977/1‑147 (MQ=255)           aGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑A‑GGGGGGGATT‑GAAACCCCGCAACGtt                                                                                                                                       >  2:827655/1‑149 (MQ=37)               tAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGGTTCGAACCCCCGATACGTTGcc                                                                                                                                    >  1:719722/1‑148 (MQ=255)                      tCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGGTTTGAAACCCCGATACGTTGCCGtata                                                                                                                               >  1:135922/1‑146 (MQ=255)                          aGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGGTTCGAACCCCCCATACGTTGCCGTCTTAACACttt                                                                                                                       >  2:184032/1‑150 (MQ=255)                                  tACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGG‑G‑G‑CGGGGGGAATCCAACCCCCCATACGCTGTCGTATAAACACTTTCCaatggtg                                                                                                              >  1:219285/1‑145 (MQ=21)                                      ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTc                                                                                                           <  2:591719/150‑1 (MQ=255)                                            gCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGGTTTGAAACCCCGATACGTTGGCGTATAAACACTTTCCCGGCGTGCTCCCTCAGc                                                                                                    >  2:563430/1‑151 (MQ=18)                                            gCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGGTTCGAACCCCCGATACGTTGCCGTATAACCACTTTCCAGGCGTGCTCCTTCAGc                                                                                                    >  1:53978/1‑151 (MQ=255)                                              ccaccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCa                                                                                                  <  1:437903/151‑1 (MQ=255)                                                  caTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCg                                                                                              <  2:801297/151‑1 (MQ=255)                                                     gTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGa                                                                                            <  2:37221/150‑1 (MQ=255)                                                             ccGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGGTTCGAACCCCCCATACGTTTCCGTATACACACCTTCCCGGCGTGCTCCCTCCGCCACCCCGGCA‑ccccccc                                                                                  >  1:27240/1‑146 (MQ=17)                                                               gTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGGGTCGAACCCCCCATACGTTTCCGCATACACACTTTCCAGGCGGGCTCCTTTCGCCACCCGGGCA‑CCCCA‑cc                                                                                 >  2:325661/1‑149 (MQ=17)                                                                 acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAAt                                                                             <  1:607310/151‑1 (MQ=255)                                                                  cagacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAAtt                                                                            <  1:538705/151‑1 (MQ=255)                                                                  cagacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAAtt                                                                            <  2:626611/151‑1 (MQ=255)                                                                  cagacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAAtt                                                                            <  2:348340/151‑1 (MQ=255)                                                                       agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGtttt                                                                       <  2:262557/151‑1 (MQ=255)                                                                         aTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGATTCGAAACCCCGATACGTTTCCGGATACAAACTTTCCAGGCGGGCTCCCTCAGCCACTCGGGAA‑CCCCA‑CCAAAAtttttt                                                                       >  1:789575/1‑149 (MQ=14)                                                                          tAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑G‑GGGGGGGGTTCGAAACCCCGCAACGGTGCCGGATACAAACCTTCCAGGCGTGGTCCTTTCACCACTCGGAGA‑CCCCC‑CCAAAAttttttt                                                                      >  1:346616/1‑148 (MQ=14)                                                                            aaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGG‑G‑A‑GGGGGGGGTTTGAACCCCCCCCACGTTGGCGGATTCACACTTTTCAGGGGGGGCCCCTCACCCCCCCCGACACCCCCC‑CCAAATtttttttt                                                                     >  2:897022/1‑147 (MQ=2)                                                                            aaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑G‑A‑GGGGGGGGTTCGAACCCCCCCTACGTTGCCCTATACACACTTTTCAGGCGGGCTCCTTCTGCCCCTCCGGCA‑CCCCC‑CCAAAATTTTTTTCt                                                                    >  2:449709/1‑149 (MQ=14)                                                                                gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGGTTTGAACCCCCCAAACGGTTCCGTATAAACACCTTTCCGGCGTGCTCCTTCAGCCACCCCGACA‑CCCCACCCAAATttttttcttccca                                                                >  1:730351/1‑143 (MQ=11)                                                                                gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCaa                                                               <  1:834790/150‑1 (MQ=255)                                                                                  gAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAA‑GGGGGGGATTTTACCCCCCGTTCGGGTG‑CTTATATATACTTTCCAGCCGCGCGCCTTCTCGCCCCCGGAAA‑CCTCA‑CCACATTGttttttttcc                                                                 >  2:807066/1‑140 (MQ=2)                                                                                    aaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCt                                                           <  2:318200/150‑1 (MQ=255)                                                                                        aaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGGTTCGAACCCCCGATACGGTGTCCTATACACACTTTCCAGGCGTGCTCCTTCAGCCCCTCCGACA‑CCCCC‑CCAAAATGTTTTTCTATCAAAACCCAttg                                                      >  2:418126/1‑149 (MQ=14)                                                                                             tCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGG‑G‑G‑GGGGGGGGTTTGGACCCCCGAACCCTTTTCCTTTACAAACTTTCCAGGCGTGCTCCCTCACACCCTCTGAAA‑CCCCC‑CCCAATTATTTTTCTTTCCAACCTCATCGgggg                                                  >  2:743200/1‑150 (MQ=2)                                                                                              cGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGGTTCGAAACCCCCATACGTTTCCCTATACACACTTTCCCAGCGTGGTCCTTCAACCCCTCGGACA‑CCCCA‑CCAAAATTTTTTTCTACCCAACCCCCTGGGGGGc                                                 >  1:196298/1‑150 (MQ=11)                                                                                               gCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGGTTCGAACCCCCCATACGTTGCCGGATACACACTTTCCCGGCGGGCTCCTTTAGCCACCCGGACA‑CCCCA‑CCCAAATTTTTTTCTACCCAACCCCCTgggggg                                                  >  1:191166/1‑148 (MQ=11)                                                                                                 aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGGTTTGGAACCCCCCTACGTTGCCGTATAAACACTTTTCAGGCGGGCTCCTTCCGCCACCCCGACA‑CCCCA‑CCAAAATTTTTTTTTACCAAAACTCCCGGGGGGC‑AAc                                             >  2:141684/1‑150 (MQ=2)                                                                                                 aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGGTTGGCACCCCCCATACGTTTCCGTATACCCACCTTCCAGGGCTGCTCCTTCAACCACCCGGGCC‑CCCCA‑ACCAATTGTTTTTCTACCCAAACTCATTGgggggga                                               >  1:209567/1‑146 (MQ=2)                                                                                                    ttCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGC‑AACggg                                          <  1:468037/150‑1 (MQ=255)                                                                                                    ttCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGC‑AACgg                                           <  1:621725/149‑1 (MQ=255)                                                                                                      cGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑G‑GGGGGGGGTTTGGACCCCCCCTACGTTGGCGTATACACACTTTTCAGGCGGGGTCCTTCAGCCACCCGGACA‑CCCCC‑CCAAAATGTTTTGTTACCCAAACTCCTGGGGGGC‑AACgggg                                         >  2:746675/1‑149 (MQ=14)                                                                                                      cGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGC‑AACGGGgc                                        <  1:636428/150‑1 (MQ=255)                                                                                                         aCGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGGTTAGAACCCCCGACAACTTTTCGTATACACACTTTCCCGGCGGGCTCCTTCAGCCCCCCGGGGA‑CCCCC‑CCCAATTTTTTTTCTTCCAAACCCCATGGGGGGCAAACggggg                                        >  2:202143/1‑147 (MQ=9)                                                                                                           gAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGC‑AACGGGGCGctact                                  <  1:707188/151‑1 (MQ=255)                                                                                                                     ttAAATATGGCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGC‑AACGGGGCGCTACTATAGGGAGt                         <  2:415818/150‑1 (MQ=255)                                                                                                                      tAAATATGGCGGT‑G‑A‑GGGGGGGGTTCGAACCCCCCATACGTTGCCGTATACACACTTTTCCGGCGTGCTCCTTCAGCCCCTCGGACA‑CCTCC‑CCAAAATGTTTTTCTACCAAACCTCATGGGTGGC‑AACGGGGCGCTACTATAAGGAGtt                        >  1:221136/1‑150 (MQ=18)                                                                                                                      tAAATATGGCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGC‑AACGGGGCGCTACTATAGGGAGtt                        <  1:444276/150‑1 (MQ=255)                                                                                                                        aaTATGGCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCG‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGC‑AACGGGGCGCTACTATAGGGAGTTGGa                     <  2:764688/151‑1 (MQ=255)                                                                                                                         atatGGCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGC‑AACGGGGCGCTACTATAGGGAGTTgg                      >  1:298443/1‑149 (MQ=255)                                                                                                                             ggCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGC‑AACGGGGCGCTACTATAGGGAGTTGGAGTaa                 <  2:740563/150‑1 (MQ=255)                                                                                                                             ggCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGC‑AACGGGGCGCTACTATAGGGAGTTGGAGTaa                 <  1:31972/150‑1 (MQ=255)                                                                                                                               cGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGC‑AACGGGGCGCTACTATAGGGAGTTGGAGTaaaa               <  2:72197/150‑1 (MQ=255)                                                                                                                               cGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGC‑AACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc              <  1:564974/151‑1 (MQ=255)                                                                                                                                              gATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGC‑AACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAgaagaa  <  2:100289/151‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                         CCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑G‑A‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGC‑AACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAA  >  NC_000913/925761‑926049

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