Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273602,0680GA53.8% 4.4 / 13.6 13Q258* (CAG→TAG) ybdKYbdK family carboxylate‑amine ligase
Reads supporting (aligned to +/- strand):  ref base G (6/0);  new base A (0/7);  total (6/7)
Fisher's exact test for biased strand distribution p-value = 5.83e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAA  >  NZ_CP009273/601979‑602127
                                                                                         |                                                           
ggCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTg                                                             >  1:311417/1‑90 (MQ=255)
                         agcagTTATTTTTTTTCTTGTTTTTAAAGCGGGCTTCCCGTCAGAACCCGGGGGGGGGTAGCTTAAATTAACCCCGCCATATTTACTGCg                                    <  1:211829/90‑1 (MQ=255)
                         agcagTTATTTTTTTTCCTGATTTTTAAGGGGGGGTTCCGTCAGTACCCGGGGGGGGGAAGCCTAAATTAATCCCGCCATATTTACTGCg                                    <  2:139/90‑1 (MQ=255)
                         agcagGTATTTTTTTTCTTTTTTTTAAAGGGGGCGTTCCGCAAGAACCCGGGGGGGGGTAGCCTAAATAACCCCCGCCATATTTACTGCg                                    <  1:182903/90‑1 (MQ=255)
                         agcagGTATTTTTTTTCTTGTTTTTTAAGGGGGCGTCCCGTCAGTACCGGGGGGGGGGTAGCCTAAATTAA‑CCCGCCATATTTACTGCGt                                   <  1:308112/90‑1 (MQ=255)
                          gcagGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGt                                   >  1:322205/1‑90 (MQ=255)
                          gcagGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGt                                   >  2:65639/1‑90 (MQ=255)
                            agGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTgg                                 >  2:342364/1‑90 (MQ=255)
                                    tttttCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGt                         >  1:469266/1‑90 (MQ=255)
                                        tCCTGTTTTTTAAGGGGGCGTCCCGTAAGAAACGGGGGGGGGGTAGCTTAAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAAc                     <  2:311417/90‑1 (MQ=255)
                                             aTGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGt                >  2:291747/1‑90 (MQ=255)
                                                        ggttttCGGTAAGTAACGGGGGGGGGGTAGCTTAAATTCCCCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCAccc     <  1:435341/86‑1 (MQ=255)
                                                        ggCGTCCCGTCAGAACCGGGGGGGGGGTCGCATAAATAAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCAccc     <  2:117490/90‑1 (MQ=255)
                                                           gTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGaa  >  2:190012/1‑90 (MQ=255)
                                                                                         |                                                           
GGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAA  >  NZ_CP009273/601979‑602127

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: ‑