Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,671,415	0	C	T	44.2%	48.2 / 37.5	52	K160K (AAG→AAA)	yhjB	putative DNA‑binding transcriptional response regulator
Reads supporting (aligned to +/- strand):  ref base C (27/2);  new base T (0/23);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 8.50e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCTGCTGGCTGTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTT‑CCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACTATTTACCGCCAACG  >  NC_000913/3671265‑3671563                                                                                                                                                         |                                                                                                                                                     ttCCTGCTGGCTGTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTT‑CCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTcc                                                                                                                                                      >  1:171774/1‑151 (MQ=255) ttCCTGCTGGCTGTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTT‑CCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTcc                                                                                                                                                      >  1:5898/1‑151 (MQ=255)            tgTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTT‑CCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTc                                                                                                                                           >  2:20238/1‑151 (MQ=255)            tgTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTT‑CCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTc                                                                                                                                           >  2:246101/1‑151 (MQ=255)                       ccTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTT‑CCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCt                                                                                                                                 >  2:45355/1‑150 (MQ=255)                       ccTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTT‑CCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCt                                                                                                                                 >  2:729612/1‑150 (MQ=255)                         tctTCAGGAGGAGGTATTTAACATCATTGCCGACTGGAAGAGATTTTTCACTT‑CCAGACGGCGAGACAGGGATTCCAGATGCGCGTTTACAGTGCCTT‑TACTGATAGTAAGCACTCTGCCGATCTCCTTATTTGATTCACCcgcagc                                                                                                                                  >  2:670951/1‑147 (MQ=255)                         tctTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTT‑CCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAAc                                                                                                                              >  1:298302/1‑151 (MQ=255)                         tctTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTT‑CCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAAc                                                                                                                              >  1:423352/1‑151 (MQ=255)                            tCAGGAGGAGATATTTAACATGATTACCGCCTGGGTGCGATTTTTCACTT‑CCAGACGGCGATACAGGGATTCCAGATACGCTTTTACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATg                                                                                                                           >  2:723022/1‑151 (MQ=255)                               ggaggaGATTTTTAACATCATTGCCGCGTGGGTGCTTTTTTTCCCTC‑CCAGACGGCAATACAGGGTTTCCAGATGCGCTTTTCCCGTTCCGT‑TACTGATATTCAGCGCTCTGCCGTTCTCTTTATTTGATTCGCCCGCCGCTAACATg                                                                                                                           <  2:768675/148‑1 (MQ=255)                               ggaggaGATTTTAAACTTCTTTGCCGCCGGGGTGCTTTTTTTCCCTC‑CCAGAGGGCAATAAGGGGATCCCAGTTCCTTTTTTCCCGTCCCGG‑TCCTGATATTCAGCGTTCTGCCGTTCTCTTTTTTTGTTTCGCCCGCCGCTAACATGGt                                                                                                                         <  2:182238/150‑1 (MQ=255)                                       tatTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTT‑CCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAAtctc                                                                                                                >  2:970791/1‑151 (MQ=255)                                       tatTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTT‑CCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAAtctc                                                                                                                >  1:92831/1‑151 (MQ=255)                                            accatcatTGCCGCCGGGGTGCGTTTTTTCACTT‑CCAGAGGGCGATACAGGGATTCCAGATGCGTTTTTCCCGTTCCGT‑TTTTTATATGCAGCGCTCGGCATACCTCTTTTTTTGATGCGCCCGCCGCTAACATGGTTAAAATCTCCCGc                                                                                                            <  1:409808/148‑1 (MQ=255)                                            aatatcatGGCCGCGGGGGTGTGTTTTTTACCTT‑CCAGAGGGCGATACAGGGATTCGAGTTGCTTTTTTCCCTTTCCGT‑TTCTGATATTCCGCTCTCGGCCGATCTCTTTATTTGATCCGCCCGCCGCTAACATGGTTAAAATCTCCCg                                                                                                             <  2:262751/146‑1 (MQ=255)                                                  agtGCCGCCGGGGTGCTTTTTTTCCCTC‑CCAGGCGGCAATAGGGGGTTTCCAGATGC‑CTTTTTCCCTTCCGGTTACTAATATCCAGCGCTCTGCCTTTCTCTTTATTTGATCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTg                                                                                                          <  1:537875/144‑1 (MQ=255)                                                   ttGCCGCCTGGGTGCGATTTTTCACTT‑CCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTggcgggc                                                                                                    >  2:420789/1‑151 (MQ=255)                                                       cgctggggTGCATTTTTACACTC‑CCAGAGGGTGAAACAGGAATCCCAGATGCTTTTTTCCCCTTCCGG‑TCCTAATATTCTCCCCTCTGCCGATCTCTTTTTTTGTTTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTggcgggc                                                                                                    <  1:306767/142‑1 (MQ=255)                                                       cgccTGGGTGCGATTTTTCACTT‑CCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCtt                                                                                                >  1:573839/1‑151 (MQ=255)                                                            gggTGCTTTTTTTCCCCC‑CCAGACGGCGATACAGGGATTCCAGATGCTTTTTTCCCGTTCCGT‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTaa                                                                                              <  1:297481/148‑1 (MQ=255)                                                             ggTGCTTTTTTTCTCTC‑CCAGCCGGCAATACAGGGATCCCAGTTCCGTTTTTCCCGTCCCGT‑TCCTAATATTCGGCGCTCGGCCGATCCCTTTTTTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGAt                                                                                          <  1:549613/151‑1 (MQ=255)                                                             ggTGCTTTTTTCCACTC‑CCAGCGGGCGACACGGGTATCCCAGTGGCTTTTTTCCCTTCCCGG‑TCCTGTTATCCAGCGCTCTGCCGTCCTCTTTTTTTGATCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACg                                                                                            <  2:1061741/149‑1 (MQ=255)                                                                 cGATTTTTCACTT‑CCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTg                                                                                       >  1:35867/1‑150 (MQ=255)                                                                  cttttttCCCTTCCCAGACGGCAAAACGGGGATTCCAGATGCTCTTTCCCCTTCCCGG‑TATTGATATCCGGCGCTCTGCCGTTCTCTTTTTTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTg                                                                                       <  1:300678/148‑1 (MQ=255)                                                                       tccACTC‑CCAGAGGGCGATACAGGGATTCCAGTTGCGTTTTTCCCGTCCCGG‑TCCTTATATTCAGCGTTCTGCCTTTCTCTTTTTTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCt                                                                                 <  2:161718/148‑1 (MQ=255)                                                                         cgttt‑CCAGCCGTCGACGCGGGGTTCCCAGATGCGCTTTTACCGTTCCGT‑TTCTGATATCCAGCGCTCTCCCGCTCCCTTTTTTTGATCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTa                                                                              <  2:436614/148‑1 (MQ=255)                                                                                  gcgggCGATACAGGGATTCCATATGCTTTTTTCCCGTCCCGT‑TCCTGTTTTCCACCGCTCTGCCTTTCTCTTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATg                                                                           <  2:36373/143‑1 (MQ=255)                                                                                    cGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCtttt                                                                    >  1:267109/1‑151 (MQ=255)                                                                                     gggaaTACAGGGATCCCAGATGCGTTTTTCCGGTCCCGT‑TACTGATTTCCGGCGCTCTGCCGTTTTCTTTTTTTGATCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATg                                                                           <  2:506352/139‑1 (MQ=255)                                                                                      gCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTcc                                                                  >  1:899414/1‑151 (MQ=255)                                                                                         aTACAGGGATTCCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGc                                                               >  2:97470/1‑151 (MQ=255)                                                                                          tagaGGTATCCCAGATGCGTTTTACCCGTTCCGG‑TATTGATTTTCAGCGTTCCGCCTTTTTTTTTATTTGTTTCCCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCg                                                                 <  1:591858/145‑1 (MQ=255)                                                                                              ggtaTTCAAGATGCTTTTTTCCCTTTCGGT‑TATTTATTTTCAGCCCTCGCCCGTTCTCTTTTTTTTTTTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGcgt                                                             <  2:376612/145‑1 (MQ=255)                                                                                              gggATTCCAGTGGCTTTTTTCCCGTCCGGG‑TACTGATATCCGGCGCTCTGCCGTTCTCTTTATTTGTTCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGcgtcgt                                                          <  2:1048616/151‑1 (MQ=255)                                                                                                  tcccAGATGCGTTTTTCCCGTCCGGG‑TCCTGATATCCAGCGCTCTGCCGATCTCTTTTTTTGATCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTcc                                                        <  1:308909/147‑1 (MQ=255)                                                                                                      aGATGCGCTTTTACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGt                                                  >  2:964166/1‑151 (MQ=255)                                                                                                             tttttTCCCTTCCGGTTTCTGATATACGGCGCTCCCCCGTCTTTTTTTTTTGTTCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGtct                                                <  2:46424/146‑1 (MQ=255)                                                                                                             cTTTTACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGAttt                                                                                        >  2:1066511/1‑106 (MQ=255)                                                                                                              ttttACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGc                                          >  1:1004676/1‑151 (MQ=255)                                                                                                                ttccctttCGGG‑TATTGATTTTCATCGCTCGCCCCATCTTTTTTTTGTTTTCCCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCa                                             <  2:275631/140‑1 (MQ=255)                                                                                                                       ttccgtttatGATATGCAGCGTTGTGCTCTTCTCTTTTTTTAATTCCCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGaaa                                     <  1:600060/139‑1 (MQ=255)                                                                                                                          gg‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACatcatc                              >  1:940274/1‑151 (MQ=255)                                                                                                                             tACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACTTGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTAAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATccc                            >  1:749540/1‑151 (MQ=255)                                                                                                                              aCTGATTTTCAGCGCTCCGCCGTTCTCTTTTTTTGTTTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATccc                            <  1:997230/150‑1 (MQ=255)                                                                                                                                    atccaccgcTCTCCCGTTCTCTTTTTTTATTCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCa                          <  1:541328/140‑1 (MQ=255)                                                                                                                                     ttCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATcgcc                     >  2:171361/1‑150 (MQ=255)                                                                                                                                            ctctGCCGTTCCTCTTATTTTATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATc                        <  1:230198/140‑1 (MQ=255)                                                                                                                                               tGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACTATTTAc          >  2:32121/1‑151 (MQ=255)                                                                                                                                               tGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACTATTTAc          >  2:955696/1‑151 (MQ=255)                                                                                                                                                 ccattctcTTTTTTTGTTTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACTATTTAc          <  1:426003/145‑1 (MQ=255)                                                                                                                                                   gATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACTATTTACCGcc      >  2:437429/1‑151 (MQ=255)                                                                                                                                                      ctcCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACTATTTACCGCCAAc   >  1:199873/1‑151 (MQ=255)                                                                                                                                                       tcCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACTATTTACCGCCAACg  >  2:322370/1‑151 (MQ=255)                                                                                                                                                         cTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACTATTTAcccc       >  1:81941/1‑142 (MQ=255)                                                                                                                                                         |                                                                                                                                                     TTCCTGCTGGCTGTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTT‑CCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGG‑TACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACTATTTACCGCCAACG  >  NC_000913/3671265‑3671563

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