Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,501,8650AC53.3% 9.7 / 15.1 15L123V (TTG→GTG) fecDFe(3+) dicitrate ABC transporter permease subunit FecD
Reads supporting (aligned to +/- strand):  ref base A (5/2);  new base C (0/8);  total (5/10)
Fisher's exact test for biased strand distribution p-value = 6.99e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAG  >  NZ_CP009273/4501776‑4501936
                                                                                         |                                                                       
cAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCaa                                                                         >  1:166017/1‑90 (MQ=255)
                    aaaGCCCCACGCGGGTGAGCCCCAGCTTCATCGGCTGGTGGGTTTTTGCCAGCATCTACAGCAATATCACCCCCGCCATGCCGCCCGCaa                                                     <  1:328535/90‑1 (MQ=255)
                     aaGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCaaa                                                    >  1:20252/1‑90 (MQ=255)
                             cGCCGGTGAGCGCCAGCTTCATCGGGTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCcagca                                            <  1:461151/90‑1 (MQ=255)
                              gCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCcagcag                                           <  1:134406/90‑1 (MQ=255)
                                  gTGGGCCCCAGTTTCACCGGCTGGTGGGTTTTGGCCAGCATTTTCAGAAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGgc                                       <  2:96019/90‑1 (MQ=255)
                                       cgcCAGCTTCATCGGCTGGTGGGTTTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCAc                                  <  1:195102/90‑1 (MQ=255)
                                         ccAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCa                                >  1:49412/1‑90 (MQ=255)
                                             cTTCATCGGCTGGGGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCAc                            <  1:172715/90‑1 (MQ=255)
                                                cATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACggg                         >  1:361951/1‑90 (MQ=255)
                                                       tggtggGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGAc                  >  2:82633/1‑90 (MQ=255)
                                                              tttttGCCAGCTTCTTCAGCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATaa           <  2:187223/88‑1 (MQ=255)
                                                                      accttctcCAACAAAACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTaga   <  2:211491/86‑1 (MQ=255)
                                                                      aGCATCTCCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTaga   <  1:429926/90‑1 (MQ=255)
                                                                       gCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTagag  <  1:57281/90‑1 (MQ=255)
                                                                       gCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTagag  >  1:193690/1‑90 (MQ=255)
                                                                                         |                                                                       
CAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAG  >  NZ_CP009273/4501776‑4501936

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