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

AACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAG  >  NZ_CP009273/3199717‑3199867
                                                                             |                                                                         
aaCTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACa                                                               <  2:805860/90‑1 (MQ=255)
  cTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAaaaaaaaa                                                             >  1:1043245/1‑90 (MQ=255)
    gTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAAATTTTTTAAACAaaaaaaaaaa                                                           >  2:870492/1‑90 (MQ=255)
           aCATTTTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGt                                                    <  1:538452/90‑1 (MQ=255)
           aCATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGt                                                    <  1:1174033/90‑1 (MQ=255)
           aCATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGt                                                    <  1:357052/90‑1 (MQ=255)
            cATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAAAAAACAAAAAATTGAggt                                                   >  1:327238/1‑88 (MQ=255)
              tatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGAAAAAAAAAAAAAAAGGAGGTgg                                                 >  1:881497/1‑89 (MQ=255)
              tatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCAACAAAAAAAAAAAAAGGAGGTgt                                                 >  1:1129693/1‑89 (MQ=255)
              tatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATTGAGGTGc                                                 >  1:382729/1‑90 (MQ=255)
              tatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGGATTTTTTTAAACAAAAAAAAAAAATAGAGGTgg                                                 >  2:387410/1‑89 (MQ=255)
              tatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGGAATTTTTTAGACAAAAAAAAAAAATGGAGGTgg                                                 >  2:1198865/1‑89 (MQ=255)
              tatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAAATTTTTTAGAAAAAAAAAAAAAATAGGGTTGc                                                 >  2:717087/1‑90 (MQ=255)
              tatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAAATTTTTTAAACAAAAAAAAAAAAAGGAGTTGc                                                 >  2:990753/1‑90 (MQ=255)
              tatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAAATTTTTCAAACAAAAAAAAAAAAAGGAGtttc                                                 >  2:1059206/1‑88 (MQ=255)
                           aaaaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGaa                                    <  1:559868/90‑1 (MQ=255)
                           aaaaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGaa                                    <  2:877683/90‑1 (MQ=255)
                           aaaaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGaa                                    <  1:1198865/90‑1 (MQ=255)
                              aaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAATTTTTTCAGACAAAAACAAAAAATGGAGGTGTCCTTGTGAGCGaaaag                                 >  1:141840/1‑88 (MQ=255)
                                  ataAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaata                             <  2:1132488/90‑1 (MQ=255)
                                      aaCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCa                         <  1:1041462/90‑1 (MQ=255)
                                          acaTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGaa                                    >  1:281030/1‑75 (MQ=255)
                                          acaTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGaa                                    <  2:281030/75‑1 (MQ=255)
                                            aTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGc                   <  2:306327/90‑1 (MQ=255)
                                            aTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGc                   <  2:382729/90‑1 (MQ=255)
                                            aTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGc                   <  1:1100982/90‑1 (MQ=255)
                                              ttttAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAg                 <  1:870492/90‑1 (MQ=255)
                                                  aGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAAAACAAAAAATGGAGTTGCCCTGATGAGCGAAAGAAAAAAGCAAAAAGCAGGGaa             >  1:1263085/1‑90 (MQ=255)
                                                     ggTAGTTTGTGGCGGTGAATTTTTTCAGACAAATAAAAAAACTGGAGGTGCCATGATGAGCGAAAGGAATAAACAACAAGCAGTGAATaa          >  1:403821/1‑90 (MQ=255)
                                                     ggTAGTTTGTGGCGGTGAATTTTTTCAGACAAATAAAAAAAATGGAGTTGCCATTATGAGCGAAAGTAATAAGCAAAAAGCAGTGAATaa          >  1:195775/1‑90 (MQ=255)
                                                     ggTAGTTTGTGGCGGGGAATTTTTTCACACAAAAAAAAAAAATGGGGGTGGCATGGTGAGGGGAAAGAAAAAACAAAAAGCAGGGAaaaa          >  2:648869/1‑90 (MQ=255)
                                                             gtgGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACag  <  1:871878/90‑1 (MQ=255)
                                                                             |                                                                         
AACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAG  >  NZ_CP009273/3199717‑3199867

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