Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,029,965	0	A	C	48.1%	64.7 / 20.1	52	L323V (TTA→GTA)	ygfT	putative oxidoreductase, Fe‑S subunit/nucleotide‑binding subunit
Reads supporting (aligned to +/- strand):  ref base A (21/6);  new base C (0/25);  total (21/31)
Fisher's exact test for biased strand distribution p-value = 1.61e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTGAATATCTCTCGCCGCTGGCTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATAGAGACTGCGCCCGAGTGATCTTTCAA  >  NC_000913/3029818‑3030112                                                                                                                                                    |                                                                                                                                                     cGGTGAATATCTCTCGCCGCTGGCTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAAccc                                                                                                                                                   >  2:687744/1‑151 (MQ=255) cGGTGAATATCTCTCGCCGCTGGCTTAATACCGTTTTATCAAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAAccc                                                                                                                                                   <  1:430602/151‑1 (MQ=255)  ggTGAATATCTCTCGCCGCTGGCTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCt                                                                                                                                                  >  1:320812/1‑151 (MQ=255)         atCTCTCGCCGCTGGCTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCT‑GCAGGc                                                                                                                                           <  1:244122/150‑1 (MQ=255)                gccgcTGGCTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCT‑GCAGGCCCAGCGcc                                                                                                                                   >  1:108286/1‑151 (MQ=255)                       gCTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCT‑GCAGGCCCAGCGCCAATCAcc                                                                                                                            <  1:838634/151‑1 (MQ=255)                       gCTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCT‑GCAGGCCCAGCGCCAATCAcc                                                                                                                            >  1:58496/1‑151 (MQ=255)                        cTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCT‑GCAGGCCCAGCGCCAATCACCg                                                                                                                           >  2:182846/1‑151 (MQ=255)                                   ttaatggAGTTTAAAAGGAGGAATGCCAAAAGTCAGCATCCCGCCAATTTCTGGATGGCAATAAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTACCCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTc                                                                                                                 <  1:849895/144‑1 (MQ=255)                                      aTCGAGTTAGAAGGGAGAAAGCCAAAAAGTCACCATCCCGCAATTTTCTGGATGGCAATAAAAGACACCGCCCTGACCCCCTGCGCGCGCCAGAATACCAGCACCCCCTACCCCG‑GCGGGCCCAGCGCCAATCACCGCCACTTTTTCGCt                                                                                                              <  1:303100/150‑1 (MQ=255)                                       tCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTAc                                                                                                            >  1:641539/1‑151 (MQ=255)                                            tttGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACgggg                                                                                                        >  2:338622/1‑150 (MQ=255)                                            tttAAAAGGAGGAATGCAAAAAGTCAGCATCCCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGACCTCCTGCGCGCGCAAGAATATCAGCACCCCCTACCCCT‑GCAGGCCCAGCCCCAATCACCGCCACTTTTTCGCTAc                                                                                                            <  2:352173/146‑1 (MQ=255)                                                aaaGGAGGAATGCCAAAAGTCAGCATCCCGCCATTTTCTGGATGGCGATCAAAGACACCGCCCTGACCTCCTGCGCGCGCAAGAATATCAGCACCCCCTACCCCG‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGaac                                                                                                     <  2:912657/149‑1 (MQ=255)                                                   ggaggaATGCAAAAAGTCAGCATCCCGCAAATTTCGGGAGGGCGATAAAAGCCACCGACCTGACCCCCTGCGCGCGCAAGAATATCAGCACCCCCTACCCCT‑GCGGGCCCAGCGCAAATCACCGCCACTTTTTCGCTACGGGGAACAACt                                                                                                 <  2:435717/150‑1 (MQ=255)                                                      ggaAGGCAAAAAGTCAGCATCCCGCAATTTTCTGGATGGCGATAAAAGACATCGCCCTGAACTCCGGCGCGCGCCAGAATATCAGCCCCCCCTACCCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTg                                                                                              <  2:863267/150‑1 (MQ=255)                                                        aaTGCCAAAAGCCACCACCCCGCCATTTTCTGGATGGCAATAAAAGCCATCGCCCTGACCCCCGGCGCGCGCAAAAATATCAGCACACCCAACCCCG‑GCAGGCCCAGCCCCAATCCCCGCCACTTTTTCGCTACGGGGAACAACttt                                                                                               <  1:491699/147‑1 (MQ=255)                                                            caaaagGTCAGCATCCCGCCATTTTCGGGAGGGAAAAAAAAGACACCGCCCTGACCCCCGGCGCGCGCAAGAATACCAGCACCCCCTACCCCT‑GCAGCCCCAGCGACAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACAt                                                                                       <  1:478275/149‑1 (MQ=255)                                                              aaagGTCAGCAACCCGCGATTTTGTGGATGGCAATCAAAGACATCGACCTGAACCCCTGCGCGCGCCAGAATATCAGCACACCCAACCCCG‑GCAGCCCCAGCGCAAAGCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCg                                                                                     <  1:954257/151‑1 (MQ=255)                                                              aaaaGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCg                                                                                     >  2:923790/1‑151 (MQ=255)                                                                       cATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGcc                                                                              >  1:853100/1‑149 (MQ=255)                                                                             ccaaaTTTCGGGATGGCGATAAAAGCCATCGCCCGGACCCCCGGCGCGCGCAAAAATATCACCACACCCAACCCCT‑GCAGCCCCAGCCCCAACCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATc                                                                      <  1:687675/148‑1 (MQ=255)                                                                              ccATTTTCGGTATGGCAACCAAAACCATCGCCCTTACCCCCGGCGCGCGCAAGAATATCAGCACCCCCAACCCCT‑GCAGCCCCGGCGCCAATCTCCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGcc                                                                          <  1:967021/146‑1 (MQ=255)                                                                                   ttCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAgc                                                                >  1:839960/1‑151 (MQ=255)                                                                                       ggATGGCGATCAAAGACATCGCCCTGAACCCCTGCGCGCGCCAAAATATCAGCACCCCCTACCCCG‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCgg                                                            <  2:378090/151‑1 (MQ=255)                                                                                       ggAGGGCGACAAAAGCCACCGACCTGACCTCCGGCGCGCGCAAGAATACCACCACCCCCTACCCCT‑GCAGCCCCAGCCCAAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGcgcg                                                             <  1:7698/150‑1 (MQ=255)                                                                                       ggAGGGCGAACAAAGACATCGCCCGGACCTCCTGCGCGCGCAAGAATATCACCACCCCCTACCCCG‑GCAGGCCCGGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGcgcg                                                             <  1:504583/150‑1 (MQ=255)                                                                                        gATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGt                                                           >  2:477899/1‑151 (MQ=255)                                                                                             cGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCgg                                                      >  1:534798/1‑151 (MQ=255)                                                                                               aTCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTg                                                    >  2:279941/1‑151 (MQ=255)                                                                                                  aaaGACATCGACCTGACCTCCTGCGCGCGCCAAAATATCAGCACCCCCTACCCCG‑GCAGCCCCAGCGCAAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGAt                                                  <  2:809440/150‑1 (MQ=255)                                                                                                    aGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGt                                                >  2:895395/1‑150 (MQ=255)                                                                                                       ctcccaCCTGACCTCCTGCGCCCCCCAAAATACCACCCCCCCCTACCCCT‑GCAGCCCCAGCCCCAATCCCCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTg                                                    <  2:833034/138‑1 (MQ=255)                                                                                                       caccgccctaacccccTGCGCGCGCAAAAAAACCAGCACCCCCTACCCCG‑GCAGGCCCAGCCCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGc                                             <  2:639788/136‑1 (MQ=255)                                                                                                       caccGACCTAACCCCCTGCGCGCGCAAGAATACCACCCCACCCTACCCCT‑GCAGCCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCg                                            <  1:792340/148‑1 (MQ=255)                                                                                                        accGACCTGACCTCCGGCGCGCGCCAAAATATCAGCCCCCCCTACCCCG‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGt                                           <  2:947470/149‑1 (MQ=255)                                                                                                          cgcccggAACTCCTGCGCGCGCCAGAATATCAGCACCCCCTACCCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTc                                         <  2:792234/145‑1 (MQ=255)                                                                                                               tGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTcc                                        >  1:537312/1‑147 (MQ=255)                                                                                                                 accTCCGGCGCGCGCCAGAATATCAGCCCCCCCTAACCCT‑GCAGCCCCAGCCCCAATCTCCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGAt                                    <  1:203598/147‑1 (MQ=255)                                                                                                                    tCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCg                               >  2:311964/1‑151 (MQ=255)                                                                                                                     ccTGCGCGCGCAAAAATATCAGCACCCCCTACCCCG‑GCGGGCCCAGCGCCAAACACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCACCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCg                               <  1:718634/150‑1 (MQ=255)                                                                                                                        gcgcgcgcGAGAATATCAGCACACCCTAACCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGTCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGa                              >  1:839589/1‑148 (MQ=255)                                                                                                                        gcgcgcgcAAGAATACCACCCCCCCCAAACCCT‑GCAGCCCCAGCGCAAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATa                            <  2:531620/150‑1 (MQ=255)                                                                                                                            gcgcCAGAATATCAGCACACCCTAACCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATAGAGAc                       >  1:910999/1‑151 (MQ=255)                                                                                                                                 aaaataccaccccccccTACCCCT‑GCAGCCCCAGCGCCAACCCCCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATAGAGACTgcg                   <  2:345673/136‑1 (MQ=255)                                                                                                                                 aGAATATCAGCACACCCTAACCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATAGAGACTgcgc                  >  1:837157/1‑151 (MQ=255)                                                                                                                                   aaTATCACCACCCCCTACCCCG‑GCAGCCCCAGCGCCAACCTCCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATAGAGACTGCGccc                <  1:712503/151‑1 (MQ=255)                                                                                                                                    atatCAGCACACCCTAACCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATAGAGACTGCGCCCg               >  2:894886/1‑151 (MQ=255)                                                                                                                                     taccAGCCCCCCCTACCCCT‑GCAGCCCCAGCCCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATAGAGACTGCGCCCg               <  2:724160/147‑1 (MQ=255)                                                                                                                                           cccccccTACCCCC‑GCGGCCCCGGCGCAAATCCCCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATAGAGACTGCGCCCGAGTGa          <  1:198903/145‑1 (MQ=255)                                                                                                                                              cccccAACCCTGGCAGCCCCAGCTCCATTCGCCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATAGAGACTGCGCCCGAGTGATCtt      <  2:641604/150‑1 (MQ=255)                                                                                                                                               cccTAACCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATAGAGACTGCGCCCGAGTGATCTTTc    >  1:753573/1‑151 (MQ=255)                                                                                                                                                   aaCCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATAGAGACTGCGCCCGAGTGATCTTTCaa  >  1:232403/1‑149 (MQ=255)                                                                                                                                                   aaCCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATAGAGACTGCGCCCGAGTGATCTTTCaa  >  2:955682/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CGGTGAATATCTCTCGCCGCTGGCTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCT‑GCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATAGAGACTGCGCCCGAGTGATCTTTCAA  >  NC_000913/3029818‑3030112

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