Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,998,178	0	G	A	41.0%	43.3 / 21.5	39	A66T (GCG→ACG)	uvrD	DNA‑dependent ATPase I and helicase II
Reads supporting (aligned to +/- strand):  ref base G (7/16);  new base A (16/0);  total (23/16)
Fisher's exact test for biased strand distribution p-value = 7.32e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACC‑TTCCACGGGCTGGC‑GCACC‑GTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACA  >  NC_000913/3998029‑3998328                                                                                                                                                       |                                                                                                                                                           aaGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC‑GTTTACCAACAAAgc                                                                                                                                                           <  1:389330/151‑1 (MQ=255)           cgcgCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCCGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGTCGGTTTCCAAACAAAgggccg                                                                                                                                                       >  1:187356/1‑141 (MQ=255)             cgcCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC‑GTTTACCAACAAAGCGGCGGCGGAGAt                                                                                                                                               <  1:102438/151‑1 (MQ=255)              gcCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGCGGCGGTTTTCAAACAAAAC‑GCCGCGGGGATg                                                                                                                                              >  2:832742/1‑151 (MQ=255)                    cagcaACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGCGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCAt                                                                                                                                        <  1:703550/151‑1 (MQ=255)                        aaCCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTa                                                                                                                                    <  2:129846/151‑1 (MQ=255)                             tCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCg                                                                                                                                 <  2:523164/149‑1 (MQ=255)                                 gTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAAc                                                                                                                           <  1:274960/151‑1 (MQ=255)                                    cTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGAC‑GTTTTTCAAAAAAAACGCGGGGGGGGTGCGGCCTCCTTTCggg                                                                                                                               >  1:462393/1‑144 (MQ=255)                                    cTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGAC‑GTTTACCAAAAAAAAGGGGGCGGGGAGGGGGCCTCGGATCGGGCaaac                                                                                                                          >  1:802541/1‑147 (MQ=255)                                          ggcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAA‑GGTTACCAACAAAACGGCGGCGGGGGTGCGGCATCGTATCGGGCAAACGATggg                                                                                                                    >  1:932797/1‑149 (MQ=255)                                           gcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCAc                                                                                                                 <  1:703353/151‑1 (MQ=255)                                               gggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGGC‑GGTTTCCCACAAAACGGCGGCGGGGAAGCGTCATCGGATCGGGCaaacgaaggg                                                                                                                    >  1:936107/1‑136 (MQ=255)                                                ggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGcc                                                                                                            <  2:659784/151‑1 (MQ=255)                                                         aaGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTa                                                                                                   <  1:951582/151‑1 (MQ=255)                                                              gcgcgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC‑GGTTTCCAACAAAACGGCGGCGGCGATGCGTCATCGGATCGGGGAACCCGTGGGGAACAGCCCGGGCGGGATgggg                                                                                              >  1:942043/1‑151 (MQ=255)                                                                         tGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAA‑GCTTTCCAACAAAGAGGCGGCGGAGATTAGTCCTCGGACCGGGCAACAGATGGGCACGAACCAGGGGGGGATGTTGGt                                                                                            >  1:114176/1‑142 (MQ=255)                                                                                   cGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCACACTCGATTATGGCGGTGAC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACC‑TTCCACGGGCTg                                                                       <  2:784335/151‑1 (MQ=255)                                                                                     ccTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACC‑TTCCACGGGCTGgc                                                                     <  2:50555/151‑1 (MQ=255)                                                                                        ggTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGAC‑GGTTACCAAAAAAAAGGCGGCGGGGGTGCGGCCTCGTATCGGGCAAATGATGGGGGCGAGCCAGGGGGGGGTGTGGGGG‑GGCACC‑CTCCCCGGGGTGgc‑g                                                                   >  2:537863/1‑149 (MQ=255)                                                                                         gTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC‑GGTTTCCCACAAAACGGGGGCGGGGATGGGGCCTCGGATCGGGGCAATGATTGGCACGAGACAGGGCGGGAAGGGGGGC‑GGGACC‑CTCCACcgg                                                                          >  1:648388/1‑142 (MQ=255)                                                                                          ttGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAA‑GTTTTCCCACAAAAAGGGGGGGGGGGGGGGGGACCGGGTGGGGGAACCGAAGGGGACGAGGCAGGGGGGGAAGTGGGGCGGGCACC‑TCCCCCggcgg                                                                        >  2:929398/1‑141 (MQ=255)                                                                                           tgatgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACC‑TTCCACGGGCTGGC‑GCAcc                                                               <  2:914024/150‑1 (MQ=255)                                                                                               gaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACC‑TTCCACGGGCTGGC‑GCACC‑Gttt                                                          <  2:526297/150‑1 (MQ=255)                                                                                                aGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGAC‑GTTTTACCACAAAACGGCCGCGGCGGTGCGTCACCCTATAGGGGAAACGCTGGGGAACAGCCCGGGCGGGATGTGGGGT‑GGCGCC‑CTCCCCCGGCTGgc‑gccc                                                                >  2:595853/1‑142 (MQ=255)                                                                                                    tGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC‑GTTTATCAACAAAACGGCGGGGGGGATGCGTCATCGTATCGGGGAAAAGATGGGGACGAGCCCGGGGGGGATGTGGGGC‑GGCACCTTTCCCCGGGC‑GGC‑GCACC‑Cttt                                                          >  1:643331/1‑145 (MQ=255)                                                                                                     ggAAAACTGCTCGCCATACTCGATTATGGCGGTGAC‑GTTTTCCAACAAAACGGCGGCGGGGAAGCGTCCTCGTATCGGGCAACTGATGGGGACGAGCCAGGGGGGGATGTGGGGC‑GGCACC‑TTCCCCCGGCGGGC‑GCCCC‑GTTttct                                                       >  1:441919/1‑147 (MQ=255)                                                                                                     ggAAAACTGCTCGCCATACTCGATTATGGCGGTGAC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACC‑TTCCACGGGCTGGC‑GCACC‑GTTTGCtgcgt                                                   <  1:460922/151‑1 (MQ=255)                                                                                                              cTCGCCATACTCGATTATGGCGGTGAC‑GTTTACCAACAAAAAGGCGGCGGGGAAGCGGCATCGTATCGGGCAACCGATGGGCACGAGCCCAGGGGGGATGTGGGGC‑GGGCCC‑CTCCACCGGCTGGG‑GCACCGGTTTGC‑GCGTGCGcccc                                            >  1:400424/1‑149 (MQ=255)                                                                                                                  ccATACTCGATTATGGCGGTGAC‑GTTTTCCCACAAAAAGGCGGCGGAGAGGCGTCATCGTTTCGGGGAACCGATGGGGACGGGCCCGGGCGGGGTGTGGGGG‑GGGACC‑CTCCCCGGGCGGGG‑GCCCC‑GGTTGTTGCGTGGGCGcc                                            >  2:770222/1‑145 (MQ=255)                                                                                                                  ccATACTCGATTATGGCGGTGAC‑GGTTTCCCACAAAACGGCGGCGGAGAAGCGGCATCCTATCCGGGCACCGATTGGCGCGGACCCGGGCGGGGTTTGGGGT‑GGGACC‑CTCCACCGGCGGGC‑GCCCC‑CGTTTCTGCGCGGGGCCCatat                                        >  2:471938/1‑149 (MQ=255)                                                                                                                        tCGATTATGGCGGTGAC‑GTTTTCCCACAAAACCGCGGCGGGGATGGGTCATTGGATCGGGGCACCGATTGGGACCAGCCAGGGGGGGGTGTGGGGC‑GGGCCC‑CTCCCCCGGCCGGCGGCCCC‑GTTT‑CTGCGGGCGCACCCCATTGACGcc                                 >  2:487746/1‑150 (MQ=255)                                                                                                                            ttATGGCGGTGAC‑GTTTACCAACAAAGCGGCGGCGGGGATGCGTCCTCTTATCGGGGAACTGATGGGGACGAGCCCGGGGGGGGTGTGGGTC‑GGCACC‑CTTCCCCGGCTGGG‑GCACC‑GTTTGCTTCGGGGGCCCCCTATGGGCGCCaa                               >  2:582211/1‑148 (MQ=255)                                                                                                                               tGGCGGTGAC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACC‑TTCCACGGGCTGGC‑GCACC‑GTTTGCTGCGTGCGCACCATATGGACGCCAATCTGcc                         >  2:208928/1‑151 (MQ=255)                                                                                                                                      gAC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACC‑TTCCACGGGCTGGC‑GCACC‑GTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGAt                  >  2:55779/1‑151 (MQ=255)                                                                                                                                          gTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACC‑TTCCACGGGCTGGC‑GCACC‑GTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTc               <  1:346858/151‑1 (MQ=255)                                                                                                                                            ttACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACC‑TTCCACGGGCTGGC‑GCACC‑GTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTcc              >  2:426278/1‑150 (MQ=255)                                                                                                                                                 aacaaAGCGGCGTCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCAGGAGCCAGGGCGGTATGTGCGTC‑GGCACC‑TTCCACCGGCTGGC‑GCACC‑GTTTGCTGCGTTCGCACCATATGGACGCCAAACTGCCGCACGATTTCCAGa           <  2:702550/148‑1 (MQ=255)                                                                                                                                                 aacaaAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACC‑TTCCACGGGCTGGC‑GCACC‑GTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATcc        >  2:780929/1‑151 (MQ=255)                                                                                                                                                  acaaAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACC‑TTCCACGGGCTGGC‑GCACC‑GTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATcc        >  2:241352/1‑150 (MQ=255)                                                                                                                                                    aaaGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACC‑TTCCACGGGCTGGC‑GCACC‑GTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCg     <  1:105158/151‑1 (MQ=255)                                                                                                                                                       gcggcggcgGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACC‑TTCCACGGGCTGGC‑GCACC‑GTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACa  <  1:518107/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                           AAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACC‑TTCCACGGGCTGGC‑GCACC‑GTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACA  >  NC_000913/3998029‑3998328

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