Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,214,1260TC52.6% 17.2 / 27.7 19I425T (ATC→ACC) ygjGputrescine aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (9/0);  new base C (0/10);  total (9/10)
Fisher's exact test for biased strand distribution p-value = 1.08e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-05
Rejected as consensus: Frequency below/above cutoff threshold.

AAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTA  >  NZ_CP009273/3214048‑3214214
                                                                              |                                                                                        
aaaTCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAc                                                                               >  2:53064/1‑90 (MQ=255)
              tttGCAAGGGAGTTGTCCCGCCGGCGCGTACTGGGGGCCGGACCGCTAAAAAACCCAAAAACGACCCGCATTGAACCGCcactgacactg                                                                 <  2:645161/90‑1 (MQ=255)
              tttCCCAGCGAGTTGTCCCGCCAGCGCGTCCGGGGGGCCGGACCGCTAAAAACCCCAAAAACGACCCGCATTGAACCGCcactgacactg                                                                 <  2:719759/90‑1 (MQ=255)
               ttGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCcactgacactga                                                                >  2:764143/1‑90 (MQ=255)
               ttGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCcactgacactga                                                                >  1:573086/1‑90 (MQ=255)
                tGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCactgacactgac                                                               >  1:701056/1‑90 (MQ=255)
                       gagaTGTCCCGCGGGGGGGTCCGGGGGGCGGAAACGCCAAAAAACGAAAAAACGCCCCGCATTAACCCGCCACTGACACTGACCATTGaa                                                        <  1:612284/90‑1 (MQ=255)
                        agaTGTCCCCCCGGCGCGTACTGGGGGCCGGAACGCTAAAAACCGCAAAAACGACCCGCATTGAACCGCCACTGACACTGACCATTGAAc                                                       <  2:255735/90‑1 (MQ=255)
                               ccGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGa                                                >  2:265120/1‑90 (MQ=255)
                                 gCCAGCGCGTACTGGGGGCCGGAACGCTAAAAAACGCAAAAACGACCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAAc                                              <  2:243022/90‑1 (MQ=255)
                                    ggcgcGTCCGGGGGGCCGAAACCCTAAAAAACCCAAAACCGACCCGAATTAAACCGCCACTGACACTGACCATTGAACAGTGTGAACTgg                                           <  1:50823/89‑1 (MQ=255)
                                    aGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTgg                                           >  1:271124/1‑90 (MQ=255)
                                               ggggCGGGACCGCTAAATAACGCAAAAACGACCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAgcg                                <  2:49490/88‑1 (MQ=255)
                                                 ggCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAgcggc                              >  2:534635/1‑90 (MQ=255)
                                                                  aCGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTggcgg             >  1:562399/1‑90 (MQ=255)
                                                                     caaaaaCGACCCGAATTGACCCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCa          <  2:125871/88‑1 (MQ=255)
                                                                     caaaaaCGACCCGAATTAACCCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCa          <  2:746037/88‑1 (MQ=255)
                                                                        aaaCGACCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGc       <  2:93136/90‑1 (MQ=255)
                                                                             aTCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTa  >  1:815547/1‑90 (MQ=255)
                                                                              |                                                                                        
AAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTA  >  NZ_CP009273/3214048‑3214214

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