Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,001,8320TG57.1% 5.7 / 10.7 14G25G (GGT→GGGyedEselenium metabolism membrane protein YedE/FdhT
Reads supporting (aligned to +/- strand):  ref base T (2/4);  new base G (8/0);  total (10/4)
Fisher's exact test for biased strand distribution p-value = 1.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.75e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAAC  >  NZ_CP009273/2001750‑2001920
                                                                                  |                                                                                        
cTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATtctct                                                                                   <  2:424341/90‑1 (MQ=255)
         tgtcatgGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTCTCtactttct                                                                          >  2:127344/1‑90 (MQ=255)
          gtcatgGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTACTTACTa                                                                         >  1:406760/1‑90 (MQ=255)
            catgGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTTTCTACTTACTAtt                                                                       >  1:131496/1‑90 (MQ=255)
                           aaCACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCa                                                        <  1:199712/90‑1 (MQ=255)
                           aaCACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCa                                                        <  1:276588/90‑1 (MQ=255)
                              acGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTCTCTACTTACTATTTTTGCATTACTGGCACCt                                                     >  1:182353/1‑90 (MQ=255)
                               cGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTCTCTACTTACTATTTTTGCATTACTGGCACCtt                                                    >  2:564720/1‑90 (MQ=255)
                                       tGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTg                                            >  1:233304/1‑90 (MQ=255)
                                         aTTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTTGCATTACTGGCACCTTTTTGGCTGTc                                          >  2:509857/1‑90 (MQ=255)
                                         aTTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTc                                          >  1:393594/1‑90 (MQ=255)
                                                   gggCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGGGAAt                                >  2:674412/1‑90 (MQ=255)
                                                               cTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGGGAATTTACCCGTTggg                    >  1:99381/1‑90 (MQ=255)
                                                                                 gTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAAc  <  2:131496/90‑1 (MQ=255)
                                                                                  |                                                                                        
CTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAAC  >  NZ_CP009273/2001750‑2001920

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