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

AAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTC  >  NC_000913/3998029‑3998324                                                                                                                                                        |                                                                                                                                                    aaGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTACCAACAAAgc                                                                                                                                                    <  2:501973/151‑1 (MQ=255)  aGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTACCAACAAAgcg                                                                                                                                                   <  2:324418/151‑1 (MQ=255)   gCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTACCAACAAAgcg                                                                                                                                                   <  2:243065/150‑1 (MQ=255)        ggCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTACCAACAAAgcggcggc                                                                                                                                              <  1:405645/150‑1 (MQ=255)                   gcagcaACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCa                                                                                                                                  <  2:49898/151‑1 (MQ=255)                          ccTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTA‑Tc                                                                                                                          <  1:673848/151‑1 (MQ=255)                             tCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTA‑TCg                                                                                                                         <  2:58738/149‑1 (MQ=255)                               tGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGGGCGTACTGGTGCAGCGTATGGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTACCAACAAAGCGGCGGCGGAGATGTGTTTTCGTA‑Tc                                                                                                                          <  2:228858/146‑1 (MQ=255)                                                ggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTA‑TCGGGCAACTGATGGGCACGAGcc                                                                                                    <  2:314012/151‑1 (MQ=255)                                                   aGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AA‑GGTTTACAAAAAAACCGCGGCGGGGATGGGGCATCGTC‑TTGGGCCACAGGTGGGGACg                                                                                                        >  1:329549/1‑144 (MQ=255)                                                            aCGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTTCCAACAAAACGGCGGCGGGGAGGCGGCATCGTA‑TCGGGCGAAAGAAGGGCGCGAGCCAGGGGGGGAtg                                                                                         >  1:49903/1‑150 (MQ=255)                                                              gcgcgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTgg                                                                                      <  1:322057/151‑1 (MQ=255)                                                                 cgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTA‑TCGGGCAAATGATGGGCACGAGCCAGGGCGGTATGTGGGTc                                                                                   <  1:87619/151‑1 (MQ=255)                                                                                    gCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTgg                                                                <  2:724015/151‑1 (MQ=255)                                                                                     ccTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGgc                                                               <  1:541284/151‑1 (MQ=255)                                                                                      cTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GGTTTCCCACAAAACCGCGGCGGGGGGGGGTCCTCCGA‑TCGGGCGACCGATGGGGACCAACCCGGGGGGGGTTTGGGGCGGCCCCCTCCCCCCGCTgg                                                                >  2:491787/1‑149 (MQ=255)                                                                                        ggTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGgcg                                                              <  1:388507/149‑1 (MQ=255)                                                                                        ggTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GGTTTCCAACAAAACGGCGGGGGGGGGGGGTCCTCGTA‑TCGGGGCACCGATGGGGCCGAGCCAGGGCGGGATGTGGGGGGGGGCCCTCCCCCGGCTGgcg                                                              >  1:16442/1‑149 (MQ=255)                                                                                          ttGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGgcg                                                              <  2:389996/147‑1 (MQ=255)                                                                                          ttGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GGTTTCCAACAAAACGGCGGCGGGGAGGCGTCATCCTA‑TCGGGGAACTGATGGGGACGAGCCCGGGGGGGGTGTGGGGGGGGACCCTCCCCCGGGTGgcggg                                                            >  1:390037/1‑147 (MQ=255)                                                                                           tgatgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTTCCAAAAAAAAGGCGGCGGGGGTGGGTCATCGTA‑TTTGGCAAACGAAGGGCACGAGCCCGGGCGGGGTGTGGGGCGGCACCTTCCCCGGGCTGGCGccccc                                                         >  1:243272/1‑150 (MQ=255)                                                                                             atgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAAG‑GTTTACC‑ACAAAACAGCGGCGGGGAGGCGTCATCGTATTCGGGC‑ACCGAAGGGCACGAGCCAAGGGGGGAGGGGGGGCGGCCCCCTTCCCCGGCGGGCGCAcccct                                                       >  1:73991/1‑148 (MQ=255)                                                                                             atgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑ACGGTTTTCCAA‑AAAACGGGGGCGGGGGGGCGTCTTCACA‑TCGGGGAACTGATGGGCCCCAGGCACGGGGGTATGTGGGGCGGGACCTTCCCCCGGCGCGCGCACCGt                                                        >  1:220646/1‑150 (MQ=255)                                                                                               gaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtg                                                     <  2:26523/151‑1 (MQ=255)                                                                                                   gTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtgctgc                                                 <  2:315878/151‑1 (MQ=255)                                                                                                     ggAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCtgcg                                                <  1:160590/150‑1 (MQ=255)                                                                                                       aaaaCTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTTCCAACAAAACGGCGGCGGGGAAGCGGCATCGTA‑TCGGGGAAATGATGGGGACGAGCCAGGGCGGGATGTGGGGCGGCCCCCTCCCCCGGCCGGCGCACCGTTTTCtgcgtg                                              >  1:149937/1‑150 (MQ=255)                                                                                                                  ccATACTCGATTATGGCGGTG‑AC‑GGTTTCCCACAAAACGGCGGCGGGGGAGGGTCATCGTA‑TCGGGGAACTGATGGGGACGAGCCAGGGCGGGGTGTGGGGGGGGACCCTCCACCGGGCGGCGCACCGGTTTCTTCGTGCGCACCcta                                     >  1:649755/1‑148 (MQ=255)                                                                                                                      aCTCGATTATGGCGGTG‑AC‑GGTTACCAACAAAAAGGCGGCGGAGAGGCGTTTTCGGA‑TCGGGGGAAAGGTGGGGACGAGCCAGGGGGGGGTGTGGGGGGGCGCCCTCCACGGGCGGGCGCACCCTGTGCTGCGGGCGCCCCatattg                                  >  2:51777/1‑145 (MQ=255)                                                                                                                            ttATGGCGGTG‑AC‑GTTTACCAACAAAACGGCGGCGGGGATGCGTCATCGTA‑TTGGGCGACTGATGGGCACGAGCCAGGGCGGGATGTGGGGCGGCACCCTCCCCCGGCCGGGGCACCGGTTGCTGCGTGGGCACCACATg                                   >  2:58718/1‑140 (MQ=255)                                                                                                                               tGGCGGTG‑AC‑GTTTACCAACAAAACGGCGGCGGAGATGCGTCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGGATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGcc                     >  1:722146/1‑151 (MQ=255)                                                                                                                                                 caacaaAGCGGCGGCGGAGATGCGTCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATc     <  2:487001/151‑1 (MQ=255)                                                                                                                                                   acaaAGCGGCGGCGGAGATGCGTCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCt   <  1:614476/151‑1 (MQ=255)                                                                                                                                                    caaAGCGGCGGCGGAGATGCGTCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTc  >  1:590643/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                    AAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTG‑AC‑GTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTC  >  NC_000913/3998029‑3998324

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