Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,206,4440TA36.0% 23.4 / 12.2 25intergenic (‑188/‑19)ttdR/ttdADNA‑binding transcriptional activator Dan/L(+)‑tartrate dehydratase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (2/14);  new base A (9/0);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 2.69e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTA  >  NC_000913/3206307‑3206582
                                                                                                                                         |                                                                                                                                          
ccATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATAc                                                                                                                                          <  2:786976/140‑1 (MQ=255)
               gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAA‑CACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGa                                                                                                                                <  2:4094/134‑1 (MQ=255)
               gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAA‑CACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGa                                                                                                                                >  1:4094/1‑134 (MQ=255)
                cgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGAGAATTTTTTCAAAAAAATAAAAAAAATGGGGTAGca                                                                                                                          >  2:558738/1‑139 (MQ=255)
                            tcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTTAAACAAAAAAAAAAAAAGGAGGTGGCCTGATTAGGGaa                                                                                                              >  2:349629/1‑140 (MQ=255)
                               tATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGCGAATTTTTTCATACAAAAAAAAAAAAAAGAGGTGCCCTGATGAGGGaaaaa                                                                                                           >  1:775793/1‑138 (MQ=255)
                                aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGATGGCATGATGTGAGAAAGta                                                                                                          >  2:366862/1‑140 (MQ=255)
                                aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAAACAAAAAAAAAAAAAAGAGGTGGCATGATGAGGGGAAgaa                                                                                                          >  2:781490/1‑138 (MQ=255)
                                 ttttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAAAAAAAAAAATGGAGATGGCATGATGAGCGAAAGgaa                                                                                                         >  2:4069/1‑140 (MQ=255)
                                 ttttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaa                                                                                                         <  1:673663/140‑1 (MQ=255)
                                  tttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAACCTGGAGTTGCCATGATGAGCGAAAGtaat                                                                                                        <  2:616838/140‑1 (MQ=255)
                                                      tGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAATGGAGGTGCCATGATGAGGGGAAATAAAAAACAAAAGGCAGGGAaaaa                                                                                    >  2:368109/1‑140 (MQ=255)
                                                       gTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAg                                                                                   <  1:558738/140‑1 (MQ=255)
                                                       gTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAg                                                                                   <  1:115763/140‑1 (MQ=255)
                                                                 tatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAgaga                                                                         <  2:576923/140‑1 (MQ=255)
                                                                  atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGGTGCCATGATGAGCGAAAAGAATAAACAAAAGGCAGTGAATAAAATGACAGAGAt                                                                        >  2:273720/1‑140 (MQ=255)
                                                                  atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAAAACAAAAAATGGAGGTGGCCTGATGAGGGGAAATAAAAAACAAAAAGCAGTGAATAAGTTGACAaagaa                                                                        >  2:687929/1‑139 (MQ=255)
                                                                          tGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCAGAAAAAAACAAAAAAAGGGGGTGTCATGATGAGCGGAAGGAAAAAACAACAAGGAGTGAATAAGATGACAGAGAGTGTTGCTa                                                                >  2:42613/1‑140 (MQ=255)
                                                                          tGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAAATGAGGTGGCATGATGAGGGGAAAAAATAAGGAAAAAGCAGGGGAAAAATTGATAGAGAGTGTTGCTa                                                                >  2:187960/1‑140 (MQ=255)
                                                                                      aaaaaCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCa                                                    <  1:229460/139‑1 (MQ=255)
                                                                                          aCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGAt                                                <  2:24100/140‑1 (MQ=255)
                                                                                              caTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCt                                            <  1:368109/140‑1 (MQ=255)
                                                                                                    tAGTGGTAGTTTGTGGCGGTGATTTTTCCAAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGACAGTGAATAAGTTGACAGAGTTTCCCGCTAACTTTACCGCCATGATTTCTACCCGa                                      <  2:585946/140‑1 (MQ=255)
                                                                                                    tAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGGTGGCATGATGAGGGAAAGGAAAAAACAAAAAGCAGTGAAAAAGTTGACAGAGATTGTTGCTAAATTTTCCGCCATGATTTTTAcccca                                      >  2:148684/1‑138 (MQ=255)
                                                                                                     aGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGAGGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGaa                                     <  1:148684/140‑1 (MQ=255)
                                                                                                        ggTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAAAGGAGGTGGCATGATGAGCGGAAGTAAAAAAAAAAAAGCAGTGAATAAATTTAAAGAGAATGTTGCTAACTTTACCGCCATTATTTCTACCCGAAtac                                  >  2:409846/1‑138 (MQ=255)
                                                                                                             tttGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAAAGGAGGTGCCCTGATGAGCGAAAAAAAAAAACAACAAGCAGGGAAAAAAATGAAAGAGATTGTCGCTAACTTTTCCCCCATGATTTCTACCCCAATGCCtgat                             >  2:353117/1‑140 (MQ=255)
                                                                                                              ttGTGGCGGTGATTTTTTCCAGAAAAATACAAAAACTGGAGTTGCCATGAGGGGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTACCTTTACCGCCATGATTTCTACCCGAATGCCtgatg                            <  1:366862/140‑1 (MQ=255)
                                                                                                                   gCGGTGAATTTTTCCAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgtg                       <  1:687929/140‑1 (MQ=255)
                                                                                                                     ggTAAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgtggt                     <  1:349629/140‑1 (MQ=255)
                                                                                                                      gTGAATTTTTTCCGACAAAAACAAAAACAGGAGTTGCCATGATGGGAGAAAATAATAAGCAACAAGCAGTGAATAAGTTGACAGAGATTGTCGCTCACTTTACCGCCATGATTTTTTACCCAATTCCTGATGACgtgggg                    >  2:62018/1‑138 (MQ=255)
                                                                                                                                        aTACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTa  >  1:371596/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                          
CCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTA  >  NC_000913/3206307‑3206582

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