Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,151,589	0	A	T	43.1%	63.7 / 28.4	51	L135F (TTA→TTT)	yegJ	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (5/24);  new base T (20/2);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 1.23e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGTTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAG  >  NC_000913/2151442‑2151731                                                                                                                                                        |                                                                                                                                               cATGTTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCTGGGATAAAATTTACGGCGGGTA‑T‑A‑CCCTTTA‑TTCtttttt                                                                                                                                              >  1:782527/1‑150 (MQ=255)         ttGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGGTT‑T‑A‑TACTTTA‑TTCTTTATTTGTCAcc                                                                                                                                       >  1:438351/1‑149 (MQ=255)                aaTCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCaaa                                                                                                                               <  2:455990/150‑1 (MQ=255)                 aTCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGCTT‑A‑A‑TCACTCA‑TCCTTTTTTTGTTTCCTATcc                                                                                                                                  >  1:512791/1‑146 (MQ=255)                  tCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATT‑T‑CAATTTA‑TACTTTTTTGGTACCCAATCAaaaa                                                                                                                              >  1:62449/1‑150 (MQ=255)                  tCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGCTT‑T‑ACCCCTTCT‑TCCTTTTTTGGTCCCCTTTCAaaaa                                                                                                                              >  2:320475/1‑150 (MQ=255)                      cAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑T‑ACATTTA‑TTCTTTTTTTGTTACCTATCCAaaaaaa                                                                                                                           >  1:618387/1‑148 (MQ=255)                       aaTGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTT‑T‑A‑CCATTGA‑ACCTTTTTTTGTGACCCCTCCaaaaaaaaaa                                                                                                                        >  2:481513/1‑150 (MQ=255)                            cccTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTT‑T‑A‑ACAATGATTCCTTTTTTGTTC‑CCTCTCCaaaaaaaaaa                                                                                                                        >  2:30034/1‑145 (MQ=255)                            cccTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCcgac                                                                                                                  <  1:438896/151‑1 (MQ=255)                            cccTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑ACATTTA‑TTCTTTTTTTGTCACCCATCCaaaaaaaaaa                                                                                                                        >  1:411924/1‑145 (MQ=255)                             ccTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑A‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCcgacg                                                                                                                 <  2:404888/151‑1 (MQ=255)                                   aTTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTa                                                                                                            <  2:619105/150‑1 (MQ=255)                                      gTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑T‑TCATTTG‑TTCTTTTTTTGTCACCCATCCCAaaaaaaaaaacgcccagcaaa                                                                                                           >  2:570504/1‑135 (MQ=255)                                      gTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAaga                                                                                                         <  2:800747/150‑1 (MQ=255)                                        tACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGc                                                                                                       <  2:932673/150‑1 (MQ=255)                                         aCAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCa                                                                                                      <  2:814567/150‑1 (MQ=255)                                             aCGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAAtt                                                                                                 <  2:482366/151‑1 (MQ=255)                                               gTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGGTA‑T‑A‑CCATTGT‑TTCTTTTTTTTTTCCCCATCCCAAAAAAAAAGCCGACGAGCTAAAAGGaaaa                                                                                                   >  2:676213/1‑147 (MQ=255)                                                tAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGGTT‑T‑T‑ACATTGG‑TCCCTTTTTTGTCGCCCTTCCCAAAAAAGAAGGCGACCAGCCAAAAGAaaaa                                                                                                   >  2:568085/1‑146 (MQ=255)                                                 aGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGtt                                                                                             <  1:544970/151‑1 (MQ=255)                                                        ggCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTaa                                                                                       <  1:486743/150‑1 (MQ=255)                                                            aGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTT‑T‑A‑ACATTTT‑TACTTTTTTTGTCACCTATTCAAAAAAAGAAAACGGGGGGCCAAGAGCAAAAATATTTTTCTAAt                                                                                      >  1:419521/1‑147 (MQ=255)                                                              gAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCg                                                                                <  1:660908/151‑1 (MQ=255)                                                                    aaaaTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTatat                                                                          <  2:143936/151‑1 (MQ=255)                                                                    aaaaTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTatat                                                                          <  2:769184/151‑1 (MQ=255)                                                                     aaaTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTAtatt                                                                         <  1:875879/151‑1 (MQ=255)                                                                          aagaaagaAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTT‑T‑A‑ACATTTA‑TCCCTTTTTTGTCACCTATCCAAAAAAAAAAAACGAACAGCTAAGAGGAAAAATATTTTGCTAATCTAACGTTAtatttt                                                                       >  2:775610/1‑147 (MQ=255)                                                                            gaaagaaGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGtt                                                                  <  1:619897/151‑1 (MQ=255)                                                                              aagaagATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTT‑T‑A‑CCATTTA‑TTCCTTATTTTGTACCTCTTCAAAAAAAGAAACCGACGAGGTAAGAGCAAAAATAGTTTTCTTATTAGCCCGTTTAtttttttttt                                                                 >  2:384834/1‑150 (MQ=255)                                                                              aagaagATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTa                                                                <  1:610139/151‑1 (MQ=255)                                                                               agaagaTATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCCTTTG‑ACCTTTATTTGTTACCCTTCCCAAAAAAAAAAACCACGGGCAAAAAGCAAAAATAATTCCGTAATAAGCCGGTATATTTAtttttt                                                                 >  2:771482/1‑149 (MQ=255)                                                                                 aagaTATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑ACAATGA‑TCCCTTATTTGTCACCCATCCAAAAGAAGAAACCGACGAGCTAAAAGAAAAAATAGTTTGCTAATCAGCCGGTATATTTTTTGTTTtaa                                                              >  2:665688/1‑150 (MQ=255)                                                                                      atTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGc                                                         <  2:113872/150‑1 (MQ=255)                                                                                      atTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCg                                                        <  2:844160/151‑1 (MQ=255)                                                                                        ttCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTT‑T‑A‑CCATTTA‑TCCTTTTTTTGTTACCCATCCAAAAAAAAAAACCGACGAGCTAAAAGCAAAAATAGTTTGCTAAACAACCGCTATATTTTTTTTTTTAAATTGCgg                                                       >  2:167481/1‑150 (MQ=255)                                                                                         tCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑T‑CCATTGA‑ATTTTTTTTTGTCAACTATTCACAAAAAAAAAACGACGAGGTAAAAGCAAAAATATTTCTCCTATCAGCCGCTATATTTTTTTTTTTAAAAGGCggg                                                      >  1:240588/1‑149 (MQ=255)                                                                                          cAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTg                                                    <  2:49818/151‑1 (MQ=255)                                                                                           aGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑TCATTTA‑TTCTTTTTTTGTCACCTATCCAAAAGAAAAAACAGACGAACTTAGAGGAAAAATATTTTGCTAATCAGCCGTTTTTTTTTTTTTTTTAACTCGCgg                                                       >  1:26122/1‑147 (MQ=255)                                                                                              cTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑TCATAGA‑TTCTTTTTTTGGCTCCTATCCAAAAAAAAAAAACGACCAGCTAAGAGCAAAATTAGTTTGCTAATCAACCGTTATTTTTTTTTTTTtaaatatcgg                                                       >  1:799416/1‑137 (MQ=255)                                                                                                ggATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑TAA‑CAATTGT‑TCCTTTTTTGGTCCCCATTCCAAAA‑AAGAAGACGAAGAGCTCAGAGGAAAAATAATTTGCTTATTAGACGATATATTTTTTGTTTtaaatagggggt                                                     >  1:6660/1‑135 (MQ=255)                                                                                                 gATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTTT‑TCCTTTTTTTGTCACCTATCCAAAAAAAAAAACCGACGAGGTAAAAGCAAAAATAATTTGCTAAACAGCCGTTATATTTTTTTTTTTAAATGGCggg                                                      >  2:421673/1‑141 (MQ=255)                                                                                                 gATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑ACACTTT‑TCCTTTTTTTGTCACCAATACCAAAAAAAAAAACGGCCAGATAAGCGCAAAAATATATCGCTCATTATCCGCTATTTTTATTTGTTATACTGGGGGGTggaatttt                                             >  2:752213/1‑144 (MQ=255)                                                                                                    gTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTTTTTGTCACCCATCCCAAAAAAGAAAGCGACGGGCTAAGAGCAAAAATATTTCGCTAATCAGCCGGTATATTTATTTTTTTAACATGCGGGTGTAttttt                                             >  2:218906/1‑148 (MQ=255)                                                                                                     tataCGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATtgttgt                                           <  2:102822/149‑1 (MQ=255)                                                                                                        aCGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTc                                       <  1:870279/150‑1 (MQ=255)                                                                                                        aCGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCa                                      <  1:481498/151‑1 (MQ=255)                                                                                                                             ggcggTTA‑T‑A‑CCATTGA‑TCTTTTTTTGGTCCCCTACCAAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCACCCTTTATATTTATTGTTAAAACGGGGGGTTGAATTGTTGTCTTCAACCGCCATGAGAAactcact                  <  2:508649/150‑1 (MQ=255)                                                                                                                                ggTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAAATGTTGTTTTCAACCGCCCTGAGAAACTCACTCTGc              >  2:703427/1‑151 (MQ=255)                                                                                                                                          ccATTGA‑TCTTTTTTTGGTCCCCTATCCAAAAGAAAAAGCCGCCGAGCTAAGAGAAAATTTATTTCGCTAACCAGCCGTTTTATTTTTTTTTTAAACGGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAg        <  2:849598/150‑1 (MQ=255)                                                                                                                                             ttGA‑TCTTTTATTGGTCCCCTATCAAAAAGAAGAAGCCGACGAGCAAAGAGAAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTAAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACCCACTCTGCCCACAgcg      <  2:688539/149‑1 (MQ=255)                                                                                                                                                     tttATGGGTCACCTATCAAAAAGAAGAAGCCGACGAGCTAAGAGAAAAATTAGTTCGTTAATCAGCCGTTATTTTTATTGTTTAAACGGGCGGTTGAATTTTTGTCTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAg  <  2:98192/146‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                               CATGTTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTA‑T‑A‑CCATTGA‑TCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAG  >  NC_000913/2151442‑2151731

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