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

TCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAG  >  NC_000913/3206294‑3206579                                                                                                                                              |                                                                                                                                                 tCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACaaa                                                                                                                                          <  1:554893/150‑1 (MQ=255)    cAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATa                                                                                                                                        >  1:290150/1‑149 (MQ=255)     aaaaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACaaaaa                                                                                                                                        >  1:911322/1‑146 (MQ=255)            tGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATAaaaaaa                                                                                                                                  >  2:515755/1‑147 (MQ=255)            tGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTg                                                                                                                               <  1:272705/150‑1 (MQ=255)               cATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGt                                                                                                                           <  1:233796/151‑1 (MQ=255)                     cTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAACGGAggttcc                                                                                                                       >  2:143375/1‑144 (MQ=255)                        ctctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAAAAAACAAAAAATGGAGtt                                                                                                                          >  2:891489/1‑143 (MQ=255)                               gcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGaa                                                                                                           <  2:162856/151‑1 (MQ=255)                                  atatTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAAAACAAAAAATGGAGTTGCCATGATGAGCGaaa                                                                                                          >  1:407053/1‑149 (MQ=255)                                   tatTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAATTGAGTTGTCCTGATTAggg                                                                                                             >  2:539230/1‑143 (MQ=255)                                    attcattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGta                                                                                                       <  2:283482/150‑1 (MQ=255)                                     ttcattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGCCAAATACAAAAAATGGAGTTGGCGTGATGGGCGGAAGaaa                                                                                                      >  2:757786/1‑150 (MQ=255)                                       cattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaataa                                                                                                   <  1:340270/151‑1 (MQ=255)                                       cattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaata                                                                                                    <  1:485859/150‑1 (MQ=255)                                           cTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGTTGCCATGATGAGCGAAAGTAATAAGCa                                                                                                >  2:805078/1‑150 (MQ=255)                                              aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACa                                                                                             <  2:30735/150‑1 (MQ=255)                                              aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAg                                                                                            <  1:685248/151‑1 (MQ=255)                                                       ttGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAATAAAAAAAATGGAGTTGGCATGGTTAGGGGAAAGAAAAAACAAAAAGCAGTg                                                                                      >  1:471690/1‑148 (MQ=255)                                                         gATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAACAAAAAATTGGGTTGCGCTTATTGACGCAATTAATAAGCaaaaaaccgggaat                                                                                   >  2:777482/1‑137 (MQ=255)                                                         gATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATa                                                                                  <  1:13780/150‑1 (MQ=255)                                                          aTTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTCCCGAAAAATAAAAAAAATTGAGGTGCCATGATGAGCGGAAGGAAAAAACAAAACGCAggggaaaa                                                                                 >  1:259672/1‑143 (MQ=255)                                                            tCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAAAAAAAAAAATGGGGTTTGCCTGATGGGCGAAAATAATAAACAACAGGCAGTGAATaa                                                                                 >  1:616592/1‑148 (MQ=255)                                                            tCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATTGAGTTGCCATGATGAGCGAAAAGAATAAACAAAAGGCAGTGAATaa                                                                                 >  1:730511/1‑148 (MQ=255)                                                                gAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAAATGGAGGTTGCATGATGAGCGGAAGTAATAAACAAAAAGCAGGGAATaa                                                                                 >  2:677449/1‑144 (MQ=255)                                                                    tGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAAACAAAAACAAAAAATGGAGTTGCCATGATGAGGGAAAGTAATAAACAACAGGCAGTGAAAAAATTGACa                                                                          >  2:124696/1‑147 (MQ=255)                                                                        aaTTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTg                                                                   <  2:595801/150‑1 (MQ=255)                                                                         aTTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAAATGGAGTTGGCCTGATGAGCGAAAGTAATAAACAACAGGGAGTGAATAAGTTGACAGAGAtt                                                                    >  1:16249/1‑148 (MQ=255)                                                                             cATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTa                                                             <  2:975149/151‑1 (MQ=255)                                                                               tatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCATACAAATACAAAAAATGGAGTTGCCATGATGAGGGGAAGAAATAAAAAAAAGGGAGTGAATAAATTTACAGAGATTGTCGCTaa                                                            >  2:611233/1‑150 (MQ=255)                                                                                atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGATAAc                                                           <  2:554184/150‑1 (MQ=255)                                                                                    aCGGTGGCAAAAAATATAAAACCGCATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAAAGATTGTCGCTATCTTTAc                                                      <  2:751167/151‑1 (MQ=255)                                                                                     cGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCAACAAATAAAAAAACCTGAGGTTTCATGATGAGGGAAAGTAAAAAGCAACAAGCAGGGAATAAGATGAAAGAGATTGTTGCTAACTTTAcc                                                     >  1:893349/1‑151 (MQ=255)                                                                                           cAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAAAAAAAAAACTGGAGGTGGCATGATGGGGGAAAGTAATAAAAAAAAGGCAGGGAATAAGATGACAGAGATTGTTGCTAACTTTatcaccctg                                               >  2:446741/1‑143 (MQ=255)                                                                                            aaaaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATAAAAAAAATGGAGTTGCCATGATGAGGGAAAGTAATAAACAAAAAGCAGGGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATg                                               >  1:175164/1‑150 (MQ=255)                                                                                                     aaaaCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGGTGCCATGATGAGGGGAAGGAATAAACAACAAGCAGTGAATAAATTGAAAGAGAATGTTGCTAACTTTACCCCCATGATTTTTAccc                                     >  2:91800/1‑151 (MQ=255)                                                                                                       aaCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGa                                   <  1:535184/151‑1 (MQ=255)                                                                                                            caTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAATACAAAAAATGGGGTTGCCCTGGTGAGGGGAAGTAATAAAGAAAAGGCAGTGAATAAGATGACAGAGATTGTCGCTAACTTTTCCGCCATGATTTTTACCCGAAAGcc                              >  2:731874/1‑151 (MQ=255)                                                                                                            caTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAACTGGAGGTGCCATGATGAGCGGAAATAATAAGAAACAAGCAGGGAATAAGGTGACAGAGAATGTTGCTAAATTTTCCGCCATGATTTTTTCCCGAATGcc                              >  1:996744/1‑151 (MQ=255)                                                                                                              tttttAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCt                             <  2:939219/150‑1 (MQ=255)                                                                                                               ttttAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGAAAAATACAAAAAATGGGGGTGCCCTGGTTAGGGGAAGTAAAAAACAAAAAGCAGGGAATAAATTGACAGAGAGTGTTGCTAAATTTTACGCCATGATTTTTACCCGAATGCCt                             >  2:898681/1‑149 (MQ=255)                                                                                                                   aGTGGTAGTTTGTGGCGGTAAATTTTCCCAGAAAAATAAAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACGGGCAGTGAATAAGTTGACAGAAATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtgat                          <  2:386731/148‑1 (MQ=255)                                                                                                                    gTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGAc                       <  1:128603/150‑1 (MQ=255)                                                                                                                     tGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAAATGGAGGTGGCCTGATGAGCGGAAGTAATAAACAACAGGCAGTGAAAAAAATGACAGAGAATTTCGCCAACTTTACCGCCATGATTTTTACCCGAATGCCCGATGACg                      >  2:488859/1‑150 (MQ=255)                                                                                                                        tattttGTGGCGGTGAATTTTTCCAGACAAATACAAAAACGGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgtgg                   <  2:189293/147‑1 (MQ=255)                                                                                                                              gtgGCGGTGAATTTTTCCAGAAAAATACAAAACCTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAATTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATa             <  2:38960/150‑1 (MQ=255)                                                                                                                               tgGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATaa            >  2:232572/1‑150 (MQ=255)                                                                                                                                   ggTAAATTTTTCCAGAAAAATACAAAAACTGGAGTTGCAATGATGAGCGAAAGTAAAAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACt         <  2:195220/149‑1 (MQ=255)                                                                                                                                      gAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAAc    <  1:494422/151‑1 (MQ=255)                                                                                                                                         tttttCCAGACAAATACAAAAACGGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAg  <  2:290143/150‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                 TCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAG  >  NC_000913/3206294‑3206579

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