Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,206,4440TA42.9% 28.8 / 15.7 28intergenic (‑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 (12/0);  total (14/14)
Fisher's exact test for biased strand distribution p-value = 5.98e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGT‑AATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAA  >  NC_000913/3206312‑3206576
                                                                                                                                    |                                                                                                                                     
cccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACaaaaa                                                                                                                                >  1:148406/1‑140 (MQ=255)
 ccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAAc                                                                                                                               <  1:789358/140‑1 (MQ=255)
    gCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATgg                                                                                                                            >  1:856119/1‑140 (MQ=255)
     ctctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAaaaaaaaaaaaaaagga                                                                                                                           >  1:465220/1‑135 (MQ=255)
      tctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAg                                                                                                                          <  1:1209526/140‑1 (MQ=255)
      tctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAg                                                                                                                          <  2:1133440/140‑1 (MQ=255)
          gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATAGAGTTgg                                                                                                                      >  2:20959/1‑139 (MQ=255)
                     attcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGTTGCCCTGATTAGCg                                                                                                           >  1:747556/1‑140 (MQ=255)
                     attcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCAGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCg                                                                                                           >  1:391041/1‑140 (MQ=255)
                      ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                                          <  1:20959/140‑1 (MQ=255)
                      ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                                          <  1:1090084/140‑1 (MQ=255)
                          tATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATAAAAAAAATTGAGTTGGCATTATTAggggaaaa                                                                                                      >  1:289938/1‑133 (MQ=255)
                            ttttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCAAAAAAAAAAAAAAAATGGAGTTGGCATGATGAGGGGAAAt‑aa                                                                                                   >  2:1170051/1‑140 (MQ=255)
                                   gTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGT‑AATAAGCaa                                                                                            <  1:959932/140‑1 (MQ=255)
                                   gTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGT‑AATAAGCaa                                                                                            <  1:594703/140‑1 (MQ=255)
                                    ttGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGT‑AATAAGCAAc                                                                                           <  2:148406/140‑1 (MQ=255)
                                              aaCTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGT‑AATAAGCAACAGGCAGTGaa                                                                                 <  1:30819/140‑1 (MQ=255)
                                               aCTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAAAAAAAAACAAAAAAAGGAGGTGCCAATATGGGGGAAAATAAAAAAGAAAAAGGCAgcaaa                                                                                 >  2:992312/1‑136 (MQ=255)
                                                  gTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAATGGAGGTGCCATGATGAGCGAAAGT‑AAAAAAAAAAAAGCAGTGAATaaa                                                                             >  2:30819/1‑139 (MQ=255)
                                                   ttAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGCATTTTTTCAGACAAAAAAAAAAAATGGAGGTGCCCTGATGGGCGGAAAT‑AATAAAAAACAAGCCGTGAATaaaa                                                                            >  2:813762/1‑138 (MQ=255)
                                                            tatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAAAACAAAAAAAGGAGGTTCCATGATGAGCGAAAGT‑AATAAAAAACAAGCAGGGAATAAGTTGACAgaga                                                                   >  2:735448/1‑140 (MQ=255)
                                                             atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCATACAAAAACAAAAAATGGAGGTGTCCTGGTGAGCGAAAGA‑AAAAAAAAAAAAGCAGAGAATAAAATTACAAAGAt                                                                  >  2:202462/1‑140 (MQ=255)
                                                             atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGGAATTTTTCAGACAAAAACAAAAAATTGAGGTTGCATGATGAGGGGAAAA‑AAAAAAAAAAAAGCAGAGAAAAAATTGACaaaaaa                                                                  >  2:493027/1‑135 (MQ=255)
                                                                    gTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGT‑AATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCt                                                           <  1:167879/140‑1 (MQ=255)
                                                                     tGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAATGGAGGTGTCCTGATGAGGGGAAGA‑AAAAAAAAAAAAGCAGAGAAAAAATTGACAGAGATTGTCGCTa                                                          >  2:343420/1‑140 (MQ=255)
                                                                     tGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAAAAAAAAAAATGGGGTTGGCCTGATGAGGGAAAAT‑AAAAAAAAAAAAGCAGGGAATAAAATGACAGAGATTGTTGCTa                                                          >  2:6160/1‑140 (MQ=255)
                                                                       gCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAAATGAGTTGGCATGATGAGGGGAAAG‑AAAAAACAAAAAGCAGTGAAAAAAGTTAAAAAGAATTTCGCTaaa                                                        >  2:574948/1‑139 (MQ=255)
                                                                                        acaTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCAAAAGT‑AATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTc                                       <  1:6160/140‑1 (MQ=255)
                                                                                                    gTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGT‑AATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGcc                           <  1:574948/140‑1 (MQ=255)
                                                                                                             ggCGGTGAATTTTTCCAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGT‑AATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgt                  <  2:652410/140‑1 (MQ=255)
                                                                                                                 gTGAATTTTTTCCAACAAAAACAAAAAAAGGAGTTGCCATGATGAGCGAAAGG‑AATAAAAAAAAAGCAGGGAAAAAAATTAAAGAGATTGTCGCTAACTTTACCCCCAAGATTTTTACCCGAAAGCCTGATGACgtgggg              >  2:1119500/1‑138 (MQ=255)
                                                                                                                             aGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGT‑AATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTaaaa  <  1:22798/140‑1 (MQ=255)
                                                                                                                                   aTACAAAAACTGGAGTTGCCATGATGAGCGAAAGT‑AATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTaaa   >  2:153231/1‑133 (MQ=255)
                                                                                                                                   aTACAAAAACTGGAGTTGCCATGATGAGCGAAAGT‑AATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTaaa   <  1:153231/133‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                     
CCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGT‑AATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAA  >  NC_000913/3206312‑3206576

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