Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	657,676	0	A	C	43.8%	68.2 / 12.2	48	V88G (GTG→GGG)	flc	fluoride efflux channel, dual topology membrane protein
Reads supporting (aligned to +/- strand):  ref base A (21/6);  new base C (0/21);  total (21/27)
Fisher's exact test for biased strand distribution p-value = 1.66e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCC‑AGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCA  >  NC_000913/657538‑657814                                                                                                                                            |                                                                                                                                           cGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCC‑AGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCg                                                                                                                                <  2:461672/151‑1 (MQ=255)   ggTTTTTTTTGGTCTTTAGTGCGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAACGACCCC‑AGAAGGCTGACGAAAACGATCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGcc                                                                                                                                 >  1:234059/1‑148 (MQ=255)              gTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCC‑AGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTgaga                                                                                                                    >  2:625087/1‑150 (MQ=255)                  ttAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCC‑AGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGgtt                                                                                                               >  2:395058/1‑151 (MQ=255)                   tAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCC‑AGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGgttg                                                                                                              >  1:542450/1‑151 (MQ=255)                   tAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCC‑AGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGgtt                                                                                                               <  1:525548/150‑1 (MQ=255)                   tAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCC‑AGAAGGATGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGTAAATGTTGAGAAGgttg                                                                                                              >  1:775178/1‑151 (MQ=255)                    aGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCC‑AGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGgttgt                                                                                                             >  1:229664/1‑151 (MQ=255)                    aGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCC‑AGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGgttgt                                                                                                             >  2:604209/1‑151 (MQ=255)                     gtgtGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCC‑AGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGgttgtt                                                                                                            <  2:296091/151‑1 (MQ=255)                           ggTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCC‑AGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGAcc                                                                                                       >  1:2841/1‑150 (MQ=255)                                    cGAAAACAGCCAGAATGCCGGTGCGGTCATGGCAAAAGACCCC‑AGAAGTTTGACGAAAACTTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAAt                                                                                             <  2:482239/151‑1 (MQ=255)                                     gAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCC‑AGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATc                                                                                            >  2:416020/1‑151 (MQ=255)                                      aaaaCAGCCAGAATCCCGGTGCGGTCATGGCAAAAGACCCC‑AGAAGGTTGACAAAAACGTTCAGTAATGCCCAGCAAAGGCGGCCCTCTTGTAACAAAAACCCCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAAt                                                                                             <  1:217413/149‑1 (MQ=255)                                          cAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCC‑AGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCggtg                                                                                       >  2:198112/1‑151 (MQ=255)                                             ccAGAATGCCAGTGCGGTCATGGCAAAAGACCCC‑AGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACCCCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCggtgg                                                                                      <  2:269202/149‑1 (MQ=255)                                             ccAGAATCCCAGTGCGTTCATGGAAAAAGCCCCA‑AGAAGGTTGACGAAAACGTTAAGAAAGCCCCAGCCAAAGCGGCCTTCTTGTAACAAAACCCCCACTTCTGCCGAAAATGTTGAGAGGGTTGTTAGACCGCCACAAAATCCGgtggtg                                                                                    <  1:762000/151‑1 (MQ=255)                                                     ccAGTGGGTTCAGGGAAAAAGCCCCA‑AAAAGGTGGACAAAAACGTTAGGAAAGGCCCAGCCAAAGCGCCCCTTTTAAAAAAAAACCCCCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATAc                                                                            <  1:627447/151‑1 (MQ=255)                                                                aTGGAAAAAGCCCCC‑AGAAGTTTAACGAAACCGTTCAGTAATGCCCAGCCAAAGCGCCCCTCTTGAAACAAAAACCCCTCTTCTGCCAAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACt                                                                  <  2:7829/150‑1 (MQ=255)                                                                aTGGAAAAAGACCCA‑AGAGGGTAGAAGAAACCTTTCATTAATGCCCGCCAAAAGCGCCCCTCTTGTAACAAAACCCCCTCTTCTCCCAAAAATGTTGAGAAGTTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCacac                                                                   <  2:289023/149‑1 (MQ=255)                                                                 tGGCAAAAGACCCC‑AGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTgg                                                                >  2:703535/1‑151 (MQ=255)                                                                  ggaaaaaGCCCCC‑AGAAGGTGGACAAAAACGTTCAGTAATGCCCCGCAAAGGCGGCCCTCTTGAAACAAAAACCCCACTTCTGCGGAAAATGTTGAGAAGTTTTTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGa                                                               <  1:329236/148‑1 (MQ=255)                                                                  ggCAAAAGCCCCC‑AGAGGTTTGACGAAAACTTCCAGTAATGCCCAGCCAATGCGCCCCTTTTTTAACAAAAACCCCACTCCCGCCGAAAATTTTAAGAAGGTTTTTGGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGa                                                               <  1:225888/151‑1 (MQ=255)                                                                  ggCAAAAGACCCC‑AGAAGGTTGACAAAACCGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAACCCCCCCTTCTGCCGAAAATGTTGAGAAGTTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTgg                                                                <  2:14428/150‑1 (MQ=255)                                                                   ggaaaaaaCCCCAAAAAGTTTGGCGACAAGGTTCAGTAATGCCCTGCCAAAGCGGCTTTCTTGTAACAAAAACACCCCTTCGGCCGAAAATGTTGAATAGGTTTTTAAACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGa                                                               <  1:747144/149‑1 (MQ=255)                                                                    cAAAAGACCCC‑AGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATc                                                             >  1:73666/1‑151 (MQ=255)                                                                     aaaaGACCCC‑AGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATc                                                             >  1:168016/1‑150 (MQ=255)                                                                             ccaggAAGGTTGACAAAAACTTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAACCCCCCCTTCTGCCGAAAATTTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCg                                                     <  1:798972/147‑1 (MQ=255)                                                                              c‑AGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGt                                                    >  2:257541/1‑150 (MQ=255)                                                                                  aaGGTTGACAAAAACTTTCATTAATGCCAAGCAAAAGCGGCCCTCTTGTAACAAAAACCCCACTTCTGCCAAAAATTTTGAGAAGGTTTTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCAt                                                 <  2:529319/150‑1 (MQ=255)                                                                                    ggTTGACAAAAACGTTCAGTAATGCCCAGCAAAAGCGGCCCTCTTGTAACAAAAACCCCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCct                                              <  2:328909/151‑1 (MQ=255)                                                                                           gAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAAc                                       >  2:748148/1‑151 (MQ=255)                                                                                           gAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAAc                                       >  1:679829/1‑151 (MQ=255)                                                                                           aaaaaCGTTCGGAAACCCCGAGCCAACGCGGCCCTCTTGTAACAAAACCCCCCCTCCGGCGGAAAATTTTGAAAAGGTTGTTAGACCGCCACAAAATCCGGTGGCGATTAATACTTTCCACACTGGATCAATGTTCGTCATCctgctg                                          <  1:395064/147‑1 (MQ=255)                                                                                            aaaaCGTTCAGTAATGCCCAGCAAAAGCGGCCCTCTTGAAAAAAAAACCCCCCTTCTGCCGAAAATGTTGAGAGGGTTTTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAAc                                       <  2:234060/150‑1 (MQ=255)                                                                                                   tCATTAATCCCCAGCCAAAGCGGCCCTTTTGTAACAAAAACACCACTTCTGCCAAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGa                                <  1:702330/150‑1 (MQ=255)                                                                                                    cAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAt                              >  1:660294/1‑151 (MQ=255)                                                                                                     attAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACCCCCCTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtc                             <  2:511630/149‑1 (MQ=255)                                                                                                     aGTAATGCCCAGCAAAAGCGGCCCTTTTGTAAAAAAACCCCCACTTCTGCCGAAAATGTTGAAAGGTTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGa                                <  2:542413/148‑1 (MQ=255)                                                                                                        aaTGCCCAGCCAAAGCGGCCCTCTTGAAACAAAAACCCCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcct                           <  2:598174/150‑1 (MQ=255)                                                                                                                 caaaaGCGGCCCTTTTTTAACAAAAACCCCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTatg                  <  2:253883/148‑1 (MQ=255)                                                                                                                       cGGCCCTCTTGTAACAAAAACCCCACTTCGGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAAt            <  2:345958/150‑1 (MQ=255)                                                                                                                           ccTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGcccc       >  1:24762/1‑151 (MQ=255)                                                                                                                           ccTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGcccc       >  2:22001/1‑151 (MQ=255)                                                                                                                               ttGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGGTTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATc   >  1:545634/1‑151 (MQ=255)                                                                                                                                 gTAACAAAACCCCCACTTCGGCCGAAAATTTTGAGAAGTTTTTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCa  <  1:81656/150‑1 (MQ=255)                                                                                                                                 gTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATc   >  1:24373/1‑149 (MQ=255)                                                                                                                                   aaaaaaaaCCCCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCaa     <  2:763475/142‑1 (MQ=255)                                                                                                                                            |                                                                                                                                           CGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCC‑AGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCA  >  NC_000913/657538‑657814

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