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

GTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTA  >  NZ_CP009273/3214038‑3214191
                                                                                |                                                                         
gttgATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATc                                                                  <  2:57182/90‑1 (MQ=255)
                        tttGCCAGCGAGTTGTCCCGCCGGCGCGTACGGGGGGCCGAAACGCTCAAAAACGCAAAACCGACCCGAATTAAACCGCcactgacactg                                          <  2:717069/90‑1 (MQ=255)
                        tttGCCAGCGAGATTTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTAAATAACGCAAAAACGATCCGCATTGAACCGCcactgacactg                                          <  2:558778/90‑1 (MQ=255)
                        ttggCCAGGGAGTTTTCCCCCCGGCGCGTACGGGGGGCCGAACCGCTAAAAAACGCAAAAACGCCCCGCATTAAACCGCcactgacactg                                          <  2:544716/87‑1 (MQ=255)
                         ttGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCcactgacactga                                         >  2:309041/1‑90 (MQ=255)
                                 gagtTGTCCCGCCGGCGCGTCCGGGGGGCCGGAACGCTCAAAAACCCAAAAACGACCCGAATTGAACCGCCACTGACACTGACCATTGaa                                 <  2:631580/90‑1 (MQ=255)
                                  agttGTCCCGCCAGCGCGTCCGGGTGGCCGGAACGCTAAAAAACCCAAAAACGACCCGCATTGAACCGCCACTGACACTGACCATTGAAc                                <  2:691795/87‑1 (MQ=255)
                                  agttGTCCCGCCAGCGCGTACTGGGGGCCGGAACGCTCAATAACGCAAAAACGACCCGCATTGAACCGCCACTGACACTGACCATTGAAc                                <  2:7768/87‑1 (MQ=255)
                                  agttGTCCCCCCGGCGGGTGCTGGGGGCCGGAACGCTAAAAACCCCAAAAACGACCCGAATGAAACCGCCACTGACACTGACCATTGAAc                                <  2:376180/87‑1 (MQ=255)
                                  agttGTCCCCCCGGCGCGTCGGGGGGGCCGGACCCCTAAAAACCCCAAAAACGACCCGCTTTAAACCGCCACTGACACTGACCATTGAAc                                <  1:65737/87‑1 (MQ=255)
                                    ctGTTCCGCCGGTGCTTCCTGGCGGCCGGACCGTTAAATACCTAAAAACCCCTCCGTGTTGACCCGCCACTGACACTGACCATTGAACAg                              <  2:611214/89‑1 (MQ=255)
                                           gCCCGCGCGTACGGGGGGCCGAACCGCCAAAAAACCAAAAAACGTCCCGCATTAAACCGCCACTGACACTGACCATTGAACAGTGTGAAc                       <  1:144675/90‑1 (MQ=255)
                                               gcgcGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGt                   >  1:190161/1‑90 (MQ=255)
                                                cgcgTACGGGTGGCCGGAACGCTAAAAAACCCAAAACCGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTg                  <  2:190161/90‑1 (MQ=255)
                                                         gtggCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAgcg         >  2:261362/1‑90 (MQ=255)
                                                           ggCGGGAACGCTAAAAAACCCAAAAACGACCCGAATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAgcggc       <  2:296522/90‑1 (MQ=255)
                                                                gAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTa  >  2:122510/1‑90 (MQ=255)
                                                                gAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTa  >  1:158212/1‑90 (MQ=255)
                                                                                |                                                                         
GTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTA  >  NZ_CP009273/3214038‑3214191

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