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

AAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTG  >  NC_000913/3206298‑3206561
                                                                                                                                          |                                                                                                                               
aaaaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCa                                                                                                                               >  2:157842/1‑140 (MQ=255)
    cATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGTAATTTTTcaaaaa                                                                                                                            >  1:758910/1‑135 (MQ=255)
                   ctctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAATGGa                                                                                                            >  2:810594/1‑140 (MQ=255)
                   ctctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTTAGACAaaaaaaaaaaaagga                                                                                                            >  2:25728/1‑135 (MQ=255)
                        gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGGTgg                                                                                                       >  2:906544/1‑139 (MQ=255)
                          gcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCa                                                                                                     <  1:751175/140‑1 (MQ=255)
                                    ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                           <  1:1108083/140‑1 (MQ=255)
                                        tATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGTTGCCATGATGAGGGaaaaa                                                                                       >  1:345070/1‑138 (MQ=255)
                                        tATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAAAAAAAAAAAAGGAGGTGGCATGATGAggggaaaa                                                                                       >  2:664874/1‑133 (MQ=255)
                                           tttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaat                                                                                    <  1:664874/140‑1 (MQ=255)
                                           tttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaat                                                                                    <  1:810594/140‑1 (MQ=255)
                                           tttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaat                                                                                    <  1:932163/140‑1 (MQ=255)
                                           tttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaat                                                                                    <  2:1183529/140‑1 (MQ=255)
                                             tCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTTGATTTTTTCCGACAAATAAAAAAAATCGAGTTGGCATGGTGAGGGAAAGtaaaaa                                                                                  >  1:759587/1‑140 (MQ=255)
                                               gcgTTGATTCACGAACTGTTAATTACATATTCACGGTTGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGTGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGc                                                                                >  1:1122183/1‑140 (MQ=255)
                                                 gTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCaa                                                                              <  2:303951/140‑1 (MQ=255)
                                                  ttGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAATGGAGGTGCCATGATGAGGGGAAGAAATAAACaaa                                                                             >  2:678779/1‑139 (MQ=255)
                                                     aTTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAgg                                                                          <  1:678779/140‑1 (MQ=255)
                                                      ttCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGc                                                                         <  2:1133223/140‑1 (MQ=255)
                                                              cTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATa                                                                <  1:157842/140‑1 (MQ=255)
                                                                    aTTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCAAACAAAAACAAAAAATGGGGTTGCCATGATGAGGGAAAAGAAAAAAAAAAAAGCAGTGGATAAgatta                                                          >  2:974697/1‑136 (MQ=255)
                                                                           atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAATAAAAAAAAAGGAGGTGCCATGATGAGCGGAAGGAAAAAACAAAAAACAGAGAAAAAAATGACAgagaa                                                   >  2:1106776/1‑139 (MQ=255)
                                                                           atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAATAGAGGTTGCATGATGAGGGAAAGTAAAAAACAAAAAGCAGAGAATAAAATGACAgagaa                                                   >  2:870587/1‑139 (MQ=255)
                                                                                 ggTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAAAAAAAAAAATGGAGTTGGCATGATGAGCGAAAGGAAAAAACAAAAAGCAAGGAAAAAATTGACAGAGAATGTCGc                                             >  2:804467/1‑140 (MQ=255)
                                                                                      cAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTTAGACAAAAAAAAAAAATGGAGTTGGCATGATGGGGGAAAATAAAAAACAAAAAGCAGAGGAAAAAGTTGCAGAGAATGTTGCTaaat                                        >  2:1005455/1‑138 (MQ=255)
                                                                                         aagatatAAACCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCAAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTa                                     <  1:71640/137‑1 (MQ=255)
                                                                                           aaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATAAAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTAcc                                   <  1:870587/140‑1 (MQ=255)
                                                                                                       caTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTAAGACAAAAACAAAAAATGGAGGTGCCATGATGAGCGGAAATAAAAAAAAACAAGCAGTGAAAAAATTTGCAAAGATTGTCGCTAAATTTACCCCCCTGAttttt                       >  2:71640/1‑138 (MQ=255)
                                                                                                       caTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCt                       <  1:974697/140‑1 (MQ=255)
                                                                                                            ttAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCg                  <  2:447202/140‑1 (MQ=255)
                                                                                                             tAGTGGTAGTTTGTGGTGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGa                 <  2:1122183/140‑1 (MQ=255)
                                                                                                              aGTGGTAGTTTGTGGGGGTAATTTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAATTTTACCGCCATGATTTCTACCCGaa                <  2:345070/140‑1 (MQ=255)
                                                                                                                      tttGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtgat        <  1:804467/140‑1 (MQ=255)
                                                                                                                            gCGGTAATTTTTTCGAGAAAAATAAAAAGTCTGGAGTTGCCATGAGGAGGAAAAGAAATAACCACCGGGCAGTGAATAAGTTGACAGAGATTGTCGCAAGTTTAACCGCCATGATTTCTACCCGAATGCCTGATGACgtg  <  2:426020/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                               
AAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTG  >  NC_000913/3206298‑3206561

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: