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

AGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAG  >  NZ_CP009273/4501777‑4501945
                                                                                        |                                                                                
aGGCTGGCCCAGCAGGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAAc                                                                                 <  1:103005/90‑1 (MQ=255)
                         ccACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCca                                                        >  1:246088/1‑90 (MQ=255)
                                     gcgcCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATg                                                                        <  1:49488/62‑1 (MQ=255)
                                     gcgcCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATg                                                                        >  2:49488/1‑62 (MQ=255)
                                          aGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCcagca                                                     <  1:226062/76‑1 (MQ=255)
                                          aGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCcagca                                                     >  2:226062/1‑76 (MQ=255)
                                                    gcgggtggGTTTTGGCCAGCATCTTCAGAAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTg                             <  2:613494/87‑1 (MQ=255)
                                                      tggtggGTTTTTGCCAGCATCTTCAGCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGAc                           <  2:540064/90‑1 (MQ=255)
                                                            gTCTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATa                     <  1:42550/90‑1 (MQ=255)
                                                                     aGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTaga            <  1:79551/90‑1 (MQ=255)
                                                                     aGCATCTTCAGCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTaga            <  2:188644/90‑1 (MQ=255)
                                                                     aGCATCTTCAGAAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTaga            <  2:547580/90‑1 (MQ=255)
                                                                          cTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGcccc       >  1:398882/1‑90 (MQ=255)
                                                                               gaaaTATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACag  <  2:229562/88‑1 (MQ=255)
                                                                               gCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACag  <  1:260511/90‑1 (MQ=255)
                                                                                        |                                                                                
AGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAG  >  NZ_CP009273/4501777‑4501945

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