Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,033,8110TG56.2% 11.9 / 17.4 16Y206D (TAT→GAT) appCcytochrome bd‑II oxidase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (3/4);  new base G (9/0);  total (12/4)
Fisher's exact test for biased strand distribution p-value = 1.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-04
Rejected as consensus: Frequency below/above cutoff threshold.

GGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTAC  >  NZ_CP009273/1033729‑1033896
                                                                                  |                                                                                     
ggTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTTACCCGGGGCATGTTTTTTTTTGCGATCAGCGCCCGGGATTTTCt                                                                                >  1:425775/1‑90 (MQ=255)
 gTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTg                                                                               <  1:459666/90‑1 (MQ=255)
              aaTTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTTTTATTGCGATCAGCGCCTGGGATTTACTTCGCGGGCGGGGgc                                                                  >  1:209690/1‑90 (MQ=255)
                   gtgCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGGATTTACTGCGCGGACGGGAGCGCAAt                                                             >  2:227035/1‑90 (MQ=255)
                          ccGTAATGGCGGGCTACGTGACCGGGGGCATGTTTTTTATGGGGATCAAGGGCGGGGTTTTACTGCGGGGCCGGGGGCGCGATGTCggca                                                      >  1:106060/1‑87 (MQ=255)
                                    gggCTACGTGACCGGGGCCATGTTTTTTATGGCGATCAGCGCCTGGGATTTTCTTCGCGGGCGGGGGCGCAATGTCGCCTTACGCTCggt                                            >  1:568748/1‑88 (MQ=255)
                                    gggCTACGTGACCGGGGCCATGTTTTTTATGGCGATCAGCGCCTGGGATTTTCTGCGCGGGCGGGGGCGCAATGTCGCATTACGCTCGtt                                            >  1:220385/1‑90 (MQ=255)
                                    gggCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGtt                                            >  1:451186/1‑90 (MQ=255)
                                    gggCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGGATTTTCTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGtt                                            >  2:523101/1‑90 (MQ=255)
                                    gggCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGGATTTACTGCGCGGGCGGGGGCGCAATGTCGCATTACGCTCGtt                                            >  1:6177/1‑90 (MQ=255)
                                    gggCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGGATTTACTGCGCGGACGGGGGCGCAATGTCGCATTACGCTCGtt                                            >  2:90631/1‑90 (MQ=255)
                                                  gggCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCg                              <  2:158274/90‑1 (MQ=255)
                                                                   gCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCg             >  2:429712/1‑90 (MQ=255)
                                                                     gATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGAt           >  2:189963/1‑90 (MQ=255)
                                                                      aTCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGAtt          <  2:6177/90‑1 (MQ=255)
                                                                              cTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTAc  <  2:568748/90‑1 (MQ=255)
                                                                                  |                                                                                     
GGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTAC  >  NZ_CP009273/1033729‑1033896

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