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

TGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTG  >  NC_000913/657561‑657820                                                                                                                    |                                                                                                                                                 tgCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAg                                                                                                                <  1:366029/150‑1 (MQ=255) tgCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAg                                                                                                                <  2:379578/150‑1 (MQ=255) tgCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGACAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAg                                                                                                                <  2:358385/150‑1 (MQ=255)    ggTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCg                                                                                                            >  1:595474/1‑151 (MQ=255)                       cAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGCGa                                                                                         >  2:453915/1‑151 (MQ=255)                              ccAGTGCGGTCATGGCAAAAGACCCCAGAAGGTGGACAAAACCGTTCAGTAATGCCCAGCAAAGGCGGCCCTTTTGAAACAAAACCCCCACTTCGGCCAAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATAc                                                                                  <  1:430995/151‑1 (MQ=255)                              ccAGTGCGGTAAGGGCAAAAACCCCCAGAAGGTTAACAAAACCGTTCAGTAATGCCCGGCAAAAGCGGCCCTCTTGTAACAAAAACCCCACTTCTGCCGAAAATTTTAGGAGGGTTGTTAGCCCGCCACAAAATCCGGTGGTGATTAATAc                                                                                  <  2:79325/151‑1 (MQ=255)                                    ggctCATGGCAAAAGACCCAAAAAGTTAGACAAAAACGTCCATAAATGCCAAGCAAAAGCGGCCCTCTTGTAAAAAAAACCCCACTTCTGCCGAAATTTTTGAAAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTcc                                                                             <  1:27964/147‑1 (MQ=255)                                          tGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGCGGTGATTAATACTTTCCACACTgg                                                                      >  2:231961/1‑151 (MQ=255)                                             cAAAAGCCCCAAAAAGTTACACAAAAACTTTCAGTAACGCCCAGCAAAAGCGGCCCTTTTGTAAAAAAAACCCCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATc                                                                   <  2:205141/151‑1 (MQ=255)                                                   gccccAGAAGGTTGACAAAAACGTTCAGTAATGCCCAACCAAAGGGGCCCTTTTGTAACAAAAACCCCCCTTCTGCCAAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATg                                                               <  2:413351/148‑1 (MQ=255)                                                   cccccAGAAGGTTAACAAAACCGTACAGTAATGCCCAGCAAAAGCGGCCCTTTTGAAAAAAAAACACCACTTCGGCCAAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGt                                                              <  1:603268/149‑1 (MQ=255)                                                    ccccAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGtt                                                             >  1:499666/1‑150 (MQ=255)                                                          aaGGTTGACGAAAACGTTCAGAAAGGCCCAGCAAAACCGGCCCTCTTGAAAAAAAAACCCCACTTCGGCCAAAAATGTTGGGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTc                                                         <  2:460825/148‑1 (MQ=255)                                                              ttGACAAAACCGTTCAGTAATCCCCAGCCAAAGCGGCCCTCTTGAAACAAAACCCCCACTTCTCCCAAAAATTTTGAAAAGGTTTTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCctgc                                                  <  2:452746/151‑1 (MQ=255)                                                                            cAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAt                                    >  2:148598/1‑151 (MQ=255)                                                                            cAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAt                                    >  1:20276/1‑151 (MQ=255)                                                                               tattGCCCAGCAAAAGCGGCCCTTTTGAAAAAAAAACCCCCCTTCTGCCAAAAATGTGGAGAAGGTTTTTAGGCCGCCGCAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcc                                  <  2:574943/147‑1 (MQ=255)                                                                                  tGCCCAGCCAAGGCGGCCCTCTTGTAACAAAACCACCACTTCTGCCGAAAATGTTGAGAAGTTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcctatt                              <  1:155355/151‑1 (MQ=255)                                                                                   ccccAGCAAAAGGGGCCCTCTTGAAACAAAACCCCCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcctattc                             <  1:157225/150‑1 (MQ=255)                                                                                    cccAGCCAAAGCGGCCCTCTTGTAACAAAACCCCCACTTCGGCCGAAAATGTTGAGAAGGTGGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcctattcc                            <  1:372071/151‑1 (MQ=255)                                                                                    cccAGCAAAAGCGGCCCTCTTGTAACAAAAACACCACTCCGGCCGAAATTGTTGAGAAGTTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcctattcc                            <  1:161861/151‑1 (MQ=255)                                                                                        gCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTatg                        >  2:206922/1‑151 (MQ=255)                                                                                         caaaaGCGGCCCTTTTGAAAAAAAAACCCCACTTCTGCCGAAATTGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcctattccta                          <  1:340237/146‑1 (MQ=255)                                                                                               cGGCCTTCTTGAAAAAAAACCCCCCCTTCTGCCGAAAATTTGGAGAGGGTTTTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAAt                  <  1:392163/150‑1 (MQ=255)                                                                                                   ccTTTTGAAAAAAAACCCCCTCTTCGGCCAAAAATTTTGAGAAGTTTGTTAGCCCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGcccc             <  2:269585/151‑1 (MQ=255)                                                                                                       ttGTAAAAAAACCCCCACTTCTGCCGAAAATGTTGAAAAGGTTGTAAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATc         <  2:55875/151‑1 (MQ=255)                                                                                                                   caccacTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTg  >  1:561536/1‑146 (MQ=255)                                                                                                                    |                                                                                                                                                 TGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTG  >  NC_000913/657561‑657820

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