Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	803,869	0	G	T	45.6%	72.9 / 25.2	46	V123L (GTG→TTG)	ybhJ	aconitase family protein
Reads supporting (aligned to +/- strand):  ref base G (1/24);  new base T (21/0);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 3.17e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGC‑GGCGGCAAAATGATCCTCGGGTCAGACAG‑CCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGC  >  NC_000913/803725‑804019                                                                                                                                                  |                                                                                                                                                         gCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTTTGcc                                                                                                                                                      >  1:12246/1‑149 (MQ=255)   tAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTTTGCCTcc                                                                                                                                                   >  1:394977/1‑150 (MQ=255)      aGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCat                                                                                                                                               <  1:547776/151‑1 (MQ=255)          tCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTTTGTCTCCCCCTAtt                                                                                                                                            >  1:295565/1‑150 (MQ=255)             ggTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTTTGCCTCCCCATATTTCgg                                                                                                                                        >  1:56365/1‑151 (MQ=255)                       gTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTTTGCCTCCGCCTATTTCGGGcatcc                                                                                                                                  >  1:524077/1‑147 (MQ=255)                           cccGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAatat                                                                                                                          <  1:380833/151‑1 (MQ=255)                                  ccGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTg                                                                                                                   <  2:598760/151‑1 (MQ=255)                                  ccGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTg                                                                                                                   <  1:515609/151‑1 (MQ=255)                                  ccGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTg                                                                                                                   <  1:695508/151‑1 (MQ=255)                                      aTGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAgat                                                                                                               <  1:149690/151‑1 (MQ=255)                                             gACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTTTGCCTCCGCATATTTCGGGCATCCCTCAATATATGCGTGGGGAGATgg                                                                                                          >  1:182068/1‑149 (MQ=255)                                              aCCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAg                                                                                                       <  2:548848/151‑1 (MQ=255)                                               ccAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTTTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAgg                                                                                                      >  1:152803/1‑151 (MQ=255)                                                cAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAggc                                                                                                     <  1:595553/151‑1 (MQ=255)                                                 aaCTGCCATAGCTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTGATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAgg                                                                                                      <  2:419524/149‑1 (MQ=255)                                                      ccATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAggc‑ggcgg                                                                                               <  1:198033/150‑1 (MQ=255)                                                      ccATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGGATTTTTTTGGCCCCGCCTATTTCGGGCATCCCTCAATATTTTCGGGGGGTGATGGGAggg‑ggcgg                                                                                               >  2:385420/1‑150 (MQ=255)                                                       cATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTTTGTCTCCGCCTATTTCGGTCATCCATCAATATATGCGGGAGAAGATGGCAggg‑ggcgg                                                                                               >  2:197881/1‑149 (MQ=255)                                                        aTAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTTTGCCACCGCAATTTTCGGTCATCCCTTCATATTTTCGTGGGATTGTGGGAggc‑gggg                                                                                                >  2:250678/1‑145 (MQ=255)                                                             tcactcTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGGATTTTTTTGTCTCCCCGTATTGTGGTGCTTCATCAATATATTGTGGGGAGGATGGCAAGG‑GGGGGCAAAATg                                                                                        >  1:707950/1‑150 (MQ=255)                                                               actcTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTTTGCCTCCCCATATTGCGGTCATCCATCAATATATGCGGGGGATGATGGGAGGC‑GGGGGCaaaa                                                                                          >  1:2929/1‑146 (MQ=255)                                                                ctctGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGGATTTTTTTGCCCCCCCCTTTTTCGGGCATCCCTCAATATATTCGTGAGATGGTGGGCGGCGGGGGGCaaaagaac                                                                                      >  2:14669/1‑146 (MQ=255)                                                                   tGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTTTGCCTCCGCCTATTGCGGTCATCCATCAATATATGCGGGGGATGATTGCAGGG‑GGCGGGAAAAAGAAcccc                                                                                  >  1:670193/1‑148 (MQ=255)                                                                    gcgcCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGC‑GGCGGCAAAATGATCCTCg                                                                                 <  1:696065/150‑1 (MQ=255)                                                                        cGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGC‑GGCGGCAAAATGATCCTCGGGTCa                                                                            <  1:341100/151‑1 (MQ=255)                                                                                        aaCGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTTTGCCTCCCCATATTGCGGGCATCCCTCAAAATATTCGTGAGATGATGGCAGGC‑GGCGGCAAAATGATCCCCGGGGCAGACAG‑CCACACCCGtt                                                           >  1:499767/1‑151 (MQ=255)                                                                                        aaCGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGC‑GGCGGCAAAATGATCCTCGGGTCAGACAG‑CCACACCCGtt                                                           <  2:410535/151‑1 (MQ=255)                                                                                        aaCGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGC‑GGCGGCAAAATGATCCTCGGGTCAGACAG‑CCACACCCGtt                                                           <  1:94320/151‑1 (MQ=255)                                                                                          cGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGC‑GGCGGCAAAATGATCCTCGGGTCAGACAG‑CCACACCCGTTAc                                                         <  2:648101/151‑1 (MQ=255)                                                                                          cGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGC‑GGCGGCAAAATGATCCTCGGGTCAGACAG‑CCACACCCGTTAc                                                         <  1:184721/151‑1 (MQ=255)                                                                                              gatgaCCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTTTGTCTCCGCATATTTCGGGCATCCCTCCAAATATGCGTGGGATGATGGGAGGG‑GGGGGGAAAAAGATTCTCCGGGTAGACCG‑GCCCACCCGGTatagg                                                      >  2:233946/1‑146 (MQ=255)                                                                                                 gaCCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGGAATTTTTTTCCTCCGCATTTTTCGGGCCTCCCTCAATATATGCGTGGGATGATGGGAGGCGGGGGGCAAAATAAT‑CTCCGGTCAGAAAA‑CCACACCCCCTACGGGGcgt                                                  >  2:30418/1‑149 (MQ=255)                                                                                                           ttGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGC‑GGCGGCAAAATGATCCTCGGGTCAGACAG‑CCACACCCGTTACGGTGCATTAGGGACaa                                         <  2:167297/150‑1 (MQ=255)                                                                                                             ggTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTTTGCCTCCGCATATTGCGGGCATCCCTCAATATATGCGGGGGATGATGGCAGGC‑GGGGGCAAAATGATCCTCGGGGCAGACAG‑CCCCACCCGGTACGGGGCATTATGGACaa                                         >  1:691386/1‑148 (MQ=255)                                                                                                               tttATCGGCGGCCCAGCGTTATGGCGGGAATTTTTTTGCCT‑CGCCTATTTCGGGCATCCATTCAAATATACGGGTGATGGGGGCAGGT‑GGCGGGGAAATGATCCTCCGGTCACAAAGCCCAACCCCGTTACgcggcttaagggcca                                          >  2:39688/1‑132 (MQ=255)                                                                                                               tttATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTTTGCCTCCGCCTATTGCGGGCATCCATCAATATATGCGTGGGATGATGGCAGGC‑GGCGGGAAAAAGATCCCCGGGGCAGACAG‑CCACACCCGCTACGGGGCATTAGGGGCaa                                         >  1:492376/1‑146 (MQ=255)                                                                                                                ttATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGC‑GGCGGCAAAATGATCCTCGGGTCAGACAG‑CCACACCCGTTACGGTGCATTAGGGACAATGGCa                                    <  1:258980/150‑1 (MQ=255)                                                                                                                  aTCGGCGGCCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGC‑GGCGGCAAAATGATCCTCGGGTCAGACAG‑CCACACCCGTTACGGTGCATTAGGGACAATGGCAGTc                                 <  1:724236/151‑1 (MQ=255)                                                                                                                     ggcggcCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGC‑GGCGGCAAAATGATCCTCGGGTCAGACAG‑CCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGt                              <  2:394974/151‑1 (MQ=255)                                                                                                                     ggcggcCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGC‑GGCGGCAAAATGATCCTCGGGTCAGACAG‑CCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGt                              <  1:339794/151‑1 (MQ=255)                                                                                                                             aGCGTTATGGC‑GGTATTTTTTTGCCTCCGCCTATTGCGGTCACCCCTCAAAATATGAGTGAGGTGATGGCAGGG‑GGGGGCAAAATGAACCTCGGGTCAGACAG‑CCACACCCGTTACGGGGTATTAGGGAGAATGGGAGTCGGGGAgggggg                      >  2:497484/1‑151 (MQ=255)                                                                                                                                gTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGC‑GGCGGCAAAATGATCCTCGGGTCAGACAG‑CCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGgcggcg                    <  2:411223/150‑1 (MQ=255)                                                                                                                                       c‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGC‑GGCGGCAAAATGATCCTCGGGTCAGACAG‑CCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTg            <  1:556349/151‑1 (MQ=255)                                                                                                                                             tttttGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGC‑GGCGGCAAAATGATCCTCGGGTCAGACAG‑CCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTaa        <  1:74307/150‑1 (MQ=255)                                                                                                                                                ttGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGC‑GGCGGCAAAATGATCCTCGGGTCAGACAG‑CCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAAc     >  1:576273/1‑150 (MQ=255)                                                                                                                                                  gtgCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGC‑GGCGGCAAAATGATCCTCGGGTCAGACAG‑CCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAgc  <  1:104945/151‑1 (MQ=255)                                                                                                                                                  gtgCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGC‑GGCGGCAAAATGATCCTCGGGTCAGACAG‑CCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAg   <  2:427249/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                         GCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGC‑GGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGC‑GGCGGCAAAATGATCCTCGGGTCAGACAG‑CCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGC  >  NC_000913/803725‑804019

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