Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,312,215	0	T	A	46.0%	49.8 / 15.6	37	D11V (GAT→GTT)	yciA	acyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base T (3/17);  new base A (17/0);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 9.31e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGC‑TAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGG  >  NC_000913/1312065‑1312355                                                                                                                                                        |                                                                                                                                              ccTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGC‑TAAAGTACGTAAAACAAGAt                                                                                                                                               <  2:430372/151‑1 (MQ=255)       aCCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGC‑TAAAGTACGTAAAACAAGATCGCCCt                                                                                                                                         <  1:572133/151‑1 (MQ=255)       aCCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGC‑TAAAGTACGTAAAACAAGATCGCCCt                                                                                                                                         <  2:138547/151‑1 (MQ=255)                   aCTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGACGGC‑TTAAGTAAGCAAAACGAAAACTCCCTggggtgcggtgt                                                                                                                             >  2:185259/1‑143 (MQ=255)                         cgACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGATGGT‑TAAAAAAAGGAAAAAAACAAAGGCCCGAGGGGCGGTATGagtt                                                                                                                        >  1:47762/1‑150 (MQ=255)                             cGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGG‑TAAAATACGTAAAAAAAAAACGCCCTGAGGGGCggttgatcggga                                                                                                                      >  2:447539/1‑137 (MQ=255)                             cGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGC‑TAAAGTACGTAAAACAAGATCGCCCTGAGGGACGtt                                                                                                                               <  1:447524/139‑1 (MQ=255)                                  gCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGC‑TAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATgg                                                                                                              <  1:360621/151‑1 (MQ=255)                                    aaTTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGC‑TAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTa                                                                                                            <  1:570898/151‑1 (MQ=255)                                    aaTTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGC‑TAAAATACGTAAAAAAAAAACGCCCCTAGGGGCGGTATGTGTTGTAGACAAgg                                                                                                              >  2:222383/1‑149 (MQ=255)                                    aaTTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCAAGGT‑TAAAGTACGGAAAACAAAAAAGACCCGAGGGAAGGTATTtgttgttaaca                                                                                                                 >  2:522241/1‑141 (MQ=255)                                        tCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGC‑TAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTaaaa                                                                                                         <  1:528209/150‑1 (MQ=255)                                        tCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGC‑TAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTaaaa                                                                                                         <  2:371496/150‑1 (MQ=255)                                               ccAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGC‑TCAAAAACGTAAAACAAGAACGCCCTGCGGGACGTTATGTGTTGTTGACATGGTAAAAACGACtt                                                                                                  >  1:92627/1‑150 (MQ=255)                                               ccAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGC‑TAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTT‑a                                                                                                <  1:411522/151‑1 (MQ=255)                                                   aaTAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACTATTGGCATTGGTATCGGCGGGCATGGC‑TAAAAGACGTAAAAAAAAAACGCCCCGAGGGGGGGTTTGGGTTGTTGAAATGGTAAAAACCACTT‑aaaa                                                                                             >  2:466848/1‑150 (MQ=255)                                                       gcgcCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCAAGGC‑CAAAATAACGAAAACAAGAAATCCCTGAGGGACGGTACTGATTGGAGACAAGGAAAAAACGACTT‑aaaaaaa                                                                                          >  1:408083/1‑146 (MQ=255)                                                          ccgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCAAGGC‑TAAAGAAAGTAAAAAAAAAAAGCCCTTAGGGGCGGTATGTGTTGTAGAAAAGGGAAAAACGACTT‑AAAAAGaaaa                                                                                       >  2:279675/1‑149 (MQ=255)                                                              cAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGC‑TAAAAAAAGTAAAACAAGAACGCCCCGAGGGGCGGTATGTGTTGTAGACATGGGAAAACCGACTT‑AAAAGGAAAATcaca                                                                                  >  2:31046/1‑150 (MQ=255)                                                                atatCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGG‑TAAAAAAAGTAAAAAAAAAACGCCCTGAGGGGCGTTATGTGTTGTTGACATGGGAAAAACGACTT‑AAAAAGAAAATaacaag                                                                                >  2:161053/1‑148 (MQ=255)                                                                atatCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGC‑TAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTT‑AAAAGGAAAATCACAGgc                                                                               <  1:269771/151‑1 (MQ=255)                                                                   tCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGC‑CAAAAAAAGTAAAAAAAAAAAGCCCCGAGGGACGTTTTGTGTTTTTGAAATGGGAAAAACGACCT‑AAAAAGAAAAAAACAAgcg                                                                              >  2:362757/1‑149 (MQ=255)                                                                    ccATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGC‑TAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTT‑AAAAGGAAAATCACAGGcgcg                                                                            <  2:489204/150‑1 (MQ=255)                                                                            gACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGC‑TAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGctacta                                                                   <  2:512728/151‑1 (MQ=255)                                                                                          accaAAGATGTCACCATTGGCATTGGTATCGGCGGCCCTGGCGCAAAGAACAT‑AAACAAAAAAGCCCTCGGGGGCGCTTTGTGTTGTTGAAAAGGGAAAAACGAATTCAAAAGAAAAAACACA‑GCGCGGTGGAACTATTAttttttttagg                                                     >  1:640030/1‑144 (MQ=255)                                                                                               aGATGTCACCATTGGCATTGGTATCGGCGGGGCTGGC‑TTAAGAACGTAAAAAAAAAACGCCCCGAGGGGCGGTATGTGTTGTAGAAAAGGGAAAAACCACCT‑TAAAAGAAAAAAACAAGCGCGGAGGTACTATTATTTTTCTTAGGGGaa                                                 >  2:541839/1‑150 (MQ=255)                                                                                                   gTCACCATTGGCATTGGTATCGGCGGGCATGGG‑TTAAAGACGTAAAAAAAAAACGCCCTGAGGGGAGGTATGTGTTGTAGAAATGGGAAAAACGACTT‑TAAAAGAAAAACAAAAGGGCGATGCTACTATTATTTTTCttggggg                                                   >  2:446818/1‑138 (MQ=255)                                                                                                     cACCATTGGCATTGGTATCGGCGGGCATGGC‑TAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAgagga                                          <  1:218069/151‑1 (MQ=255)                                                                                                       ccATTGGCATTGGTATCGGCGGGCACGGC‑TCAAAAACGCAAAACAAGAATGCCCTCCGGGGCGTTATTTGTTGTTGAAAAGGTAAAAACGACTT‑AAAAAGAAAATCAAAAGCGCGGTGGGACTATGATTTTTCtttggggaaaaaagg                                           >  2:617098/1‑134 (MQ=255)                                                                                                         attggcattggTATCGGCGGGCAAGGC‑TAAAATAAGTAAAACAAAAACACCCTGAGGGGCGGTTTGTGTTGTAGAAAAGGTAAAAACCACTT‑AAAAAAAAAAAAACAGGGGCGGTGCTACTATTATTTTTCTTAGGAGaaaaaa                                             >  1:26029/1‑140 (MQ=255)                                                                                                            ggcattggTATCGGCGGGCATGGC‑TAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACa                                    <  2:572265/150‑1 (MQ=255)                                                                                                                  ggTATCGGCGGGCATGGC‑TAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCa                              <  2:90355/150‑1 (MQ=255)                                                                                                                  ggTATCGGCGGGCATGGC‑TAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCa                              <  2:205301/150‑1 (MQ=255)                                                                                                                   gTATCGGCGGGCATGGC‑TAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGc                            <  2:454593/151‑1 (MQ=255)                                                                                                                              cATGGC‑TAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACt                 >  1:158281/1‑151 (MQ=255)                                                                                                                                       aaGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCt         <  2:364324/151‑1 (MQ=255)                                                                                                                                              tAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAgg  >  2:189957/1‑151 (MQ=255)                                                                                                                                              tAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAgg  >  1:576748/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                              CCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGC‑TAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGG  >  NC_000913/1312065‑1312355

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