Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,033,5560CG52.9% 13.0 / 24.6 17L200V (CTG→GTG) msrPprotein‑methionine‑sulfoxide reductase catalytic subunit MsrP
Reads supporting (aligned to +/- strand):  ref base C (0/8);  new base G (9/0);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 4.11e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-04
Rejected as consensus: Frequency below/above cutoff threshold.

GCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAA  >  NZ_CP009273/2033473‑2033638
                                                                                   |                                                                                  
gCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAat                                                                              <  2:256635/90‑1 (MQ=255)
  aaTGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGTGGAaaat                                                                            >  1:398481/1‑90 (MQ=255)
   aTGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGTGGAaaata                                                                           >  1:375712/1‑89 (MQ=255)
   aTGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGTGGAAAATc                                                                           >  1:232174/1‑90 (MQ=255)
        aaGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGTGAAAAAACCtttt                                                                      >  2:1772/1‑88 (MQ=255)
             tGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGGGAAAAAATCTTATTTCgg                                                                 >  1:294247/1‑89 (MQ=255)
                    ttCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGGAAAATCCTTATGTTGAAGGAATg                                                          >  1:484800/1‑90 (MQ=255)
                    ttCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGTGAAAAATCCTTATGTCGAAGGATTg                                                          >  1:139630/1‑90 (MQ=255)
                             aTTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTcgac                                                 <  1:74050/90‑1 (MQ=255)
                                       cGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATTCTTATTTCGAAGGAATGCGTCTTGACGAAGCAATga                                       >  1:348295/1‑89 (MQ=255)
                                         gAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCAt                                     <  2:219772/90‑1 (MQ=255)
                                                tGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTca                              <  1:282604/90‑1 (MQ=255)
                                                      gccagCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACtga                        <  2:398481/90‑1 (MQ=255)
                                                                ccGCTTTATCGGCGGTGGGCTGAAATGTCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAgg              <  1:344482/90‑1 (MQ=255)
                                                                ccGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAgg              <  2:60123/90‑1 (MQ=255)
                                                                ccGCTTTATCGGCGGGGGGGTGGAAAATCCTTTTGTCGAAGGAGTGCGTCTCGACGAAGGAATGCATCCGCCCACACTGATGAACGGAgg              >  2:206349/1‑90 (MQ=255)
                                                                            cggTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCaa  <  1:189806/90‑1 (MQ=255)
                                                                                   |                                                                                  
GCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAA  >  NZ_CP009273/2033473‑2033638

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