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

TTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTC‑GCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTA  >  NC_000913/3206293‑3206582                                                                                                                                               |                                                                                                                                                     ttCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGaaaa                                                                                                                                               >  2:211293/1‑150 (MQ=255)    gCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGAaaaaaa                                                                                                                                            >  1:544769/1‑148 (MQ=255)      aaaaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACaa                                                                                                                                         <  2:801091/151‑1 (MQ=255)          cATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTAGTAGTTTGTGGCGGTGAATTTTTTCCGACAaaaaaaaaaa                                                                                                                                      >  1:405866/1‑150 (MQ=255)                 aTTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAAAAAAAAAAAAATGGAGt                                                                                                                               >  1:853882/1‑150 (MQ=255)                 aTTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGtt                                                                                                                              <  1:724582/151‑1 (MQ=255)                     ccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAAAAAATAAAAAAAAAAGAGTTGcc                                                                                                                           >  1:103290/1‑150 (MQ=255)                        gCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatga                                                                                                                       <  2:836721/151‑1 (MQ=255)                         ctctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAATACAAAAACTGGAGGTGcc                                                                                                                           >  1:791501/1‑146 (MQ=255)                             ggCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatgatg                                                                                                                     <  1:298492/148‑1 (MQ=255)                              gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCg                                                                                                                 <  2:624061/151‑1 (MQ=255)                                gcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAAAGGAGTTGGCATGATGAGCg                                                                                                                 >  1:624479/1‑149 (MQ=255)                                     attcattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAACGGAGTTGGCATGGTGAggggaa                                                                                                              >  2:538626/1‑142 (MQ=255)                                        cattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaataa                                                                                                       <  2:686478/151‑1 (MQ=255)                                              tATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCaa                                                                                                   <  2:372351/149‑1 (MQ=255)                                               aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGTTTCCATGATGAGCGGAAGGAATAAACaa                                                                                                   >  1:519191/1‑148 (MQ=255)                                               aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAg                                                                                                <  2:351083/151‑1 (MQ=255)                                                 tttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAgg                                                                                               <  1:37132/150‑1 (MQ=255)                                                    cgcgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGt                                                                                           <  1:838376/151‑1 (MQ=255)                                                     gcgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTg                                                                                          <  2:543582/151‑1 (MQ=255)                                                          gATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATa                                                                                      <  1:380712/150‑1 (MQ=255)                                                           aTTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAg                                                                                    <  2:876581/151‑1 (MQ=255)                                                              cACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAACAAAAAATGGAGTTTCCCAGGTTAGCGGAAATAATAAAAAAAAGGCAGGGAAAAAGtt                                                                                  >  1:628393/1‑150 (MQ=255)                                                              cACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAAAAAAAACAAAAAATGGAGTTGTCATGGTGAGGGGAAGTAATAAACAACAAGCAGGGAATaa                                                                                     >  2:717473/1‑147 (MQ=255)                                                                cGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGAc                                                                               <  2:590820/151‑1 (MQ=255)                                                                 gAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAAAGGAGTTGCCCTGATGAGCGAAAAAAATAAACAAAAAGCAGGGAATAAGtt                                                                                  >  1:450471/1‑147 (MQ=255)                                                                        tAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCAGACAAAAAAAAAAACAGGAGTTGCCCTTATGAGCGAAAATAATAAAAAAAAAGCAGGGAAAAAGGTGAAAGAGAtt                                                                        >  1:511938/1‑150 (MQ=255)                                                                          aTTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGAGAATTTTTTCCAACAAAAAAAAAAAAAGGGGTTGCCAAGATGAGCGGAAGTAATAAACAACAAGCAGTGAAAAAATTGACAGAGAtt                                                                        >  1:798673/1‑148 (MQ=255)                                                                          aTTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTAGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTc                                                                     <  2:405815/151‑1 (MQ=255)                                                                            tACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAAACAAAAACAAAAAATGGAGTTTCCATGATGAGGGAAAAAAAAAAAAAAAAGGGAGGGAATAAAATGAAAGAGATTGtgg                                                                    >  2:9515/1‑148 (MQ=255)                                                                             aCATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTC‑Gc                                                                  <  2:255956/150‑1 (MQ=255)                                                                               atatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAAAACAAAAAATGGAGTTGCCATGATGAGGGGAAGGAATAAACAAAAGGCAGGGAATAAATTTACAGAGATTGTC‑GCt                                                                 >  1:455038/1‑149 (MQ=255)                                                                                tatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGGGTTGCCATGATGAGCGGAAGAAATAAACAAAAAGCAGTGAATAAGTTGACAGAGATTGTC‑GCTaa                                                               >  1:109681/1‑150 (MQ=255)                                                                                tatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAAAACAAAAAATGGAGGTGGCATGATGAGGGAAAGGAATAAACAAAAGGCAGGGAATAAGTTGACAGAGATTTTC‑GCTaa                                                               >  1:532567/1‑150 (MQ=255)                                                                                 atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGTTGCCCTGATGAGCGAAAGGAATAAACAAAAAGCAGTGAATAAATTGACAGAGATTGTC‑GCTaa                                                               >  1:596582/1‑149 (MQ=255)                                                                                   tCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTC‑GCTAACtt                                                            <  2:821279/150‑1 (MQ=255)                                                                                         tGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAAAAAAAAAAAAATGGAGGTGTCCTTCTGAGGGAAAATAAAAATAAAAAAGCAGTGGAAAAATTGATAGAGATTGTT‑GCTAACTTTACCGcc                                                     >  2:490198/1‑151 (MQ=255)                                                                                            cAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAACAAAAACTGGAGGTGGCCTGATGAGCGGAAGAAAAAAACAACAGGCAGTGAAAAAATTGACAGAGAATGTCGGCAAATTTTACC‑CCATg                                                  >  1:833437/1‑151 (MQ=255)                                                                                             aaaaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGCGAATTTTTTCAGACAAAAACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAAAAGGCAGAGAATAAATTGACAGAGATTGTT‑GCTAAATTTAACGccc                                                    >  2:778441/1‑147 (MQ=255)                                                                                                     taAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAAAAAAAAACAAAAAATGGGGTTTGCATGATGAGGGAAAAGAATAAACAACAAGGAGTGAATAAAATGACAGAGATTTTC‑GCTAACTTTAACCCCATGAttt                                              >  2:819139/1‑146 (MQ=255)                                                                                                       aacccACATTTTTAGTGTTAGTTTGTGGGTTTAAATTTTACCAGACAAATAAAAAAACTGAAGTTGCCAGGATGAGGAAAAGTAATAAGCAACGGGCAGTGAAAAGGTTGACAGATATTGTC‑GTTATCTTTACCGCCATGATTTCTACCCg                                       <  2:4824/148‑1 (MQ=255)                                                                                                        aaCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAAATGGGGTTGCCATGATGAGCGGAAGAAATAAACAACAGGCAGGGAATAAGTTGACAGAGATTGTC‑GCTAACTTTACCGCCCTGAttttt                                            >  1:185562/1‑143 (MQ=255)                                                                                                            acaTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTCCCAGAAAAATACAAAAACTGGAGTGGCCATGATGAGCGAAAGTAATAAGCAACAGGCATTGAATAAGTTGACAGAGATTGTC‑GTTACCTTTACCGCCATGATTTCTACCCGAAt                                    <  2:29994/149‑1 (MQ=255)                                                                                                               tttttAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAACTGGGGGTGCCATGATGAGCGGAAGTAAAAAGCAAAAGGGAGTGAATAAAATGACAGAGATTGTT‑GCTAACTTTACCGCCCTGATTTTTTCCCCAAAGcc                                 >  1:828458/1‑149 (MQ=255)                                                                                                               tttttAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGTTGCCATGATGAGCGAAAGTAATAAACAACAGGCAGTGAATAAGTTGACAGAGATTGTT‑GCTAAATTTACCGCCATGATTTCTTCCCCAATGcc                                 >  1:211497/1‑149 (MQ=255)                                                                                                                  ttAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGGGCAAAAGAAATAAGCAACAGGCAGTAAATAAGTGGACAGAGATTGTC‑GCTAACTTTACCGCCATGATTTCTACCCGAATGCCtgat                             <  2:509978/150‑1 (MQ=255)                                                                                                                     gTGGTAGTTTGGGGCGGTGAATTTTTCCAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAAAAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTC‑CCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGAc                          <  2:277801/150‑1 (MQ=255)                                                                                                                       ggTATTTTGTGGGGGTGAATTTTTCCAGAAAAATAAAAAACCGGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTC‑GTTACCTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgtg                       <  2:33954/151‑1 (MQ=255)                                                                                                                            tttGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTC‑GCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGtggtgg                   <  2:596524/150‑1 (MQ=255)                                                                                                                                         aTTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTC‑GCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAAc       >  1:516343/1‑149 (MQ=255)                                                                                                                                             ttCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTC‑GCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTa  >  2:447446/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                                     TTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTC‑GCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTA  >  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: ‑