Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,436,6360AC68.8% ‑1.8 / 11.2 16V274G (GTA→GGA) secYpreprotein translocase subunit SecY
Reads supporting (aligned to +/- strand):  ref base A (5/0);  new base C (0/11);  total (5/11)
Fisher's exact test for biased strand distribution p-value = 2.29e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGT  >  NZ_CP009273/3436558‑3436714
                                                                              |                                                                              
cAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCata                                                                     >  2:66503/1‑90 (MQ=255)
          cccccccAACCCAGGGCGCGTTGGGCGCCGGAAAAAAAATAATACTGGAAGGAAAGATGGCGGGAATACCCCCCGCCATATTCACTTTCa                                                           <  2:379252/88‑1 (MQ=255)
          cccccccAACCCAGGCCGCGAGGGCCGCCGGAAAAAAAAAAATCCGGGAAGAAAAGTTTGCCGGGATCCCCCCCGCCATATTCACTTTCa                                                           <  1:80747/88‑1 (MQ=255)
           gCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAg                                                          >  2:35216/1‑90 (MQ=255)
            cccccGACCCAGGACGGGGTGGCCGCGGGAAAAAAAATAATACTGGAAGCGAAGATGGCGGGGTTTCCCCCCGCCATATTCACTTTCAGc                                                         <  2:350107/90‑1 (MQ=255)
                 gAACCATGACGCGAGGGTCGCCGGAAAAAAAATAATACTGGAAGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTaa                                                    <  2:195893/90‑1 (MQ=255)
                      aTGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtg                                               >  2:206792/1‑90 (MQ=255)
                           ggggtGGCCGCCGGAAAAAAAATAATACTGAAAGAAAAGTTGGCCGGATTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctc                                          <  1:99518/86‑1 (MQ=255)
                           gggATGGTCCCCGGAAAAAAAATAATCCGGGAAGAGAAGTTTGCCGGTTTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctc                                          <  1:161210/88‑1 (MQ=255)
                           gggATGGCCGCGGGAAAAAAAAAAATACTGAAAGGAAAGTTTGCGGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGttc                                          <  2:299072/88‑1 (MQ=255)
                           gcgcTGGTCGCCGGGAAGAAAAAAATACTGGGAGAAAAGGTTGCGGGCATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctc                                          <  1:459894/90‑1 (MQ=255)
                               tGGCCGCGGGAAACAAAATAATACTGGAAGAGAAGTTTGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg                                      <  2:211270/90‑1 (MQ=255)
                                   cGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgcagc                                  >  2:459894/1‑90 (MQ=255)
                                                          aaGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGAc           >  1:536114/1‑90 (MQ=255)
                                                                   ttGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt  <  2:78358/90‑1 (MQ=255)
                                                                   tggCCGGAATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt  <  2:300359/88‑1 (MQ=255)
                                                                              |                                                                              
CAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGT  >  NZ_CP009273/3436558‑3436714

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