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

CAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGAC‑AAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCT  >  NC_000913/3206297‑3206581                                                                                                                                           |                                                                                                                                                    cAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGAC‑AAATACa                                                                                                                                        <  1:457590/151‑1 (MQ=255)  aaaaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGAC‑aaa                                                                                                                                            >  2:359048/1‑146 (MQ=255)    aaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTATCCAGCC‑AAAAACaaa                                                                                                                                      >  1:173590/1‑150 (MQ=255)      cATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCAAC‑Aaaaaaaaaaa                                                                                                                                    >  1:581731/1‑150 (MQ=255)          gCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGAC‑AAATACAAAAACTGGa                                                                                                                               <  2:923262/151‑1 (MQ=255)             aTTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGAA‑AAATACaaaaaaaag                                                                                                                                >  2:507825/1‑143 (MQ=255)             aTTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGAC‑AAATACAAAAACTGGAGtt                                                                                                                            <  1:947360/151‑1 (MQ=255)                    gCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAA‑AAATAAAAAAAATGGAGGTGcc                                                                                                                         >  2:187586/1‑147 (MQ=255)                     ctctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGTATTTTTTTAGAC‑AAAAACAAAAAATTGAGTTggcgagg                                                                                                                     >  1:445828/1‑144 (MQ=255)                               atatTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGAC‑AAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGt                                                                                                          <  2:547266/151‑1 (MQ=255)                                   tcattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAA‑AAAAACAAAAAAAGGAGTTGTCATGATGAGCGGAAGaaa                                                                                                        >  1:836230/1‑149 (MQ=255)                                   tcattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGCGGATTTTTTTACACAAAATACAAAAA‑TAAAGTTTGCATGTAGAGTGAAAGtaa                                                                                                        >  1:806759/1‑149 (MQ=255)                                    cattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCAAC‑AAATAAAAAAAATGGAGTTGCCATGATGAGGGGAAGaaa                                                                                                        >  2:752077/1‑148 (MQ=255)                                    cattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGAC‑AAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaataa                                                                                                     <  2:939248/151‑1 (MQ=255)                                       tcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGAC‑AAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGc                                                                                                   <  2:229762/150‑1 (MQ=255)                                           aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAAA‑AAAAAAAAAAAAAGGAGTTTGCCTGATTAGCGAAAGGAAAAAACa                                                                                                  >  2:456178/1‑147 (MQ=255)                                             tttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCATAA‑AAAAACAAAAAACTGAGGTGCGCTGGAGAGGGAACGTAATAAAAAACa                                                                                               >  1:282790/1‑148 (MQ=255)                                              ttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGAC‑AAAAAAAAAAAACGGAGGTGGCATGATGAGCGAAAAGAATAAacaaca                                                                                               >  2:552201/1‑147 (MQ=255)                                              ttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGAC‑AAATACAAAAACTGGAGTCGCCATGATGAGCGAAAGTAATAAGCAACAGGc                                                                                            <  1:33761/150‑1 (MQ=255)                                               tCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGAC‑AAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCa                                                                                           <  2:322012/150‑1 (MQ=255)                                               tCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGAC‑AAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCa                                                                                           <  2:927694/150‑1 (MQ=255)                                                   gTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCTGAC‑AAATAAAAAAAATTGAGGTGCCATGTTGAGAGAAAGGAATAAAAAACAAGCAGTg                                                                                        >  1:744581/1‑149 (MQ=255)                                                         tCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGAC‑AAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGt                                                                                 <  2:204983/150‑1 (MQ=255)                                                          cACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGGAATTTTTCTTAA‑AAAAAAAAAAAATGTAGTTTCCCTGATAGGCGAAAGAAATAAGAAAAAAGACATGAATAAGt                                                                                 >  1:53720/1‑149 (MQ=255)                                                          cACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTCTGTGGCGGTGAATTTTTCCAGAC‑AAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTg                                                                               <  2:928391/151‑1 (MQ=255)                                                                  gTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGAC‑AAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACagag                                                                         <  2:913927/149‑1 (MQ=255)                                                                         aCATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTTAGAC‑AAAAAAAAAAAATGGAGTTGGCATGGTGAGCGAAAGGAATAAAAAACAAGCAGGGAATAAAATGGAAGAGAATGttg                                                                  >  1:89011/1‑147 (MQ=255)                                                                              ttCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTTAAATTTTTTCCAC‑AAAAAAAAAAAATGGAGTTTGTCTGATGAGGGAAAAGAAAAAACAAAAAGAAGGGAATAAATTGTCAGAGAATTTCGCt                                                                >  1:354058/1‑146 (MQ=255)                                                                               tCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGAC‑AAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACtt                                                           <  1:700018/150‑1 (MQ=255)                                                                                        cAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGAC‑AAAAACAAAAACTGGAGTTGGCATGATGAGGGAAAGGAATAAAAAAAAAGCAGGGAATAAGATGAAAGAGATTTTCGCTAACTTTACCGCCa                                                   >  2:147985/1‑149 (MQ=255)                                                                                         aaaaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGAC‑AAATAAAAAAACTGGAGTTGGCATGATGAGCGAAAATAATAAACAAAAGGAAGTGAATAAATTGACAGAGATTGTCGCTAACTTTACCGCCCTg                                                 >  1:323142/1‑150 (MQ=255)                                                                                              atatAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGAATTTTTCACAA‑AAAAAAAAAAAATGGAGTTGTCCTGATGAGGGGAAGTAATAAAAAAAAGACAGGGAAAAAAATGACAGAGATTGTTGATAAATTTACCGCCATGAttttt                                           >  1:554376/1‑149 (MQ=255)                                                                                                  aaaaCCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGAC‑AAAAAAAAAAAATGGAGGTGGCATGATGGGGGAAAGTAATAAACAACAAGCAGGGAAAAAAATTACAGAGATTGTCGCCTACTTTACCGCCCTGATTTTTAacc                                       >  2:547028/1‑151 (MQ=255)                                                                                                      ccACTTTTTTAGTGGTATTTTGTGGCGGTGAATTTTTCAAGAA‑AAATAAAAAAACTGGAGTTGCCATGATGAGCGAAATAAATAACCAACGGGCAGTAAATAAGTTGACAAAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAAt                                   <  1:552168/151‑1 (MQ=255)                                                                                                            ttttAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAC‑AAAAACAAAAAATGGAGGTGCCCTGATGAGCGGAAGGAATAAACAACAAGCAGTGAAAAAGTTGACAGAGAATGTTGCTAACTTTTCCGCCCTGATTTTTACCCCAAAGcc                                >  2:130214/1‑148 (MQ=255)                                                                                                              ttAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAC‑AAAAACAAAAAAAGGAGGTGGCATGATGAGCGAAAAGAATAAACAAAAGGGAGGGAATAAAGTGAAAGAGATTGTCGCTAAATTTTTCCCCACGATTTTTTCCCGAATGcc                                >  2:561238/1‑146 (MQ=255)                                                                                                                aGTGGTAGTTTGTGGCGGTGAATTTTTTCCGAC‑AAAAAAAAAAAACGGAGGTGGCCTGATTAGGGAAAGGAAAAAACAAAAAGCAGGGAATAAAATGAAAGAGATTGTCGCTAACTTTACCGCCCTGATTTTTTCCCCAAAGcc                                >  2:525913/1‑144 (MQ=255)                                                                                                                  tGGTAGTTTGTGGCGGTGAATTTTTTCCGAC‑AAAAACAAAAACTGGGGGTGGCCTGATGAGCGGAAGTAATAAACAACAGGCAGTGGATAAATTGACAGAGATTGTCGCTAACTTTACCGCCCTGATTTTTACCCGAATGCCCTATGACgt                       >  2:157563/1‑151 (MQ=255)                                                                                                                        tttGTGGCGGTGATTTTTTCCAGAC‑AAATACAAAAACTGGATTTGCCATGTTGAGCGAAAGTAATAAGCAACAGGCAGTGAATAATTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgtggt                    <  2:279790/148‑1 (MQ=255)                                                                                                                        tttGTGGCGGTGAATTTTTCCAGAC‑AAATAAAAAAACTGGATTGGCCAGGATGGGGGAAAGTAAAAAGCAACAGGCAGTGAAAAAGTGGAAAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACg                        <  2:500455/144‑1 (MQ=255)                                                                                                                        tttGTGGCGGTGAATTTTTCCAAAC‑AAATACAAAAACTGAATTTCCCATAATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAATTGGAAAGAGTTTTTCGCTAATTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgtg                      <  1:260947/146‑1 (MQ=255)                                                                                                                        tttGTGGCGGTAAATTTTTCCAGAC‑AAATACAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgtggt                    <  2:667771/148‑1 (MQ=255)                                                                                                                           gtgGCGGTGAATTTTTTCCGAC‑AAAAAAAAAAAATGGAGGTGGCCTGATGGGGGGAAGTAATAAACAAAAAGCAGTGAATAAGGTGACAGAGAATGTCGCTAACTTTACCGCCATGATTTTTAACCGAAAGCCTGATGACGtggggg                  >  1:476673/1‑147 (MQ=255)                                                                                                                             gggggTGAATTTTTCAAGAA‑AAATAAAAAAACTGGAGTGCCCATGATGGGCAAAAGTAATAAGCAACAGGCAGTAAATAATTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATaaa             <  1:175587/147‑1 (MQ=255)                                                                                                                             ggCGGTGAATTTTTCCAGAC‑AAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAAc            >  2:698372/1‑151 (MQ=255)                                                                                                                              ggggtaatttttcCCAAAA‑AAATACAAAAACGGAAGTTGCCAGGATGAGCGAAAGTAATAAGCAACAGGCAGTAAATAATTTGACAGAGATTTTCGCTATCTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACt           <  1:862325/143‑1 (MQ=255)                                                                                                                               cGGTGAATTTTTCCAGAC‑AAGTACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACt           >  1:513459/1‑150 (MQ=255)                                                                                                                                   gAATTTTTCCAGAC‑AAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAAc      >  1:530795/1‑151 (MQ=255)                                                                                                                                       ttttCCAGAC‑AAATACAAAAACTGAAGTTGCCATGATGAGCGAAAGCAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCt  >  1:785884/1‑151 (MQ=255)                                                                                                                                        tttCCAAAA‑AAATAAAAAAACGGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATGGTCTCTAATTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCt  <  1:873291/150‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                    CAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGAC‑AAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCT  >  NC_000913/3206297‑3206581

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