Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,620,6020GC53.8% 2.5 / 13.4 13R346G (CGG→GGG) dcppeptidyl‑dipeptidase Dcp
Reads supporting (aligned to +/- strand):  ref base G (5/1);  new base C (7/0);  total (12/1)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.70e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CACCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGA  >  NZ_CP009273/1620530‑1620670
                                                                        |                                                                    
cacCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGGTGGGCCTCTTCAAGAACAAATTTTTCCCCCCGTACCTCTTCCGCAt                                                     >  1:779040/1‑90 (MQ=255)
cacCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTTTCCCCCCGTACCTGTTCCGCAt                                                     >  1:19845/1‑90 (MQ=255)
cacCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTTTCCCCCCGTACCTGTTCCGCAt                                                     >  1:489612/1‑90 (MQ=255)
                cGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCCTCAAGATCAAATTTTTCCCCCCGTACCTTTTCCGCATAAAATGCCCCGTcccc                                     >  2:629920/1‑89 (MQ=255)
                 gtttttAATTAAAAATATGGTTTGAGCTGCCCCTCAAAAAAATAAATTTCCCCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCAc                                    <  2:433921/87‑1 (MQ=255)
                 gtgtTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCAc                                    >  2:1231935/1‑90 (MQ=255)
                 gtgtTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCCCCGTACCTGTTCCGCATAAAATGCCCAGTCCCAc                                    >  2:1025202/1‑90 (MQ=255)
                  tgtTTAATTCAAAATATGGCTTGGGCTGCGCCTCATCAAGATCAAATTTTTCCCCCCGTACCTGTTCCGCATAAAATTCCCAGTCCCACg                                   >  1:308382/1‑90 (MQ=255)
                           cAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTTTCCCCCCGTACCTGTTCCGCCTAAAAAGCCCCGTCCCCCCGCTgcccgc                          >  1:718012/1‑90 (MQ=255)
                                  tGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAAccc                   >  2:1197425/1‑90 (MQ=255)
                                        gAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCctg             >  1:49767/1‑90 (MQ=255)
                                             gcgcCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTgctgct        >  1:925447/1‑90 (MQ=255)
                                                   catcaAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGa  >  2:796253/1‑90 (MQ=255)
                                                                        |                                                                    
CACCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGA  >  NZ_CP009273/1620530‑1620670

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