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

GTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAG  >  NC_000913/3206289‑3206585                                                                                                                                                    |                                                                                                                                                       gtATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCCaa                                                                                                                                                     >  1:242779/1‑150 (MQ=255) gtATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTTCaaa                                                                                                                                                    >  1:367834/1‑149 (MQ=255) gtATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTTCAGa                                                                                                                                                    >  2:595163/1‑151 (MQ=255) gtATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTTCAGa                                                                                                                                                    >  2:595153/1‑151 (MQ=255) gtATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGACTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTTCaaa                                                                                                                                                    >  2:296772/1‑149 (MQ=255)           aaaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTTCAGACAAATACaaa                                                                                                                                          >  2:41418/1‑151 (MQ=255)              cATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAAc                                                                                                                                       <  2:149575/151‑1 (MQ=255)                       tCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTTCAGACAAATACAAAAAATGGAGGTgg                                                                                                                              >  1:418413/1‑150 (MQ=255)                       tCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTTCAGACAAATAAAAAAAATGGGGtttg                                                                                                                              >  1:491775/1‑149 (MQ=255)                       tCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AAATTTTTCAGACAAATAAAAAAAATgggggtgg                                                                                                                              >  1:611816/1‑145 (MQ=255)                         ccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCa                                                                                                                            <  2:568544/151‑1 (MQ=255)                              tctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatgatg                                                                                                                       <  2:301799/151‑1 (MQ=255)                                       atatTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTTCCAACAAAAAAAAAAAATGGAGGTGGCATGATGAGAGaaaag                                                                                                              >  1:43377/1‑149 (MQ=255)                                       atatTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTTCAGACAAATACAAAAAATGGAGTTGTCATGATGAGCGAAAGt                                                                                                              >  1:453983/1‑151 (MQ=255)                                           tcattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaata                                                                                                          <  1:148483/151‑1 (MQ=255)                                              ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTTCAAACAAATAAAAAAAATGGAGTTGCCATGATTAGCGAAAGTAATAAg                                                                                                        >  1:241258/1‑150 (MQ=255)                                                     tttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGc                                                                                                <  1:41668/151‑1 (MQ=255)                                                              gATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTTCAAACAAAAACAAAAAATGGAGTTGGCATGATGAGCGAAAGGAAAAAACAACAAACAAAGAATaa                                                                                       >  1:129464/1‑151 (MQ=255)                                                                ttCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGt                                                                                     <  2:476473/151‑1 (MQ=255)                                                                         tGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTTCAAACAAAAAAAAAAACTTGAGTTGCCCTGATGAGGGAAAGTAATTAGAAAAAGGCAGTGAATAAGTTGACAgaga                                                                            >  2:828267/1‑151 (MQ=255)                                                                          gTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGGAAAAAACAAAAGGCAGTGAATAAATTGACAGAGAt                                                                           >  2:864931/1‑151 (MQ=255)                                                                          gTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AAATTTTTCAGACAAAAACAAAAAATGGAGTTGCCATGATGAGCGAAAATAAAAAGAAAAAAGCAATGAATAAATTGACAGAGAt                                                                           >  1:766232/1‑151 (MQ=255)                                                                          gTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AAATTTTTCAAAAAAAAACAAAAACTGGAGTTGTCATGATGGGGGAAAAGAATAAACAAAAAGCAGAGAAAAAATTGACAGAGAt                                                                           >  1:576356/1‑151 (MQ=255)                                                                            tAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTTCCGACAAATAAAAAAAATTGGGTTTCCATGATGAGCGGAAAGAATAAACAACCAGCAGTGGAAAAATTGACAGAGATTg                                                                         >  2:328156/1‑151 (MQ=255)                                                                              aTTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTTCAAACAAAAACAAAAAATGGAGTTGGCATGATGAGAGAAAGTAATAAACAAAAAGCAATGAAAAAATTGACAAAGAATGtt                                                                       >  1:203244/1‑150 (MQ=255)                                                                                 aCATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTTCAGACAAAAACAAAAAATGGAGGTTGCATGATGAGCGAAAGTAATAAACAACAAGCAGGGAATAAGTTGACAGAGATTGTCg                                                                      >  1:612078/1‑149 (MQ=255)                                                                                   atatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTTCAGACAAATAAAAAAAATGGAGTTGCCATGATGAGCGAAAGTAATAAACAACAGGCAGGGAATAAGTTGACAGAGATTGTCGCTaa                                                                  >  2:364508/1‑151 (MQ=255)                                                                                   atatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTTCAGACAAAAACAAAAAATGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAAGCAGTGAATAAGTTGACAGAGATTGTCGCTaa                                                                  >  2:647659/1‑151 (MQ=255)                                                                                     atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACt                                                                <  2:195648/151‑1 (MQ=255)                                                                                             tGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGc                                                         <  2:367836/150‑1 (MQ=255)                                                                                             tGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGc                                                         <  1:448186/150‑1 (MQ=255)                                                                                                   aaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGAtt                                                  <  1:248041/151‑1 (MQ=255)                                                                                                                acaTTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGc                                     <  1:457239/151‑1 (MQ=255)                                                                                                                     tttAGTGGTAGTTTGTGGCGGTG‑AATTTTTTCAGACAAATACAAAAAATGGAGTTGCCATGATGAGCGGAAGGAATAAACAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAAATTTACCGCCATGATTTCTACCCGAATGCCtgat                                >  2:101231/1‑151 (MQ=255)                                                                                                                       tAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtgatga                              <  1:178327/151‑1 (MQ=255)                                                                                                                         gTGGTAGTTTGTGGCGGTG‑AATTTTTTCAGACAAATACAAAAAATGGAGGTGGCATGATGAGCGAAAGTAATAAACAAAAAGCAGTGAAAAAATTGACAGAGATTGTTGCTAACTTTACCGCCATGATTTTTACCCGAATGCCTGATGACg                            >  1:533574/1‑151 (MQ=255)                                                                                                                                      gCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAAc                <  2:649313/150‑1 (MQ=255)                                                                                                                                            aaaTTTTTCAAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAAc          <  1:701366/150‑1 (MQ=255)                                                                                                                                             aaTTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAAc          <  1:491614/150‑1 (MQ=255)                                                                                                                                                    ccAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAg  <  2:508/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                       GTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAG  >  NC_000913/3206289‑3206585

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