Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,386,862	0	G	A	44.6%	60.5 / 39.6	57	intergenic (+31/‑72)	yfbL/yfbM	putative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base G (3/28);  new base A (25/0);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.18e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACC  >  NC_000913/2386714‑2387007                                                                                                                                                     |                                                                                                                                                    cGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCg                                                                                                                                                   <  1:229331/151‑1 (MQ=255)        tATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑c                                                                                                                                            <  1:433250/150‑1 (MQ=255)             taataaACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAG‑CACaatt                                                                                                                                      >  1:125963/1‑151 (MQ=255)                  aaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGAAAAG‑CACaattttttt                                                                                                                                 >  2:583736/1‑151 (MQ=255)                    cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGAAATG‑CACaatttttttt                                                                                                                                >  2:226848/1‑150 (MQ=255)                     aaTACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTcc                                                                                                                              <  1:869781/151‑1 (MQ=255)                            cTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAG‑CACATTTTTTTTTTCCTCATCt                                                                                                                       >  2:11018/1‑151 (MQ=255)                              tGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTg                                                                                                                      <  1:91648/150‑1 (MQ=255)                                     tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACAAG‑CACATTTTTTTTTTCCTCATCTGtttt                                                                                                                  >  2:566155/1‑147 (MQ=255)                                       aTACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATg                                                                                                            <  2:365369/151‑1 (MQ=255)                                           cGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATG‑CACATTTTATTTTTCCTCATCTGTTTTTAGATGCg                                                                                                          >  1:159726/1‑149 (MQ=255)                                               gacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAG‑CACATTTTTTTTTTCCTCATCTGTTTTTAGATGCGaaaaa                                                                                                     >  1:629968/1‑148 (MQ=255)                                                      tGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAAtt                                                                                              <  2:767585/150‑1 (MQ=255)                                                       gAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATG‑CACATTTTTTTTTTCCCCCTCTTTTTTTAAATGGGGAAAATGAAAtt                                                                                              >  1:265139/1‑149 (MQ=255)                                                       gAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACATG‑CACATTTTTTTTTTCCTCATCTGTTTTTAAATGCGAAAAAAGaaa                                                                                                >  1:303080/1‑147 (MQ=255)                                                       gAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑a                                                                                           <  1:68712/151‑1 (MQ=255)                                                         aTTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AAc                                                                                         <  2:246849/151‑1 (MQ=255)                                                          ttATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAA‑CGACAAG‑CACATTTTTTTTTTCCTCATCTGTTTTTAGATGCGAAACATGAAAATC‑aa                                                                                          >  1:452203/1‑148 (MQ=255)                                                              cAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACggag                                                                                     <  2:913002/150‑1 (MQ=255)                                                                gAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACATG‑CACATTTTATTTTTCCTCATCTGTTTTTAGATGCGGAAAATGAAATTT‑AAaggagg                                                                                    >  1:35506/1‑149 (MQ=255)                                                                   aaTGGCTCAGGTAGTGGATGGTGTTATGACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGaaaaa                                                                               <  1:488934/151‑1 (MQ=255)                                                                   aaTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGaaaaa                                                                               <  1:444520/151‑1 (MQ=255)                                                                     tGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACATG‑CACAATTTTTTTTTCCTCCTCTGTTTTTAGATGCGAAAAATGAAATTT‑TACGGAGGaaaaa                                                                               >  1:638272/1‑149 (MQ=255)                                                                      ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAaaata                                                                             <  2:663584/150‑1 (MQ=255)                                                                        cTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATG‑CACATTTTTTTTTTCCTCCTCTGTTTTTAAATGCGAAAAAAGAAAATT‑AAACGAGGGaaaaaaaa                                                                           >  1:402753/1‑150 (MQ=255)                                                                          cAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAAAATAAATgg                                                                        <  2:28022/151‑1 (MQ=255)                                                                              tAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATG‑CACATTTTTTTTTTCCCCCTCTGTTTTTAAATGCGAAAAATGAAATTT‑AACGGAGGAAAAATAAATgggga                                                                     >  1:885392/1‑148 (MQ=255)                                                                               aGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCTAAAAAACGACAAG‑GAAAATTTATTTTTCCCCCTCTTGTTTAAGATTAGGAACATCAAAATCTAAAGGGGGAAAAAAAAAggg                                                                        >  2:73726/1‑147 (MQ=255)                                                                               aGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAG‑AACATTTTTTTTTTTCTCCTCTTTTTTTAAATGCGGAAAAAGAAAATT‑CACAGGGGAAAAAAAAATggg                                                                       >  2:762612/1‑147 (MQ=255)                                                                                gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAAAATAAATGGGAATGa                                                                  <  1:696276/151‑1 (MQ=255)                                                                                  ggATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAAAATAAATGGGAATGAtt                                                                <  2:755741/151‑1 (MQ=255)                                                                                       gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTa                                                           <  1:900055/151‑1 (MQ=255)                                                                                        tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCAGAAAAAAAGAAATG‑CAAATTTTATTTTTCCTCCTCTGTTTTAAAAAACGAAAAAAGAAAATT‑AAAGGAGGGAAAATAAATGGGAATGATTGGCTAt                                                          >  2:105497/1‑151 (MQ=255)                                                                                          ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATG‑CACATTTTTTTTTTCCTCCTCTGTTTTTAGATGCGAAACAAGAAATTT‑AACGGAGGGAAAAAAAATGGGGATGATTGGGTAttt                                                        >  1:898321/1‑151 (MQ=255)                                                                                           tataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAAAAAAAACGACAAG‑AACAATTTTTTTTTTCTCATCTGTTTTTAAAAGGGAAAAATGAAAATC‑AACGGAGGGAAAAAAAATGGGGATTATTGGCTAtttt                                                       >  2:926108/1‑151 (MQ=255)                                                                                            ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTAtttt                                                       <  2:69180/150‑1 (MQ=255)                                                                                                  ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGa                                                  <  1:355070/149‑1 (MQ=255)                                                                                                   tGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGaaa                                                <  1:336219/150‑1 (MQ=255)                                                                                                   tGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGaaa                                                <  2:402616/150‑1 (MQ=255)                                                                                                   tGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGaaa                                                <  1:384141/150‑1 (MQ=255)                                                                                                    gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAG‑CACATTTTTTTTTTCCTCCTCTTTTTTTAGAAGCGAAAAAAGAAATTC‑AAAGGAGGAAAAAAAAATTGGAAAGATTGGCTATTTTTTAGaaa                                                >  2:809935/1‑149 (MQ=255)                                                                                                               gTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACATG‑CACATTTTTTTTTTCCCCATCTGTTTTAAAAAGGGGAACATGAAATTT‑AAAGGGGGAAAAAAAAATGGGAATGATTTGCTATTTTTTAGAAAATTAATCAGaaa                                    >  1:373965/1‑150 (MQ=255)                                                                                                                   ataataAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATG‑CACATTTTTTTTTTCCTCATCTGTTTTTAGATTCGAAAAATGAAATTT‑AACGGGGGAAAAATAAATGGGGATGATTGGCTATTTTTCAGAAAATTACTCAGAAAAGa                                 >  1:113830/1‑149 (MQ=255)                                                                                                                     aataaAAATTTCTCGAGGCGGGTGCAGAAAAAACGAAAAG‑CACAATTTTTTTTTCCTCATCTGTTTTTAGATGCGGAAAATGAAAATT‑AACGGAGGGAAAATAAATGGGGATGATTGGGTATTTTTTAGAAATTGACTCAGaaaaaaat                               >  2:15994/1‑145 (MQ=255)                                                                                                                     aataaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTa                              <  2:125952/150‑1 (MQ=255)                                                                                                                     aataaAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCCAAATTTTTTTTTCCC‑CCTCTGTTTTTAGATGCGAAAAAAGAAAATT‑AACAGGGGGAAAAAAAATGGGGATGATTGGGTATTTTTTAGAAAATGAAACAGaaaaaaata                              >  2:460336/1‑145 (MQ=255)                                                                                                                         aaaaTTTCTCGAGGCGGGTCCAGAAAAAACGACAAG‑CACATTTTTTTTTTCCTTTTCTGTTTTTAGATGCGAAAAATGAAAATT‑AAAGGGGGAAAAAAAAATGGGAAAGATTTGCTATTTTTTAGAAAATTACTCAGAAAAAATTaa                             >  2:306614/1‑147 (MQ=255)                                                                                                                           aaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACt                       <  2:662986/151‑1 (MQ=255)                                                                                                                            aTTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTg                      <  2:428082/151‑1 (MQ=255)                                                                                                                                  cGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTg                  <  2:503941/149‑1 (MQ=255)                                                                                                                                  cGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAg                <  2:208450/151‑1 (MQ=255)                                                                                                                                   gAGGCGGGTCCCGAAAAAACGACATG‑CACATTTTATTTTTCCTCATCTGTTTTAAGATGCGAAAAATGAAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTTCAGAAAATTACTCAGAAAAGATTAAACAACTGCtt                   >  2:45006/1‑147 (MQ=255)                                                                                                                                         ggTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTg         >  2:334027/1‑151 (MQ=255)                                                                                                                                         ggTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTg         <  1:583792/151‑1 (MQ=255)                                                                                                                                          gTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGa        <  2:4828/151‑1 (MQ=255)                                                                                                                                              aGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaaa     >  2:516515/1‑150 (MQ=255)                                                                                                                                                aaaaaGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAAcc  >  2:380902/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATG‑CACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACC  >  NC_000913/2386714‑2387007

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