Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,390,9090AC57.1% 6.3 / 11.7 14G82G (GGT→GGGrngribonuclease G
Reads supporting (aligned to +/- strand):  ref base A (5/1);  new base C (0/8);  total (5/9)
Fisher's exact test for biased strand distribution p-value = 3.00e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATAT  >  NZ_CP009273/3390834‑3390979
                                                                           |                                                                      
ccccTGCACCATCAGATCTTCCCCCTGACGAACCAGTTCCGAGATGTCGCGCCCCGTGAATTGCTTTTGTTCTTCCCCCGCCACACATTc                                                          <  1:239241/88‑1 (MQ=255)
  ccTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCgg                                                        >  2:181442/1‑90 (MQ=255)
    tGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGgtg                                                      >  1:96745/1‑90 (MQ=255)
       cccATCAGTTCTTCCCCCTGACGACCCATTCCCGAGATGTCGCGCCCCGTAATTTGTTTTTTTTCTCCCCCCGCCACACATTCGGTGTGc                                                   <  2:225188/89‑1 (MQ=255)
          aTCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGc                                                >  2:123268/1‑90 (MQ=255)
            cAGATTTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCCCCGTGAATTGTTTTTTTTCTTCACCCGCCACACATTCGGTGTGCGGCat                                              <  1:26658/90‑1 (MQ=255)
                                    tccggagaTGCCCCGCCCCGTGATTTTTTTTTTTTCTTCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAAcgc                      <  2:243414/86‑1 (MQ=255)
                                    tcccGAGATGTCGCGCCCCGTGAATTGCTTTTGTTCTTCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAAcgc                      <  1:104460/88‑1 (MQ=255)
                                    tcccGAGATGTCGCGCACCGTGAATTGCTTTTTTTCTTCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAAcgc                      <  1:260307/88‑1 (MQ=255)
                                    tcccGAGATGTCGCCCCCCGTAATTTTTTTTTTTTCTCCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAAcgc                      <  1:623089/88‑1 (MQ=255)
                                    tcccGAGAGGTCGCCCCCCGTAATTTTTTTTTTTTCTCCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAAcgc                      <  2:177630/88‑1 (MQ=255)
                                             gTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATc             >  1:536729/1‑90 (MQ=255)
                                             gTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATc             >  1:655861/1‑90 (MQ=255)
                                                        tGATTTGTTTTTTTTCTTCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAatat  <  1:641299/90‑1 (MQ=255)
                                                                           |                                                                      
CACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATAT  >  NZ_CP009273/3390834‑3390979

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