Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,206,4350CT39.3% 25.1 / 18.6 28intergenic (‑179/‑28)ttdR/ttdADNA‑binding transcriptional activator Dan/L(+)‑tartrate dehydratase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (0/17);  new base T (11/0);  total (11/17)
Fisher's exact test for biased strand distribution p-value = 4.66e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGA  >  NC_000913/3206298‑3206566
                                                                                                                                         |                                                                                                                                    
aaaaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCa                                                                                                                                    >  2:433375/1‑140 (MQ=255)
             tCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACaaaa                                                                                                                       <  2:903919/140‑1 (MQ=255)
                      tGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCACACAAATAAAAAAAAAGGAggt                                                                                                             >  1:408226/1‑138 (MQ=255)
                        gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGc                                                                                                           <  1:829405/140‑1 (MQ=255)
                        gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAATACAAAAAATGGAGGTgg                                                                                                           >  2:811784/1‑139 (MQ=255)
                             gtatTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatga                                                                                                      <  1:408797/139‑1 (MQ=255)
                             atatTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGTCatga                                                                                                      <  1:1009530/140‑1 (MQ=255)
                             atatTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAAAAAAAAAAATAGAGTTGGCatta                                                                                                      >  2:823842/1‑138 (MQ=255)
                                ttcattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATAAAAAAAAAGGAGTTGGCgtgatgg                                                                                                   >  1:246845/1‑139 (MQ=255)
                                  cattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAGAAAAACAAAAACCGGGGATTCGATCATGAGc                                                                                                 >  1:57628/1‑140 (MQ=255)
                                    ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                               <  2:977256/140‑1 (MQ=255)
                                      cTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGaaa                                                                                             <  1:734347/140‑1 (MQ=255)
                                        tATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAAAAAAAAAAAAGGAGGTGCCCTGATGAGCGaaaaa                                                                                           >  1:651378/1‑138 (MQ=255)
                                         aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTTAGACAAAAACAAAAAATGGAGTTGCCATGATGAGGGAAAAta                                                                                          >  2:685694/1‑140 (MQ=255)
                                         aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGAAAAATAAAAAAAATAGGGGTGCCATTATGAGGGGAAGta                                                                                          >  1:191060/1‑140 (MQ=255)
                                         aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGGTGCCATGATGAGCGaaaaaa                                                                                          >  1:1055292/1‑137 (MQ=255)
                                         aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAAAAAAAAAAAAAAATGGAGTTGCCATGATGAGCGAAAggt                                                                                          >  2:408797/1‑138 (MQ=255)
                                               gcgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCACACAAATACAAAAACTGGGGTTGTCATGATGAGCGAAAGgaataaaa                                                                                    >  1:309768/1‑138 (MQ=255)
                                                  ttGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAAc                                                                                 <  2:25523/140‑1 (MQ=255)
                                                     aTTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAgg                                                                              <  2:18040/140‑1 (MQ=255)
                                                     aTTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAgg                                                                              <  2:329422/140‑1 (MQ=255)
                                                      ttCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGc                                                                             <  2:993963/140‑1 (MQ=255)
                                                            aaCTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGaa                                                                       <  1:433375/140‑1 (MQ=255)
                                                            aaCTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGaa                                                                       <  1:390510/140‑1 (MQ=255)
                                                                 ttAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGGTGCCATGATGAGCGAAAAAAAAAAACAACAAGCAGTGAATaaat                                                                  >  1:71332/1‑138 (MQ=255)
                                                                         atatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAAAAAAAAAAAAATGGGGGTGGCATGATGAGGGAAAATAAAAAACAACAAGGAGGGAAAAAATTTAAagag                                                          >  2:686534/1‑140 (MQ=255)
                                                                             tCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAAAAAAAAAAAAATGGAGGTGCCATGATGAGCGGAAGTAATAAACAAAAAGCAGTGAATAAAGTGACAGAGAttt                                                      >  2:307465/1‑139 (MQ=255)
                                                                                cGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCAGACAAAAACAAAAAATGGGGGTTGCCTGATGAGCGGAAGAAATAAAAAAAAAGCAAGGAAAAAAATGACAGAGATTGttg                                                   >  2:122406/1‑138 (MQ=255)
                                                                                                             tAGTGGTAGTTTGGGGCGGTGAATTTTTCCAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGa                      <  1:276371/140‑1 (MQ=255)
                                                                                                                    aGTTTGTGGCGGTGATTTTTTCCAGAAAAATAAAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtg               <  2:1055292/140‑1 (MQ=255)
                                                                                                                       ttGTGGGGGTGAATTTTTCCAGAAAAAAACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACGGGAAGTGAATAAGTTGACGGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtgatg            <  2:552671/140‑1 (MQ=255)
                                                                                                                        tgtgGCGGTAAATTTTTCCAAAAAAATAAAAAACTGGAGTTACCATGATGAGCAAAAGAAATAAGCAACAGTCAGTAAATAAGTTGACAGAGATTTTCGCTTACTTTACCGCCATGATTTCTACCCGAATGCCTGATGAc          <  2:426137/140‑1 (MQ=255)
                                                                                                                            gCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgtg       <  1:823842/140‑1 (MQ=255)
                                                                                                                            gCGGTGAATTTTTCCAGAAAAATACAAAAACGGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGTTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgtg       <  2:722140/140‑1 (MQ=255)
                                                                                                                                 aaaTTTTTCAAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGTTAATTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGa  <  2:209792/139‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                    
AAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGA  >  NC_000913/3206298‑3206566

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: