Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,231,094	0	T	C	48.1%	57.1 / 17.3	52	F736S (TTC→TCC)	ygjK	alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base T (18/9);  new base C (0/25);  total (18/34)
Fisher's exact test for biased strand distribution p-value = 1.21e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGG‑TTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAA  >  NC_000913/3230948‑3231236                                                                                                                                                    |                                                                                                                                                 gggAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGG‑TTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCg                                                                                                                                               <  2:1006749/150‑1 (MQ=255)     aCGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGG‑TTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACg                                                                                                                                          <  2:191781/151‑1 (MQ=255)      cggcggCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGG‑TTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGc                                                                                                                                         <  1:504780/151‑1 (MQ=255)        gcggcgTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGG‑TTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCa                                                                                                                                       <  2:645738/151‑1 (MQ=255)                   cAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGG‑TTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGG‑TT‑AAcc                                                                                                                          <  2:205252/151‑1 (MQ=255)                        cGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGG‑TTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGG‑TT‑AAccgccg                                                                                                                      >  2:876579/1‑150 (MQ=255)                                     tGATTTCTACTGGGGGGGGCGGTTTGGGGGGGATCTTTTCTGTTTTGGTCGAAAAGGGAGGGTGCTTTCCGT‑TTATCGTGATGACGCCCTGAGGTGGGGGGATATTTTCTCCCGACACGCAAAGGGG‑TT‑AACCGCCGATGGACCAATTc                                                                                                          <  1:347912/149‑1 (MQ=255)                                         atCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGG‑TTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGaa                                                                                                     >  2:132325/1‑150 (MQ=255)                                          tCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGG‑TTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGaaa                                                                                                    >  2:654986/1‑150 (MQ=255)                                                ggggggggggtgtgtGGGTGGTTCATTTTTGGTTGGTTGGGAAGGGGATGGAGTTTTGCGT‑TTATGGTGAGGATGCCCTGAAGGTGGCGGTTACTTTCTCCCGCCACGCAAAGGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTa                                                                                                <  1:613519/134‑1 (MQ=255)                                                ggggggggCGCGTAGGGGTGGATCTGTTCTGTTTTGGTCTGAAAGGGATGGAGCGTTACGG‑TTATCGTGATGAGCCCCTGAAGGGGGCGGATACTTTCTTCGGGCCCGAAAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACaa                                                                                             <  1:525244/146‑1 (MQ=255)                                                    tGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGG‑TTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCg                                                                                         >  2:342171/1‑151 (MQ=255)                                                      ggCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGG‑TTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCt                                                                                       >  1:262730/1‑151 (MQ=255)                                                        cgcgTATGGGTGGATCAATTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGG‑TTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTg                                                                                      >  1:94279/1‑150 (MQ=255)                                                         gcgTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGG‑TTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGAc                                                                                    >  2:998326/1‑151 (MQ=255)                                                               ggggggATCAGTTCTGTTTTGGTCGGAAGGGGAGGGGGCTTTGCGG‑TTATCGTGATGATGCCCTGAAGGGGGCGGATACGTTCTCCCGGCACGCAAAGGGG‑TT‑AGCCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAgg                                                                                 <  1:329096/148‑1 (MQ=255)                                                                  tGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGG‑TTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGcac                                                                            >  1:728357/1‑150 (MQ=255)                                                                   gattcgttttgggTTTGTTGAGAAAGGGAGGGGGCTTTGTGT‑TTTTCGTAATGACCCCCGGAAGCGGGCGTAT‑CCTTCCTCCCGCCCCCCCAAGGG‑TTAAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGcac                                                                            <  1:361995/138‑1 (MQ=255)                                                                        ggtgtgggTGTGGTCGAAAAGGGAGGGGGGGTTCCGT‑TTATCGTTTTGAGGCCCGGATGCGGGCGAATACTTCCTCCCGCCCCGCCAAGGGT‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGgcg                                                                               <  1:425195/135‑1 (MQ=255)                                                                         gTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGG‑TTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGgc                                                                    >  1:250718/1‑151 (MQ=255)                                                                          ttttGGTTTGGTCGGAAGGGAATGGCG‑TTTATGGTTTATCGTGATGATGCCTTAAAGCTGGGGAATAGTTTCTTCCGCCACCCAAAGGGG‑TA‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAgg                                                                     <  2:24957/146‑1 (MQ=255)                                                                           tCTGGTTGGATCAGAAAGGGATGGTGTTTTCGGT‑TTATGGTGAGGACGCCGAGAAGCGGGAGAATCTTTTCTCCCGGCCCACCAAAGGT‑TA‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGgcg                                                                   <  1:648564/150‑1 (MQ=255)                                                                                tttGGTCTGAAGGGGAGGGAGCGTTCCGG‑TTATCGTGATGATGCCCTGAAGCGGGCGGATCCGTTCTCCCGGCCCGCCAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCa                                                              <  2:121426/150‑1 (MQ=255)                                                                                 ttgatataaaaGGGATGGGGCTTTACGT‑TTACCGTAAGGATGCCGGGAAGCGGGCGAATACTTTCTCCCGGCACGCCAAGGGT‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCaa                                                             <  2:667911/142‑1 (MQ=255)                                                                                   ggtgagAAAGGGATGGGGCTTTGGGT‑TTGTCGTGATGATGCCCGGAAGCGGGGGGATACTTTCTCCCGGCACGCAAAAGGT‑TA‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCaa                                                             <  2:1065926/143‑1 (MQ=255)                                                                                         aaggggAGGGGGCGTTCCGG‑TTATGGTGATGTTCCCCTGAGGCGGGCGGATCCTTTCTCCCGGCAC‑CCAAAGGGTTT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTcc                                                     <  2:373231/147‑1 (MQ=255)                                                                                           aGGGATGGAGCGTTACGG‑TTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTgg                                                  >  2:1079368/1‑151 (MQ=255)                                                                                               aTGGAGCGTTACGG‑TTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTg                                              >  2:165852/1‑151 (MQ=255)                                                                                                  gggCGTTACGT‑TTATCGTGATGATGCCCTGAAGCGGGGGGATACGTTCTCCCGGCCCGCAAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCcg                                           <  2:33656/149‑1 (MQ=255)                                                                                                      gTTACGG‑TTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCgcgc                                        >  1:645226/1‑150 (MQ=255)                                                                                                          cGG‑TTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTg                                   >  2:972402/1‑151 (MQ=255)                                                                                                          acgtttttctggatgttGCCCTGAGGCGGGCGGATACTTTCTCCCGGCACGCAAAAGGT‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCAtt                                     <  2:692224/139‑1 (MQ=255)                                                                                                            g‑TTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGt                                  >  1:946194/1‑150 (MQ=255)                                                                                                            atttATCTTGATTGCGCCCTGAAGCTGGCGGTTTTTTTCTCCCGGCACCCCAAGGGT‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCAttt                                    <  1:95555/147‑1 (MQ=255)                                                                                                              ttATGGTGATGATCCCCTGAGGCTGGCGGATCCTTTCTCCCGGC‑CCCCAAAGGGTTT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCAttt                                    <  2:362800/147‑1 (MQ=255)                                                                                                              ttATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCa                                                              >  2:1010185/1‑121 (MQ=255)                                                                                                                aTCGTAATGATGCCCTGAGGCTGGCGGATACTTTCTCCCGGCAC‑CCAAAGGGTTT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGtata                               <  2:609726/150‑1 (MQ=255)                                                                                                                   gTGATGATCCCGGGAGGGTGGCGGTTACGTTCTCCCGGCACGCAAAAGGT‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATAtgtt                           <  2:533128/151‑1 (MQ=255)                                                                                                                      acgatgCCCTGAGGCGGGCGGATCCGTTCTCCGGACACGCCAAGGGG‑TA‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGtat                                <  1:741717/141‑1 (MQ=255)                                                                                                                         atgCCAGGAAGCTGGCGGTTCTGTCCTCCCGGCACGCAAAGGGT‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATAtgttgt                         <  1:120260/147‑1 (MQ=255)                                                                                                                           gCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACg                   >  1:1053679/1‑151 (MQ=255)                                                                                                                           gCCCTGAAGCTGGCGGATACGTTCTTCCGCCACGCCAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACg                   >  1:277796/1‑151 (MQ=255)                                                                                                                            cccTGAGGCGGGCGGTTACTTTCTCCCGGCCCGCAAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAAc                    <  2:701869/149‑1 (MQ=255)                                                                                                                            cccGGAGGCGGGCGGTTCCTTTCTCCGGCCCCGCAAAAGGT‑TA‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATAtgtt                           <  2:533156/142‑1 (MQ=255)                                                                                                                            cccGGAAGGTGGGGGATACGTTCTCCCGCCACGCAAAGGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGtata                      <  1:726412/147‑1 (MQ=255)                                                                                                                             ccgaagggggggggATCCTTTCTCCGGCCACGCAAAGGGG‑TA‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGtat                       <  2:23567/133‑1 (MQ=255)                                                                                                                             ccgaaggggggggaaTACTTCCTCCCGGCCCTCAAAGGGT‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGta                        <  1:716445/130‑1 (MQ=255)                                                                                                                                gacgccggcggcTCCGTTCTTCGGGCACCCCAAAGGA‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGAtt                <  1:1074082/143‑1 (MQ=255)                                                                                                                                 aaGCGGGCGAATCCTTTCTCCCGCCACGCAAAGGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGAttt               <  2:135911/149‑1 (MQ=255)                                                                                                                                    ttggtgcattttttctccCGCCCCGCCAAAGTA‑AG‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGTCAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGAtttttt            <  1:510460/137‑1 (MQ=255)                                                                                                                                     gggCGGTTCCGTTTTCCCGGCAC‑CCAAAGGGTTT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCg         <  1:269187/150‑1 (MQ=255)                                                                                                                                       gCGGATACTTTCTCCCGCCACGCAAAGGGT‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTa       <  2:606003/151‑1 (MQ=255)                                                                                                                                           aTACGTTCTTCCGGCACGCCAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGc    >  2:1038128/1‑150 (MQ=255)                                                                                                                                                ttctccCGGCGCGCCAAGGGT‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCaa  <  2:142540/147‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 GGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGG‑TTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGG‑TT‑AACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAA  >  NC_000913/3230948‑3231236

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