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

GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGA  >  NC_000913/925751‑926048                                                                                                                                                       |                                                                                                                                                         gAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGAt                                                                                                                                                         <  2:117622/151‑1 (MQ=255)  aGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGAtt                                                                                                                                                        <  1:1059465/151‑1 (MQ=255)      cACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGa                                                                                                                                                     <  2:259152/150‑1 (MQ=255)          gCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAAcc                                                                                                                                                  <  1:434545/149‑1 (MQ=255)               aCAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGAt                                                                                                                                            <  2:823394/150‑1 (MQ=255)                      gCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGcc                                                                                                                                    <  2:667426/151‑1 (MQ=255)                             aaTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑aca                                                                                                                            <  2:556426/151‑1 (MQ=255)                                 cacaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACAC‑ACt                                                                                                                       <  2:14056/151‑1 (MQ=255)                                   caGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAACCCCCCCTACGTTGCCGGAT‑ACAC‑ACctt                                                                                                                     >  1:387481/1‑151 (MQ=25)                                   caGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACAC‑ACttt                                                                                                                     <  1:390023/151‑1 (MQ=255)                                      gTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGTTTGAACCCCCCACACGTTGCCGTAT‑ACAC‑ACttt                                                                                                                     >  1:970739/1‑148 (MQ=255)                                          cATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACAC‑ACTTTCCAGGCg                                                                                                              <  1:981213/151‑1 (MQ=255)                                          cATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACAC‑ACTTTCCAGGCg                                                                                                              <  1:354463/151‑1 (MQ=255)                                             aCTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCGA‑GGGGGGGGTTTGGACCCCCCCAACGGTGCCGGAT‑AAAC‑ACCTTTCCGGCggg                                                                                                            >  2:789650/1‑148 (MQ=21)                                              cTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGTTTGAACCCCCCATACGTTTTCGCAT‑ACAA‑ACTTTCCAGGCGTg                                                                                                            >  1:503915/1‑149 (MQ=21)                                                 cACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACAC‑ACTTTCCAGGCGTGCTCCt                                                                                                       <  1:534378/151‑1 (MQ=255)                                                 cACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGGGTTGGAACCCCC‑ATACGTTTCCGTAT‑ACAC‑AACTTCCCGGCGGGCTcc                                                                                                        >  2:770816/1‑150 (MQ=21)                                                     cGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACAC‑ACTTTCCAGGCGTGCTCCTTCAg                                                                                                   <  1:536283/151‑1 (MQ=255)                                                     cGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACAC‑ACTTTCCAGGCGTGCTCCTTCAg                                                                                                   <  1:643237/151‑1 (MQ=255)                                                             aTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTGGAACCCCCGGTACGTTGCCGTAT‑ACAC‑ACTTTTCAGGCGGGGTCCCTCCGCCAccccg                                                                                           >  1:262035/1‑147 (MQ=18)                                                                 tCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCGG‑GGGGGGGATTCGAAACCCCCCTACGTTGCCCGAT‑ACAC‑CCTTTTCCCGCGTGCTCCTTCAGCCCCCCGGac                                                                                         >  1:529658/1‑149 (MQ=17)                                                                         gTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCCAACCCCCCGCACGGTGCCCTGT‑AAAC‑AATTTTCCGGCGTGGTCCTTCAGCCCCCCGGGAA‑CCACA‑cc                                                                               >  1:231215/1‑149 (MQ=17)                                                                             agacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGATTCGAACCCCCCATACGGTTCCGAAT‑ACAC‑ACTTTTCACGGGGGCTCCCTCCCCCACCCGGGCA‑CCACC‑CCaaa                                                                            >  2:269682/1‑148 (MQ=14)                                                                              gacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTg                                                                         <  2:927478/150‑1 (MQ=255)                                                                              gacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGt                                                                        <  1:375180/151‑1 (MQ=255)                                                                                      atatacGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGAGTTGGAACCCCGGTAAGTTGTGGGAA‑AATCAAATTTCCAGG‑GGGGTCCTTTAGCCCCCCGGACC‑CCCCACCCAAATtttttttc                                                                   >  2:339798/6‑147 (MQ=2)                                                                                        aTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGT‑GGGGGGGGGTCTAACCCCCCAAACGATGTCGGAT‑ACAC‑CATTTCTAGGCGCGCCCCTTCCACCACCCGGACA‑CCCCA‑CCCAAATGttttttttcc                                                               >  1:673753/1‑144 (MQ=9)                                                                                        aTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGGACCCCCCATACGTTGCCGTAT‑ACAC‑ACTTTTCCCGCGGGCTCCTTCCGCCACCCGGACA‑CCCCA‑CCCAAATGTTTTGTTAcc                                                               >  1:149974/1‑150 (MQ=14)                                                                                           cGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGGTTGAACCCCCGATACGTTGCCGTAT‑AAAC‑ACCTTTCCGGCGGGCTCCTTCCGCCACCCCGGCA‑CCCCA‑CCAAAATTTTTTGCTACCaaa                                                            >  2:539868/1‑150 (MQ=14)                                                                                                aaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAACCCCCCCTCACTTTTCGTAT‑AAAA‑ACTTTTCAAGGGGGCGCCTTTCGCCCCCCCGGCA‑ACCCA‑CCAAAATTTTTTTTTACCaaa                                                            >  2:265387/1‑145 (MQ=2)                                                                                                aaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGAACCCCCCCAACACTGCCGCAA‑AAAA‑CCATTCCAGGGGTGCCCCCTCCGCCCACCGGACT‑CCTCC‑CCAAATTGTTTTTTTTCCCAACCtta                                                       >  1:322820/1‑148 (MQ=2)                                                                                                     gCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGCAACCCCCCTACGTTTCCGCAT‑ACAA‑ACTTTTCCGGCGGGCTCCTTCAGCCCCCCCGGCACCCCCC‑CCAAATTTTTTTCTTCCCAAACC‑CATggggg                                                 >  1:730391/1‑149 (MQ=2)                                                                                                       tCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTTG‑GGGGGGGGGTTGTACCCCCCCTACGTTGTCGCGT‑ACAA‑ACCTTCCACGGGAGCGCCTTCACCCACCCCGGGA‑CCACA‑CCCAAATTTTTTTCTACCCAACCCCATgggggg                                                >  1:555816/1‑150 (MQ=9)                                                                                                       tCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTCCGAACCCCC‑CTACCTTTTCGTATAAAAA‑ACTTTCCAGGGGGGGCCCTCCGCCCCCCCGGA‑A‑CCCCC‑CCAAAATGTTTTTCTTCCAAAACTCATtgg                                                   >  2:198486/1‑147 (MQ=2)                                                                                                        cGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCCTACGTTGCCGGAT‑ACAC‑ACTTTTCACGGGTGCTCCCTCCGCCACTCCGGCA‑CCCCA‑CCCAAATTTTTTTTTACCAAACCCCCTgggggg                                                >  1:25/1‑149 (MQ=2)                                                                                                          cACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCGG‑GGGGGGGGGTTGAAACCCCGATACGTTTGCGTAT‑AAAA‑ACCTTTCCGGCGGGGTCCCTCACCCCCCCGGGCA‑CCCCC‑CCAAAATATTTTGCTTCCCAACCCCATGGGTggg                                               >  2:1020063/1‑147 (MQ=9)                                                                                                           aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAAACCCCCATACGGTGCCGGAT‑ACAC‑ACTTTTCCGGCGGGCGCCCTTAGCCCCTCCGACA‑CCCCA‑CCCAAATGTTTTTCTACCCAACCTCATGGGGGGGaa                                             >  1:424748/1‑149 (MQ=9)                                                                                                             tttCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCCTACGTTGGCGTAT‑ACAC‑ACTTTTCCGGGGTGCTCCCTCACCCACCCCGGCA‑CCCCA‑CCCAAATTTTTTTTTACCCAACCCCATTGGGGGGAACggg                                         >  1:1044759/1‑151 (MQ=2)                                                                                                             tttCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCATACGTTGCCGGAT‑AAAC‑ACTTTTCCGGCGGGCTCCCTCAACCCCCCCGACA‑CCCCC‑CCAAAATGTTTTTTTTCCAAACCCCATGGGTGGCAAAggg                                         >  1:741023/1‑151 (MQ=9)                                                                                                                    cGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGAACCCCCCCTACGTTTGCGTAT‑AAAC‑ACTTTTCCGGCGTGCTCCCTTCGCCCCCCGGGCA‑CCCCC‑CCCAAATGTTTTGTTACCAAACCTCCTCGGGGGGAACGGGgggcca                                   >  1:596774/1‑148 (MQ=9)                                                                                                                       gCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGGGTTGAAACCCCCCTACGTTGTCGTAT‑AAAA‑ACTTTCCAGGGGGGGGCCTTCAGCCACCCGGACA‑CCTCC‑CCAAAATTTTTTTTTTCCCAACCCCATGGGGGGGAAAGGGgggcca                                   >  2:637185/1‑145 (MQ=2)                                                                                                                          cttcttTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACtata                              <  2:780876/149‑1 (MQ=255)                                                                                                                            tcttTAAATATGGCGGTGA‑GGGGGGGGGTTGAACCCCCCATACGTTTTCGTAT‑ACAC‑ACTTTTCCGGGGGGCTCCTTCACCCCCCCGGGCA‑CCCCA‑CCCAAATGTTTTTTTACCCAAACCCATGGGGGGGAAAGGGGGGCGACTATAggg                           >  2:332939/1‑150 (MQ=9)                                                                                                                            tcttTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAggg                           <  1:68012/150‑1 (MQ=255)                                                                                                                               ttAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGa                          <  2:104397/148‑1 (MQ=255)                                                                                                                               ttAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGt                        <  2:1064980/150‑1 (MQ=255)                                                                                                                                tAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGtt                       <  2:212073/150‑1 (MQ=255)                                                                                                                                         cGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc             <  2:616037/151‑1 (MQ=255)                                                                                                                                           gTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACgg           <  2:419007/151‑1 (MQ=255)                                                                                                                                            tGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGt          >  2:710850/1‑151 (MQ=255)                                                                                                                                                   ggggATTCGAACCCCCGATACGTTGCCGTAT‑ACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaa    <  1:156782/151‑1 (MQ=255)                                                                                                                                                    gggATTCGAACCCCCGCTACGTTGCCGTAT‑ACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaag   <  2:606773/151‑1 (MQ=255)                                                                                                                                                      gATTCGAACCCCCGATACGTTGCCGTAT‑ACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAagaaga  >  2:1022049/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                         GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGA  >  NC_000913/925751‑926048

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