Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,994,5290AC58.8% 4.8 / 14.8 17T289P (ACC→CCC) xdhAxanthine dehydrogenase molybdenum‑binding subunit XdhA
Reads supporting (aligned to +/- strand):  ref base A (6/1);  new base C (0/10);  total (6/11)
Fisher's exact test for biased strand distribution p-value = 5.66e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.36e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTC  >  NZ_CP009273/2994440‑2994607
                                                                                            |                                                                              
agCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTa                                                                                >  1:101630/1‑90 (MQ=255)
agCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTa                                                                                >  1:256646/1‑90 (MQ=255)
agCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTa                                                                                >  2:108733/1‑90 (MQ=255)
agCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTa                                                                                >  1:536210/1‑90 (MQ=255)
     aTGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGc                                                                           >  2:352233/1‑90 (MQ=255)
      gggcgccATCCCGGACCAGCAAGCTTGGCGGCCTCCGGGTAAATTTTCCCTTTAGCCGGAAAGGGTTTTCCCCCCCACCCCGTCCCCGcc                                                                          <  2:101630/89‑1 (MQ=255)
           gcATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCt                                                                     >  1:2299/1‑90 (MQ=255)
                        gCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACg                                                        >  1:192923/1‑90 (MQ=255)
                               tGGCGGCTTTCGGTTAAATTTTCCCCTTAGCGGTGAAGGGTTTTTCCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAAt                                                 <  2:256646/90‑1 (MQ=255)
                                             gtaaaaTTTCCCCTAACCGGTGAAGATTTTTCCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcgc                                    <  1:112328/88‑1 (MQ=255)
                                                    ttCCCTTTAGCTGTAAAGTGTTTTCCCCCCCACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGaa                             <  1:39272/90‑1 (MQ=255)
                                                           ttAGCCGTGAAGAGTTTTCCCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGa                       <  1:195293/90‑1 (MQ=255)
                                                             ggCCGGGAAGAGTTTTCCCCCCCCCCCCGCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGaaa                     <  2:536210/89‑1 (MQ=255)
                                                             cgCCGGGATGGGTATTCCCCCGCACCCCGCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGaaa                     <  2:524669/89‑1 (MQ=255)
                                                             aGCCGTGAAGAGTTTTCCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGaaa                     <  1:232662/90‑1 (MQ=255)
                                                             aGCCGTAAAGGGTTTTCCCCCCCACCCCGCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGaaa                     <  1:108733/90‑1 (MQ=255)
                                                                  gTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGtt                 >  2:157610/1‑90 (MQ=255)
                                                                   tGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTta                >  1:288941/1‑90 (MQ=255)
                                                                         gtgtTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCt          >  1:456283/1‑90 (MQ=255)
                                                                                 ccGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTc  <  1:34501/89‑1 (MQ=255)
                                                                                 ccGCACCCCGCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTc  <  2:456283/89‑1 (MQ=255)
                                                                                            |                                                                              
AGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTC  >  NZ_CP009273/2994440‑2994607

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