Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,336,4920TG53.3% 10.9 / 19.0 15T218P (ACC→CCC) fumBclass I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (5/2);  new base G (8/0);  total (13/2)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-04
Rejected as consensus: Frequency below/above cutoff threshold.

TGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGC  >  NZ_CP009273/4336408‑4336581
                                                                                    |                                                                                          
tGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCCGTACCGAGGGGACGCa                                                                                       >  1:204011/1‑90 (MQ=255)
               tggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGGTGCAGTACCCAGGGGACGCATTTTTTCTACGAggg                                                                       >  1:388508/1‑89 (MQ=255)
                                    cAATCACAAACGCGATATGGTACGGCGGGGAGGGGGCAGTAACGAGGGGACGCAATTTTTCGGCGAGGGAGGTTTTTTGTTTGCCggggg                                                  >  2:123598/1‑88 (MQ=255)
                                            aaCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGGGGGGACGCATTTTTTCGACGAGGGAGGTTTTTTGTTTTCCGGGGGGCACCAgg                                          >  1:631807/1‑90 (MQ=255)
                                                 gATATGGTACGGCGGGCGGGGAGCAGGACCGAGGGGACGCACTTTCTCGAGCAGGGAAGTTTTTAGTTTTGCCGGAGTGGGCAGGGGttt                                     >  1:338225/1‑90 (MQ=255)
                                                  atatGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTGCCGGGGGGCAGCAGGGCTTTg                                    >  1:270757/1‑90 (MQ=255)
                                                            gcggGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTGCCGGGGGGCAGCAGGGCTTTTGGTTCCTGGt                          >  1:545402/1‑90 (MQ=255)
                                                            gcggGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTCAGTTTGCCGGGGGGCAGCAGGGCTTTGGTTTTCTGGt                          >  1:173159/1‑90 (MQ=255)
                                                                             ccGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTttgttg         >  1:455525/1‑90 (MQ=255)
                                                                               gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTcagca                                            <  1:654581/53‑1 (MQ=255)
                                                                               gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTcagca                                            >  2:654581/1‑53 (MQ=255)
                                                                               gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTgg                           >  1:493032/1‑70 (MQ=255)
                                                                               gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTgg                           <  2:493032/70‑1 (MQ=255)
                                                                               gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGc       >  1:371106/1‑90 (MQ=255)
                                                                                 gggTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCag     >  1:389672/1‑90 (MQ=255)
                                                                                    tACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAgc  >  1:378471/1‑90 (MQ=255)
                                                                                    tACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAgc  >  2:241539/1‑90 (MQ=255)
                                                                                    |                                                                                          
TGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGC  >  NZ_CP009273/4336408‑4336581

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