Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,151,589	0	A	T	40.0%	68.5 / 24.4	50	L135F (TTA→TTT)	yegJ	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (3/27);  new base T (20/0);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 3.76e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

CATGTTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAA  >  NC_000913/2151442‑2151737                                                                                                                                                      |                                                                                                                                                        cATGTTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTg                                                                                                                                                      <  1:1053143/151‑1 (MQ=255)  aTGTTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTTT‑TACATTTATTCtttttt                                                                                                                                                       >  2:309886/1‑149 (MQ=255)  aTGTTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTgg                                                                                                                                                     <  1:940538/151‑1 (MQ=255)    gTTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTTT‑ACCATTGATTCTTTATTTGt                                                                                                                                                    >  2:582725/1‑150 (MQ=255)       tATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑AACATTTATCCttttttttttacc                                                                                                                                                >  2:904570/1‑144 (MQ=255)          tGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTTTTCCTTTTTTTGTTAcc                                                                                                                                                >  1:204719/1‑148 (MQ=255)            gCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTT‑ACCATT‑ATTCTTTTTTTGTTAcc                                                                                                                                                >  1:1010359/1‑146 (MQ=255)              gTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ATAATTTATTCTTTTTTTGTTAc                                                                                                                                                 >  1:222170/1‑143 (MQ=255)                  tCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ATCAATGTTTCTTTTTTTGTTACCTA‑T‑CCa‑aaaa                                                                                                                                   >  1:273344/1‑149 (MQ=255)                      cAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑AACCTTTTTCCTTTTTTTGTCTACTA‑T‑CCA‑aaaaaa                                                                                                                                 >  1:574288/1‑148 (MQ=255)                       aaTGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTTT‑ACCATTTATTCTTTTTTTGTCACCTA‑T‑CCa‑aaaaaaaaa                                                                                                                              >  2:444289/1‑150 (MQ=255)                        aTGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGc                                                                                                                            <  2:1195585/151‑1 (MQ=255)                            cccTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑GTATAACAATTAATCTTTTTTTGG‑GACCTC‑T‑CCACAaagaaaa                                                                                                                               >  1:1105938/1‑143 (MQ=255)                             ccTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCcgacg                                                                                                                       <  2:1009634/151‑1 (MQ=255)                                   aTTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACGAGCTa                                                                                                                  <  2:789130/150‑1 (MQ=255)                                   aTTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑AACATTTATTCTTTTTTTGGTACCTATT‑CCA‑aaaaaaaa                                                                                                                               >  2:15559/1‑136 (MQ=255)                                      gTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTTT‑AACCTTGATTCTTTTTTTGTCAACTA‑T‑CCA‑AAAAAAGAAAGCCACGAGCTaa                                                                                                                 >  2:1062195/1‑148 (MQ=255)                                        tACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACGAGCTAAGAGc                                                                                                             <  2:6616/150‑1 (MQ=255)                                           aaaCGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑ATAT‑ACCGTTTATGCTTTTTTGGTGACCTA‑T‑TCACAAAGAA‑AAAACGACCAGCTaataaca                                                                                                            >  1:904468/1‑143 (MQ=255)                                                  gAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTc                                                                                                  <  1:998465/151‑1 (MQ=255)                                                     tatGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑AACATTTATTCTTTTTTTGTTACCTA‑T‑CCA‑AAAAAAAAAGACGACGAGCTAAAAGCAAAAATAGTTTGCt                                                                                               >  2:355007/1‑151 (MQ=255)                                                      atGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCt                                                                                               <  2:1243091/150‑1 (MQ=255)                                                        ggCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTTT‑AACCTTTATCCTTTTTTTGTTACCTA‑T‑TCC‑AAAAAAAAAACCCACGAGCTAAAAGGAAAATTAGTTTGCTaa                                                                                             >  2:1172994/1‑150 (MQ=255)                                                        ggCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAAt                                                                                            <  1:601335/151‑1 (MQ=255)                                                          ccAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑AACATTTATTCTTTTTTTGTTACCCA‑T‑ACC‑AAAAAAAAAACCGACCAGCTAAGAGGAAAAATAAtt                                                                                                   >  1:260421/1‑142 (MQ=255)                                                           cAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTTATCCTTTTTTTGTTACCCA‑T‑CCA‑AAAAAAAAAAACGACGAGCTAAAAGAAAAAATAATTTGCt                                                                                               >  2:1075494/1‑145 (MQ=255)                                                            aGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGc                                                                                        <  2:253516/151‑1 (MQ=255)                                                                       aTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑CTAT‑AACATTTATCCTTTTTTTTTCACCCA‑TCCCA‑AAAAAAAAAGCC‑ACGGAGGAAGAACAAAAATAGTTTGCCAATCAGTCGGTAtattt                                                                              >  2:305184/1‑150 (MQ=255)                                                                         cAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTAtatttatt                                                                           <  1:719890/151‑1 (MQ=255)                                                                         cAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTAtatttatt                                                                           <  2:95157/151‑1 (MQ=255)                                                                         cAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTAtatttatt                                                                           <  1:597654/151‑1 (MQ=255)                                                                          aagaaagaAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTGTATTTATTg                                                                          <  1:1158840/151‑1 (MQ=255)                                                                          aagaaagaAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTg                                                                          <  2:1166812/151‑1 (MQ=255)                                                                            gaaagaaGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCCTTTATCCTTTTTTTGTTACCCA‑T‑CCA‑AAAAAAAAAACCGACGAGCTAAAAAGAAAAATAATTTGCTAATCAGACGTTTTATTTATtttt                                                                        >  2:321149/1‑151 (MQ=255)                                                                              aagaagATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTa                                                                      <  2:433543/151‑1 (MQ=255)                                                                              aagaagATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑AACAATGATTCTTTTTTTGGGACCTA‑T‑ACA‑AAAAAAAAAAACGAAGAGGTAAGAGGAAAAATAGTTCGCTAATAAGACGGTATATTTATTGttt                                                                       >  2:1151047/1‑150 (MQ=255)                                                                               agaagaTATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTaa                                                                     <  1:1155164/151‑1 (MQ=255)                                                                                  agaTATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACt                                                                  <  2:962022/151‑1 (MQ=255)                                                                                    atatTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑GTTT‑AACATTTATCCTTTTTTTGTTACCTA‑T‑CCC‑AAAAAAAAAACCGAACGGCTAAGAGAAAAAATAATTTGCTAATCAGACGCTATATTTTTTGTTTtaa                                                                    >  2:438982/1‑147 (MQ=255)                                                                                       tttCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTTATCCTTTTTTTGTCACCTA‑T‑CCA‑AAAAAAAAAAGCCACGAGCTAAGAGCAAAAATAATTTTCTAATTAGCCCGTATATTTTTTTTTTTAAATTGCgg                                                             >  2:1240039/1‑151 (MQ=255)                                                                                                ggATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACTAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAAttg                                                     <  1:1045325/150‑1 (MQ=255)                                                                                                ggATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAAtt                                                      <  1:1105587/149‑1 (MQ=255)                                                                                                   tGTATACGCGCGGGGATAAAATTTACGGCGG‑TTTT‑AACATTTATCCTTTTTTTGTTACCTA‑T‑TCA‑AAAAAAGAAGCCGACGAGCTAAGAGCAAAAATAGTTTGCTAATCCACCGTTATATTTTTTTTTTTAAATGGCGGGTGTAttttt                                                   >  2:268780/1‑149 (MQ=255)                                                                                                    gTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATtgttgt                                                 <  1:261410/150‑1 (MQ=255)                                                                                                             cgGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGc                                       <  2:668443/151‑1 (MQ=255)                                                                                                                    aaaaTTTACGGCGT‑TTAT‑CCCATTGATCTTTTATGGTTCCCCTA‑T‑CAA‑AAAAAAAAAGCCGACGAGCAAAGAGAAAAATTATTTCGCTAATCAGCCGTTATATTTATTTTTAAAACGGGCGGTTGAATTGTTGTTTACAACCGCCATgag                                 <  1:411690/150‑1 (MQ=255)                                                                                                                                 gtttAT‑ACCATTGACCCTTTATTGGTCACCTA‑T‑CAA‑AAAGAAGAAGCCGACGAGATAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTTTTAAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTg                     <  2:712200/148‑1 (MQ=255)                                                                                                                                    tat‑aCCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGccc                  >  1:435076/1‑150 (MQ=255)                                                                                                                                              gttCCTTTATTGGTCCCCTA‑C‑CAA‑AAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGTAAATCAGCCTTTATATTTATTGTTTAACCGGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGcc          <  2:605159/147‑1 (MQ=255)                                                                                                                                                 cttttATGGGTCCCCTA‑C‑CAA‑AAAGAAGAAGCCGACGAGCAAAGAGCAAAATTAGTCCGTTAACCACCCTTTATTTTTTTTTTTTAAACTGGCGGTTGAATTGTTTTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAg        <  1:825500/146‑1 (MQ=255)                                                                                                                                                    ttATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGCTGAATTGTTGTCTTCCACCGCCATGAGAAACTCACTCTGCCCACAACGCCAggggca   >  2:1268101/1‑145 (MQ=255)                                                                                                                                                     tATTGGTCCCCTA‑T‑CAA‑AAAAAAGAAGCCACCGAGCTAAGAGAAAAATTAGTTCGCTAATCAGCCTTTATTTTTATTTTTTAAACGGGGGTTTGAATTTTTGTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCaa  <  2:908112/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                        CATGTTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGG‑TTAT‑ACCATTGATCCTTTATTGGTCACCTA‑T‑CCA‑AAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAA  >  NC_000913/2151442‑2151737

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