Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,105,8450CA54.5% 2.2 / 15.4 11E188D (GAG→GATgspCtype II secretion system protein GspC
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base A (3/3);  total (6/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGCACACTGACAGCAGGCTCAGCAACAGCTTGCTTTGCTTCGTCACTGACGGCTTTTTTGTTGGTCACAGCAACAGTGGAACGTTCCTCTTCTGCCAGACTCAGACGCTCTATTTTTCCCTGATAACGCAGCATCACATGGTCGCGGTTGATTTCCTCAATCACCGCGTTGTG  >  NZ_CP009273/3105759‑3105931
                                                                                      |                                                                                      
ggCACACTGACAGCAGGCTCAGCAACAGCTTGCTTTGCTTCGTCACTGACGGCTTTTTTGTTGGTCACAGCAACAGTGGAACGTTCctct                                                                                     >  2:168217/1‑90 (MQ=255)
          cagcagGCTCAGCAACAGCTTGCTTTGCTTCGTCACTGACGGCTTTTTTGTTGGTCACAGCAACAGTGGAACGTTCCTCTTCTGCCAGAc                                                                           >  1:115216/1‑90 (MQ=255)
              agGCTCAGCAACAGCTTGCTTTGCTTCGTCACTGACGGCTTTTTTGTTGGTCACAGCAACAGTGGAACGTTcc                                                                                        >  1:143117/1‑73 (MQ=255)
              agGCTCAGCAACAGCTTGCTTTGCTTCGTCACTGACGGCTTTTTTGTTGGTCACAGCAACAGTGGAACGTTcc                                                                                        <  2:143117/73‑1 (MQ=255)
                                   tgcttCGTCACTGACGGCTTTTTTGTTGGTCACAGCAACAGTGGAACGTTCATCTTCTGCCAGACTCAGACGCTCTATTTTTCCCTGATa                                                  >  2:207290/1‑90 (MQ=255)
                                   tgcttCGTCACTGACGGCTTTTTTGTTGGTCACAGCAACAGTGGAACGTTCATCTTCTGCCAGACTCAGACGCTCTATTTTTCCCTGATa                                                  >  2:288413/1‑90 (MQ=255)
                                   tgcttCGTCACTGACGGCTTTTTTGTTGGTCACAGCAACAGTGGAACGTTCATCTTCTGCCAGACTCAGACGCTCTATTTTTCCCTGATa                                                  >  2:92945/1‑90 (MQ=255)
                                       tCGTCACTGACGGCTTTTTTGTTGGTCACAGCAACAGTGGAACGTTCCTCTTCTGCCAGACTCAGACGCTCTATTTTTCCCTGATAACgc                                              <  1:85685/90‑1 (MQ=255)
                                                             tgGTCACAGCAACAGTGGAACGTTCATCTTCTGCCAGACTCAGACGCTCTATTTTTCCCTGATAACGCAGCATCACATGGTCGCGGTTGa                        <  1:207290/90‑1 (MQ=255)
                                                             tgGTCACAGCAACAGTGGAACGTTCATCTTCTGCCAGACTCAGACGCTCTATTTTTCCCTGATAACGCAGCATCACATGGTCGCGGTTGa                        <  1:288413/90‑1 (MQ=255)
                                                             tgGTCACAGCAACAGTGGAACGTTCATCTTCTGCCAGACTCAGACGCTCTATTTTTCCCTGATAACGCAGCATCACATGGTCGCGGTTGa                        <  1:92945/90‑1 (MQ=255)
                                                                             tgAACGTTCCTCTTCTGCCAGACTCAGACGCTCTATTTTTCCCTGATAACGCAGCATCACATGGTCGCGGTTGATTTCCTCAATCACcgc        >  1:8992/2‑90 (MQ=255)
                                                                             ggAACGTTCCTCTTCTGCCAGACTCAGACGCTCTATTTTTCCCTGATAACGCAGCATCACATGGTCGCGGTTGATTTCCTCAATCACcgc        >  1:230222/1‑90 (MQ=255)
                                                                                   ttCCTCTTCTGCCAGACTCAGACGCTCTATTTTTCCCTGATAACGCAGCATCACATGGTCGCGGTTGATTTCCTCAATCACCGCGTtgtg  <  1:168217/90‑1 (MQ=255)
                                                                                      |                                                                                      
GGCACACTGACAGCAGGCTCAGCAACAGCTTGCTTTGCTTCGTCACTGACGGCTTTTTTGTTGGTCACAGCAACAGTGGAACGTTCCTCTTCTGCCAGACTCAGACGCTCTATTTTTCCCTGATAACGCAGCATCACATGGTCGCGGTTGATTTCCTCAATCACCGCGTTGTG  >  NZ_CP009273/3105759‑3105931

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