Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,157,271	0	T	G	44.4%	88.8 / 11.0	54	G165G (GGT→GGG)	ycfH	putative DNase
Reads supporting (aligned to +/- strand):  ref base T (5/25);  new base G (24/0);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 9.09e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCG  >  NC_000913/1157123‑1157420                                                                                                                                                     |                                                                                                                                                      atatCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTaa                                                                                                                                                     <  2:302278/151‑1 (MQ=255)   atCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGaaa                                                                                                                                                    >  2:620975/1‑150 (MQ=255)        gATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACt                                                                                                                                               >  2:1015400/1‑150 (MQ=255)          ttGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGaaa                                                                                                                                                    >  2:414093/1‑143 (MQ=255)              tCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTGGAtct                                                                                                                                         >  2:593258/1‑150 (MQ=255)                  gaactgaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCAGGTAAATTACTGGATCTCGGa                                                                                                                                     <  1:412224/150‑1 (MQ=255)                      tgaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGTCAATACTGGGTCTCCGGtttt                                                                                                                                 >  2:97827/1‑150 (MQ=255)                         acAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTAATGGGTCTCCGGtttt                                                                                                                                 >  1:672316/1‑147 (MQ=255)                           aaGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACAtct                                                                                                                           <  1:558346/151‑1 (MQ=255)                            aGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACAtctc                                                                                                                          <  1:389359/151‑1 (MQ=255)                             gCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTcc                                                                                                                         <  2:587003/151‑1 (MQ=255)                                ggTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCttt                                                                                                                      <  1:728308/151‑1 (MQ=255)                                    aTCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGAGCTTGGGTTTTACAACTCCtttt                                                                                                                     >  1:442850/1‑148 (MQ=255)                                        tCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTAATGGGATTCGGGTTTTACATCTCCCTTTTCCGCAt                                                                                                              >  2:22602/1‑151 (MQ=255)                                         ccATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCAtt                                                                                                             <  2:768140/151‑1 (MQ=255)                                           aTACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATtgt                                                                                                           <  2:641739/151‑1 (MQ=255)                                             aCCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGGTCTCGGGTTTTACATCTCCTTTTTCGGCATtgt                                                                                                           >  2:917684/1‑149 (MQ=255)                                               ccGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGATCTCGGGTTTTTTATCTCCtttttccgcctttttgcc                                                                                                       >  1:980954/1‑137 (MQ=255)                                                cGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCt                                                                                                      <  1:1015233/151‑1 (MQ=255)                                                      gCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATTTTGGATCTCGGGTTTTAAAACTCCTTTTTCCGCATTGTGGCCtttc                                                                                                   >  2:332519/1‑146 (MQ=255)                                                       cccGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTa                                                                                                <  1:352526/150‑1 (MQ=255)                                                                 tACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGATCTCGGGTTTTATATCTCCTTTTTCCGCCTTGTGACCTTCCGTAATGCGGAGc                                                                                       >  1:984556/1‑149 (MQ=255)                                                                       ggCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTgcg                                                                                <  1:110020/150‑1 (MQ=255)                                                                          aaTTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGAt                                                                            <  1:610399/151‑1 (MQ=255)                                                                            ttCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGGAATTACTTGGTTTCGGGTTTTTTATATCCTTTTTCGGCATTTTGACCCTCCGGAATGCGGAGCAAATGCGCGGtg                                                                           >  2:712296/1‑150 (MQ=255)                                                                                gcgcGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGAtgctg                                                                        <  2:224027/149‑1 (MQ=255)                                                                                    gaagaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGATCTTGGGTTTTACATATCCCTTTTCGGCATTGTGACCTTCCGTTATGCGGAGCAACCGCGCGATGCTGcgcggt                                                                  >  2:111743/1‑149 (MQ=255)                                                                                     aagaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGAGCTCGGATTTTACATCTCCTTTTTCCGCATTGTGACCTTCCCTAAAGCGGAGCAACTGCGCGATGCTgcg                                                                      >  2:153349/1‑146 (MQ=255)                                                                                     aagaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTTGGTCTCGGGTTTTAAAACCCCTTTTTCGGCATTTTGGGCTTCCGTAATTCGGGGGAACTGCGCGGTGCTGcgcg                                                                    >  1:208555/1‑148 (MQ=255)                                                                                     aagaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATTGGTCTCGGGTTTTATATCTCCTTTTTCCGCCTTGTGACCTTCCGTAATTCGGGGGAACTGCGCGATGCTGCGCGtt                                                                  >  1:587664/1‑150 (MQ=255)                                                                                      agaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTAt                                                                <  1:168725/151‑1 (MQ=255)                                                                                          aaaGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTcc                                                            <  2:593474/151‑1 (MQ=255)                                                                                          aaaGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGAATTTAAATCTCCTTTTTCCGCATTGGGACCTTCCGTTATGCGGGGGAACATCGCGGTGCTGCGCGGTAttt                                                              >  2:444104/1‑147 (MQ=255)                                                                                             gTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTGAAATAATCGGTTTCGGGTTTTAAAACTCCTTTTTCCGCCTTGTGAACCTTCGTAATGCGGGGGAACTGCGGGGTGCcgcggggttttttcccc                                                         >  2:845503/1‑138 (MQ=255)                                                                                             gTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTccccc                                                         <  2:710663/151‑1 (MQ=255)                                                                                             gTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGGTCTCGGGTTTTAAATCTCCTTTTTCCGCATTGTGACCTTCCCTAATGCGGGGCAACTGCGCGCTGCTGCGCGTTTTGTtcccc                                                         >  2:94586/1‑151 (MQ=255)                                                                                               gACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGGAATTACTGGGACTCGGGTTTTACATCTCCTTTTTCGGCACTGTGGCCTTCCGTAATTCGGGGCAACTTCGCGGTGCTGCGCGTTTTGTtcccc                                                         >  1:373336/1‑149 (MQ=255)                                                                                                     ttGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTAATGGATCTCGGATTTTAAATCTCCTTTTTCGGGATTGTGAACTTCCGTAATGCGGAGCAACTGCGCGGTGCTTCGCGGTATTTCCCCCTGGAGCgg                                                 >  1:792681/1‑151 (MQ=255)                                                                                                        cGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGGTCTCGGATTTTTAATCTCCTTTTCCGGCCTTGTGACCTTCCGTAATGCGGAGGAACTGCGCGATGCTGCGCGGTTTGTCCCCCCGGGTCgg                                                 >  1:205370/1‑148 (MQ=255)                                                                                                         ggTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGATCTCGGGTTTTACATCTCCTTTTTCGGCATTGTGACCTTCCGTAATGCGGGGGAACTGCGCGATGCTGCGCGGTATGTCCCCCCGGATCgggt                                               >  2:788512/1‑147 (MQ=255)                                                                                                            ggCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACtg                                           <  1:712175/150‑1 (MQ=255)                                                                                                                   tacACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTTCTGGATCTCGGATTTTACATCTCCTTTTTCCGCATTTTGACCTTCCGTAATGCGGGGGAACTGCGCGATGCTGCGCGGTATGTCCCCCTTGGTCGGGTTCTGGGGGgaa                                    >  1:90396/1‑150 (MQ=255)                                                                                                                       cTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGGTCTCGGGTTTTACATCTTCTTTTCCCGCCTTGGGGCCCTCCGTAATGAGGGGGAACTGCGCGGTGCTGCGCGGTTTGTTCCCCTGGCTGGGGTTCTTGTGGgaa                                    >  1:975698/1‑146 (MQ=255)                                                                                                                               cAGAGGACAGAGAAACGGCGGGGAAAATACTGGATCTCGGATTTTATATCTCCTTTTTCGGCATTGTGACCTTTCGTAATGCGGAGCAACTGCGCGCTGCTGCGCGTTATGTCCCCCCGGGTCGGGTACTGGGGGGAACTGACTccccct                        >  1:86099/1‑148 (MQ=255)                                                                                                                                agagGACAGAGAAACGGCGGGGAAAATACTGGGTTTCGGATTTTAAATCTCCTTTTCCGGCATTGTGACCTTCCGTAATTGGGGGCAAAATCGCGGTGCTGCCCGGTATTTTCCCCCGGGTCGGGTTCTTGGGGGAACTGACTccccct                        >  1:767809/1‑147 (MQ=255)                                                                                                                                    gACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTg                  <  1:388799/151‑1 (MQ=255)                                                                                                                                    gACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTg                  <  1:917529/151‑1 (MQ=255)                                                                                                                                    gACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTg                  <  2:916225/151‑1 (MQ=255)                                                                                                                                    gACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTg                  <  2:235899/151‑1 (MQ=255)                                                                                                                                           aaaCGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCgg            <  1:750059/150‑1 (MQ=255)                                                                                                                                            aaCGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCTTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTa          <  1:488908/151‑1 (MQ=255)                                                                                                                                            aaCGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTa          <  1:1059551/151‑1 (MQ=255)                                                                                                                                                 cggGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCa     <  2:775983/151‑1 (MQ=255)                                                                                                                                                  gggTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCAt    >  2:685057/1‑151 (MQ=255)                                                                                                                                                     tAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCg  >  2:736892/1‑150 (MQ=255)                                                                                                                                                     tAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCg  >  2:957459/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                      ATATCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCG  >  NC_000913/1157123‑1157420

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