Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,206,435	0	C	T	46.5%	33.4 / 23.2	43	intergenic (‑179/‑28)	ttdR/ttdA	transcriptional activator of ttdABT/L‑tartrate dehydratase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (4/19);  new base T (20/0);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 1.65e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAA  >  NC_000913/3206294‑3206576                                                                                                                                               |                                                                                                                                              tCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAAATTTTTCATACaaa                                                                                                                                       >  1:87145/1‑150 (MQ=255)  cGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTAGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAAATTTTCCAGACaaa                                                                                                                                       >  1:641653/1‑149 (MQ=255)    cAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTCCAGACAAATACa                                                                                                                                   <  2:266183/151‑1 (MQ=255)        aCATTGCGATTCCCTGCTCTGGCGCAATAGTCATTCTTATTTTCGCGTTGATTCACGAACTTTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTGTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTCCAGACAAATACaaaaa                                                                                                                               <  1:149478/151‑1 (MQ=255)        aCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTCCAGACAAATACaaaaa                                                                                                                               <  1:253261/151‑1 (MQ=255)           ttGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAAATTTTTTAGACAAATAaaaaaaa                                                                                                                              >  1:134629/1‑148 (MQ=255)               cATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTCCAGACAAATACAAAAACTGGAGt                                                                                                                        <  2:646354/151‑1 (MQ=255)                       gCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCat                                                                                                                  <  2:405288/149‑1 (MQ=255)                          ctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTTCCGACAAAAATAAAAAATTGAGTTGGCatgttg                                                                                                              >  2:708417/1‑150 (MQ=255)                                cAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGaaa                                                                                                       <  1:609748/151‑1 (MQ=255)                                  atatTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGC‑GGAATTTTTCACACAAATACAAAAAATAGAGATGACaaaattggt                                                                                                           >  2:489647/1‑137 (MQ=255)                                       cattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaata                                                                                                 <  2:703932/150‑1 (MQ=255)                                        attcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAATTATAAAACCACATTTTTAGTGGTAGTTTGGGGCGGG‑GAAATTTTTTACAAAAAATCAAAAAATTGGGTTGCCCTGTTGAggggaaaaaaaaa                                                                                                >  1:531316/1‑138 (MQ=255)                                          tcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTTCAGACAAAAACAAAAAACGGGGATGGCATGATGAGCGAAAGgaata                                                                                                 >  2:483087/1‑147 (MQ=255)                                                tttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAgg                                                                                        <  2:831188/150‑1 (MQ=255)                                                   cgcgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGt                                                                                    <  1:623318/151‑1 (MQ=255)                                                   cgcgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGt                                                                                    <  2:729380/151‑1 (MQ=255)                                                          aTTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTTAATACAAATAAAAAAAATGGAGTTGGCATGAAGGGCGTAAGGAATAAGCAACAGGCAGGGAATa                                                                               >  1:864132/1‑149 (MQ=255)                                                           ttCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTTCAGACAAAAACAAAAAATGGAGGTGCCATGATGAGCGGAAACAATAAACAAAAAGCAGTGAaaaa                                                                              >  2:428908/1‑149 (MQ=255)                                                                                atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACt                                                       <  2:280331/151‑1 (MQ=255)                                                                                   cACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAAATTTTTCCGACAAAAACAAAAAATGGAGGTGGCATTATGAGCGGAAGAAATAAACAAAAAGCAGTGGATAAATTGACAGAGATTGTCGCTAACttt                                                     >  2:112975/1‑150 (MQ=255)                                                                                     cGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTAcc                                                  <  1:816418/151‑1 (MQ=255)                                                                                        tGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGG‑GAAATTTTTTAGAAAAATACTAAAAATAGAGTTGTCATGATGGGCGAAAAGAAAAAGAAACAAAAAATGGAAAAAATTACAGAGAATTTCGCTAACATTACcccc                                               >  1:490576/1‑151 (MQ=255)                                                                                         ggCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GGATTTTTTCACACAAAAACAAAAAACTGAGGTGCCATGATGGGGGGAAAGAATAAACAAAAGGCAGGGGAAAAGTTGAAAGAGAATTTTGCTCACttt                                                     >  2:76709/1‑144 (MQ=255)                                                                                              aaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAAATTTTTCAGACAAATACAAAAAATGGAGTTGGCATGATGAGGGAAAGTAATAAAAAAAAGGAAGAGAATAAGTTGACAGAGAATGTTGCTAAATTTAACCCCACGAtt                                         >  1:702158/1‑151 (MQ=255)                                                                                               aaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTTCAGACAAAAACAAAAAATGGAGGTGCCATGATGAGCGAAAGGAATAAGCAACAGGCAGTGAATAAATTGACAGAGATTGTCGCTAAATTTACCGCCATGAttt                                        >  2:312569/1‑151 (MQ=255)                                                                                                          cacaTTTTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTTCAGACAAATAAAAAAAATGGAGTTGCCATGATGAGCGAAAGTAAAAAAAAAAAGGCAGGGAAAAAATTTACAGAGATTGTTGCTAAATTTACCCCCCTGATTTTTAacc                                  >  2:571669/1‑146 (MQ=255)                                                                                                           acaTTTTTAGTGGTAGTTTGTGGCGGTGGAATTTTCCTGGACAAAAAAAAAAATTGG‑GTTGGCATGGTGTGGGGAAGAAATAAACAAAAAGAAGGGGAAAAATTTAAAGAGATTGTCGCTAAACTTAcacccaacggtttt                                      >  1:838676/1‑128 (MQ=255)                                                                                                           acaTTTTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTTCAGAACAATATAAAAAATGGAGTTGTCCTGATGAGAGAAAGTAAAAAGCAACAGACAGAGAATAAGATGACAGAGATTGTCGCTAACTTTAACGCCATGATTTTTAcccca                                >  1:136919/1‑145 (MQ=255)                                                                                                             aTTTTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTTCCGACAAATACAAAAAATGGAGGTGCCATGATGAGCGGAAAAAATAAACAAAAGGCAGGGAAAAAGATGAAAGAGATTGTCGCTAACTTTACCGCCCTGATTTTTACCCGAAAGcc                           >  2:853105/1‑150 (MQ=255)                                                                                                              tttttAGTGGTAGTTTGTGGCGGG‑GAATTTTTTCAGACAAATACAAAAAATGGCGGTTGCATGATGGGGGGAAGGAAAAAAAAAAAGGCCGGGAATAAATTGTCAGAGATTGTCGCTTACTTTAACCCCATGATTTTTACCCGCATGcc                           >  2:692713/1‑149 (MQ=255)                                                                                                              tttttAGTGGTAGTTTGTGGCGGG‑GAAATTTTTCCGACAAATATAAAAAATTGGGTTGGCAAGATGAGGGAAAGTAATAAACAACAAGCAGTGAATAAGATTGCAGAGAATGTCGCTCACTTTATCGCCCTGAATTTTACCCgca                               >  1:449173/1‑143 (MQ=255)                                                                                                                      ggTAGTTTGTGGCGGT‑GAATTTTTTCAGACAAATACAAAAAATGGAGTTGCCCTGATGAGCGAAAGGAATAAGAAACAGGCAGGGAATAAGTTGACAGAGATTGTCGCTAAATTTACCGCCATGATTTCTACCCGAATGCCTGATGACgtg                 >  2:238297/1‑151 (MQ=255)                                                                                                                      ggTAGTTTGTGGCGGT‑GAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgtg                 <  2:97515/151‑1 (MQ=255)                                                                                                                       gTAGTTTGTGGCGGT‑GAATTTTTTCAGACAAATACAAAAACATGAGTTTGCATGATGAGGGAAAGGAATAAACAAAAAGCAGGGAATAAATTGAAAGAGAATGTCGCTAACTTTAACACCCTGATTTTTAACCCAAAGCCCGGTGACgtgg                >  1:547738/1‑151 (MQ=255)                                                                                                                        tAGTTTGTGGCGGT‑GATTTTTTCCAGAAAAATACAAAAACTGGAGTTCCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGTTTGTCGCTAATTTTACCGCCATGATTTCTACCCGAATGCCtgatg                      <  2:363379/144‑1 (MQ=255)                                                                                                                              ggggCGGT‑GAATTTTTCCAGAAAAATAAAAAACCTGGATTTGCAATGTTGAGCGAAAGAAATAAGCAACAGGCAGTGAATAAGTTGACAGAAATTGTCGTTACCTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGAt           <  2:402550/147‑1 (MQ=255)                                                                                                                                ggCGGT‑GAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGTTTGTCGTTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAAc       <  2:456705/151‑1 (MQ=255)                                                                                                                                 gCGGT‑GAATTTTCCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATGGTCGCTACCTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACt      <  1:529902/151‑1 (MQ=255)                                                                                                                                  cGGT‑GAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACt      >  2:505652/1‑150 (MQ=255)                                                                                                                                       gAATTTTACTAGACAAATACAAAACCTGGAGTTGCCATGATGAGCGAACATAGTAAGCAACAGGCAGAGAATAAGTGGACAGATATTGTCGCTAACTTTACCGCCATGATTTCTATGCGAATGCCTGATGAGGCGGTGGAtaaaataaaa  >  1:740850/1‑150 (MQ=255)                                                                                                                                        aaTTTTTCCAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTaaaa  <  2:349927/149‑1 (MQ=255)                                                                                                                                        aaTTTTTCCAGAAAAATAAAAAAACTGGAGTGGCCATGATGAGCGAAAGAAAAAAGCAACAGGAAGTAAATAAGTGGAAAGGGTTGGTCGCTAATTTTACCGCCATGATTTCTACCAGAATGCCTGATGACGTGGTGGATAAACTaaaa  <  1:880043/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                              TCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGT‑GAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAA  >  NC_000913/3206294‑3206576

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