Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	803,869	0	G	T	40.0%	67.1 / 31.8	55	V123L (GTG→TTG)	ybhJ	aconitase family protein
Reads supporting (aligned to +/- strand):  ref base G (5/28);  new base T (22/0);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 6.21e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAAC  >  NC_000913/803720‑804016                                                                                                                                                      |                                                                                                                                                     aGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTtg                                                                                                                                                      <  2:375974/150‑1 (MQ=255)      gCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTG‑cccc                                                                                                                                               >  1:1065601/1‑149 (MQ=255)         aaaGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTcc                                                                                                                                              <  1:227180/149‑1 (MQ=255)         aaaGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGc                                                                                                                                            <  1:135697/151‑1 (MQ=255)          aaGCGACCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCCGCGGCACTATTAACGGTGATGACCATGATTTTGTTTTATCGGCGGCCCAGCGTTATGTCGGTAATTTtgt                                                                                                                                                     >  2:822447/1‑142 (MQ=255)             cGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTG‑ccctcgcct                                                                                                                                          >  2:78784/1‑142 (MQ=255)                   gTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTGTTTTTTTT‑CCCCCGccaattttcgg                                                                                                                                  >  1:690024/1‑141 (MQ=255)                   gTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTG‑GCCCCGCCTATTGCgg                                                                                                                                   >  1:526159/1‑150 (MQ=255)                   gTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTG‑CCCCCGCCTAtt                                                                                                                                       >  2:865326/1‑146 (MQ=255)                     aTGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTc                                                                                                                                 <  2:136079/150‑1 (MQ=255)                       ggAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTG‑GCCCCCCCTAtttggggcat                                                                                                                               >  1:122138/1‑142 (MQ=255)                        gAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTcatcc                                                                                                                             <  1:557918/151‑1 (MQ=255)                           cGTTTCCCGCTGCCGAATGTGCTGACCAGCTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTcatccat                                                                                                                           <  2:77943/150‑1 (MQ=255)                             tttCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTG‑CCTCCCCCTATTTCgggcgttcct                                                                                                                           >  1:193091/1‑140 (MQ=255)                                 ccGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTCATCCATCAatata                                                                                                                    <  2:840640/151‑1 (MQ=255)                                 ccGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTG‑TCCCCCCCTATTTCgggtacgcctccaaat                                                                                                                     >  2:624613/1‑136 (MQ=255)                                  cGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTT‑CCTCCCCATTTTTGGGGCATCCCTCaa                                                                                                                        >  1:217228/1‑146 (MQ=255)                                      gCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGt                                                                                                               <  2:484886/151‑1 (MQ=255)                                          tATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTG‑CCTCCGCATATTGCGGGCATCCATCAATATATGCGTgag                                                                                                            >  1:543066/1‑150 (MQ=255)                                              tgCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCCGCGTTATGGCGGGAATTTTTTT‑CCTCCGCTTAttttggttcacca                                                                                                                            >  2:925205/1‑119 (MQ=255)                                               gCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAgatgat                                                                                                       <  2:50098/150‑1 (MQ=255)                                                    ccAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAg                                                                                                  <  1:1006231/150‑1 (MQ=255)                                                     cAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTGTTTTTTTGTCCCCCGC‑TATTGGTGTCAACCCCCCATATATTTGGGGGGTGATGGGAggg                                                                                                >  1:200402/1‑150 (MQ=255)                                                     cAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAggc                                                                                                <  2:543047/151‑1 (MQ=255)                                                      aaCTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAg                                                                                                  <  1:236087/148‑1 (MQ=255)                                                      aaCTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTAATTTTTTG‑GCTTCGCCTATTTCGGGCATCCATCAATATTTGTGGGGGATGATGGGAgg                                                                                                 >  1:1130368/1‑149 (MQ=255)                                                      aaCTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGTTTTTTTTG‑TCTCCGCCTATTTCGGGCCTCCCTCAATATATGCGGGGGGAGATGGGAgg                                                                                                 >  2:848531/1‑149 (MQ=255)                                                       aCTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAggcgg                                                                                              <  1:522732/151‑1 (MQ=255)                                                         tGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTG‑GCCCCCCCTATTTCGGGCCTTCATCCATATATTCGTGGGGTGATggggggggg                                                                                              >  2:992414/1‑142 (MQ=255)                                                          gCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTG‑TCCCCCCCTATTGGGGGGCTCCCTCCATATTTTCGGGGGGTGATGGCAggggg                                                                                              >  2:713944/1‑148 (MQ=255)                                                          gCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAggcggcgg                                                                                           <  1:1080468/151‑1 (MQ=255)                                                             aTAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTG‑CCTCCCCCTATTTCGGTCATCCCTCCATATATGCGTGGGATGATGGCAggcggcgg                                                                                           >  1:475864/1‑148 (MQ=255)                                                             aTAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCaa                                                                                        <  1:1087498/151‑1 (MQ=255)                                                               aaCTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCaaaa                                                                                      <  1:629049/151‑1 (MQ=255)                                                                actcactcTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTG‑CCCCCGCCTATTTCGGGGCTCCCTCCATATATTCGTGAGATGAAGGCAGGGGGCGGCcaaa                                                                                      >  1:832249/1‑150 (MQ=255)                                                                actcactcTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAAt                                                                                     <  1:81525/151‑1 (MQ=255)                                                                   cactcTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGa                                                                                   <  1:4128/150‑1 (MQ=255)                                                                     ctctGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTG‑CCTCCGCCTATTGCGGGCATCCATCCAAATATGCGTGGGGTGATGGCAGGCGGGGGCAAAATGATcc                                                                                >  1:937225/1‑151 (MQ=255)                                                                         gcgcCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCgg                                                                            <  1:504874/151‑1 (MQ=255)                                                                             cGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCa                                                                        <  1:536396/151‑1 (MQ=255)                                                                                cggcggCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGGATTTTTTG‑CCGGCCCGCTTTACGGCGGGCCACCACTATATGCGTGGGACGTGGGTGGGGGGGGGGAAAAAAATCCCCGGGTCAGa                                                                      >  2:41597/1‑150 (MQ=255)                                                                                       aCTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTG‑CCTCCGCCTATTTCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGAACCTCCGGGCAGACAGcc                                                                 >  1:76850/1‑148 (MQ=255)                                                                                            tAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCg                                                          <  2:175326/150‑1 (MQ=255)                                                                                                        ccATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTTTTTTTTTT‑CCTCCACCTATTGCGGTGCTCCTTCCCCATATGCGTGGGGTGGTGGGGGGCCGCGGCTAAATGATCCTCCCGCCCTACAGCCCCACCCCTTTCGGTGCtttt                                             >  1:187788/1‑150 (MQ=255)                                                                                                         cATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGGCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGAACCTCGGGGCAGACAGCCACACCCGTTACGGTGCAtt                                              >  1:309705/1‑149 (MQ=255)                                                                                                                ttGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCCCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGGGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACaa                                      <  2:138510/150‑1 (MQ=255)                                                                                                                         cggcggCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCg                             <  2:278125/150‑1 (MQ=255)                                                                                                                         cggccgCCCAACGGTATGGCGGTATTTTTGTG‑CCTCCGAAGATGGCGGCGATCCAGCAATCTAGGCTTGAGAGGATGGAAGGCGGCGGCAAAAAAATCGCGGGGTAAGACAGCCGCAGCCGGTACGGGGCAGTAGGGACAACGGcg                                 <  2:711941/146‑2 (MQ=255)                                                                                                                            cggcCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGTGACAATGGCAGTCGGTGa                         <  1:575498/151‑1 (MQ=255)                                                                                                                                ccAGCGTTATGGCGGTATTTTTGTT‑CCTCCCCATATTTCGGGCCTCCATCAATATATGCGGGAGATGATGGCAGGGGGCGGCAAAATGATCCTCCGGGCAGCAAGCCCCACCCGGTACGGGGCCTTAGGGACAATGGCCGTCgggggggg                      >  2:562489/1‑150 (MQ=255)                                                                                                                                          ggCGGTATTTTTGTG‑CCTCCGCATATTGCCGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAg            <  1:643332/150‑1 (MQ=255)                                                                                                                                            cGGTATTTTTGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTg         >  2:2032/1‑151 (MQ=255)                                                                                                                                             ggTATTTTTGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCAGTGAGTTgg        <  1:639822/151‑1 (MQ=255)                                                                                                                                                   tttGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTaaaa   <  2:680817/150‑1 (MQ=255)                                                                                                                                                    ttGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAAc  >  1:396588/1‑150 (MQ=255)                                                                                                                                                    ttGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAAc  <  1:924882/150‑1 (MQ=255)                                                                                                                                                    ttGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAAc  >  2:965075/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                     AGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTG‑CCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAAC  >  NC_000913/803720‑804016

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