Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,065,598	0	G	A	42.9%	56.8 / 21.3	50	L57L (CTG→TTG)	yihO	putative sulphoquinovose importer
Reads supporting (aligned to +/- strand):  ref base G (22/6);  new base A (0/21);  total (22/28)
Fisher's exact test for biased strand distribution p-value = 7.74e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTC  >  NC_000913/4065449‑4065734                                                                                                                                                          |                                                                                                                                         gATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGa                                                                                                                                         <  2:206651/151‑1 (MQ=255)    aaaCTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGaaat                                                                                                                                      <  1:366534/151‑2 (MQ=255)    aaaCTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGaa                                                                                                                                        >  2:404961/1‑149 (MQ=255)    aaaCTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAg                                                                                                                                      <  2:291379/151‑1 (MQ=255)     aaCTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGaa                                                                                                                                        >  2:320258/1‑148 (MQ=255)     aaCTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGa                                                                                                                                         >  2:544995/1‑147 (MQ=255)       cTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAGAt                                                                                                                                    >  2:582221/1‑150 (MQ=255)         ggAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAGATaa                                                                                                                                  >  2:682886/1‑150 (MQ=255)           aGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAGATAATcc                                                                                                                               <  2:343630/151‑1 (MQ=255)             cGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAGATAATccc                                                                                                                              >  1:397338/1‑150 (MQ=255)             cGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAGATAATCCCg                                                                                                                             >  1:856791/1‑151 (MQ=255)                            cgccgGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGc                                                                                                              >  1:77504/1‑151 (MQ=255)                                     cGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAAt                                                                                                     >  2:460331/1‑151 (MQ=255)                                           atatAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCg                                                                                               >  2:397559/1‑151 (MQ=255)                                           ataaAAAAAAAAGGGACAAAATTTGCTTTTAGGCCAAATATTTTC‑CCGCGAGTCGAGAAAAAACCCGGTGAGCATAGCGGAAAACGCGGAAAAAAA‑CTC‑CGCCCCA‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATc                                                                                                <  1:780011/150‑1 (MQ=255)                                             aaaaaataaaGCAACTAAATTGGTTTTTA‑GCCCATTTTTTTTCCCGCGAGTGGGAAAAAAAGCGGGTGAGCATATCGATAACGGCGTAAAAAAA‑CTC‑CGCCACG‑AAAAAGAAACCCCCGCCATAGTAAGCAGGCATCCCTAATTCATCgg                                                                                              <  2:214362/148‑1 (MQ=255)                                             aaaaaaaaaaGGGACAAAATTTGCCTTTAGGCCCAATATTTTT‑CCGCGAGTCGAGAAAAAATCCGGTGAGCATATCGGTAAACGCGGAAAAAAA‑CTT‑CGCCACC‑AAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCg                                                                                               <  2:317078/147‑1 (MQ=255)                                               aaaaaaaaGGAACAAAATTGGCCTTTAGCCCAAATATTTTC‑CCGGGAGTCGATAAAAAATCCGGTGAGCATATCGGTAAACGGGGAAAAAAA‑CTT‑CGCCACC‑AAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtat                                                                                           <  2:574155/151‑1 (MQ=255)                                                 aaaaagggaaagAAATTCGCCTTTAGGCCCATTTTTTTC‑CAGCGAGGGGAAAAAAAACCGGGGGAGCTAATGGAAAAACG‑GGTAAAAAA‑TCT‑GCCCACCAAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATaa                                                                                         <  2:18059/140‑1 (MQ=255)                                                  ataaaGGGACGAAATTTGCCTTTAGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATaa                                                                                         >  2:278919/1‑150 (MQ=255)                                                       gggACGAAATTTGCCTTTAGGCCAAATTTTTTC‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAC‑CTT‑CGCCACC‑AAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTg                                                                                   <  2:765463/151‑1 (MQ=255)                                                         gACGAAATTTGCCTTTAGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTg                                                                                   >  1:21878/1‑149 (MQ=255)                                                            gcaaTTTGCCTTTAGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGAATCCCTAATTCATCGGTATAAAATTTgagg                                                                                >  1:635249/3‑149 (MQ=255)                                                               attgcacTTTATGCCAAATATTTTT‑CCGCGAGGCGATAAAAAAGCAGGTGACCTCACGGGAAGACGCGGAAAAAAA‑CAT‑CCCCACA‑AAAAAAATATCCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATa                                                                           <  1:758981/145‑1 (MQ=255)                                                                  tGCCTTTAGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATAc                                                                          >  1:631725/1‑149 (MQ=255)                                                                    cttttAGGCCAAATATTTTC‑CCGCGAGCATAGAAAAAATCCGGTGAGCAAACCGTAAACCGCGGAAAAAAT‑CTC‑CCCCACA‑AAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACa                                                                         <  1:700903/146‑1 (MQ=255)                                                                    ccTTTAGGCCAAATATTTTC‑CCGCGAGTCGAGAAAAAACCCGGTGAGAATATCGGTAAACGCGGAAAAAAC‑CTT‑CGCCACC‑AAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACaga                                                                       <  2:12899/150‑1 (MQ=255)                                                                       ttAGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTAc                                                                   >  1:34318/1‑151 (MQ=255)                                                                         aGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCg                                                                 >  1:471920/1‑151 (MQ=255)                                                                            cccAATTTTTTC‑CCGCGAGTCGAGTAAAAACCCGGTGAGCATATCGGTAAACGAGTAAAAAAA‑CTT‑CGCCACC‑AAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATg                                                              <  2:140506/151‑1 (MQ=255)                                                                            cccAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATg                                                              >  1:584657/1‑151 (MQ=255)                                                                               aaTATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGc                                                             >  2:445318/1‑149 (MQ=255)                                                                                   ttttc‑CCGCGAGTGGAAAAAAAATCCGGGAGGAATACCGGTAAACG‑GGTAAAAAATCTT‑CGCCACC‑GAAAAGATATTCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga                                                          <  2:273717/148‑1 (MQ=255)                                                                                     ttt‑CCGCGAGCCGAGAAAAAACCGGGTGAACATTCGGGAAAACGCGAAAAAAAA‑CTC‑CCCCACA‑AAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCAt                                                      <  1:703279/150‑1 (MQ=255)                                                                                         ccGCGAGTGAAGAAAAAAGCGGGTGAGCATATCGTAAAACGCGG‑AAAAAA‑CTTCCGCCACA‑AAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTaa                                                   <  1:766743/150‑1 (MQ=255)                                                                                                   agaaaaaaCCGGGTGAGCATATGGGAAAACGCGAAAAAAAC‑CTC‑CCCCACA‑AAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGaa                                            <  1:280577/144‑1 (MQ=255)                                                                                                     aaaaaaCCCGGTGAACATATCGAAAAACGGGAAAAAAAC‑CTC‑CACCACC‑AAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACg                                      <  1:609742/150‑1 (MQ=255)                                                                                                      aaaaaTACGGGGAGCATATCGAAAAACGGGATAAAAAA‑CTC‑CGCCACA‑AAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGAc                                       <  1:807882/149‑1 (MQ=255)                                                                                                      aaaaaCCCGGTGAGCAAACGGGTAAACGCGGAAAAAAC‑CTC‑CCCCACA‑AAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGaa                                            <  2:524126/144‑1 (MQ=255)                                                                                                         aaTCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGc                                  >  2:325585/1‑151 (MQ=255)                                                                                                           tCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGccc                                >  1:498144/1‑151 (MQ=255)                                                                                                              ggggAGCATATCGGAAAACGCGGTAAAAAA‑CTC‑CGCCACA‑AAAAAAATATCCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccat                              <  1:678883/147‑1 (MQ=255)                                                                                                              ggTGAGCATATCGTTAAACGCGAAAAAAAA‑CTT‑CCCCCCC‑AAAAAGATACTCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccata                             <  2:24992/151‑1 (MQ=255)                                                                                                                  agaattccGGAAAACGGGGAAAAACA‑TTC‑CCCCTCA‑AAAAAAAAAAACCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccataccat                         <  1:491342/144‑1 (MQ=255)                                                                                                                     atatCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAgg                      >  1:834729/1‑151 (MQ=255)                                                                                                                      tatCGGAAAACGCGGAAAAAAC‑CTC‑CGCCACA‑AGAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccataccata                        <  1:365853/148‑1 (MQ=255)                                                                                                                             aacacgtgTAAAAAT‑CTC‑CCCCACA‑GTAAAAATAGTCCCGCCATAGTAAGCAGGCAGCCCTGATTCATCGGTATAAAAGTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGACCGGACGTCGCCCATACCATAGGCGAttt                <  1:235258/142‑1 (MQ=255)                                                                                                                                caagatAAAAAA‑CTC‑CGCCTCA‑AAAAAGAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGAttt                <  2:580790/141‑1 (MQ=255)                                                                                                                                    aaaaaaaa‑CTT‑CCCCACA‑AAAAAAATAAACCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCg         <  1:461635/147‑1 (MQ=255)                                                                                                                                         aaa‑CTT‑CGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTc  >  1:828454/1‑151 (MQ=255)                                                                                                                                                          |                                                                                                                                         GATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTT‑CCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTT‑CGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTC  >  NC_000913/4065449‑4065734

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