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

TTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTA  >  NC_000913/3206293‑3206582                                                                                                                                               |                                                                                                                                                    ttCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACaa                                                                                                                                              <  1:680964/150‑1 (MQ=255)   cGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATa                                                                                                                                           <  1:447405/151‑1 (MQ=255)        aaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACaaaa                                                                                                                                      <  2:190422/151‑1 (MQ=255)        aaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACaaaa                                                                                                                                      <  2:189819/151‑1 (MQ=255)            ttGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGACCTGTTAATTACATATTCACGGTGTCAAAAATTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGCGAAAAATTTCAGACAAAAATAATAACTg                                                                                                                                  >  1:71540/1‑151 (MQ=255)              gCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGa                                                                                                                                <  2:73733/151‑1 (MQ=255)              gCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGa                                                                                                                                <  2:550477/151‑1 (MQ=255)                 aTTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAATAAAAAAAATgggggt                                                                                                                             >  2:602872/1‑147 (MQ=255)                    cccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGc                                                                                                                           <  1:574264/150‑1 (MQ=255)                     ccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCa                                                                                                                         <  2:480291/151‑1 (MQ=255)                       tGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATgggggtgttat                                                                                                                        >  2:162797/1‑141 (MQ=255)                                     attcattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGAATTTTCCAGACAAAAAAAAAAAAAGGGGTTGCCATGGTGAggggaa                                                                                                             >  2:252730/1‑142 (MQ=255)                                     attcattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaa                                                                                                         <  2:392189/151‑1 (MQ=255)                                        cattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATTGGGGTGCCATGGTGAGGGAAAGtaaaa                                                                                                       >  1:264639/1‑148 (MQ=255)                                        cattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaataa                                                                                                      <  2:102719/151‑1 (MQ=255)                                          ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAATAAAAAAACTGGAGTTGCCATGATGAGGGAAAGtaaaaa                                                                                                      >  2:393394/1‑149 (MQ=255)                                                 tttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGCGAAATAAAAAAAACGGAGTTTGCACGTTGAGCGGAAGAAATAAGaaa                                                                                                  >  2:533792/1‑146 (MQ=255)                                                 tttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAgg                                                                                              <  2:53843/150‑1 (MQ=255)                                                   tCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGGTGCCATTATGAGGGGAAGGAATAAACaa                                                                                                  >  2:82568/1‑144 (MQ=255)                                                      cgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGa                                                                                        <  1:85267/151‑1 (MQ=255)                                                           aTTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGGTGCCATGATGAGCGGAAGGAAAAAGCAACAAGCAggg                                                                                         >  1:519750/1‑143 (MQ=255)                                                                 gAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGTTTCCCTGATGAGCGAAAGTAATAAGAAACAAGGAGGGAATAAGATGAc                                                                              >  2:384803/1‑150 (MQ=255)                                                                  aaCTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACag                                                                            <  1:137249/151‑1 (MQ=255)                                                                        tAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAAACAAAAACAAAAAATGGGGGTGGCATGATGGGGGGAAGTAATAAACAAAAGGGAGTGAATAAGATGACAAAGAtt                                                                       >  1:662077/1‑150 (MQ=255)                                                                         aaTTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGt                                                                     <  2:214940/151‑1 (MQ=255)                                                                          aTTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTc                                                                    <  2:571545/151‑1 (MQ=255)                                                                          aTTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTc                                                                    <  2:662125/151‑1 (MQ=255)                                                                              cATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTa                                                                <  2:383434/151‑1 (MQ=255)                                                                                 atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAAAAAACAAAAACTTGAGGTGCCATGGTGAGGGGAAGGAATAAAAAAAAAGCAGGGAATAAGTTGAAAGAGATTGTCGCTaa                                                               >  2:159779/1‑149 (MQ=255)                                                                                    cACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACttt                                                           <  2:267998/150‑1 (MQ=255)                                                                                             aaaaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAAACAAAAAAAAAAACTGGAGTTGTCATGATGAGGGGAAATAATAAAAAAAAGGCAGGGAATAAGTTGACAGAGATTGTTGCTAACtttttcccc                                                     >  1:342238/1‑141 (MQ=255)                                                                                                  atatAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGTTGCCATGATGAGGGGAAGTAATAAGCAACCAGCAGGGAAAAAATTGACAGAGATTTTCGCTAACTTTACCCCCATGAttttt                                            >  2:453468/1‑149 (MQ=255)                                                                                                    ataAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAACAGGGGTTGCCATGATGAGCGAAAGTAATAAGCAAAAAGCAGGGAAAAAATTGACAGAGAATGTCGCTAACTTTACCGCCCTGAttt                                              >  2:357968/1‑147 (MQ=255)                                                                                                    ataAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAAATGGAGTTGCCATTATGAGCGAAAGGAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGAttttt                                            >  2:408217/1‑147 (MQ=255)                                                                                                         aCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACCAAAACAAAAAATGGGGTTGGCATTATGAGGGGAAGGAATAAACAAAAAGGAGGGAATAAGTTGACAGAGATTGTTGCTAAATTTTACGCCATGATTTTTTccc                                        >  1:198158/1‑148 (MQ=255)                                                                                                                 tttAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACCAATACAAAAAATGGAGGTTGCATGATGAGGGGAAGGAATAAACAACAGGCAGTGAAAAAATTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTTTACCCCAAAGACtg                               >  1:231188/1‑149 (MQ=255)                                                                                                                 tttAGTGGTAGTTTGGGGCGGTGAATTTTTCAAGACAAATAAAAAAACGGGAGTTGCCATGATGAGCGAAAGTAAAAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAATTTTACCGCCATGATTTCTACCCGAATGCCtga                              <  2:573377/150‑1 (MQ=255)                                                                                                                     gTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAACTGGAGGTGGCATGATGAGGGGAAGGAATAAGGAACAGGCAGTGAAAAAATTGACAGAGGTTTTTGTTAAATTTAACGCCACGATTTCTACCCGAATGCCTGATGACg                         >  1:363567/1‑151 (MQ=255)                                                                                                                           gTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGGGAAAGTAAAAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTACCTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgtg                       <  2:486307/147‑1 (MQ=255)                                                                                                                                 ggCGGTGAATTTTTCCAGACAAATACAAAAACGGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATaaa              <  2:198179/150‑1 (MQ=255)                                                                                                                                             ttCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTa  >  2:176240/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                                    TTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTA  >  NC_000913/3206293‑3206582

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: ‑