Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,147,733	0	T	G	43.4%	83.3 / 12.0	54	G20G (GGT→GGG)	yegD	Hsp70 chaperone family protein
Reads supporting (aligned to +/- strand):  ref base T (2/28);  new base G (23/0);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 4.81e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTAC  >  NC_000913/2147583‑2147875                                                                                                                                                         |                                                                                                                                               gACGCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGt                                                                                                                                                <  1:14731/151‑1 (MQ=255)  aCGCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTa                                                                                                                                               <  1:902665/151‑1 (MQ=255)   cGCTGCGCTCTTCTTTCATATTCCGGTTGTCGCGACTGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCTTGTTTATTGGTTTTGATTACGGTACAGCAATCTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTaa                                                                                                                                              <  2:714578/151‑1 (MQ=255)      tGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GGGGG‑GCGGGaaa                                                                                                                                             >  1:3967/1‑149 (MQ=255)        cgcTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGc                                                                                                                                         <  1:702290/151‑1 (MQ=255)        cgcTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGGTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGc                                                                                                                                         <  2:883881/151‑1 (MQ=255)              gctttCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTg                                                                                                                                    <  1:885506/150‑1 (MQ=255)               ctttCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GTGTG‑AGGGGAAAACGCCttt                                                                                                                                     >  2:524884/1‑148 (MQ=255)                  tCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTaa                                                                                                                                <  1:425060/150‑1 (MQ=255)                    atatTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGGAAAACGCATTTGCTtaaa                                                                                                                              >  1:256916/1‑150 (MQ=255)                       ttCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGGAAAACGCCCTTGCTaaaaa                                                                                                                             >  2:175741/1‑148 (MQ=255)                        tCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATgg                                                                                                                          <  2:578155/150‑1 (MQ=255)                          cGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTTAT‑GCGTG‑ACGGGAAAACCCAATTTCTAAAAAtagga                                                                                                                        >  2:643382/1‑146 (MQ=255)                           ggTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑GCGGGGAACCCCCTTTTCTTAAAAAGGaa                                                                                                                        >  2:797773/1‑149 (MQ=255)                            gTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAAc                                                                                                                     <  2:79186/151‑1 (MQ=255)                                       ggCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACgct                                                                                                          <  2:850470/151‑1 (MQ=255)                                              tttCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGGG‑ACGGGGAACCGCATTTGCTTAAAATGGAAAAAGACAGCACGCTGCTGcc                                                                                                    >  2:359395/1‑150 (MQ=255)                                                     aaaGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATg                                                                                             <  2:428920/150‑1 (MQ=255)                                                          cACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTGCT‑GCGGG‑ACGGGAAACCGCCCTTTCTCAAAAAGGGAAAAGAAAGCACGCCGCCGCCTTCAATGCttt                                                                                         >  2:391937/1‑149 (MQ=255)                                                            cGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTgcgc                                                                                     <  1:459764/151‑1 (MQ=255)                                                                  ttcttcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGGGGTG‑A‑GGGAAAACGCATTTTCTTAAAAAGGGAAACGACAGCACGCTGCTGCCCTCAATGCTTTTCGCGcc                                                                                  >  2:383393/1‑148 (MQ=255)                                                                   tcttcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAAc                                                                               <  2:276940/150‑1 (MQ=255)                                                                   tcttcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGGAAAACGCATTTTCTAAAAAAGGAAAAAGACAACACGCTGCTGCCTTCAATGCTTTTCGCGcc                                                                                  >  1:55190/1‑147 (MQ=255)                                                                   tcttcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGACAT‑GCGTT‑TCGGGAAAACACATTTTTTTAAAAAGGAAAAACGCAGCACGCTGCTTCCTTTCATGCttttgggccc                                                                                  >  2:459434/1‑140 (MQ=255)                                                                     ttcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑GCGGGAAAACGCATTTGCTAAAAAATGGAAACGGAAGCACGCTGCCGGCTTCCACTCTTTGGGCGcc                                                                                  >  2:870622/1‑145 (MQ=255)                                                                      tcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAAcg                                                                              <  2:561277/148‑1 (MQ=255)                                                                            gCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCg                                                                     <  1:848311/151‑1 (MQ=255)                                                                            gCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCg                                                                     <  2:721519/151‑1 (MQ=255)                                                                               tGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGGAAAACGCCTTTTCTAAAAAAGGAAAAAGACAGCACGCTGCTGCCTTCAATGCTTTTCGCGCCAACGCGGGAAGCgg                                                                    >  1:384113/1‑149 (MQ=255)                                                                               tGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGGAAAACGCATTTTCTAAAAAAGGAAAAAGACAGCACGCGGCTGCCTTCAATGCTTTTCGCGCCAACGCGGGAAGCgg                                                                    >  1:272951/1‑149 (MQ=255)                                                                                ggggATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGGAAAACGCCTTTTCTAAAAAAAGGAAACGACAGCCCGCCGCTGCCTTCCATGCTTTTCGCGCCCACGCGGGAAACgg                                                                    >  2:767754/1‑148 (MQ=255)                                                                                   gATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑GCGGGAAAACCCATTTTCTAAAAATGGAAAAAGACAGCACGCTGCCGCCTTCCATGCTTTTCGCGCCAACGcgggaagcgggaag                                                                >  2:388707/1‑149 (MQ=255)                                                                                    attattTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGGAAAACGCCTTTGCTAAAAAAGGAAAAAGACAGCACGCCGCTGCCTTCAATGCTTTGCGCGCCAACGCGCGAAGCGGGaa                                                                 >  1:331863/1‑147 (MQ=255)                                                                                       attTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATg                                                          <  2:639617/151‑1 (MQ=255)                                                                                        tttCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGGG‑ACCGGAAAACGCATTTGCTTAAAATGGAAAAAGACAGCACGCTGCCGCCTTCAATGGTTTTTGCGCCCACGCGTGAAACGGGAAGCgga                                                            >  1:671032/1‑146 (MQ=255)                                                                                                  aTTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGcc                                               <  1:524995/151‑1 (MQ=255)                                                                                                        tttGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCcatcat                                          <  1:52776/150‑1 (MQ=255)                                                                                                         ttGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAACCCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCatg                                         <  2:2979/150‑1 (MQ=255)                                                                                                         ttGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCatg                                         <  1:72539/150‑1 (MQ=255)                                                                                                            aTTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGGG‑ACGGGAAAACGCATTTTCTAAAAAAGGGAAAAAAAAACACGCTGCTGCCTTCAAAGCTTTTCGCGCCAAAGCGGGAAGCGGGAAGGGAAAGGGTGTACCGCCAACatgaag                                      >  1:253037/1‑148 (MQ=255)                                                                                                            aTTACGGTACAGCAAACTGTTCAGTGGCGGCCAT‑GGTTGCATGGAAAAAC‑CACTTATTAAAAAAACAAAACGACAGCAACCGGCTGGCCTCCATGCTTTTCGCGCCAAAACGTTACGCGGAAAGTGGCTGGCTTTATCCCCCTCatgc                                        >  2:555858/1‑147 (MQ=255)                                                                                                                cGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCg                                 <  2:321243/151‑1 (MQ=255)                                                                                                                 ggTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCgg                                <  2:465561/151‑1 (MQ=255)                                                                                                                    aCAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGGAAAACGCATTTTCTAAAAAAGGAAAACGACAGCACGCCGCTGCCTTCAATGCTTTTCGCGCCAACGCGTGAAGCGGGAAGCGGATGGCTGTACCCCCCTCATGATGTTTCgg                                >  1:371598/1‑148 (MQ=255)                                                                                                                      aGCAAACTGTTCAGTGGCGGTCCT‑GCGGG‑AAGGGAAAACGCACTTACTAAAAATGGAAAAAGACAGCACGCTGCTGCCTTCAATGGTTTTCCCGCCAACGCGTGAAGAGGGAAACGAAAGGCCGTACCCCCATCATTGTGtttc                                  >  2:702906/1‑142 (MQ=255)                                                                                                                       gCAAACTGTTCAGTGGCGGTCAT‑GCGGG‑ACGGGAAACCGCATTTGCTAAAAATGGAAAACGACAACACGCCGCTGCCTTCAATGCTTTTCGCGCCAACGCGGGAAGCGGGAAGCGAATGGCTGTACCGCCCTCATGATGTTCCGGCAgac                           >  1:641554/1‑150 (MQ=255)                                                                                                                        cAAACTGTTCAGTCGCGGTCAT‑GCGTG‑ACGGTAAACCGCATTGGCTAAAAAGTGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATTGCTGTACCTCCATCAGGATGTTCCGGCAgacga                         <  2:417877/151‑1 (MQ=255)                                                                                                                           aCTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGAcgacga                      <  2:288282/151‑1 (MQ=255)                                                                                                                            cTGTTCAGTGGCGGTCAT‑GCGTG‑GCGGGAAACCCCCTTTTCTAAAAAAGGGAAAAGACAGCACGCCGCTGCCTTTAAATCTTTGCGCGCCAAAGCGGGAAGCGGGTAGCGAATGGCTGTACCCCCCTCCTGTTGTTCCGGGAg                             >  2:143999/1‑143 (MQ=255)                                                                                                                              gTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAAACGCATTTGCTAAAAATGGGAAACGACAGCACGCTGCTGCCTTCAAAGCTTTGCGCGCCAACGCGGGAAGCGGGAAGCGAATGGCTGTACCGCCATCCTGATGTTTCGGCAgacga                         >  1:115395/1‑145 (MQ=255)                                                                                                                                     ggCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGc              >  1:138160/1‑149 (MQ=255)                                                                                                                                     ggCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAg            <  1:609544/151‑1 (MQ=255)                                                                                                                                     ggCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAg            <  1:412996/151‑1 (MQ=255)                                                                                                                                        ggTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGgcg         <  2:519459/151‑1 (MQ=255)                                                                                                                                                cGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTAc  <  2:3967/151‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                               GACGCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAT‑GCGTG‑ACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTAC  >  NC_000913/2147583‑2147875

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