Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,214,1180CA62.5% 1.9 / 18.0 16A422A (GCC→GCAygjGputrescine aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (6/0);  new base A (0/10);  total (6/10)
Fisher's exact test for biased strand distribution p-value = 1.25e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAG  >  NZ_CP009273/3214034‑3214193
                                                                                     |                                                                           
gTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCaaacc                                                                        >  2:556282/1‑88 (MQ=255)
                            tttCCCAGCGAGATGTCCCGCCCGCCCGTACGGGGGCCCGAACCGCTCAAAACCGCAAAACCGACCCGCATTAACCCGCcactgacactg                                            <  2:9762/90‑1 (MQ=255)
                             ttGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCcactgacactga                                           >  1:592931/1‑90 (MQ=255)
                                 cAGGGAGTTGTCCCGCCGGCGGGTACGGGGGGCCGGACCGCTAAAAAACGCAAAAACGACCCGTATTGAACCGCCACTGACACTGACCAt                                       <  1:646062/90‑1 (MQ=255)
                                     gagaTTTCCCGCCAGCGCGTACGGGGGGCCGGAACGCTAAAAAACGCAAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGaa                                   <  1:276323/90‑1 (MQ=255)
                                      agttGTCCCGCCAGCGCGTCCTGGTGGCCGGAACCCTCAATAACCCAAAAACGACCCGCATTAAACCGCCACTGACACTGACCATTGAAc                                  <  1:17368/87‑1 (MQ=255)
                                      agaTTTCCCCCCGGCGCGTCCGGGTGGCCGGACCCCAAAAAACCCCAAAAACGACCCGCATTGAACCGCCACTGACACTGACCATTGAAc                                  <  1:508904/90‑1 (MQ=255)
                                      agaTGTCCCGCCGGCGCGTCCTGGGGGCCGGAACCCTAAAAAACGCAAAAACGACCCGCATTGAACCGCCACTGACACTGACCATTGAAc                                  <  1:19623/90‑1 (MQ=255)
                                        aTGTCCCGCCGGGGCGTCGGGGGGCCCGAACCGCTAAAAAACCCAAAACCGCCCCGCATTAACCCGCCACTGACACTGACCATTGAACAg                                <  2:592931/90‑1 (MQ=255)
                                        aTGTCCCGCCAGCGCGTACTGGGGGCCGGAACGCTAAAAAACGCAAAAACGACCCGCATTGAACCGCCACTGACACTGACCATTGAACAg                                <  1:160547/90‑1 (MQ=255)
                                               gCCAGCGCGTGCTGGTGGCGGGAACGCTAAATAACGCAAAAACGACCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAAc                         <  1:183883/90‑1 (MQ=255)
                                                  ggcgcGTCCGGGGGGCGGGACCGCCAAAAACCGCAAAAACGACCCGCATTAACCCGCCACTGACACTGACCATTGAACAGTGTGAACTgg                      <  1:377262/89‑1 (MQ=255)
                                                               ggCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAgcggc         >  2:107468/1‑90 (MQ=255)
                                                                    gAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTa    >  2:235061/1‑90 (MQ=255)
                                                                    gAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGAGATCAAAGCGGGGCGTa    >  2:58818/1‑90 (MQ=255)
                                                                      aCGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAg  >  1:429301/1‑90 (MQ=255)
                                                                                     |                                                                           
GTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAG  >  NZ_CP009273/3214034‑3214193

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