Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,151,608	0	G	A	43.8%	42.3 / 29.3	50	E142K (GAA→AAA)	yegJ	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (10/17);  new base A (16/5);  total (28/22)
Fisher's exact test for biased strand distribution p-value = 9.41e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTG‑A‑TCCTTTATTGGTCACCT‑A‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATC  >  NC_000913/2151470‑2151752                                                                                                                                               |                                                                                                                                                 ccTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTACCATTT‑T‑TTCTTTTTTTGGCTCCC‑A‑TCCCAaaaaaaaaaacga                                                                                                                                        >  1:96681/1‑144 (MQ=255)    ggTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTG‑A‑TCCTTTATTGGTCACCT‑A‑TCCAAAAGAAGAAGCcgacga                                                                                                                                     <  1:476069/149‑1 (MQ=255)      tATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGGTTTAACCTTT‑A‑TTCTTTATTTGTTGCCC‑A‑TCCAAAAAAAGAAGAcggcga                                                                                                                                     >  2:257897/1‑147 (MQ=255)          gTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTG‑A‑TCCTTTATTGGTCACCT‑A‑TCCAAAAGAAGAAGCCGACGAGCTAagag                                                                                                                             <  1:642630/151‑1 (MQ=255)          gTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGGTTTACCATTG‑A‑TTCTTTTTTTGTTACCT‑A‑TCCCAAAAAAGAAGCCCACGAggtca                                                                                                                                >  1:339765/1‑144 (MQ=255)             aCAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTAACATTT‑T‑TTCTTTTTTTGTGACCT‑ATTCCAAAAAAAAAAGcggcgggctaa                                                                                                                                 >  1:619747/1‑135 (MQ=255)                    tAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGCTATATCAATGTA‑TCCTTTTTTGGCAGCCT‑A‑TCAAAAAAAAGAA‑CCGGCGCGCTAAAAGGaaaa                                                                                                                        >  2:318239/1‑146 (MQ=255)                     aGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTG‑A‑TCCTTTATTGGTCACCT‑A‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGtt                                                                                                                  <  2:707441/151‑1 (MQ=255)                          atGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTG‑A‑TCCTTTATTGGTCACCT‑A‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCt                                                                                                              <  2:634802/150‑1 (MQ=255)                          atGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGATATACCATTT‑ATTCCTTATTTGGTAACAA‑A‑ACCAAAA‑AAGAAAGCGGCGGGCTTAGAAAAAAAATAGtt                                                                                                                  >  2:327511/1‑146 (MQ=255)                           tGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTAACATTT‑A‑TCCTTTTTTTGGTACCTAA‑TCCAAAA‑AAGAAAGCGACGAGCCAAAAGCAAAAATAGt                                                                                                                   >  2:585855/1‑144 (MQ=255)                           tGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCAATT‑A‑TTCTTTTTTTGTTACCC‑A‑TCCAAAAAAAGAAAGCGACCAGCTAAAAGCAAAAATAATTTGCTa                                                                                                             >  2:225007/1‑150 (MQ=255)                            ggCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTG‑A‑TCCTTTATTGGTCACCT‑A‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTaa                                                                                                            <  1:253757/150‑1 (MQ=255)                            ggCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTG‑A‑TCCTTTATTGGTCACCT‑A‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTaa                                                                                                            <  2:87629/150‑1 (MQ=255)                                aGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTAACATTT‑A‑TTCTTTTTTTGTTACCT‑C‑TCCAAAAAAAAAAAACGCCGCGCTAAAAGGAAAAATAGTTTGCTTATAAg                                                                                                        >  1:45747/1‑150 (MQ=255)                                aGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAAAAATT‑A‑TTCTTTATTTGGTACCC‑A‑TCCAAAAAAAAAAAACGGCCAGCGAAAAGCAAAAATAATTTGCt                                                                                                              >  2:235452/1‑144 (MQ=255)                                      ggAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTT‑A‑TCCTTTTTTTGGGTCCC‑A‑TCCCAAAAAAAAAAACGACGAGCTAAGAGCAAAAATAGTTTGCTAAACAAcc                                                                                                      >  2:884956/1‑146 (MQ=255)                                        aaaaTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTTTCATTT‑T‑TCCTTTTTTTGGTTCCT‑A‑TCCCAAAAAAAAAAACGGCGAGCTAAGAGGAAAAATAAATTTCTAATTATCCGtttttt                                                                                               >  1:319848/1‑147 (MQ=255)                                        aaaaTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTG‑A‑TCCTTTATTGGTCACCT‑A‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTatat                                                                                               <  2:312026/151‑1 (MQ=255)                                         aaaTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACATTT‑A‑TCCTTTTTTTGTTTCCT‑A‑TTCAAAAAAAGAAGCCGACCAGCCAAAAGCAAAATTATTTTTCTAATCAGCCGTTAtatt                                                                                              >  1:523002/1‑151 (MQ=255)                                         aaaTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACATTG‑T‑TTCCTTTTTGGTCACCC‑A‑TTCAAAAAAAAAAAGCGACGAGCGAAAAGGAAAAATAAATTGCTaaataaccggtaaat                                                                                               >  2:909749/1‑137 (MQ=255)                                           aTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATATCATTT‑A‑TTCTTTTTTTGGCACCT‑A‑TCCAAAAAAAAAAGACGACGAGCTAAGAGGAAAAATAATTTGCTAATTAACCCCTAtattt                                                                                             >  2:23992/1‑150 (MQ=255)                                           aTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTG‑A‑TCCTTTATTGGTCACCT‑A‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTAtattta                                                                                            <  1:608385/151‑1 (MQ=255)                                           aTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTG‑A‑TCCTTTATTGGTCACCT‑A‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTAtattta                                                                                            <  2:846596/151‑1 (MQ=255)                                                gaaagaaGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTG‑A‑TCCTTTATTGGTCACCT‑A‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGtt                                                                                       <  1:232008/151‑1 (MQ=255)                                                 aaagaaGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTT‑A‑TACCTTTTTGGTGTCCC‑A‑TTCAAAAGAAAAAGCCGACGAGCTAAGAGCAAAAATAATTCTCTTATCAACCCGTATAtttttttttt                                                                                      >  2:294672/1‑151 (MQ=255)                                                  aagaagATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTG‑A‑TTCTTTATTTGTCACCT‑A‑TCCAAAAAAAAAAACCGACGGGCTAAAAGCAAAAATAAGTTGCTAATTAGCCCGTATAtttttttttt                                                                                      >  1:137758/1‑150 (MQ=255)                                                    gaagaTATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTG‑A‑TCCTTTATTGGTCACCT‑A‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTaaa                                                                                   <  2:761620/151‑1 (MQ=255)                                                       gaTATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGCTATACCATTG‑A‑TTCTTTATTTGTCCCCT‑A‑TTCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTTATCAGCCGTTATATTTATTGTTTTAAATg                                                                                >  2:617850/1‑151 (MQ=255)                                                                     gATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTG‑A‑TCCTTTATTTGTCACCC‑A‑TCCAAAAAAAGAAACCGACGAGCTAAAAGCAAAAATAGTTCGCTAATCAACCGTTATATTTTTTTTTTTAAATGGCGGGTggaatttt                                                                  >  1:619242/1‑144 (MQ=255)                                                                        gTATACGCGCGGGGATAAAATTTACGGCGGTTATAACATTT‑A‑TTCTTTTTTGGGCACCT‑A‑TCCAAAAAAAGAAACCGACGAGCTAAGAGCAAAAATAATTCGCTAAACAAGCGGTATATTTTTTTTTTTAAATGGGGGGTggaattttt                                                                 >  1:217981/1‑141 (MQ=255)                                                                            aCGCGCGGGGATAAAATTTACGGCGGTTTTACCATTT‑A‑TCCTTTTTTTGTGACCC‑A‑TCCAAAAAAAGAAGCCCACGAGCTAAAAGAAAAAATAATTTGCTAATAAACCGGTATATTTTTTTTTTTAAATTGCGGGTTTAttgttttttt                                                             >  2:805008/1‑144 (MQ=255)                                                                            aCGCGCGGGGATAAAATTTACGGCGGTTATACCATTG‑A‑TCCTTTATTGGTCACCT‑A‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCa                                                           <  2:527937/151‑1 (MQ=255)                                                                                    gaaaaaaaTTTACGGCGGTTATACCTTTG‑A‑TCTTTTATTGGTCACTT‑A‑CCAAAAAGAAGAAGCCGACGAGCAAAGAGAAAAATTAGTTCGTTAATCAGCCGTTATATTTATTTTTTAAACGGGCGGTTGAATTGTTGTCTTCAACCGCCAt                                                   <  2:477963/147‑1 (MQ=255)                                                                                            tttACGGCGTTTATCCCTTTA‑T‑TCTTTTATTGGTCCCCT‑A‑TCCAAAAAAAGAAGCCGACGAGCTAAGAACAAAATTATTTCGCTAATCAGCCTTTATTTTTATTGTTAAAGCTGGGGGTTGATTTTTTGTCTTCAACCGCCATGAGAAac                                            <  2:729733/150‑1 (MQ=255)                                                                                                 ggcggTTATCCCTTTG‑A‑TCTTTTATTGGTCCCTT‑A‑CCAAAAAAAAAAAGCCGACGAGCTAAGAAAAAAATTATTCCTCTAATCAGCCGTTATTTTTATTGTTTAACCTGGCGGTTAATTTGTTGTTTTCAACCGCCATGAGAAactc                                          <  2:45743/147‑1 (MQ=255)                                                                                                    gtttATACCATTG‑A‑CCCTTTATTGGTCCCCT‑A‑CCAAAAAGAAGAAGCCGACGAGCAAAGAGAAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTAAAACTGGCGGTTGAATTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGc                                   <  1:909842/149‑1 (MQ=255)                                                                                                    ggTTATACCATTG‑A‑TCCTTTATTGGTCACCT‑A‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTg                                    >  1:433199/1‑150 (MQ=255)                                                                                                    ggTTATACCATTG‑A‑TCCTTTATTGGTCACCT‑A‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTg                                    >  2:155547/1‑150 (MQ=255)                                                                                                      ttATCCCTTTG‑A‑TCTTTTTTTGTTCCCCT‑A‑CCAAAAAAAAGAAGCCGACGAGCAAAGAGAAAAATTAGTTCTTTAATCACCCTTTATTTTTATTGTTAAACCTGGCGTTTGATTTTTTTTCTTCAACCGCCATGAGAAACTCACTCTg                                    <  2:377583/148‑1 (MQ=255)                                                                                                      ttATACCATTG‑A‑TCCTTTATTGGTCACCT‑A‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTGTTCAAACGCCATGGAAAAATTACTCTGccc                                 >  2:488209/1‑151 (MQ=255)                                                                                                           ccATTG‑A‑TCTTTTATTGTTCACCT‑A‑TCAAAAAAAAAAGGCCAAGGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCTTTATATTTATTGTTTAAACTGGCGTTTAAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAg                             <  1:86849/150‑1 (MQ=255)                                                                                                                      cTTTATTGGTCACCT‑A‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCCGTGcc                   >  1:51805/1‑151 (MQ=255)                                                                                                                       tttATTGGTCACCT‑A‑CCCAAAAAAAGAAGCCGACGAGCAAAGAGCAAAATTAGTTCGCAAACCAGCCTTTATATTTATTGTTTAAACTGGGGGTTGAATTGTTGTCTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTg                     <  1:478308/148‑1 (MQ=255)                                                                                                                       tttATTGGCCCCCT‑T‑TCAAAAAGAAGAAGCCGAGGAGCAAAGAGAAAAATTAGTTCGCTAATCAGCCTTTATATTTATTGTTAAAACGGGGGTTTAAATTGTTGTTTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGt                      <  1:8389/147‑1 (MQ=255)                                                                                                                        ttATTGGTCACCT‑A‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCCGTGCCa                  >  1:573955/1‑150 (MQ=255)                                                                                                                          aTTGGTCACCT‑A‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCCCAGCGCCAGTGCCaa                 >  2:641756/1‑149 (MQ=255)                                                                                                                              gTCACCT‑A‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCCACGtttt           >  2:949846/1‑151 (MQ=255)                                                                                                                                  ccT‑A‑CCAAAAAGAAGAAGCCGACGAGCTAAGAAAAAAATTAGTTCGCTAATCAGCCTTTATATTTATTGTTAAAACTGGCGGTTAATTTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatc        <  2:619685/150‑1 (MQ=255)                                                                                                                                        tCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCCGTGCCAACGCTTTatcatcat   >  1:358561/1‑151 (MQ=255)                                                                                                                                           aaaaGAAGAAGCCGAGGAGATAAGAGCAAAATTAGTCCGCTAATCAGCCGTTTTATTTATTTTTTAAATGGGCGGTTGAATTGTTGTCTCAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatcatcatc  <  1:257920/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                 CCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTG‑A‑TCCTTTATTGGTCACCT‑A‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATC  >  NC_000913/2151470‑2151752

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