Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,331,790	0	T	G	43.5%	79.4 / 11.3	46	V430G (GTT→GGT)	proP	proline/glycine betaine transporter
Reads supporting (aligned to +/- strand):  ref base T (3/23);  new base G (20/0);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 6.32e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACCGGCGGCGTCAGATATACA‑GGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCA  >  NC_000913/4331641‑4331937                                                                                                                                                      |                                                                                                                                                        cccGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGtggtt                                                                                                                                                        <  2:343343/151‑1 (MQ=255)     aCGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTggttggtt                                                                                                                                                    <  1:190722/151‑1 (MQ=255)     aCGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTggttggtt                                                                                                                                                    <  2:595667/151‑1 (MQ=255)        cATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTaa                                                                                                                                                 <  1:112556/151‑1 (MQ=255)          tatCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGGGGGTGGGTTAAt                                                                                                                                                >  2:779511/1‑150 (MQ=255)              cGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGCGGTTGGTTTAATCAcc                                                                                                                                            >  1:488709/1‑150 (MQ=255)               gTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGGTGGGTTTATctac                                                                                                                                            >  1:840771/1‑146 (MQ=255)                  aCAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGGGGGTGGGTTTATTACCcgcgg                                                                                                                                       >  2:320748/1‑150 (MQ=255)                      cgcgcTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGGTGGGTTTATCACCCGCGG‑AAcc                                                                                                                                  >  1:410813/1‑151 (MQ=255)                           tGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGGTGGTTTTATTACCggggcgaacacgt                                                                                                                               >  1:846641/1‑139 (MQ=255)                            ggcggcggcATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGaaa                                                                                                                            <  1:425880/151‑1 (MQ=255)                                gcggcATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGGTGGGTTTATTACCCGCGG‑AAACATGaaa                                                                                                                            >  1:276759/1‑147 (MQ=255)                                 cggcATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGGGGGTGGGTTTATTACCCGCCG‑AAcccataaa                                                                                                                            >  2:821117/1‑140 (MQ=255)                                       ttAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGGTGGGTTAATCACCCGCGG‑AAACATGaaaaa                                                                                                                          >  1:538944/1‑140 (MQ=255)                                             tttCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGGTGGGTTTATTAACCGCGG‑AAACATGAAAGAG‑AAGGCaaa                                                                                                                >  1:166043/1‑145 (MQ=255)                                                 ggTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGGTTTATCACCGGCGG‑TACCATTGAAAAG‑GCGGCaaaac                                                                                                              >  1:687883/1‑141 (MQ=255)                                                 ggTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGGTGGTTTTATTACCCGCGG‑AACCCTGAAAAAG‑AAGGCAAAAAGTcc                                                                                                          >  2:655076/1‑147 (MQ=255)                                                          tGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGc                                                                                             <  2:21610/150‑1 (MQ=255)                                                           gCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGa                                                                                           <  2:505963/151‑1 (MQ=255)                                                            ccGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGGGGGTGGGGTAATCACCGGCGA‑AAACAAAAAAAAGAAGGGCAAAACGTCCG‑tgaacaggggg                                                                                             >  2:405591/1‑139 (MQ=255)                                                                   gACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACcggc                                                                                    <  1:460444/150‑1 (MQ=255)                                                                               ggcggcCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACCGGCGGCGTCAGatat                                                                        <  2:951018/150‑1 (MQ=255)                                                                                    ccTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGGTTTATCACCCGCGT‑AAACCTGAAAGAG‑ACGGGAAAACCTCCG‑TTGGAAGGGGGGGCACCCGCGGGGGCCGAtataaaag                                                                   >  1:460654/1‑146 (MQ=255)                                                                                     ctggctggTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGGTGGTTTTATAACCCGCGG‑TACCATGAAAAAG‑AAGGCAAAACGCTCC‑TTGTAAGGGGGGGCAcc                                                                                       >  2:681104/1‑129 (MQ=255)                                                                                       ggctggTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACCGGCGGCGTCAGATATACA‑GGaa                                                                <  2:644576/149‑1 (MQ=255)                                                                                         ctggTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGGTGGGTTAATTACCGGCGG‑AACCATGAAAAAG‑AAGGGAAAACGTCCG‑TTGAAAAGGGGGACCCCGGCGGGGGCAGATATAAA‑AGAAGCg                                                             >  1:201547/1‑150 (MQ=255)                                                                                              cGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACCGGCGGCGTCAGATATACA‑GGAAGCGAAGGa                                                        <  1:790158/150‑1 (MQ=255)                                                                                               gAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGAGGGTGGGTTTATTAACGGCGT‑AAACAAGAGAAAA‑AAGGGAAAACGTGCG‑TTGGAAAGGGGGGCACCCGCGGCGGGAGata                                                                         >  2:409926/1‑133 (MQ=255)                                                                                                  aGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGGTTGTTTAATAACCCGCGG‑AAACCTGAAAAAA‑AAGGCAAATCCTCCG‑TTGAAAGGGGGGACACCCGCGGCGTCAGATAAAAA‑AGAAACGaaaggaa                                                      >  1:5302/1‑143 (MQ=255)                                                                                                      cGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACCGGCGGCGTCAGATATACA‑GGAAGCGAAGGAAATtctc                                                 <  1:447454/149‑1 (MQ=255)                                                                                                      cGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACCGGCGGCGTCAGATATACA‑GGAAGCGAAGGAAATTCtcgt                                               <  1:629649/151‑1 (MQ=255)                                                                                                      cGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACCGGCGGCGTCAGATATACA‑GGAAGCGAAGGAAATTCtcg                                                <  2:646960/150‑1 (MQ=255)                                                                                                         aGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACCGGCGGCGTCAGATATACA‑GGAAGCGAAGGAAATTCTCGTCGa                                            <  1:490608/151‑1 (MQ=255)                                                                                                           aaTCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGGTGGGTTTATCACCCGCGG‑GACCATGAAAGAG‑ACGGCAAAACGGCCG‑CTGGAAGGGGGGGAACCCGCGGCGCCCGAAATAAA‑AGAAGCGGAGGAAAATTtcctcgcg                                           >  1:805561/1‑148 (MQ=255)                                                                                                           aaTCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGGGGGGTTTATAACCCCCGG‑AAACCCGAAAGAA‑AGGGGGAAAAGGCCC‑CTGAGAGGAGGGGCACCGGCGGCGGCCGAGATACA‑GGAAGGgaaagaaa                                                      >  1:307302/1‑139 (MQ=255)                                                                                                            aTCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGGTGGGGTTATTACCCCCGTAAACCAT‑AAAAAA‑ACGGGGAAAAGGCCG‑TTGTAAAGGGGGGCACCCGCGGGGGGAGAGATAAA‑AGAAGGGGAGGAAAATAttgtcgcg                                           >  2:64072/1‑147 (MQ=255)                                                                                                             tCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGGTGGGTTTATTACCGGCGG‑AAACAAGAAAAAG‑ACGGCCAAAAGCCCGTTTGAAAGG‑GCGAGACCCGCGGCGGCAGATATAAAGGGAACGGAGGGAAAttttgtcctg                                            >  1:764218/1‑140 (MQ=255)                                                                                                                gatgatGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACCGGCGGCGTCAGATATACA‑GGAAGCGAAGGAAATTCTCGTCGAGCAtt                                       <  2:401503/149‑1 (MQ=255)                                                                                                                gatgatGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACCGGCGGCGTCAGATATACA‑GGAAGCGAAGGAAATTCTCGTCGAGCATTa                                      <  1:411470/150‑1 (MQ=255)                                                                                                                gatgatGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACCGGCGGCGTCAGATATACA‑GGAAGCGAAGGAAATTCTCGTCGAGCATTAc                                     <  2:276753/151‑1 (MQ=255)                                                                                                                gatgatGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACCGGCGGCGTCAGATATACA‑GGAAGCGAAGGAAATTCTCGTCGAGCATTAc                                     <  1:214524/151‑1 (MQ=255)                                                                                                                gatgatGCCTGCCTATTACCTGATGGTAGTGGCGGTGGGTGGGTTTATTACCCGCGG‑AACCATTAAAAAA‑ACGGCAAAACGTCCGTTTGAAAGG‑GCGGCACCCGCGGCGGCCGATATAAA‑AGAAAGGAAGGAAAATTTCGTCGGGCCtt                                       >  1:421904/1‑149 (MQ=255)                                                                                                                          gcctATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACCGGCGGCGTCAGATATACA‑GGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAatat                             <  2:225274/149‑1 (MQ=255)                                                                                                                          gcctATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACCGGCGGCGTCAGATATACA‑GGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATc                            <  1:277234/150‑1 (MQ=255)                                                                                                                          gcctATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACCGGCGGCGTCAGATATACA‑GGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCg                           <  1:780705/151‑1 (MQ=255)                                                                                                                              aTTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACCGGCGGCGTCAGATATACA‑GGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCa                       <  1:98827/151‑1 (MQ=255)                                                                                                                                           aGTGGCGATGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACCGGCGGCGTCAGATATACA‑GGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCgatga           <  1:94873/150‑1 (MQ=255)                                                                                                                                              ggCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACCGGCGGCGTCAGATATACA‑GGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGatat        <  1:311964/150‑1 (MQ=255)                                                                                                                                               gCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACCGGCGGCGTCAGATATACA‑GGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTg      <  2:840751/151‑1 (MQ=255)                                                                                                                                                    ggttggttTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACCGGCGGCGTCAGATATACA‑GGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCa  <  1:862510/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                        CCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGT‑AACCATGAAAGAG‑ACGGCAAATCGTCCG‑TTGAAAGGTGCGACACCGGCGGCGTCAGATATACA‑GGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCA  >  NC_000913/4331641‑4331937

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