Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,698,2980GT57.9% 13.1 / 28.0 19H90N (CAT→AAT) ydgJoxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (8/0);  new base T (0/11);  total (8/11)
Fisher's exact test for biased strand distribution p-value = 1.32e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-05
Rejected as consensus: Frequency below/above cutoff threshold.

CAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGG  >  NZ_CP009273/1698212‑1698379
                                                                                       |                                                                                 
cagCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGttt                                                                                >  1:1029134/1‑90 (MQ=255)
         ccccAGGCTTTTTGCCATCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTc                                                                       >  1:987004/1‑90 (MQ=255)
               gttttttGCCGGCCCACCCAGCTCCCGCGTTTGTGACGGTGCCAGGGTAAGGGTTTTATCAACGCCCCCTTTTTTCCCCGCCTCAAgcgc                                                                 <  1:408695/88‑1 (MQ=255)
               gttttttCCCGGCGCTCCCACCTCCCGCGTTTGTGACAGTGTCAGGGTAAGGGTTTTACAACCGCCCACTTTTTTACCCGCCTCAAgcgc                                                                 <  1:449407/88‑1 (MQ=255)
                   tttGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCt                                                             >  2:934752/1‑90 (MQ=255)
                         aGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCt                                                       >  2:816094/1‑90 (MQ=255)
                              aTCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACgg                                                  >  2:633607/1‑90 (MQ=255)
                                 cgcctctCGCGTTTGTGACAGTGTCAGGGTAAGGGGTTTATCACCGCCCACTTTTTTCCCCGCCTCAAGCGCCGCTTTGGCTAACGGGaa                                               <  1:403237/87‑1 (MQ=255)
                                             ttGTGACAGTGCCCGGGTAAGGGTTTTAACAACCCCCACTTTTTTCCCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGtt                                   <  2:987004/90‑1 (MQ=255)
                                              tgtgACGGTGTCAGGGAAAGGGTTTTACAAACGACCACTTTTTTCCCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTg                                  <  1:197230/90‑1 (MQ=255)
                                              tgggACGGTGTCAGGGAAAGGGTTTTACCAACGCCCAATTTTTTCCCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTg                                  <  1:591447/87‑1 (MQ=255)
                                              tgggACAGTGCCGGGGAAAGGGTTTTAACAACGACCATTTTTTTCCCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTg                                  <  2:993081/87‑1 (MQ=255)
                                                cgACGGTGCCACGATAAGGGTTTTACCAACGTCCATATTTTTCCCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTggg                                <  1:633607/89‑1 (MQ=255)
                                                     gtgtCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAg                           >  1:590897/1‑90 (MQ=255)
                                                            ggaaaaGGTTTTACAAACGACCATTTTTTTCCCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGAc                    <  2:1029134/87‑1 (MQ=255)
                                                                aaGGGTTTTACCAACAACCATTTTTTTCCCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAAt                 <  2:590897/90‑1 (MQ=255)
                                                                      tttATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTc           >  1:340901/1‑90 (MQ=255)
                                                                        tATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTa         >  1:914499/1‑90 (MQ=255)
                                                                               ccccAAATTTTTCCCCGCCTCAAGCGCCGCTTTGGATAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTggg  <  1:402787/88‑1 (MQ=255)
                                                                                       |                                                                                 
CAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGG  >  NZ_CP009273/1698212‑1698379

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