Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	775,709	0	T	C	41.5%	51.9 / 11.8	41	I186T (ATC→ACC)	tolQ	membrane spanning protein in TolA‑TolQ‑TolR complex
Reads supporting (aligned to +/- strand):  ref base T (18/6);  new base C (0/17);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 6.89e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTG‑CCGCTATCCCGGCAGTTATGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGAAGAGTTTACCGCGATTCTGCACCGCCAGGCGTTTACCGTTAGCGAGAGCAACAAGGGGTAAGCCATGGCC  >  NC_000913/775560‑775854                                                                                                                                                       |                                                                                                                                                   cAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTG‑CCGCTATc                                                                                                                                                   >  1:136883/1‑151 (MQ=255)    cccGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTG‑CCGCTATCCCg                                                                                                                                                <  2:104077/151‑1 (MQ=255)      cGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTG‑CCGCTATCCCGGc                                                                                                                                              >  2:138777/1‑151 (MQ=255)                   gTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTG‑CCGCTATCCCGGCAGTTATGGCC‑TAc                                                                                                                                >  1:84597/1‑151 (MQ=255)                       ggTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTG‑CCGCTATCCCGGCAGTTATGGCC‑TACAAcc                                                                                                                            <  2:14293/151‑1 (MQ=255)                           cGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTG‑CCGCTATCCCGGCAGTTATGGCC‑TACAACCGCCt                                                                                                                        <  2:545270/151‑1 (MQ=255)                                 gggggATCATGCACGCTTTTACCGCCCTGGGGGGGTAAAAACAAGCACCACTGCAAAGGGTTCCCCCGGGTACCGAAGAGGCGTTAATGGCGCCTCAATTTGGTTTGTTGG‑CCTCTACCCCGGCAGTTATGGCC‑TACAACCGCCTCAACCa                                                                                                                  <  1:167687/151‑1 (MQ=255)                                 gggggATCATCCCCGCTTTTATCGCCCGCGGGGCGGTAAAACAGGCAACCCTGCAAATGGTTGCCCCCGGTACCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTTTTGG‑CCGCTACCCCGGCAGTTATGGCC‑TACAACCGCCTCAAcc                                                                                                                   <  1:231456/150‑1 (MQ=255)                                        cATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTG‑CCGCTATCCCGGCAGTTAAGGCC‑TACAACCGCCTCAACCAGcgcg                                                                                                             >  1:174489/1‑149 (MQ=255)                                             ccGCTTTTATCCCCCGGGGGGGGGAAAAACAAGAAACTCGGCAAAGGTTGGCCCCCGTTCTCGCAGAAGGGTTAATCGCGCCCGAAATTGGTCTTTTGC‑CCGCTACCCCGGCAGTTAGGCCC‑TACAACCGCCTCAACCAgcgc                                                                                                              <  1:267740/142‑1 (MQ=255)                                                  tttATCGCCCGGGGGGGGGAAAAAAAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACGCCATTGGGTCTGTTGG‑CCGCTACCCCGGCAGTTAGGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTg                                                                                                  <  1:5325/150‑1 (MQ=255)                                                  tttACCGCCCTCGGGGCGGAAAAACAAGAAACACTGAAAATGGTTGCGCCCGGTACCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTGG‑CCGCTACCCCGGCAGTTATGGCC‑TACAACCGCCTCAACCAGCGCGTaaac                                                                                                        <  2:487273/144‑1 (MQ=255)                                                  tttACCCCCGTCGGGGCGGAAAAACAACCAACCCGGCAAGTGGTGGCCCCCGGTTCGGCAGAGGCGTTGATTGCGACGGAATTTGTTTTGTTTC‑CCGTAACCCCGCTAGTTATGGCC‑TAAACCCGCCTCAACCAGCGCGTAAACAAACt                                                                                                   <  1:509355/149‑1 (MQ=255)                                                       cGCCCTCGGGGCGGTAAAACCAGCAACACTGCCAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGGCTGCAATTGGTCTGTTTG‑CCGCTATCCCGGCAGTTATGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTg                                                                                            >  1:681201/1‑151 (MQ=255)                                                         cccTGGGGGCGGAAAAACAACCACCACTGAAAATGGTTGCCCCCGTTATCGAAGAAGCTTTGATTGCGACGGCATTTGTTCTGTTGG‑CCGCTACCCCGGCAGTTATGCCC‑TAAAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTg                                                                                            <  1:443560/149‑1 (MQ=255)                                                              ggggggGAAAAAAAAGCAACACTGCAATTGGTGGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGAATTTGGTCTGTTTG‑CCGCTACCCCGGCAGTTATGGCT‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTAc                                                                                      <  1:266295/150‑1 (MQ=255)                                                                ggCGGAAAAAAAAGAACCACTCCAAATGGTTGCCCCCGGTACCGAAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTC‑CCGCTACCCCGGCAGTTATGCCC‑TAAAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGAc                                                                                   <  2:95391/151‑1 (MQ=255)                                                                      aaaaaaAGCACCCCTGCAAATGGTTGCGCCCGTTATCGAAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTG‑CCGCTACCCCGGCAGTTATGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACttt                                                                             <  2:354578/151‑1 (MQ=255)                                                                        aaaaaccAACACTGAAAAGGGTGGCCCCCGGTATCGAAGAGGCGTTGATTGCGACTGAAATTGGTCTGTTTC‑CCGCTACCCCGGAAGTTAGGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTAt                                                                           <  1:532124/145‑1 (MQ=255)                                                                         aCAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTG‑CCGCTATCCCGGCAGTTATGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATg                                                                          >  1:303120/1‑151 (MQ=255)                                                                            accAACACGGCAAATGGTTGCGCCCGGTACCGCAGAAGCGTTGATTGCGCCTGCATTTGGTCTGTTTG‑CCGCTACCCCGGCAGTTATGGCC‑TAAAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGaa                                                                       <  2:84617/149‑1 (MQ=255)                                                                              cAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTG‑CCGCTATCCCGGCAGTTATGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGa                                                                        >  2:450995/1‑148 (MQ=255)                                                                                   cTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTG‑CCGCTATCCCGGCAGTTATGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGAAGAGttt                                                                 >  1:188894/1‑150 (MQ=255)                                                                                    tGCAAATGGTTGCGCCCGGTATCGCAGAGGCTTTGATTGCGACTGCATTTGGTCTGTTTG‑CCGCTATCCCGGCAGTTATGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGAAGAGttt                                                                 <  1:377535/149‑1 (MQ=255)                                                                                       aattGGTTGCGCCCGGTATCAAAGAGGGGTTGATTGCGACTGC‑ATTGGTTTTTTTGCCCGCTACCCGGGCAGTTATGCCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGAAGAGTTTAcc                                                              <  1:586669/146‑1 (MQ=255)                                                                                       aaaTGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTG‑CCGCTATCCCGGCAGTTATGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGAAGAGTTTACcgc                                                            >  1:128131/1‑151 (MQ=255)                                                                                        aaTGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTG‑CCGCTATCCCGGCAGTTATGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGAAGAGTTTACcgcg                                                           >  2:518221/1‑151 (MQ=255)                                                                                         aTGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTG‑CCGCTATCCCGGCAGTTATGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGAAGAGTTTACcgcg                                                           >  2:500190/1‑150 (MQ=255)                                                                                                              gaggCGTTGTTTGCGACTACATTTGTTCTGTTTG‑CCGCTCTCCCGGCA‑TTAAGGCCTTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGAAGAGTTTACCGCGATTCTGCACCGCCAGGCGttt                                      <  1:494637/147‑1 (MQ=255)                                                                                                                   gTTGATTGCGACTGCAATTGGTCTGTTTG‑CCGCTATCCCGGCAGTTATGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGAAGAGTTTACCGCGATTCTGCACCGCCAGGCGTTTACCGtt                                >  2:169867/1‑151 (MQ=255)                                                                                                                             aCTGCAATTGGTCTGTTTG‑CCGCTATCCCGGCAGTTATGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGAAGAGTTTACCGCGATTCTGCACCGCCAGGCGTTTACCGTTAGCGAGAGc                       >  2:443/1‑150 (MQ=255)                                                                                                                                      gttCTGTTTG‑CCCCCTCCCCGGCATTAGTGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGAAGAGTTTACCGCGATTCTGCACCGCCAGGCGTTTACCGTTAGCGAGAGCAACAAgggg              <  2:72772/148‑1 (MQ=255)                                                                                                                                          tGTTTG‑CCGCTATCCCGGCAGTTATGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGAAGAGTTTACCGCGATTCTGCACCGCCAGGCGTTTACCGTTAGCGAGAGCAACAAGGGGTAAGc         >  2:641355/1‑151 (MQ=255)                                                                                                                                            tttG‑CCGCTATCCCGGCAGTTATGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGAAGAGTTTACCGCGATTCTGCACCGCCAGGCGTTTACCGTTAGCGAGAGCAACAAGGGGTAAGcc        >  2:34234/1‑150 (MQ=255)                                                                                                                                            tttG‑CCGCTATCCCGGCAGTTATGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGAAGAGTTTACCGCGATTCTGCACCGCCAGGCGTTTACCGTTAGCGAGAGCAACAAGGGGTAAGcc        >  2:366422/1‑150 (MQ=255)                                                                                                                                            tttG‑CCGCTACCCCGGCAGTTAGGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGAAGAGTTTACCGCGATTCTGCACCGCCAGGCGTTTACCGTTAGCGAGAGCAACAAGGGGTaa           <  1:343126/147‑1 (MQ=255)                                                                                                                                                 ccGTTATCCGGGCATTTATGCCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGAAGAGTTTACCGCGATTCTGCACCGCCAGGCGTTTACCGTTAGCGAGAGCAACAAGGGGTAAGCCATgg    <  2:298422/150‑1 (MQ=255)                                                                                                                                                 ccGCTACCCCGGCAGTTATGGCC‑TAAAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGAAGAGTTTACCGCGATTCTGCACCGCCAGGCGTTTACCGTTAGCGAGAGCAACAAGGGGTAAGCCATGGc   <  2:307733/151‑1 (MQ=255)                                                                                                                                                 ccGCTACCCCGCCATTTATGCCC‑TAAAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGAAGAGTTTACCGCGATTCTGCACCGCCAGGCGTTTACCGTTAGCGAGAGCAACAAGGGGTAAGCCATg     <  2:153409/149‑1 (MQ=255)                                                                                                                                                   gCTATCCCGGCAGTTATGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGAAGAGTTTACCGCGATTCTGCACCGCCAGGCGTTTACCGTTAGCGAGAGCAACAAGGGGTAAGCCATGGcc  >  1:634872/1‑150 (MQ=255)                                                                                                                                                     tATCCCGGCAGTTATGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGAAGAGTTTACCGCGATTCTGCACCGCCAGGCGTTTACCGTTAGCGAGAGCAACAAGGGGTAAGCCATGGcc  >  1:182759/1‑148 (MQ=255)                                                                                                                                                       |                                                                                                                                                   CAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTG‑CCGCTATCCCGGCAGTTATGGCC‑TACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAACTTTATGGAAGAGTTTACCGCGATTCTGCACCGCCAGGCGTTTACCGTTAGCGAGAGCAACAAGGGGTAAGCCATGGCC  >  NC_000913/775560‑775854

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