Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,312,200	0	G	A	46.9%	50.3 / 25.8	49	T16I (ACT→ATT)	yciA	acyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base G (5/21);  new base A (23/0);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 2.87e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATC‑GCCCTGAGGGACGTTATGTGTTGTAGAC‑ATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCT  >  NC_000913/1312050‑1312348                                                                                                                                                       |                                                                                                                                                       cGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAg                                                                                                                                                        <  2:658246/151‑1 (MQ=255) cGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAg                                                                                                                                                        <  1:81815/151‑1 (MQ=255)       aaaGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGATGGCTAAAgaaa                                                                                                                                                     >  1:649664/1‑145 (MQ=255)         aGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGTATTGGTATCGGCGGGCATGGCTAAAGTACGTaaa                                                                                                                                                <  1:155495/151‑1 (MQ=255)             aTTCCTTCAACCCGAACAGTCACTACGCGTCTGTGGGCAATTGCATGTCCCACAATAGCGCCGCCAATGTCCATTTGTGACACTAACCAACCACAAAAGATGTCAACATTGGCATTGGTATCGGCGGGCATGGCTAATGTACGTAAAACaa                                                                                                                                            <  1:222049/151‑1 (MQ=255)                  ttCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCGAAAAAACGGAAAACaaaa                                                                                                                                          >  1:82210/1‑146 (MQ=255)                   tCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATC‑g                                                                                                                                      <  1:224939/150‑1 (MQ=255)                   tCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAATAAGGAAAAAAAGATC‑Gc                                                                                                                                     >  2:5017/1‑151 (MQ=255)                    cAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCGATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATC‑Gc                                                                                                                                     <  2:401142/150‑1 (MQ=255)                      aCCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGAACGTAAAAAAAAAAC‑Gccc                                                                                                                                   >  1:868761/1‑150 (MQ=255)                      aCCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAACGTAAAAAAAAAAC‑Gccc                                                                                                                                   >  2:900778/1‑150 (MQ=255)                       cccGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGATGGCTAAAATAAGTAAAAAAAGAAA‑Gccc                                                                                                                                   >  1:715047/1‑149 (MQ=255)                       cccGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATC‑GCCCTg                                                                                                                                 <  1:854081/151‑1 (MQ=255)                            acaGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAGACGTAAAACAAAAAC‑ACCCCGAggg                                                                                                                             >  1:36832/1‑150 (MQ=255)                            acaGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAAAGTAAAAAAAAATC‑GCCCTGAggg                                                                                                                             >  2:920234/1‑150 (MQ=255)                               gTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCAAGGGTAAAATACGTAAAAAAAAAAC‑GCCCCGggggggg                                                                                                                          >  2:282763/1‑147 (MQ=255)                                 cACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATC‑GCCCTGAGGGACGtt                                                                                                                        <  2:364521/150‑1 (MQ=255)                                 cACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATC‑GCCCTGAGGGACGtt                                                                                                                        <  2:49366/150‑1 (MQ=255)                                         gACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCCGGGGTAAAAAAAAGAAAACAAAAAC‑GGCCTCAGGGACGGTATTtgtggt                                                                                                               >  1:702914/1‑151 (MQ=255)                                           ccGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAACGTAAAAAAAAAAC‑GCCCTGAGGGGGGGTAtgtggt                                                                                                                 >  1:704733/1‑145 (MQ=255)                                                   aaTTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGGTAAAAGACGTAAAAAAAAAAA‑GCCCCGGGGGGCGGTATGTGTTGTAGAC‑ATgg                                                                                                      >  2:337563/1‑149 (MQ=255)                                                    aTTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATC‑GCCCTGAGGGACGTTATGTGTTGTAGAC‑ATGGTaa                                                                                                   <  2:587964/151‑1 (MQ=255)                                                      ttCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATC‑GCCCTGAGGGACGTTATGTGTTGTAGAC‑ATGGTaaa                                                                                                  <  2:261205/150‑1 (MQ=255)                                                           ttGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGAATGGCTAAAATACGTAAAACAAAAAA‑GCACCGCGGGGAGGGATTTGTTGTAGAC‑ATGGTgaaaacaac                                                                                            >  1:911767/1‑146 (MQ=255)                                                             gCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGAGC‑GCCCTGAGGGACGTTATGTGTTGTAGAC‑ATGGTAAAACCGACtt                                                                                          <  2:459611/151‑1 (MQ=255)                                                                aGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATC‑GCCCTGAGGGACGTTATGTGTTGTAGAC‑ATGGTAAAACCGACTTaa                                                                                        <  2:193864/150‑1 (MQ=255)                                                                  aaTAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATC‑GCCCTGAGGGACGTTATGTGTTGTAGAC‑ATGGTAAAACCGACTTAAAAg                                                                                     <  1:243833/151‑1 (MQ=255)                                                                        gccgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAATACGTAAAAAAAGAAC‑GCCCTGCGGGGGGGTATGTGTTGTAGAC‑ATGGGAAAAACCACTTAAAAAGaaaa                                                                                >  1:200185/1‑150 (MQ=255)                                                                         ccgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGATGGGCAAAAAAAGGAAAAAAAAAAC‑GCCCCGAGGGGGGGTTTGTGTTGTTGAA‑ATGGTAAAAAAGACTTTAAAAGaaaa                                                                                >  1:491718/1‑149 (MQ=255)                                                                          cgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAAAGTAAAACAAAATC‑GCCCCGAGGGGCGGTATGTGTTGTAGAA‑ATGGGAAAAACGACTTAAAAAGaaaa                                                                                >  2:493397/1‑148 (MQ=255)                                                                            ccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTCAAAAAAGGAAAAAAAAAAC‑GCCCCCAGGGGGGGTATGTGTTTTAGAA‑AAGGGAAAAACCACTTAAAAAGaaaa                                                                                >  1:345903/1‑146 (MQ=255)                                                                            ccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATC‑GCCCTGAGGGACGTTATGTGTTGTAGAC‑ATGGTAAAACCGACTTAAAAGGAAAATcaca                                                                           <  1:590697/151‑1 (MQ=255)                                                                                   ccATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATC‑GCCCTGAGGGACGTTATGTGTTGTAGAC‑ATGGTAAAACCGACTTAAAAGGAAAATCACAGGcgcg                                                                     <  2:641578/150‑1 (MQ=255)                                                                                   ccATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAATACGTAAAAAAAAAAC‑GCCCTGAGGGGCGTTATGTGTTGTAGAA‑AAGGTAAAAACGACTTAAAAAGAAAAAAACAAgcg                                                                       >  2:819462/1‑148 (MQ=255)                                                                                             cATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCCAAAAAAAGTAAAAAAAAAAC‑GCCCTTCGGGGGGGTTTGTGGTGTAGAC‑ATGGGAAAAACCACTTTAAAAGAAAAAAACAAGGGCGGTGCTACTAt                                                           >  1:437726/1‑150 (MQ=255)                                                                                                tAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTCAAATACGTAAAAAAAAATC‑GCCCCGTGGGGAGGTATGTGTTGTTGAA‑ATGGGAAAAACGCCTTTAAAAGAAAATCACAAGGGCGATGCTTCTattatt                                                       >  2:507607/1‑151 (MQ=255)                                                                                                    caaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCCCCGCCAAAATAACTAAAAAAAGATCGGCCCTG‑GGGGCGGTCTGTGTTGTAGAT‑ATGGGAAAAACCAATTAAAAAGAAAAAAACAGGGGCGATGgtacttttaattt                                                     >  2:82061/1‑141 (MQ=255)                                                                                                    caaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGAATTGCTACAAAACGCAAAAAAAAAAA‑ACCCTTAAGGACGGTCTGTGTTGTTGACAAAGGGAAAACC‑ACCTAAAAAGAAAAACACAGGCGGGGTGGTACTATAGtttt                                                     >  1:607148/1‑149 (MQ=255)                                                                                                     aaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAATACGCAAAAAAAAATC‑GCCCTGAGGGGCGCTATGTGTTGTAGAC‑AAGGTAAAAACAACTTAAAAAGGAAATCACAGGCGCGATGCTACTATTAtttttc                                                   >  2:7140/1‑148 (MQ=255)                                                                                                       ccaccaAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGGAAGGAAAAAAAAATC‑GCCCTGGGGGGGGGTATGTGTTGTAGAA‑AAGGGAAAAACGACTTCAAAAGAAAAACAAACGGGCGGGGGTACTATTAttttt                                                    >  1:187905/1‑146 (MQ=255)                                                                                                                 tGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAATACGTAAAACAAAAAC‑GCCCTGAGGGGCGGTATGTGTTGTAGAC‑ATGGTAAAAACGACTTTAAAAGAAAAACACAGGCGCGATGCTACTATGATTTTTCTGAGGAGaa                                          >  2:444280/1‑147 (MQ=255)                                                                                                                  gTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATC‑GCCCTGAGGGACGTTATGTGTTGTAGAC‑ATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACaga                                      <  1:53343/150‑1 (MQ=255)                                                                                                                    cACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATC‑GCCCTGAGGGACGTTATGTGTTGTAGAC‑ATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAgagga                                   <  1:5017/151‑1 (MQ=255)                                                                                                                             cattggTATCGGCGGGCATGGCTAAAGGACGTAAAAAAAAAAC‑GCCCTGAGGGGGGGTATTTGTTGTAGAA‑ATGGGAAAAACGACTTAAAAAGAAAAAAACAGGCGCGATGCTACTATGATTTTTCTTTGGGGAAAAGAGGAGATaaa                             >  1:7215/1‑146 (MQ=255)                                                                                                                             cattggTATCGGCGGGCATGGCTAAAAGAAGGAAAACAAAAAC‑GCCCTGAGGGGCGGTATGTGTTGTTGAC‑ATGGTAAAAACGAATTAAAAAGAAAAAAACAAGCGCGGTGCTCATATGATTTTTCTGAGGAGAACAGAGGAGATACa                             >  1:402622/1‑148 (MQ=255)                                                                                                                                      cggcggGCATGGCTAAAGTACGTAAAACAAGATC‑GCCCTGAGGGACGTTATGTGTTGTAGAC‑ATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCg                 >  2:589730/1‑151 (MQ=255)                                                                                                                                         cggGCATGGCTAAAGTACGTAAAACAAGATC‑GCCCTGAGGGACGTTATGTGTTGTAGAC‑ATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAg              <  2:734486/151‑1 (MQ=255)                                                                                                                                         cggGCATGGCTAAAGTACGTAAAACAAGATC‑GCCCTGAGGGACGTTATGTGTTGTAGAC‑ATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAg              <  1:790440/151‑1 (MQ=255)                                                                                                                                          gggCATGGCTAAAGTACGTAAAACAAGATC‑GCCCTGAGGGACGTTATGTGTTGTAGAC‑ATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAgg             <  1:603939/151‑1 (MQ=255)                                                                                                                                           ggCATGGCTAAAGTACGTAAAACAAGATC‑GCCCTGAGGGACGTTATGTGTTGTAGAC‑ATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGa            <  2:715128/151‑1 (MQ=255)                                                                                                                                               tGGCTAAAGTACGTAAAACAAGATC‑GCCCTGAGGGACGTTATGTGTTGTAGAC‑ATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTgg        >  2:811169/1‑151 (MQ=255)                                                                                                                                                     aaGTACGTAAAACAAGATC‑GCCCTGAGGGACGTTATGTGTTGTAGAC‑ATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCt  <  2:52534/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                       CGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATC‑GCCCTGAGGGACGTTATGTGTTGTAGAC‑ATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCT  >  NC_000913/1312050‑1312348

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