Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,199,7940CT55.0% 5.3 / 24.2 20intergenic (‑179/‑28)ttdR/ttdADNA‑binding transcriptional activator TtdR/L(+)‑tartrate dehydratase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (1/8);  new base T (11/0);  total (12/8)
Fisher's exact test for biased strand distribution p-value = 7.14e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTG  >  NZ_CP009273/3199674‑3199920
                                                                                                                        |                                                                                                                               
tGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCCGACAAATAaaaaaaat                                                                                                              >  1:497544/1‑137 (MQ=255)
tGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTTAGAcaaaaacaaaaaat                                                                                                              >  2:539307/1‑137 (MQ=255)
  ctctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAAATgg                                                                                                            >  2:166636/1‑139 (MQ=255)
   tctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAAAAAAAAAAAAATg                                                                                                             >  1:436007/1‑137 (MQ=255)
       gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAg                                                                                                          <  1:302751/136‑1 (MQ=255)
          cAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGCCAAATACAAAAAATGGAGTTGGCa                                                                                                    >  1:228964/1‑139 (MQ=255)
                     cTTATTTTCGCGTTTATTCACGAACTGTTAATTACATATTCACGGTTGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                          <  1:186164/138‑1 (MQ=255)
                     cTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                          <  1:396792/138‑1 (MQ=255)
                                     ttCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAgg                                                                         <  2:113187/139‑1 (MQ=255)
                                        aCGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCa                                                                       <  1:539307/138‑1 (MQ=255)
                                                        atatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGCATTTTTTCCGACAAATATCAAAAACAGAGTTGCCATCCTGAGCGAAAGTAATAAAAAACCAGCAGCGAATAAGttttga                                                        >  1:65953/1‑133 (MQ=255)
                                                          atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAATAAAAAAAATGGAGTTGGCATGGTGAGCGAAAGTAAAAAGCAAAAGGCAGTGAATAAGGTGaac                                                        >  2:531177/1‑133 (MQ=255)
                                                               cGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGTTGCCATGGTGAGGGAAAGAAATAAAAAAAAAGCAGTGAAAAAATTGACAAAGATTGtg                                               >  2:192241/1‑138 (MQ=255)
                                                                    gCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCacac                                                                                                                            >  2:360974/1‑58 (MQ=255)
                                                                    gCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGAc                                                                                                                            <  1:360974/58‑1 (MQ=255)
                                                                     cAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTTAGACAAAAAAAAAAACTGGGGGTGGCCTTATGGGCGAAAGTAAAAAAAAACAAGGAATGAAAAAATTGACAGAGATTGTCGCTAAc                                         >  2:366787/1‑139 (MQ=255)
                                                                                      caTTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCAGACAAAAACAAAAACTGGAGGTGCCCTGATGAGCGGAAGAAATAAACAACAGGCAGTGAATAAGTTGAAAGAGATTGTCGCTAAATTTACCCCCATGAtttt                        >  1:422435/1‑138 (MQ=255)
                                                                                      caTTTTTAGTGGTAGTTTGTGGCGGGGAAATTTTTCCGACAAATACAAAAACAGGAGGTGCCACGATGAGGGGCAGCAATAAACAAAAAGCAGCGAATAAGTTGACAGAGATTGTCGCTAACTTTTCCGCCCTGAtttt                        >  1:124218/1‑138 (MQ=255)
                                                                                                     tttGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtga         <  1:614051/139‑1 (MQ=255)
                                                                                                            cGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTACCTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgtg  <  1:155811/139‑1 (MQ=255)
                                                                                                                        |                                                                                                                               
TGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTG  >  NZ_CP009273/3199674‑3199920

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