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

ATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAG  >  NC_000913/3206291‑3206585                                                                                                                                                 |                                                                                                                                                       aTTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGAc                                                                                                                                                   <  2:1069072/150‑1 (MQ=255)    tCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACaaa                                                                                                                                                <  1:923659/150‑1 (MQ=255)      gCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATAc                                                                                                                                             <  2:447997/151‑1 (MQ=255)       cAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAaaaaaa                                                                                                                                              >  2:291547/1‑147 (MQ=255)        aaaaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCTGAaaaaaa                                                                                                                                              >  2:833041/1‑146 (MQ=255)               tGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTg                                                                                                                                     <  1:850422/150‑1 (MQ=255)                   aTTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATgg                                                                                                                                    >  1:395443/1‑147 (MQ=255)                      cccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCCAACAaaaaaaaaaaacag                                                                                                                                    >  1:647284/1‑140 (MQ=255)                         tGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGAATTTTTCAGACAAAAACAAAAAATGGAGtttgcct                                                                                                                           >  1:814156/1‑145 (MQ=255)                          gCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCTAGACAAATACAAAAACTGGAGTTGTCatg                                                                                                                          <  2:305510/150‑1 (MQ=255)                            tctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCAGACAAAAAAAAAAAATGGGGGTTGCatgat                                                                                                                        >  1:987363/1‑150 (MQ=255)                               ggCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAg                                                                                                                     <  1:1027226/150‑1 (MQ=255)                               ggCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGc                                                                                                                    <  2:598867/151‑1 (MQ=255)                                  gcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAATAAAAAAAATGGAGGTGGCCTGATGAggg                                                                                                                   >  2:243338/1‑147 (MQ=255)                                        ttcattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTTAGACAAATAAAAAAAATGGAGTTGCCATGATGAGGGAAAGgaa                                                                                                            >  1:1006600/1‑150 (MQ=255)                                          cattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaata                                                                                                          <  2:24013/150‑1 (MQ=255)                                                 aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCa                                                                                                      <  2:6932/147‑1 (MQ=255)                                                      cgcgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCAGACAAAAACAAAAAATGGAGTTGCCATGATGAGCGGAAGTAATAAACAAAAAGCAg                                                                                              >  1:636905/1‑150 (MQ=255)                                                      cgcgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAATTTTTTCAAAAAAAAACAAAAACTTGGGTTTGCATTATTGGGGAAAGgaataaaaaaaaa                                                                                                  >  2:77383/1‑139 (MQ=255)                                                        cgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGTTGCCATGGTGAGCGGAAGTAAAAAACAACAAGCAggg                                                                                            >  1:266612/1‑148 (MQ=255)                                                         gTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGaa                                                                                          <  2:8377/151‑1 (MQ=255)                                                                    aaCTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACag                                                                               <  2:60041/151‑1 (MQ=255)                                                                    aaCTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACa                                                                                <  2:676510/150‑1 (MQ=255)                                                                     aCTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAATTTTTTTACACAAAAAAAAAAAAAGGAGTTGCCATGATTAGGGAAAATAATAAACAAAAAGCAGTGAATaa                                                                                       >  1:348501/1‑142 (MQ=255)                                                                        gTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGAt                                                                           <  1:625068/151‑1 (MQ=255)                                                                             ttACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATTGGGTTGCCATGATGAGGGGAAAGAATAAACAAAAGGGAGGGAATAAGATGACAGAGATTGt                                                                        >  1:34139/1‑149 (MQ=255)                                                                             ttACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGAtt                                                                          <  1:639682/147‑1 (MQ=255)                                                                                   atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAAc                                                                 <  1:1052392/150‑1 (MQ=255)                                                                                   atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACt                                                                <  2:84818/151‑1 (MQ=255)                                                                                          gTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAAAAAAAAAAATGGAGGTGTCCTGATGAGCGGAAGTAAAAAACAAAAAGCAGGGGATAAATTGACAGAGATTGTTGCTAACTTTAcccc                                                         >  1:440303/1‑149 (MQ=255)                                                                                              cAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTTAGAAAAAAAAAAAAAATGTAGGTTGCATGGTTTGCGGAAGGAATAAAAAACAGGCAAGGAAAAAATTGAAGGAGAGTGTCTCTaaactttctgcc                                                        >  2:461061/1‑138 (MQ=255)                                                                                               aaaaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCAGACAAAAAAAAAAAAAAGGGGTGCCCTTGTGGGGGGAAGAAATAAACCACAAGCAGTGAATAAATTGACAGAGAATTTCGCTAAATTTACCg                                                          >  1:364584/1‑145 (MQ=255)                                                                                                aaaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACGGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGAt                                                   <  1:976547/151‑1 (MQ=255)                                                                                                   aaTATAAAACCACATTTTTAGTGGTAGTTTGGGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTc                                                <  1:327892/151‑1 (MQ=255)                                                                                                    atatAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATTGGGGTTCCATGATGAGCGGAAAGAAAAAACAAAAAGCAGTGAAAAAATTGGCAGAGATTTTCGCTAAATTTAACCCCATGAttttt                                               >  1:1015682/1‑149 (MQ=255)                                                                                                     tataAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATAAAAAAAATGGGGTTGGCCTGGTGAGCGGAAAGAATAAAAAAAAGGCAGGGGATAAATTGACAGAGATTGTCGCCAAATTTACCCCCCTGAttttt                                               >  1:555043/1‑148 (MQ=255)                                                                                                       taAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAAATGGAGATGCCATGATGAGCGAAAGGAATAAACAAAAAGCAGTGAATAAGTTGAAAGAGATTGTCGCTAACTTTACCGCCATGATTTCt                                               >  1:660797/1‑148 (MQ=255)                                                                                                        aaaaCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAATTGAGGTGGCATGATGAGGGGAAGAAATAAGCAACAGGCAGTGAATAAAATGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTAccc                                           >  2:447200/1‑151 (MQ=255)                                                                                                            ccACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATAAAAAAAATTGAGTTGCCATGATGAGCGAAAGAAATAAACAAAAAGCAGGGAATAAATTGACAGAGATTGTTGCTAACTTTTCCGCCATGATTTTTTCCCGaa                                        >  1:23180/1‑150 (MQ=255)                                                                                                                 tttttAGTGGTAGTTTGTGGCGGTAAATTTTACAAGAAAAAAACAAAAACAGGATTTGCCATGAGGAGCGAAAGTAATAAGCAAAAGGTGGTGAATAAGTTGACAGATTTTGTCGTTAACTTTACCGCCATGATTTCTACCCGAATGCCt                                   <  2:643682/150‑1 (MQ=255)                                                                                                                     tAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAAATGGGGTTGCCCTGATGAGCGAAAGTAATAAACAACAAGCAGTGAATAAATTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTTTACCCCAAAGCCtgatg                               >  1:173908/1‑150 (MQ=255)                                                                                                                        tGGTAGTTTGTGGCGGGGAAATTTTACCCAAAAATACAAAAACTGGAGGTGCCATGATTGGCGGAAGTAATAAGCAACAGGCAGTGAAAAAATTGATAGAGAATGTCGCTTAATTTTCCGCCATGATTTTTAACCGAAAGcc                                    >  2:165991/1‑142 (MQ=255)                                                                                                                         ggTAGTTTGTGGGGGTAAATTTTTCAAGAAAAAAAAAAAAACTGGAGTTGCCATAATGAGCGAAAGTAATAAGCAACAGGCATTGAATAAGTTAACAGAGATTGGCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtgatga                              <  2:718115/147‑1 (MQ=255)                                                                                                                         ggTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGTTTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGAc                             <  1:113710/148‑1 (MQ=255)                                                                                                                         ggTAGTTTGTGGCGGTGAATTTTTCCAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtgatga                              <  1:833048/147‑1 (MQ=255)                                                                                                                            aGTTTGTGGCGGTAAATTTTTCCAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgtgg                         <  1:434119/149‑1 (MQ=255)                                                                                                                                   gggggTGATTTTTCCCAGACAAGAACAAAAACTGGAGTTGCCACGGTGAGCAAAAGTAATAAGCAACAGGTAGTGAATATGTGGACAGAGATTGGCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAAc                <  2:617295/148‑1 (MQ=255)                                                                                                                                      ggTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTa              >  1:422534/1‑150 (MQ=255)                                                                                                                                          aaTTTTCCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTaaaa           <  1:132475/149‑1 (MQ=255)                                                                                                                                            tttttCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAg        <  1:184357/150‑1 (MQ=255)                                                                                                                                              tttCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCGGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCt      <  1:908958/150‑1 (MQ=255)                                                                                                                                              tttCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCt      >  2:818799/1‑150 (MQ=255)                                                                                                                                              tttCAAGAAAAATACAAAAAGTGGATTTGCCATGATGAGCGAAATTAATAAGCAAGGGGCAGTGAATAAGTGGACGGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCt      <  2:731653/150‑1 (MQ=255)                                                                                                                                                 ccAGAAAAATACAAAAACTGGGTTTGCAATGAGGAGCAAAAGTAATAAGCAACAGGCGTTGAATAAGTGGACAGAGTTTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAg  <  2:740112/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                       ATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAG  >  NC_000913/3206291‑3206585

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