Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,199,7940CT48.4% 11.8 / 20.6 31intergenic (‑179/‑28)ttdR/ttdADNA‑binding transcriptional activator TtdR/L(+)‑tartrate dehydratase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (2/14);  new base T (15/0);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 5.09e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGA  >  NZ_CP009273/3199659‑3199925
                                                                                                                                       |                                                                                                                                   
aaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGAGAATTTTTTCa                                                                                                                                   >  1:552022/1‑138 (MQ=255)
               tGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAATACaaaaaat                                                                                                                   >  2:456670/1‑137 (MQ=255)
                      gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAg                                                                                                               <  2:158307/136‑1 (MQ=255)
                      gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGAAAAATAAAAAAACAGGGGTTg                                                                                                            >  1:240295/1‑139 (MQ=255)
                              ttcattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatgat                                                                                                     <  2:277782/138‑1 (MQ=255)
                                      tATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAAATGGAGTTGTCCTGGTGAGCGGAAg                                                                                            >  1:609856/1‑139 (MQ=255)
                                             gcgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGTATTTTTCCCTACCAATACAAAAACTTAGGTTTCCATGATTGGCGAAAGTAATAAg                                                                                     >  1:418843/1‑139 (MQ=255)
                                                  gATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGTTGCCATGATGAGGGAAAGGAATAAGaaaa                                                                                 >  2:555225/1‑137 (MQ=255)
                                                    ttCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAgg                                                                              <  2:252112/139‑1 (MQ=255)
                                                    ttCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAgg                                                                              <  1:456670/139‑1 (MQ=255)
                                                      cACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGc                                                                             <  1:655529/138‑1 (MQ=255)
                                                                         atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATAAAAAAACTGGGGTTGGCATGATGGGCGGAAGAAATAAGCAACAAGCAGTGAATAAATTGACa                                                             >  1:629456/1‑135 (MQ=255)
                                                                         atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGGTGCCATGATGAGCGAAAGTAATAAACAAAAGGCAGTGAAAAAATTGACa                                                             >  2:28547/1‑135 (MQ=255)
                                                                         atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAAAATAAAAAAAAAGGGGGTTGCATGGTGAGGGAAAGGAAAAAAAAAAAAGGAATGAATAAGATGTCaga                                                           >  2:302238/1‑137 (MQ=255)
                                                                         atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGCGAATTTTTTCAAACAAAAAAAAAAACAGGGGTTGCCCTGATGAGCGGAAGAAATAAGAAAAAAGCAGAGAAAAAGTTGaac                                                             >  1:531159/1‑133 (MQ=255)
                                                                                 tGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATAAAAAAACTGGGGGTGTCATGATGAGCGAAAATAATAAGCAACAAGCAGAGAATAAATTGACAGAGATTGTCGc                                                  >  2:63579/1‑138 (MQ=255)
                                                                                 tGGCAAAAAATATAAAACCACATTTTTAGTGGGAGTTTGTGGCGGTGAAATTTTTCAGACAAAAACAAAAACTGGAGTTGCGATGATGGGCGGAAGGAATAAGAACAAAGCAGAGAATAAATTGAAAGAGAgtgtggg                                                  >  1:581624/1‑135 (MQ=255)
                                                                                             taAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAAATTTTTCAGACCAAAAAAAAAAATAGGGGTGCCATGATGAGCGGAAGGTAAAAGGAAAAAGGAGGGAAAAAGGTGACAGAGATTGTCGCTTAATTTACCGcc                                     >  1:395324/1‑139 (MQ=255)
                                                                                               aaaCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATAAAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCa                                    <  1:45795/138‑1 (MQ=255)
                                                                                                    acaTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGTTGCCATGATGAGCGGAAGTAATAAACAACAAGCAGGGAATAAATTGACAGAGATTGTCGCTTACTTTACCGCCCTGAttt                              >  2:651097/1‑139 (MQ=255)
                                                                                                    acaTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTCCCAGAAAAATACAAAAACTGGAGTTGCCAGGATGAGCGAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGAttt                              <  2:293438/139‑1 (MQ=255)
                                                                                                     caTTTTTAGTGGTAGTTTGTGGCGGTGGAATTTTTCAAACAAAAAAAAAAAATGGAGTTGTCCTGGTGGGCGAAAAGAATAAAAAACAAGCAGTGAAAAAAATGCCAGAGGTTGTTGATAAACTTAACCCCCTGAtttt                             >  2:592272/1‑138 (MQ=255)
                                                                                                      aTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCt                            <  1:555225/139‑1 (MQ=255)
                                                                                                           tAGGGGTATTTTTGGGCGTTAAATTTTTCCAAAAAAAAAAAAAAACTGGAGTTGCAATGATGAGCAAAAGAAATAAGCAACAGGCAGTAAATAAGTTAACAGAGATTGTCGTTACCTTTACCGCCATGATTTCTAccc                        <  2:326422/138‑1 (MQ=255)
                                                                                                              tgttggTTTGGGGCGGTAAATTTTCCCAAAAAAATACAAAAACTGGAGTAGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACATAGATGGTCGCTAGCTTTACCGCCATGATTTCTACCCGAAt                    <  1:592272/134‑1 (MQ=255)
                                                                                                                gTAGTTTGTGGCGGTGAATTTTTTCCGAAAAAAACAAAAAATGGAGGTGGCCTGATGAGGGGAAGTAAAAAAAAACAGGCAGAGAAAAAAATTACAGAGGTTGTCGCTAAAATTAACGCCATGATTTTTACCCCAAAGc                  >  2:160513/1‑139 (MQ=255)
                                                                                                                  aGTTTGTGGCGGTGAATTTTTTCAGACAAATAAAAAAAATGGAGTTGCCCTGATGAGGGAAAGTAATAAACAACAAGCAATGAATAAATTGACAAAGATTGTCGCTAACTTTACCGCCCTGATTTCTACCCCaaatcct                >  1:465132/1‑134 (MQ=255)
                                                                                                                    tttGGGGCGGGGAATTTTTCCAGACAAATACAAAAACTGGAGTTCCCATAATGAGCGAAAGAAATAACCAACAGGCAGTGAAAAAGTTGACAGAGATGGCCGCTACCTTAGCCGCCATGATTTCTACCCGAATGCCtga              <  2:629456/139‑1 (MQ=255)
                                                                                                                           cGGTGAATTTTACAATAAAAATAAAAAAACTGGAGTTGCCATGATGAGCGAAAGAAAAAAGCACCGGGCAGTAAATAATTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgtg       <  2:485151/139‑1 (MQ=255)
                                                                                                                                aaTTTTTCCAGACAAATACAAAAACTGAAGTTGCCATTATGAGCGAAATAAATAAGAAACAGGCAGTGAATAAGTTAACAGAGATTGTCGTTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGa  <  2:124486/139‑1 (MQ=255)
                                                                                                                                aaTTTTTCCAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTGAATAAGTGGACAGAGATTGTCGCTATTTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGa  <  2:609856/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                   
AACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGA  >  NZ_CP009273/3199659‑3199925

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