Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,151,592	0	G	T	41.1%	89.3 / 39.6	56	L136F (TTG→TTT)	yegJ	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (4/29);  new base T (23/0);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 5.54e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGTTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTT  >  NC_000913/2151442‑2151741                                                                                                                                                        |                                                                                                                                                      cATGTTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTg                                                                                                                                                       <  1:324074/151‑1 (MQ=255)         ttGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑AACCTTTATTCCTTATTTGTCAc                                                                                                                                                  >  2:94059/1‑148 (MQ=255)           ggCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATc                                                                                                                                             <  1:609863/151‑1 (MQ=255)           ggCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATc                                                                                                                                             <  2:935663/151‑1 (MQ=255)            gCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGGTAT‑TACAATTTTCCTTTTTTTGTTACCt                                                                                                                                                >  1:119466/1‑147 (MQ=255)             cGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCa                                                                                                                                           <  2:89926/151‑1 (MQ=255)               tAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑AACCTTGATTCTTTTTTTGTTAcc                                                                                                                                                 >  1:497259/1‑143 (MQ=255)                aaTCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCaaaa                                                                                                                                        <  2:958767/151‑1 (MQ=255)                aaTCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCaaaa                                                                                                                                        <  1:451762/151‑1 (MQ=255)                  tCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGGTTT‑AACATTTATTCTTTTTTTGTCGCCAATCCaaaa                                                                                                                                        >  1:38430/1‑149 (MQ=255)                   cAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAaaga                                                                                                                                      <  2:761481/150‑1 (MQ=255)                      cAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAaagaagaa                                                                                                                                  <  2:1024451/151‑1 (MQ=255)                         tGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ATCATTTATTCTTTTTTTGTTACCTCTTCAAaaaaagaa                                                                                                                                  >  2:116790/1‑148 (MQ=255)                           aCCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCcga                                                                                                                             <  2:641395/151‑1 (MQ=255)                           aCCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑AACATTTATCCTTTTTTTGTCTCCTTTTCaaaaaaaaaa                                                                                                                                  >  1:1009013/1‑146 (MQ=255)                                  tATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑AACATTTATCCTTTTTTTGTTACCTATCCAAAAAAAGAAACCGACGAg                                                                                                                         >  2:962352/1‑148 (MQ=255)                                         aCAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCaa                                                                                                               <  1:714868/151‑1 (MQ=255)                                         aCAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCa                                                                                                                <  2:726994/150‑1 (MQ=255)                                           aaaCGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑AACATTGTTTCTTTTTTTGTCCCCTCTCCAAAAGAAGAAAACGACGAGCTAAAAGaaaa                                                                                                              >  1:1214703/1‑150 (MQ=255)                                                tAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTT‑AACATTTATCCTTTTTTTGTCGCCTTTCCCAAAAAAAAAAACGGCCAGCTAAGAGGaaaa                                                                                                             >  2:1210271/1‑146 (MQ=255)                                                  gAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTc                                                                                                      <  1:1251498/151‑1 (MQ=255)                                                    atatGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACCATT‑ATACTTTTTTTGTTCCCTATTCCAAAAAAGAAAACGCCGAGCTAAAAGCAAAAATAAtt                                                                                                       >  2:997638/1‑148 (MQ=255)                                                      atGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCt                                                                                                   <  2:528969/150‑1 (MQ=255)                                                       tGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑AACAATGTTTCTTTTTTTGTTACCCAACCAAAAAAAAAAACCCACGAGAAAAAAGAAAAAATAATTTGCt                                                                                                   >  2:1087147/1‑149 (MQ=255)                                                          ccAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGGTTT‑ACCATTTATCCTTTATTTGTTACCCATCCAAAAAAAAAAAACGACGGGCTAAGAGGAAAAATAATTTGCTaa                                                                                                 >  1:800514/1‑148 (MQ=255)                                                           cAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTATTTCGCTAATCAg                                                                                             <  1:146610/151‑1 (MQ=255)                                                            aGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGc                                                                                            <  1:441687/151‑1 (MQ=255)                                                            aGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑AACATTTATCCTTTTTTTGTGACCTATTCCAAAAAAAAAAACGACGAGCTAAGAGAAAAAATatatttct                                                                                                   >  1:170814/1‑137 (MQ=255)                                                               aGTGGAAAATCACGAAAGAAGATATGTCAGACTGGAGGTATACGCGCGGGGATAAAAATTACGGAGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGATGAACGCGACGAGCTAAGCGCAAAATTAGTTCGCTTAACAGTCg                                                                                          <  1:1045515/150‑1 (MQ=255)                                                               aGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑TTCCTTTTTCCTTTTTTTGTCACCTTTCCAAAAAAAGAAACCGACAAGCTAAAAGGAAAATAAGTTTGCATATCa                                                                                              >  2:178787/1‑146 (MQ=255)                                                                gTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTT‑ATCATTGTTCCTTTTTTTGTTACCTATCCAAAAAAAAAAAACGACGAGGTAAGAGCAAAATTAATTTTCTAAtaaacc                                                                                           >  1:655181/1‑143 (MQ=255)                                                                 tGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTT‑TTCATTTATCCTTTTTTTGGGACACCTCCAAAAAAAAAAGCCCCCGAGCTCAGAGCAAAAATAGTTCTCTAATCAAcccgt                                                                                        >  2:699280/1‑147 (MQ=255)                                                                    aaaaTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTatat                                                                                    <  2:1066546/151‑1 (MQ=255)                                                                      aaTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTAtattt                                                                                  <  2:188431/151‑1 (MQ=255)                                                                       aTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCAGTTAT‑ACCACTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAACCGGCGAGCTAAGAGCAAAATTAATTCGCTAATCAGCCGTTAtatttt                                                                                 <  1:904228/151‑2 (MQ=255)                                                                         cAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTT‑TACATTGATTCTTTATTTGTGCCCTATTCAAAAAAAAAAAACGAAGAGGTAAGAGCAAAAATAAATTCTTAATCAAACGGTAtatttttt                                                                               >  1:148164/1‑151 (MQ=255)                                                                         cAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTAtatttatt                                                                               <  2:136407/151‑1 (MQ=255)                                                                          aagaaagaAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ATCATTTATTCTTTTTTTGTTACCTATTCAAAAAAAAAAAACGGACAGCTAAAAGGAAAAATAATTTGCTAATTAACCGGTAtatttttt                                                                               >  2:891624/1‑150 (MQ=255)                                                                           agaaagaaGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGt                                                                             <  2:497260/151‑1 (MQ=255)                                                                             aaagaaGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTT‑TACAATTATTCTTTATTTGTCACACATCCAAAAAAAAAAAGCGGCGGGACAAAAGGAAAATAAGTTTGCTAATAAACCGTTTTAtttttttttt                                                                           >  2:175053/1‑151 (MQ=255)                                                                               agaagaTATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTT‑AACAATTATTCTTTTTTTGGTGCCCATTCAAAAAAAAAAAACGACGAGGTAAAAGCAAAAATAATTTTCTTAAAAACCCCTATATTTAtttttt                                                                           >  1:353762/1‑149 (MQ=255)                                                                                gaagaTATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTaaa                                                                        <  2:1115469/151‑1 (MQ=255)                                                                                   gaTATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTg                                                                     <  2:381310/151‑1 (MQ=255)                                                                                     tatTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑TACAATTATTCTTTATTTGTGCCCCATCCAAAAAAAAAAAACCACGAGGTAAAAGCAAAAATAATTTGCTTATCAACCCGTATATTTTTTTTTTtaa                                                                        >  2:312953/1‑146 (MQ=255)                                                                                      atTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCg                                                                  <  2:8611/151‑1 (MQ=255)                                                                                       tttCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTT‑TACATTATTTCTTTCTTTGTTACCTTTCCCAAAAAAAAAAGCGACGCGGCAAGAGCAAAAATAATTCGTTAATCAAGCGATAAATTTTTTTTTTTAAATg                                                                     >  1:750828/1‑147 (MQ=255)                                                                                        ttCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTT‑ACCAATTAACCCTTTTTTGTTACCCATCCAAAAAAAGAAACCGCCCAGATAAGAACAAAAATATTTTGCTTAACAACCGGTTTATTTTTTTTTTTAAATAGCgg                                                                 >  2:487323/1‑150 (MQ=255)                                                                                          cAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGtt                                                               <  2:986132/150‑1 (MQ=255)                                                                                          cAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTg                                                              <  2:1129917/151‑1 (MQ=255)                                                                                            gACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGaa                                                            <  1:127263/151‑1 (MQ=255)                                                                                                               gggATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCa                                         <  2:211852/151‑1 (MQ=255)                                                                                                                  aaaaatTTTAGGGCGTTTAT‑ACCATTGATCTTTTATGGGTCCCCTATCAAAAAGAAGAAGCCGACGAGCAAAGAGAAAATTTAGTTCGTTAACCAGCCTTTATATTTATTGTTTAAACGGGCGGTTGATTTGTTGTCTTCAACCGCCATg                                       <  1:559083/148‑1 (MQ=255)                                                                                                                               cggTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAactcactc                           >  2:301885/1‑149 (MQ=255)                                                                                                                                 gTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCCTGAGAAAATCACTCTGcc                       >  2:1207670/1‑151 (MQ=255)                                                                                                                                  ttAT‑CCAATTGATCTTTTATTGGCCACCTACCAAAAAGAAAAACCCGACAAGCAAAGAGAAAATTAAGTTCGCTAACCAGCCTTTATATTTATTTTTAAAAGTGGGGGTTAAATTGTTTTCTTCAACCGCCATGAGAAACTCACTCTg                         <  1:712088/148‑1 (MQ=255)                                                                                                                                       aCCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCcac                    >  1:1109753/1‑149 (MQ=255)                                                                                                                                                       tGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGtt  >  1:1008608/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                      CATGTTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTT  >  NC_000913/2151442‑2151741

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