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

AGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAA  >  NC_000913/3206286‑3206575                                                                                                                                                      |                                                                                                                                              aGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTGTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTcc                                                                                                                                              <  2:209072/151‑1 (MQ=255)     tATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGa                                                                                                                                           <  2:753304/150‑1 (MQ=255)         tCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACaaa                                                                                                                                       <  2:57505/150‑1 (MQ=255)          cGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTTTGGCGGTGAATTTTTCCAGACAAATa                                                                                                                                     <  1:691394/151‑1 (MQ=255)                aCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGacaaaaacaaaa                                                                                                                                >  1:350468/1‑150 (MQ=255)                    tGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTg                                                                                                                            <  1:795760/150‑1 (MQ=255)                     gCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGAAAAAAACaaaaa                                                                                                                               >  2:683115/1‑146 (MQ=255)                     gCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTgg                                                                                                                           <  1:275372/150‑1 (MQ=255)                      ccATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑g                                                                                                                        <  1:140541/151‑1 (MQ=255)                        aTTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑Gtt                                                                                                                      <  2:822329/151‑1 (MQ=255)                               gCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTTTATTTTTTAAGACCAAAAAAAAAACTTGA‑GGTGGCa                                                                                                                  >  2:121575/1‑148 (MQ=255)                                 tctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAAACAAAAAAAAAAAATGGA‑GGTggcctg                                                                                                                >  1:399196/1‑143 (MQ=255)                                    ggCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAAACAAATAaaaaaaattgggtttgcat                                                                                                                  >  2:58783/1‑131 (MQ=255)                                        cAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAAAAAAAACaaaaa                                                                                                                               >  2:506337/1‑127 (MQ=255)                                         aaTATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTTAGAAAAAAAAAAAAAATGGA‑GTTGTCATGATGAgggga                                                                                                       >  2:768219/1‑145 (MQ=255)                                              tcattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGtaat                                                                                                 <  1:323905/150‑1 (MQ=255)                                                      aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGCGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGCGCATTTTTCCCGCCAAAAATAAAAAACGTA‑GTTGTCTTGATAAGCGAAataaataaaaa                                                                                            >  1:702652/1‑144 (MQ=255)                                                          tCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCa                                                                                     <  1:595856/150‑1 (MQ=255)                                                           cgcgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGt                                                                                   <  2:828318/151‑1 (MQ=255)                                                            gcgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAACAAAAAATGGA‑GGTGCCAATGGGAGCTCAAGTGATAAACAAAAGGccgg                                                                                   >  2:467136/1‑147 (MQ=255)                                                                tGATTCACGAACTGTTAATTACATATTCACGGTGGCACAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAAATTTTACAGACAAATAAAAAAACTGGA‑GTTTCCATGATGAGCGAAAGCAATAAACAACAAGGAGGGAAt                                                                               >  2:555915/1‑150 (MQ=255)                                                                    tCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAAACAAAAACAAAAACTTGA‑GTTGTCACGATGAGTGTAAGGAATAAAAAACATGCAGTGAAAAAGtt                                                                          >  2:461345/1‑151 (MQ=255)                                                                     cACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTg                                                                         <  1:320952/151‑1 (MQ=255)                                                                     cACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGAGAATTTTTTCAGACAAAAACAAAAAATGGA‑GTGGGCATGGTGAGCGGAAGGAAAAAGCAACAGGCAGGGAaaaa                                                                             >  2:347539/1‑147 (MQ=255)                                                                          aCTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAAACCAAAAAAAAAAAATGA‑GTTTGCCTGATGAGGGGAAGTAATAAACAAAAAGCAGGGAaaaaagttaa                                                                       >  2:61701/1‑142 (MQ=255)                                                                              ttAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAAAATAAAAAAAATGGAGGTTGCC‑TGATGAGGGGAAGGAATAAACAAAAAGCAGGGGATAAATTGAAAGAGAtt                                                                >  1:447528/1‑151 (MQ=255)                                                                                    aCATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGTTAAATAAAAAAAATTGA‑GGTGCCCTGATGGGGGGAAAAAAAAAACAACAGGCAGTTGATAAATTggaaaagaatatcccc                                                          >  1:254717/1‑136 (MQ=255)                                                                                     cATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTa                                                         <  2:812490/151‑1 (MQ=255)                                                                                           cACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCCACAAAAACAAAAACTGGA‑GGTGGCCTGATTAGGGGAAGAAAAAAGCAACAGGCAGGGGAAAAGTTGACAGAGATTGGCGCTAACttt                                                    >  1:244608/1‑150 (MQ=255)                                                                                           cACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTa                                                   <  1:182568/151‑1 (MQ=255)                                                                                           cACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAATTTTTTCAAACAAAAACAAAAACTGGG‑GTTGGGTAGGTTAGAGGAAGTAATAAAAAAAAAGGAGGGAATAAGTTGAAAGAGGTTGTGGGTAAAttt                                                    >  2:386291/1‑150 (MQ=255)                                                                                             cGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAAAACAAAAACTGGA‑GTTGGCATGATGAGCGAAAGAAATAAACAACAAGCCGGGAATAAGTTGACAGAGATTGTTGCTAACTTTTcc                                                 >  1:474767/1‑151 (MQ=255)                                                                                                 ggCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGAAAAATAAAAAAAATGGA‑GTTGCCAAGATGAGGGAAAGTAATAAACAACAGGCAGTGAATAAGGTGACAGAGATTGTCggtcaattttccccc                                              >  2:64888/1‑136 (MQ=255)                                                                                                     aaaaaTATAAACCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGAt                                         <  1:784946/151‑1 (MQ=255)                                                                                                               aaCCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCACACAAATACAAAAAATGGA‑GTTGGGCTGGTGAGGGGAAGAAATAAAAAACAAGCAGGGAATAAGGTGACAGAAATTGTAGCTGACTTTAACCCCCTGATTTTTAccc                                 >  2:368758/1‑149 (MQ=255)                                                                                                                                   tttGTGGCGGTGAATTTTTCCAGAAAAATACAAAAACTGGA‑GTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTACCTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGtggtgg            <  2:546875/150‑1 (MQ=255)                                                                                                                                       tgGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATaa        >  1:478380/1‑150 (MQ=255)                                                                                                                                        ggCGGTGAATTTTTCCAAACAAATACAAAAACGGGA‑TTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATGGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATaaa       <  2:616127/150‑1 (MQ=255)                                                                                                                                            gTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTaa   >  1:52482/1‑150 (MQ=255)                                                                                                                                             tGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTaaa  <  2:703814/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                              AGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAA  >  NC_000913/3206286‑3206575

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