Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,199,8030TA54.5% 12.3 / 31.1 22intergenic (‑188/‑19)ttdR/ttdADNA‑binding transcriptional activator TtdR/L(+)‑tartrate dehydratase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (1/9);  new base A (12/0);  total (13/9)
Fisher's exact test for biased strand distribution p-value = 2.01e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-05
Rejected as consensus: Frequency below/above cutoff threshold.

TATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAA  >  NZ_CP009273/3199731‑3199880
                                                                         |                                                                             
tatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGGATTTTTTCAGACAAAAACAAAAAAAGGAGGTgg                                                              >  2:3279/1‑89 (MQ=255)
tatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAAATTTTTCAAAAAAAAAAAAAAAAAGGAGGTGc                                                              >  1:812241/1‑90 (MQ=255)
  ttCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAAAAAACAAAAAAAGGAGGTGCCa                                                            >  2:129387/1‑90 (MQ=255)
     aCGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatga                                                         <  1:174615/90‑1 (MQ=255)
        gTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAAAAAAAAAAAAAAGGAGGTGGCatgttta                                                      >  2:911767/1‑86 (MQ=255)
          ggCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAAAGGAGGTGGCATTATGAgg                                                    >  1:178165/1‑89 (MQ=255)
          ggCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGc                                                    <  1:3264/90‑1 (MQ=255)
                    ataAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaata                                          <  1:479025/90‑1 (MQ=255)
                    ataAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaata                                          <  2:812241/90‑1 (MQ=255)
                    ataAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAAATTTTTTAAAAAAAAAAAAAAAAAGGGGTTTGCATGATGGGGgaaagaaaaa                                          >  1:883234/1‑88 (MQ=255)
                           cacaTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACa                                   <  2:421142/90‑1 (MQ=255)
                            acaTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGTTGCCATGATGAGCGGAAGGAATAAGCAAAAg                                  >  2:254482/1‑90 (MQ=255)
                            acaTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAAAGGAGTTGCCATGATGAGGGGAAGGAATAAGCAACaa                                  >  2:403124/1‑89 (MQ=255)
                               tttttAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCa                               <  1:556908/90‑1 (MQ=255)
                               tttttAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCa                               <  1:826901/90‑1 (MQ=255)
                                     gTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAAt                         <  2:183538/90‑1 (MQ=255)
                                       ggTAGTTTGTGGCGGTGGAATTTTTCCGACAAAAACAAAAAATGGGGGTGGCCTGGTGGGGGGAAGTAATAAAGAAAAAGCAGGGAaaaa                       >  1:79632/1‑90 (MQ=255)
                                       ggTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGGTGCCATGATGAGGGAAAATAAAAAAAAAAAAGCAGTGAaaaa                       >  1:292364/1‑90 (MQ=255)
                                       ggTAGTTTGTGGCGGGGGAATTTTTTAGACAAAAACAAAAAATGGGGGTGGCTTGGTGGGGGGAAGGAATAAACAAAAGACAGGGAaaaa                       >  1:894341/1‑90 (MQ=255)
                                         tAGTTTGTGGCGGTGAAATTTTTCAAACAAAAAAAAAAAATGGGGGTGCCATGATGAGCGAAAAAAAAAAACAAAAAGCAGGGAaaaaat                     >  2:665354/1‑88 (MQ=255)
                                                   cGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGAt           <  1:403124/90‑1 (MQ=255)
                                                             tttCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTaa  >  2:720289/1‑90 (MQ=255)
                                                                         |                                                                             
TATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAA  >  NZ_CP009273/3199731‑3199880

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: