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

AGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCT‑AACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAA  >  NC_000913/3206286‑3206576                                                                                                                                                       |                                                                                                                                               aGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTcc                                                                                                                                               <  2:609126/151‑1 (MQ=255)      aTTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTCCAGAc                                                                                                                                           <  1:292854/150‑1 (MQ=255)            cAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTCCAGACAAATAc                                                                                                                                     <  2:69527/150‑1 (MQ=255)                   ttGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTCCAGTCCaaaaaaaaaa                                                                                                                                >  2:384483/1‑148 (MQ=255)                    tGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTTTAGAAAAAAAAAAAAAATg                                                                                                                             >  2:658952/1‑150 (MQ=255)                    tGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTCCAGACAAATACAAAAACTg                                                                                                                             <  2:300674/150‑1 (MQ=255)                       cATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑ATTTTTCAGACAAAAACAAAAAATgg                                                                                                                            >  2:287239/1‑148 (MQ=255)                            ccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTCCAGACAAATACAAAAACTGGAGTTGCCa                                                                                                                    <  2:141355/151‑1 (MQ=255)                               gCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTCCAGACAAATACAAAAACTGGAGTTGCCatg                                                                                                                  <  2:383927/150‑1 (MQ=255)                                      cgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCTAGA‑AAATACAAAAACTGGAGTTgtcgaggt                                                                                                                >  1:190946/1‑138 (MQ=255)                                       gcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTTCAAACCAAAAAAAAAAAATGAGGTGCCCTGATGAggg                                                                                                           >  2:702078/1‑147 (MQ=255)                                            attcattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGATGAA‑TTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaa                                                                                                    <  1:69691/151‑1 (MQ=255)                                                attcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTTCAAACAAAAAAAAAAAATGGAGGTGGCATGATTGGGGGAAGtaataa                                                                                                 >  2:265250/1‑150 (MQ=255)                                                      aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑ATTTTTCAAACAAATAAAAAAAAATGAGGTGTCATTATGAGGGGAAGGAAAAAACaaaa                                                                                           >  2:655232/1‑148 (MQ=255)                                                            gcgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAG‑TTTTTTTAGAAAAAAACAAAAAATTGAGGTGTCATGATGGGGGGAAATAATAAGCAAAAggacgt                                                                                     >  2:641513/1‑146 (MQ=255)                                                             cgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTTCAGACAAAAACAAAAAATGGAGTTTCCATGATGAGCGAAAGAAATAAAAAACAAGCAGTg                                                                                    >  2:421462/1‑150 (MQ=255)                                                                tGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAAt                                                                                 <  1:288521/150‑1 (MQ=255)                                                                  aTTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTTCCGACAAATACAAAAACTGGAGTTGCCCTGATGAGCGGAAGTAATAAACAACAGGCAGTGAATaa                                                                               >  1:4038/1‑150 (MQ=255)                                                                     cACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGtt                                                                            <  1:500413/150‑1 (MQ=255)                                                                     cACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGtt                                                                            <  2:374975/150‑1 (MQ=255)                                                                          aCTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACaga                                                                      <  1:84774/151‑1 (MQ=255)                                                                             gTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACagag                                                                     <  2:360001/149‑1 (MQ=255)                                                                                          tCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCT‑AACtt                                                      <  1:710993/150‑1 (MQ=255)                                                                                           cACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCT‑AACTTTa                                                    <  1:597743/151‑1 (MQ=255)                                                                                           cACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAA‑TTTTTTCAAAAAAAAAAAAAAAATGGAGGTGCCCTGATGAGGGAAAGTAATAAACAAAAAGCAGAGAATAAAGTGAAAGAGAATGTCGGT‑AAAttt                                                     >  2:423583/1‑150 (MQ=255)                                                                                                tGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAA‑ATTTTTCAAACAAAAAAAAAAAATGGAGTTGTCATGATGAGGGAAAGGAATAAACAACAGGCAGAGAAAAAATTTACAAAGATTGTCGCT‑AACttttacccc                                               >  2:268805/1‑145 (MQ=255)                                                                                                    aaaaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTCCAGAAAAAAACAAAAAAAGGAGTTGCCATGATGGGCGAAAGAAATAAGCAACAGGCAGGGAATAAGATGACAGAGAATGTCGCT‑AACTTTACcccc                                               >  1:647212/1‑147 (MQ=255)                                                                                                          tataAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑ATTTTTCAGACAAAAACAAAAAACGGGGGTGGCCTGATGAGCGGAAAGAATAAACAAAAAGCAGTGAATAAATTGACAGAGATTTTCGCT‑AACTTTAACCCCCTGAttttt                                      >  1:671521/1‑148 (MQ=255)                                                                                                           ataAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTTCAGACAAAAACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAAGCAGTGAATAAGTTGACAGAGATTGTCGCT‑AACTTTACCGCCCTGAttttt                                      >  1:664487/1‑147 (MQ=255)                                                                                                           ataAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTTCAGACAAAAAAAAAAAATGGAGTTGCCATGATGAGCGGAAGTAAAAAGCAAAAGGCCGTGAAAAAATTGACAGAGATTGTCGCT‑AAATTTTACGCCATGAttttt                                      >  1:75450/1‑147 (MQ=255)                                                                                                            taAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑ATTTTTCAGACAAATACAAAAAATGGAGTTGCCCTGATGAGGGAAAGTAATAAACAAAAGGCAGTGAATAAGTTGACAGAGATTGTTGCT‑AACTTTACCGCCATGAttttt                                      >  1:243868/1‑146 (MQ=255)                                                                                                             aaaaCCACATTTTTAGTGGTAGTTTGTGGCGGTAAA‑TTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCT‑AACTTTACCGCCATGATTTCTAcc                                   <  1:172264/150‑1 (MQ=255)                                                                                                              aaaCCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTTCCTACAAAAAAAAAAAAATGAGGTGGCATGATTGGGGGCAGTAATAAGCAAAAAGGAGGGGATAAATTTGCAAAGATTGGTGCTAAATTTTACCCCCA‑GATTTCTAccc                                  >  1:617717/1‑150 (MQ=255)                                                                                                                 ccACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTCCAGACAAATACAAAAACGGGAGTGGCCAGGATGAGCGAAAGTAATAAGCAACAGGCAGTAAATAAGTGGACAGAGATTGTCGCT‑AACTTTACCGCCATGATTTCTACCCGa                                <  2:688847/149‑1 (MQ=255)                                                                                                                   acaTTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTTCAGACAAAATTAAAAAATGGAGTTGGCATGTGGAGGGGAAAGAAAAAGAAAAAAGAAGAGAAAAAGTTGACAGAGATTTTCGCT‑CACTTTACACCCCTGATTTCTTCCCGAAAGc                            >  2:98754/1‑151 (MQ=255)                                                                                                                       ttttAGTGGTAGTTTGTGGCGGGGAA‑TTTTTTCCGACAAAAAAAAAAAATGGAGGTGGCATGATGGGGGGAAGTAATAAACAAAAGGCAGGGAAAAAATTGAACGGGAATGTTGGT‑AAATTTTCCGCCCAGATTTTTACCCGAAAGCCtg                         >  2:302564/1‑150 (MQ=255)                                                                                                                                    ttGTGGGGGTGAA‑TTTTTCCAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGTT‑AACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGa            <  2:543194/150‑1 (MQ=255)                                                                                                                                    ttGTGGCGGTGAA‑TTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCT‑AAC‑TTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGAt           <  1:98820/150‑1 (MQ=255)                                                                                                                                       tgGCGGTGAA‑TTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCT‑AACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATaa         >  2:71931/1‑150 (MQ=255)                                                                                                                                               aa‑TTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCT‑AACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTaaaa  <  1:37186/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                               AGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAA‑TTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCT‑AACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAA  >  NC_000913/3206286‑3206576

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